% 編寫求解空載建立電壓微分方程的函數(shù)(sy_ge_nodamp_rduanlu_sol)

function sy_ge_nodamp_rduanlu_sol 
%設(shè)定時間和初始值
tspan = [0  10];
y0 = [0;0;0];
% 方程求解
[t,y]=ode113('sy_ge_nodamp_rduanlu_ode',tspan,y0);
%繪制定子電流和勵磁電流變化規(guī)律
subplot(2,1,1)
plot(t,y(:,3))

xlabel('Time[s]')
ylabel('Ifd[A]')
hold on
subplot(2,1,2)
w=377;
ia=y(:,1).*sin(w*t)+y(:,2).*cos(w*t);
plot(t,ia)
xlabel('Time[s]')
ylabel('Ia[A]')