File:ACPower-en.svg

import numpy as np
import matplotlib.pyplot as plt


def current(i_hat, omega, theta_u, theta_i, r):
    return i_hat*np.cos(omega*r + theta_u - theta_i)

def voltage(u_hat, omega, r):
    return u_hat*np.cos(omega*r)

def power(P, omega, r):
    return P + P*np.cos(2*omega*r)

# parameteres

omega = theta_u = theta_i = P = 1
u_hat = 3
i_hat = 2

# x axes

rad = np.linspace(0, 2*np.pi, 100)
deg = np.linspace(0, 360, 100)

# y axes

v = voltage(u_hat, omega, rad)
c = current(i_hat, omega, theta_u, theta_i, rad)
p = power(P, omega, rad)
a = np.mean(p)*np.ones(100)

# plot values

for y in (v, c, p, a):
    plt.plot(deg, y, linewidth=2)

plt.xlabel("Degrees", fontsize=10)
plt.ylabel("Volt (V) / Ampere (A) / Watt (W) / Watt (W)", fontsize=10)
plt.legend(["Voltage", "Current", "Power", "Average"], loc=4, fontsize=15)
plt.savefig("ACPower-en.svg")
plt.show()