sh_ge_se_ex_ode.m

包里有直流发电机simulink模型

% 编写建立自激过程微分方程的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);