%編寫求解空在建立電壓微分方程的函數（sy_ge_damp_noload_sol)
function sy_ge_damp_noload_sol
% 設定時間和初始值
tspan=[0 10];
y0=[0;0;0;0;];
% 方程求解
[t,y]=ode113('tongbufadianji',tspan,y0);
% 空載建立電壓過程中勵磁電流的變化規律
ifd=y(:,3);
plot(t,ifd)
%下面指定縱橫軸標簽
xlabel('time[s]')
ylabel('ifd[A]')