% 編寫起動過程微分方程的M—函數
% 將該M—函數定義為sh_motor_vary_ode(shunt_ motor_self_ode)
function dydt=sh_motor_vary_ode(t,y)
% 下面輸入電機基本數據：
uf=250;u=220;IaN=633;nN=617;TLN=3275;Lf=9;Rf=12;Ra=0.012;La=0.00035;Laf=0.18;J=30;Bm=0;
% 下面列寫標準形式微分方程：
dydt=[uf/Lf-Rf/Lf*y(1);u/La-Ra/La*y(2)-Laf/La*y(1)*y(3);Laf/J*y(1)*y(2)-Bm/J*y(3)-TLN/J];