% 編寫求解自激過程微分方程的腳本函數
% 將該腳本函數定義為sh_ge_se_ex_so(shunt_generator_self_excited_solver)
[t,iff]=ode23(@sh_ge_se_ex_ode,[0 10],[0])
% 指定全局變量（注意全局變量必須在同時使用的多個函數中同時指定）
global   a1 a2 a3 a4 Rf
% 繪圖（這里采用多子圖法同時繪制ea=f(if)、uf=f(if)、ea=f(t)、if=f(t)等多條曲線）
% 首先繪制ea=f(if)和uf=f(if)曲線
subplot(3,1,1),
ea=a1*iff.^3+a2*iff.^2+a3*iff.^1+a4;
plot(iff,ea)
hold on
uf=Rf.*iff;
plot(iff,uf,'r'); % 'r'表示用紅顏色繪制曲線
% 下面指定縱橫軸標簽
xlabel('If[A]')
ylabel('Ea\Uf[V]')
% 下面繪制if=f(t)曲線
subplot(3,1,2),
plot(t,iff);
xlabel('Time[T]')
ylabel('If[A]')
% 下面繪制ea=f(t)曲線
subplot(3,1,3),
plot(t,(a1*iff.^3+a2*iff.^2+a3*iff.^1+a4));
xlabel('Time[T]')
ylabel('ea[V]')
% 下面計算勵磁電流和感應電勢的數值，然后通過指令窗口讀取二者的穩態值
iff,a1*iff.^3+a2*iff.^2+a3*iff.^1+a4