sh_ge_start_so.m

包里有直流发电机simulink模型

% 编写求解起动过程微分方程的脚本函数
% 将该脚本函数定义为sh_ge_start_so(shunt_motor_self_excited_solver)
uf=25;IaN=633;w=61;TLN=3275;Lf=9;Rf=12;
[t,y]=ode23(@sh_ge_start_ode,[0:.0001:.025],[uf/Rf;0])
% 首先绘制if=f(t)曲线
subplot(2,1,1),
plot(t,y(:,1))
xlabel('Time[T]')
ylabel('If[A]')
grid on
% 下面绘制ia=f(t)曲线
subplot(2,1,2),
plot(t,-y(:,2)),
xlabel('Time[T]')
ylabel('Ia[A]')
grid on