RLC电路matlab仿真m文件

L=0.5;C=0.02;%设定RLC元件的参数满足R≤2L/C即R≤225=10
uc0=1;iL0=0;%设定储能元件上的初值uc（0+）=uc（0-）=Us=1ViL（0+）=iL（0-）=0A

for R=2:2:10%设定R的变化范围为2Ω4Ω6Ω8Ω10Ω
alpha=R/（2*L）;w0=sqrt（1/（L*C））;%定义衰减因子α=R2L和电路谐振角频率ω0=1LC

dt=0.001;t=0:dt:1;%设置时间步进长度为0.001秒仿真时间从0～1秒 

%方法1?微分方程法
wd=sqrt（w0^2-alpha^2）;%设定阻尼振荡角频率ω0=ω02-α2

theta=atan（wd/alpha）;%设定θ=arctanωd/αuc1=（w0/wd）3uc03exp（-alpha3t）.3sin（wd3t+theta）;?????%uc（t）的解

iL1=-uc0/（wd*L）*exp（-alpha*t）.*sin（wd*t）;%iL（t）的解
figure（1）plot（tuc0）hold onfigure（2）plot（tiL1）hold on
%方法2?拉普拉斯变换法
num=[uc0R/L*uc0+iL0/C];%由Uc（s）的分子系数构成的数组
den=[12*alphaw0^2];%由Uc（s）的分母系数构成的数组
[rpk]=