% 編寫他勵發電機空載特性分析腳本函數
% 將該函數定義為se_ge_no_load(seperate_generator_no_load)
% 下面輸入電機基本數據：
n1=750;nN=1000;IfN=2.5;
% 下面輸入750轉／分鐘時的空載特性試驗數據(Ifdata是勵磁電流,Eadata是感應電勢)：
Ifdata=[ 0, 0.4, 1.0, 1.6,	2.0, 2.5, 2.6, 3.0,3.6,4.4	]; 
Eadata=[5, 33, 78, 120,  150, 176, 180, 193.5,206,225 ]; 
% 下面繪制空載特性曲線(原始數據)：
plot(Ifdata,Eadata,'*');
hold on
% 下面進行三次樣條插值：
If=0:.01:3.5;% 進行勵磁電流參數重新設置：
Ean1=spline(Ifdata,Eadata,If);% 根據樣點數據求If對應的樣條插值
% 下面進行空載特性曲線重構：
plot(If,Ean1)
%計算額定轉速時的空載特性
EanN=Ean1*nN/n1;
plot(If,EanN,'r')
hold on
%計算空載額定電壓
If=IfN;
EaN=spline(Ifdata,Eadata,If)*nN/n1
%在勵磁電流工作點作一條垂直于電流軸的虛線與空載特性曲線相交，交點對應的電勢即電機的額定工作電勢
plot([IfN,IfN],[0,300],'-.')
hold on
xlabel('If[A]')
ylabel('Ea[V]')