dydt3.m

四阶龙格库塔求解电机启动的仿真程序 ！

function Yd=dydt3(t,Y)
a=Y(1);b=Y(2);c=Y(3);d=Y(4);w=Y(5);o=Y(6);
uds=0.5*cos(t-o);uqs=0.5*sin(t-o); 
a1=(uds+0.057445*c+(2.21284*b+2.12102*d)*w-0.06574*a)/0.21257;
b1=(uqs+0.057445*d-(2.21284*a+2.12102*c)*w-0.06574*b)/0.21257;
%c1=(-0.06355*c-uds-(2.21284*b+2.12102*d)*w+0.06574*a)/0.22538;
c1=-(2.12102*a1+0.060913*c)/2.24906;
%d1=(-0.06355*d-uqs+(2.21284*a+2.12102*c)*w+0.06574*b)/0.22538;
d1=-(2.12102*b1+0.060913*d)/2.24906;
w1=(2.12102*(b*c-a*d)-0.007012*w)/16.39;
o1=w;
Yd=[a1;b1;c1;d1;w1;o1];
end