% 編寫建立自激過(guò)程微分方程的M—函數(shù)
% 將該M—函數(shù)定義為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  
% 下面輸入電機(jī)基本數(shù)據(jù)：
Lf=22;Rf=22.8;
% 下面輸入空載特性試驗(yàn)數(shù)據(jù)(Ifdata是勵(lì)磁電流的實(shí)驗(yàn)數(shù)據(jù),Eadata是感應(yīng)電勢(shì)的實(shí)驗(yàn)數(shù)據(jù))：
Ifdata=[ 0,   2, 3.2, 4.8,	5.5, 6.5, 8.2, 11.7	]; 
Eadata=[10, 100, 150, 200,  220, 244, 266,  281 ]; 
% 下面進(jìn)行曲線擬合：
a=polyfit(Ifdata,Eadata,3);% 求取空載特性的多項(xiàng)式擬合系數(shù)a(i)
% 下面列寫標(biāo)準(zhǔn)形式微分方程并將感應(yīng)電勢(shì)用勵(lì)磁電流的擬合多項(xiàng)式函數(shù)表示：
% 這里的a(1)、a(2)、a(3)、a(4)是指上一步的多項(xiàng)式擬合系數(shù)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);