% 編寫建立兩相空載短路微分方程的M—函數
% 將該M—函數定義為tongbufadianji
function dydt=tongbufadianji(t,y)
% 下面輸入電機基本數據： 
r=.26;Rfd=0.13;Rkd=0.0224;Rkq=0.02;Ufd=17.8876;w=377;Ll=1.14e-3;Lmd=13.7e-3;Lmq=11.0e-3;
Llfd=2.1e-3;Llkd=1.4e-3;Llkq=1e-3;
Lfd=Llfd+Lmd;Lkd=Llkd+Lmd;Lkq=Llkq+Lmq;Ld=Lmd+Ll;Lq=Lmq+Ll;
%--------------------------------------------------------------------------------------------------------------------------------------   
%電感矩陣 
 L=[  
        Lfd,          Lmd,         0,            -Lmd*sin(w*t);
        Lmd,          Lkd,         0,            -Lmd*sin(w*t);
        0,            0,           Lkq,          -Lmq*cos(w*t);
        Lmd*sin(w*t), Lmd*sin(w*t),Lmq*cos(w*t), -(Ld+Lq)/2+(Ld-Lq)/2*cos(2*w*t)] ;
G=[ 
         0,       0,        0,  -Lmd*cos(w*t);
         0,       0,        0,  -Lmd*cos(w*t);
         0,       0,        0,  Lmq*sin(w*t);
     Lmd*cos(w*t), Lmd*cos(w*t),-Lmq*sin(w*t), -(Ld-Lq)*sin(2*w*t)];
 %--------------------------------------------------------------------------------------------------------------------------------------   
% 電阻矩陣
 R=[ 
       Rfd,  0,   0,  0+0+0;
       0,    Rkd, 0,  0+0+0;
       0,    0,   Rkq,0+0+0;
      0,    0,   0,   -r];
 %--------------------------------------------------------------------------------------------------------------------------------------   
 %電壓向量
Ualphabeta=[17.8876,0,0,0]';
 %微分方程
dydt=L\(Ualphabeta-w*G*y-R*y);





% 下面列寫標準形式微分方程并將感應電勢用勵磁電流的擬合多項式函數表示：
% 這里的a(1)、a(2)、a(3)、a(4)是指上一步的多項式擬合系數a(i)
%dydt=[1/Lf*(a(1)*iff^3+a(2)*iff^2+a(3)*iff^1+a(4))-Rf/Lf*iff];
% 給全局變量賦值
%a1=a(1);a2=a(2);a3=a(3);a4=a(4);