（2.1.2有阻尼三相突然短路中突然短路微分方程的函数程序）

电机仿真,含阻尼微分方程的M函数程序和仿真图

%编写求解空载建立电压微分方程的函数（sy_ge_damp_3duanlu_sol）
function sy_ge_damp_3duanlu_sol(t,y)
%设定时间和初始值
tspan=[0 6];
y0=[20;0;0;0;0];
%方程求解
[t,y]=ode113('sy_ge_damp_3duanlu_ode',tspan,y0);
%绘制定子电流和励磁电流变化规律
subplot(3，1，1)；
plot(t,y(:,1));
%下面指定纵横轴标签
xlabel('Time[s]')
ylabel('Id[A]')
hold on
subplot(3,1,2);
plot(t,y(:,3));
%下面指定纵横轴标签
xlabel('Time[s]')
ylabel('Id[A]')
hold on
w=314
ia=y(:,1).*sin(w*t)+y(:,2).*cos(w*t);
subplot(3,1,3);
plot(t,ia)
%下面指定纵横轴标签
xlabel('Time[s]')
ylabel('Id[A]')