% 編寫異步發電機dqo坐標下的局部數學模型M—函數
% 將該M—函數定義為asy_ge,其中x為電流(id1,iq1,id2,iq2)的一階導數,y=(id1,iq1,id2,iq2,ud1,uq1,w)
function x=asy_ge(y)
% 下面輸入電機基本數據：
LR0=0.01391719745223;LS0=0.01391719745223;RR=2.16;RS=2.27;
% 下面計算互感LM與激磁電流之間的關系(考慮磁路飽和作用)：
mm=((y(1)+y(3))^2+(y(2)+y(4))^2)^.5;
LM=(0.2209268+0.1829319*mm-0.1096410*mm^2+0.0257472*mm^3-0.0028836*mm^4);
LR=LM+LR0;LS=LM+LS0;
%下面計算系數矩陣A和B
A=1/(LS*LR-LM^2)*[-RS*LR,      y(7)*LM^2,   RR*LM,       y(7)*LM*LR;
                  -y(7)*LM^2, -RS*LR,      -y(7)*LM*LR,  RR*LM;
                   RS*LM,     -y(7)*LM*LS, -RR*LS,      -y(7)*LS*LR;
                   y(7)*LM*LS, RS*LM,       y(7)*LS*LR, -RR*LS];
B=1/(LS*LR-LM^2)*[LR,0;0,LR;-LM,0;0,-LM];
%下面計算x與y之間關系
x=[(A*[y(1);y(2);y(3);y(4)]+B*[y(5);y(6)])];