MATLAB程序--RLC振荡电路分析

figure(1)
L=1;C=10^(-6);R=200;
uc0=10;iL0=0;
alpha=R/(2*L);w0=sqrt(1/(L*C));
w=sqrt(w0^2-alpha^2);
dt=0.00000001;t=0:dt:0.01;
A=uc0*w0/w;
a=atan(alpha/w);
uct=A*exp(-alpha*t).*cos(w*t-a);
uLt=A*exp(-alpha*t).*cos(w*t+a);
iLt=A*exp(-alpha*t).*cos(w*t-pi/2)/w0*L;
subplot(3,1,1),plot(t,uct),title('Uc(t) R<2000'),xlabel('time'),ylabel('Uc(t)');
subplot(3,1,2),plot(t,iLt),xlabel('time'),ylabel('iL(t)');
subplot(3,1,3),plot(t,uLt),xlabel('time'),ylabel('uL(t)');
figure(2)
L=1;C=10^(-6);R=4000;
uc0=10;iL0=0;
alpha=R/(2*L);w0=sqrt(1/(L*C));
p1=-alpha+sqrt(alpha^2-w0^2);
p2=-alpha-sqrt(alpha^2-w0^2);
dt=0.00000001;t=0:dt:0.01;
A=uc0*p2/(p2-p1);
B=-uc0*p1/(p2-p1);
uct=A*exp(p1*t)+B*exp(p2*t);
uLt=B*exp(p1*t)+A*exp(p2*t);
iLt=-(A/p2*exp(p1*t)+B/p1*exp(p2*t))/L;
subplot(3,1,1),plot(t,uct),title('Uc(t) R>2000'),xlabel('time'),ylabel('Uc(t)');
subplot(3,1,2),plot(t,iLt),xlabel('time'),ylabel('iL(t)');
subplot(3,1,3),plot(t,uLt),xlabel('time'),ylabel('uL(t)');
figure(3)