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]=