sm phasor diagram L=length(t) Ifd=iff(L) Ifq=0 Efd=0 Efq=Lmd*iff(L) Id=i(L,1) Iq=i(L,2) I=Id+j*Iq Vd=v(L,1) Vq=v(L,2) V=Vd+j*Vq VILq=-(Ra+Lq*j)*I Vqp=V+VILq Vqpd=real(Vqp) Vqpq=imag(Vqp) VLd=vL(L,1) VLq=vL(L,2) - sm phasor di

从初值出发，先计算y对时间t的导数，得到y的变化速度；求出x对t的导数，得到x随时间的变化速度。然后时间前进一小步h=0.002，同时x和y也都发生一个小的变化，即在原来数值基础上加上时间步长h和变化速度的乘积。如果x到达1或者-1，强制让y发生变化。这个过程不断迭代，得到结果。 -completed chart