Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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,i t)
- 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
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement