【文章內(nèi)容簡介】
x2*z1。 i1=x2*du39。 y=[,]。% part 2 t= 。 for i=1:5 b(i)=y(i)。 end dy=z2*(y).39。+i1。 tt(1)=t。 di(1)=y(1)+di0。 qi(1)=y(2)+qi0。 fi(1)=y(3)+fi0。 zdi(1)=y(4)。 zqi(1)=y(5)。 ai(1)=di(1)*cos(t+cita0)qi(1)*sin(t+cita0)。 tm(1)=qi(1)*(xd*di(1)+xaf*fi(1)+xazd*zdi(1))di(1)*(xq*qi(1)+ xazq*zqi(1))。 % while (tt(n)p)% k=0。 for n=2:10 a(1)=h/2。 a(2)=a(1)。 a(3)=h。 a(4)=h。 for k=1:3 for i=1:5 c(i)=b(i)+a(k)*dy(i)。 y(i)=y(i)+a(k+1)*dy(i)/。 end dy=z2*(c).39。+i1。 end for i=1:5 y(i)=y(i)+h*dy(i)/。 end t=t+h。 for i=1:5 b(i)=y(i)。 end dy=z2*(y).39。+i1。 % k=k+1。% part 3 tt(n)=t。 di(n)=y(1)+di0。 qi(n)=y(2)+qi0。 fi(n)=y(3)+fi0。 zdi(n)=y(4)。 zqi(n)=y(5)。 ai(n)=di(n)*cos(t+cita0)qi(n)*sin(t+cit