w = 2 * pi * 50
R = 5
C1 = 0.0005
C2 = 0.0009
C3 = 0.0011
m = 50
n = 100
U = 1
T = R*C1
t = 0:T/n:m*2*C3
A = asin(w/(sqrt((1/T)ˆ2 + wˆ2))) ;
it = U/R *(1/(sqrt(1+(w*T)ˆ2))) * (-exp(-t/T) * sin(A) + cos(w*t-A)*w*T);
plot(t,it)
hold on
T = R*C2
t = 0:T/n:m*2*C3
A = asin(w/(sqrt((1/T)ˆ2 + wˆ2)));
it2 = U/R * (1/(sqrt(1+(w*T)ˆ2))) * (-exp(-t/T) * sin(A) + cos(w*t - A) *w*T);
plot (t,it2,'r - - ')
hold on
T = R*C3
t = 0:T/n:m*2*C3
A = asin(w/(sqrt((1/T)ˆ2 + wˆ2)));
it3 = U/R * (1/(sqrt(1+(w*T)ˆ2)))*(-exp(-t/T) * sin(A) + cos(w*t - A)*w*T);
plot(t, it3, 'g: ')
legend('C1 : 5 nmu F' ,'C2 : 9 nmu F' ,'C3 : 11 nmu F' ,1);
xlabel('t');
ylabel('i(t)');
hold off