硕527 马彦昭

留数,其格式为:其中num和den分别为分子、分母多项式系数组成的数组。进而写出:这样就无需求出其显式,使得程序特别简明。II、Matlab程序:clear,closeall,pactUs=72;R1=5000;R2=7500;L=2.5;C=1e-7;%给出原始数据R=R1+R2;dt=0.00001;t=0:dt:0.01;%设定时间数组uc0=Us;iL0=0;%算出初值iL0及uc0Va0=uc0+iL0*R1%算出初值Va0disp('dVa/dt初值(初值计算)=');dVa=-iL0/C-R1*(uc0-(R1+R2)*iL0)/L;%算出初值dVa/dtdisp(dVa)num=[uc0,R/L*uc0+iL0/C];%uc(s)的分子系数多项式den=[1,R/L,1/L/C];%uc(s)的分母系数多项式[r,s,k]=residue(num,den);%=r(1)*exp(s(1)*t)+r(2)*exp(s(2)*t);%求时域函数iLn=C*diff(ucn)/dt;%求导得到电流iLniLn=[iL0,iLn];Va=ucn+iLn*R1;%计算Va的值dVa=diff(Va)/dt;%通过差分计算初值dVa/dtdisp('dVa/dt初值(差分计算)=');disp(dVa(2));subplot(3,1,1),plot(t,ucn),grid%绘图xlabel('t'),ylabel('ucn')subplot(3,1,2),plot(t,iLn),gridxlabel('t'),ylabel('iLn')subplot(3,1,3),plot(t,Va),gridxlabel('t'),ylabel('Va')III、程序运行结果:通过差分求得的dva/dt初始值为-1.4113e+005,这与通过初始条件计算的制-1.44e+5差不多。