h,a33);Рk=(k11+2*k22+2*k33+k44)/6;Рw=w(1)+h*k;Рfunction y=pw(a,c)Рfi=1;B=0.1;J=0.002;Рy=(fi*a-B*c)/J;Рfunction y=f11(w,iq,Lq)Рfi=1;Rq=0.5;Рy=(-w*fi-Rq*iq)/Lq;РrungeonРfunction [iq,w]=rungeon(iq0,w0,Lq,Udc)Рh=0.00005;Рw(1)=w0;iq(1)=iq0;Рk1=f1(w(1),iq0,Lq,Udc);a1=iq(1)+h*k1/2;Рk2=f1(w(1)+0.5*h,a1,Lq,Udc);a2=iq(1)+h*k2/2;Рk3=f1(w(1)+0.5*h,a2,Lq,Udc);a3=iq(1)+h*k3;Рk4=f1(w(1)+h,a3,Lq,Udc);Рk=(k1+2*k2+2*k3+k4)/6;Рiq=iq(1)+h*k;Рk11=pw(iq,w(1));a11=w(1)+h*k11/2;Рk22=pw(iq+0.5*h,a11);a22=w(1)+h*k22/2;Рk33=pw(iq+0.5*h,a22);a33=w(1)+h*k33;Рk44=pw(iq+h,a33);Рk=(k11+2*k22+2*k33+k44)/6;Рw=w(1)+h*k;Рfunction y=pw(a,c)Рfi=1;B=0.1;J=0.002;Рy=(fi*a-B*c)/J;Рfunction y=f1(w,iq,Lq,Udc)Рfi=1;Rq=0.5;Udc=200;Рy=(Udc-w*fi-Rq*iq)/Lq; Р五、仿真结果及其分析Р1.仿真结果Р图 1plot(t,w,t,i,iq)Р图 2plot(t,pwm)Р图 3plot(t,w,t,iq,t,pwm*200)
全文预览
