【正文】 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=0。 ff=50。 nx=13。 fn=18000。h=1/ff/80。x=[。1。0。0。0。0。0。0。]。 z=zeros(nx,fn)。 z(:,1)=x。 ak=zeros(nx,1)。ad=[h/2。 h/2。 h。 h]。 for fi=2:fn ak=RLPSSfun(x,nx,t)。 b=x。 for fj=1:3 x=z(:,fi1)+ak.*ad(fj)。 b=b+ak.*ad(fj+1)/3。 ak=RLPSSfun(x,nx,t)。 end x=b+h*ak/6。 z(:,fi)=x。 t=t+h。 end Angles=zeros(1,fn)。Angles(1,:)=z(4,:)。 plot(Angles39。) %Angles=zeros(2,fn)。%Angles(1,:)=z(1,:)。%Angles(2,:)=z(2,:)。 %dA1=Angles(1,:)39。%dA2=Angles(2,:)39。 %save(39。d:\39。,39。dA139。,39。ASCII39。)%save(39。d:\39。,39。dA239。,39。ASCII39。) 4. RLPSSfunfunction ak = RLPSSfun(x,nx,t) ak=zeros(nx,1)。Te=。Ke=100。Xd=。Xq=。Xd1=。Xq1=。r=。Td=。Tq=。Tj=。D=0。Re=。Xe=。Re1=。Xe1=。Kg=20。Ts=。Tch=。Trh=8。a1=1/3。M0=。K0=3。T1=。T2=。T5=。T6=。a=。b=。c=。d=。%a=。b=。c=。d=。V=。Vr=1。Wb=。Wr=1。Ug=Wrx(3)。Ur=Wrx(3)。C4=。C2=。C1=。 if t0 amp。amp。 t= Id=((1+a*rb*Xd1)*(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4)))(a*Xd1+b*r)*(b*x(1)+a*x(2)+c*V*cos(x(4))+d*V*sin(x(4))))/((1+a*rb*Xd1)*(1+a*rb*Xd1)+(a*Xd1+b*r)*(a*Xd1+b*r))。 Iq=(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4))(1+a*rb*Xd1)*Id)/(a*Xd1+b*r)。elseif tamp。amp。t= Id=((r+Re1)*(x(1)+V*sin(x(4)))(Xd1+Xe1)*(x(2)V*cos(x(4))))/((r+Re1)*(r+Re1)+(Xd1+Xe1)*(Xd1+Xe1))。 Iq=(x(1)+V*sin(x(4))(r+Re1)*Id)/(Xd1+Xe1)。 else Id=((r+Re)*(x(1)+V*sin(x(4)))(Xd1+Xe)*(x(2)V*cos(x(4))))/((r+Re)*(r+Re)+(Xd1+Xe)*(Xd1+Xe))。 Iq=(x(1)+V*sin(x(4))(r+Re)*Id)/(Xd1+Xe)。end Vd=x(1)r*IdXd1*Iq。Vq=x(2)r*Iq+Xd1*Id。v2=Vd*Vd+Vq*Vq。Vt=sqrt(v2)。 if x(13)C4 x(13)=C4。elseif x(13)C4 x(13)=C4。end if x(8)+M0C2 x(8)=C2M0。elseif x(8)+M0C1 x(8)=(C1M0)。end ak(1)=(x(1)(XqXq1)*Iq)/Tq。ak(2)=(x(13)x(2)+(XdXd1)*Id)/Td。ak(3)=(x(10)+a1*x(9)D*x(3)(x(1)*Id+x(2)*Iq(Xq1Xd1)*Id*Iq))/Tj。ak(4)=(x(3)1)*Wb。ak(5)=((1T5/T6)*Ugx(5))/T6。ak(6)=((1T5/T6)*(x(5)+T5/T6*Ug)x(6))/T6。ak(7)=((1T5/T6)*(x(6)+T5/T6*(x(5)+T5/T6*Ug))x(7))/T6。ak(8)=(Kg*(x(7)+T5/T6*(x(6)+T5/T6*(x(5)+T5/T6*Ug)))x(8))/Ts。ak(9)=(x(8)+M0x(9))/Tch。ak(10)=((1a1)*x(9)x(10))/Trh。ak(11)=(K0*(1T1/T2)*Urx(11))/T2。ak(12)=((1T1/T2)*(x(11)K0*T1/T2*Ur)x(12))/T2。ak(13)=(Ke*(VrVt+x(12)+T1/T2*(x(11)K0*T1/T2*Ur))x(13)+)/Te。說明：同樣RLPSS和RLPSSfun在同一個文件夾中。5．myfunfunction fit=myfun(pop) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=0。 ff=50。 nx=13。 fn=18000。h=1/ff/80。x=[。1。0。0。0。0。0。0。]。 z=zeros(nx,fn)。 z(:,1)=x。 ak=zeros(nx,1)。ad=[h/2。 h/2。 h。 h]。 for fi=2:fn ak=myfunfun(x,nx,t,pop)。 b=x。 for fj=1:3 x=z(:,fi1)+ak.*ad(fj)。 b=b+ak.*ad(fj+1)/3。 ak=myfunfun(x,nx,t,pop)。 end x=b+h*ak/6。 z(:,fi)=x。 t=t+h。 end %Angles=zeros(1,fn)。%Angles(1,:)=z(4,:)。 %plot(Angles39。) %Angles=zeros(2,fn)。%Angles(1,:)=z(1,:)。%Angles(2,:)=z(2,:)。 %dA1=Angles(1,:)39。%dA2=Angles(2,:)39。 %save(39。d:\39。,39。dA139。,39。ASCII39。)%save(39。d:\39。,39。dA239。,39。ASCII39。)fit=[z(4,:)]*[(z(4,:))]39。5. myfunfunfunction ak = myfunfun(x,nx,t,pop) ak=zeros(nx,1)。Te=。Ke=100。Xd=。Xq=。Xd1=。Xq1=。r=。Td=。Tq=。Tj=。D=0。Re=。Xe=。Re1=。Xe1=。Kg=20。Ts=。Tch=。Trh=8。a1=1/3。M0=。T1=。T2=。T5=。T6=。a=。b=。c=。d=。%a=。b=。c=。d=。V=。Vr=1。Wb=。Wr=1。Ug=Wrx(3)。Ur=Wrx(3)。C4=。C2=。C1=。 if t0 amp。amp。 t= Id=((1+a*rb*Xd1)*(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4)))(a*Xd1+b*r)*(b*x(1)+a*x(2)+c*V*cos(x(4))+d*V*sin(x(4))))/((1+a*rb*Xd1)*(1+a*rb*Xd1)+(a*Xd1+b*r)*(a*Xd1+b*r))。 Iq=(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4))(1+a*rb*Xd1)*Id)/(a*Xd1+b*r)。elseif tamp。amp。t= Id=((r+Re1)*(x(1)+V*sin(x(4)))(Xd1+Xe1)*(x(2)V*cos(x(4))))/((r+Re1)*(r+Re1)+(Xd1+Xe1)*(Xd1+Xe1))。 Iq=(x(1)+V*sin(x(4))(r+Re1)*Id)/(Xd1+Xe1)。 else Id=((r+Re)*(x(1)+V*sin(x(4)))(Xd1+Xe)*(x(2)V*cos(x(4))))/((r+Re)*(r+Re)+(Xd1+Xe)*(Xd1+Xe))。 Iq=(x(1)+V*sin(x(4))(r+Re)*Id)/(Xd1+Xe)。end Vd=x(1)r*IdXd1*Iq。Vq=x(2)r*Iq+Xd1*Id。v2=Vd*Vd+Vq*Vq。Vt=sqrt(v2)。 if x(13)C4 x(13)=C4。elseif x(13)C4 x(13)=C4。end if x(8)+M0C2 x(8)=C2M0。elseif x(8)+M0C1 x(8)=(C1M0)。end ak(1)=(x(1)(XqXq1)*Iq)/Tq。ak(2)=(x(13)x(2)+(XdXd1)*Id)/Td。ak(3)=(x(10)+a1*x(9)D*x(3)(x(1)*Id+x(2)*Iq(Xq1Xd1)*Id*Iq))/Tj。ak(4)=(x(3)1)*Wb。ak(5)=((1T5/T6)*Ugx(5))/T6。ak(6)=((1T5/T6)*(x(5)+T5/T6*Ug)x(6))/T6。ak(7)=((1T5/T6)*(x(6)+T5/T6*(x(5)+T5/T6*Ug))x(7))/T6。ak(8)=(Kg*(x(7)+T5/T6*(x(6)+T5/T6*(x(5)+T5/T6*Ug)))x(8))/Ts。ak(9)=(x(8)+M0x(9))/Tch。ak(10)=((1a1)*x(9)x(10))/Trh。ak(11)=(pop*(1T1/T2)*Urx(11))/T2。ak(12)=((1T1/T2)*(x(11)pop*T1/T2*Ur)x(12))/T2。ak(13)=(Ke*(VrVt+x(12)+T1/T2*(x(11)pop*T1/T2*Ur))x(13)+)/Te。同樣myfun和myfunfun在一個文件夾中。參考文獻致 謝