sh_motor_stop_so.m

来自「经典的50例电机仿真」· M 代码 · 共 22 行

% 编写求解起动过程微分方程的脚本函数
% 将该脚本函数定义为sh_motor_stop_so(shunt_motor_self_excited_solver)
[t,y]=ode23tb(@sh_motor_stop_ode,[0:.001:.1],[ 20.83333333328  ; 873.64064640326 ;  55.87178121392])
% 首先绘制if=f(t)曲线
subplot(3,1,1),
plot(t,y(:,1))
xlabel('Time[T]')
ylabel('If[A]')
grid on
% 下面绘制ia=f(t)曲线
subplot(3,1,2),
plot(t,y(:,2)),
xlabel('Time[T]')
ylabel('Ia[A]')
grid on
% 下面绘制w=f(t)曲线
subplot(3,1,3),
plot(t,y(:,3))
xlabel('Time[T]')
ylabel('w[s/s]')
grid on