sy_ge_nodamp_3duanlu_sol.m

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

% 编写求解空载建立电压微分方程的函数(sy_ge_nodamp_3duanlu_sol)
function sy_ge_nodamp_3duanlu_sol 
%设定时间和初始值
tspan = [0  10];
y0 = [0;0;10];
% 方程求解
[t,y]=ode23tb('sy_ge_nodamp_3duanlu_ode',tspan,y0);
%绘制定子电流和励磁电流变化规律
subplot(3,1,1);
plot(t,y(:,3));
xlabel('Time[s]')
ylabel('Id[A]')
hold on
w=377;
ia=y(:,1).*sin(w*t)+y(:,2).*cos(w*t);
subplot(3,1,2);
plot(t,ia)
xlabel('Time[s]')
ylabel('Ia[A]')
hold on
subplot(3,1,3);
plot(t,y(:,3));
xlabel('Time[s]')
ylabel('Ifd[A]')