% 編寫他勵發(fā)電機(jī)空載特性分析腳本函數(shù)
% 將該函數(shù)定義為se_ge_no_load(seperate_generator_no_load)
% 下面輸入電機(jī)基本數(shù)據(jù)：
n1=750;nN=1000;IfN=2.5;
% 下面輸入750轉(zhuǎn)／分鐘時的空載特性試驗數(shù)據(jù)(Ifdata是勵磁電流,Eadata是感應(yīng)電勢)：
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 ]; 
% 下面繪制空載特性曲線(原始數(shù)據(jù))：
plot(Ifdata,Eadata,'*');
hold on
% 下面進(jìn)行三次樣條插值：
If=0:.01:3.5;% 進(jìn)行勵磁電流參數(shù)重新設(shè)置：
Ean1=spline(Ifdata,Eadata,If);% 根據(jù)樣點數(shù)據(jù)求If對應(yīng)的樣條插值
% 下面進(jìn)行空載特性曲線重構(gòu)：
plot(If,Ean1)
%計算額定轉(zhuǎn)速時的空載特性
EanN=Ean1*nN/n1;
plot(If,EanN,'r')
hold on
%計算空載額定電壓
If=IfN;
EaN=spline(Ifdata,Eadata,If)*nN/n1
%在勵磁電流工作點作一條垂直于電流軸的虛線與空載特性曲線相交，交點對應(yīng)的電勢即電機(jī)的額定工作電勢
plot([IfN,IfN],[0,300],'-.')
hold on
xlabel('If[A]')
ylabel('Ea[V]')