% 直流發(fā)電機(jī)外特性分析
% 將該函數(shù)定義為dc_ge_load(dc_generator_load)
%--------------------------------------------------------------------------
% 下面輸入電機(jī)基本數(shù)據(jù)：
nN=1000;IfN=2.5;Rf=293/12;Ra=1.8;k=.1;E0=271
% 下面輸入750轉(zhuǎn)／分鐘時(shí)的空載特性試驗(yàn)數(shù)據(jù)(Ifdata是勵(lì)磁電流,Eadata是感應(yīng)電勢(shì))：
 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 ]; 
 %-------------------------------------------------------------------------
 % 下面進(jìn)行空載特性曲線擬合：
p=polyfit(Ifdata,Eadata,3);
If=0:.01:12;% 進(jìn)行勵(lì)磁電流參數(shù)重新設(shè)置：
Ea=polyval(p,If);% 計(jì)算對(duì)應(yīng)于擬合曲線的感應(yīng)電勢(shì)
%--------------------------------------------------------------------------
% 計(jì)算端電壓和負(fù)載電流：
U=If*Rf;
Ia=(Ea-U)/Ra;
I=Ia-If;
% 計(jì)算并勵(lì)發(fā)電機(jī)外特性：
 plot(I,U,'k')
 hold on
 axis([0,100,0,300])
 %-------------------------------------------------------------------------
% 計(jì)算串勵(lì)發(fā)電機(jī)外特性
UU=2*(Ea-If*(Rf+Ra));
plot(6*If,UU,'b')
hold on
% 計(jì)算他勵(lì)發(fā)電機(jī)外特性
UUU=(E0-If*Ra-If.^2*k);
plot(6*If,UUU,'r')
xlabel('I[A]')
ylabel('U[V]')
%--------------------------------------------------------------------------