% 編寫求解空載建立電壓微分方程的函數(sy_ge_damp_rduanlu_sol)
function sy_ge_damp_rduanlu_sol 
%設定時間和初始值
tspan = [0 10];
y0 = [0;0;0;0;0];
% 方程求解
[t,y]=ode113('sy_ge_damp_rduanlu_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('Ifd[A]')
hold on
subplot(3,1,3);
plot(t,y(:,4))
xlabel('Time[s]')
ylabel('Ikd[A]')