% 編寫建立自激過程微分方程的M—函數
% 將該M—函數定義為sh_ge_se_ex_ode(shunt_generator_self_excited_ode)
function dydt=sh_ge_se_ex_ode(t,iff)
% 首先指定全局變量
global   a1 a2 a3 a4 Rf  
% 下面輸入電機基本數據：
Lf=22;Rf=22.8;
% 下面輸入空載特性試驗數據(Ifdata是勵磁電流的實驗數據,Eadata是感應電勢的實驗數據)：
Ifdata=[ 0,   2, 3.2, 4.8,	5.5, 6.5, 8.2, 11.7	]; 
Eadata=[10, 100, 150, 200,  220, 244, 266,  281 ]; 
% 下面進行曲線擬合：
a=polyfit(Ifdata,Eadata,3);% 求取空載特性的多項式擬合系數a(i)
% 下面列寫標準形式微分方程并將感應電勢用勵磁電流的擬合多項式函數表示：
% 這里的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);