%編寫求解空載建立電壓微分方程的函數(fangzhen)
function fangzhen 
%設定初始值
tspan = [0  1];
y0 = [17.8876/0.13;0;0;0];
%解微分方程
[t,y]=ode23tb('tongbufadianji',tspan,y0);
%求解各電流
subplot(3,1,1)
plot(t,y(:,1))
xlabel('Time[s]')
ylabel('Ifd[A]')
hold on 
subplot(3,1,2)
plot(t,y(:,2))
xlabel('Time[s]')
ylabel('Ikd[A]')
hold on
subplot(3,1,3)
plot(t,y(:,4))
xlabel('Time[s]')
ylabel('Ibeta0[A]')