% 直流發電機調整特性分析
% 將該函數定義為dc_ge_set(dc_generator_set)
%--------------------------------------------------------------------------
% 下面輸入電機基本數據：
UN=250;IfN=2.5;Rf=293/12;Ra=1.8;k=.1;
% 下面輸入750轉／分鐘時的空載特性試驗數據(Ifdata是勵磁電流,Eadata是感應電勢)：
 Ifdata=[ 0, 0.4, 1.0, 1.6,	2.0,  3.0, 8.5  12]; 
 Eadata=[15, 43, 88, 130,  150,  193.5,   290, 291 ]; 
 %-------------------------------------------------------------------------
 % 下面進行空載特性曲線擬合：
p=polyfit(Eadata,Ifdata,3);
%--------------------------------------------------------------------------
% 他勵發電機調整特性
I=0:0.01:50;
Ea=UN+I*Ra;
If=polyval(p,Ea);% 計算對應于擬合曲線的勵磁電流
plot(I,If,'r')
hold on
xlabel('I[A]')
ylabel('If[A]')
axis([0,90,0,60])
%--------------------------------------------------------------------------
% 并勵發電機調整特性：
Iab=0:.01:100;
Eab=UN+Iab*Ra;
Ifb=polyval(p,Eab);
Ib=Iab-Ifb;
plot(Ib,Ifb)
hold on
 %-------------------------------------------------------------------------