% 直流發電機轉矩特性分析
% 將該函數定義為dc_mo_tor(dc_motor_torque)
%--------------------------------------------------------------------------
% 下面輸入電機基本數據：
nN=1000;IfN=2.5;Rf=293/12;Ra=1.8;k=.1;E0=271
% 下面輸入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(Ifdata,Eadata,3);
If=0:.01:12;% 進行勵磁電流參數重新設置：
Ea=polyval(p,If);% 計算對應于擬合曲線的感應電勢
%--------------------------------------------------------------------------
% 計算端電壓和負載電流：
U=If*Rf;
Ia=(Ea-U)/Ra;
I=Ia-If;
% 計算并勵發電機外特性：
 plot(I,U,'k')
 hold on
 axis([0,100,0,300])
 %-------------------------------------------------------------------------
% 計算串勵發電機外特性
UU=2*(Ea-If*(Rf+Ra));
plot(6*If,UU,'b')
hold on
% 計算他勵發電機外特性
UUU=(E0-If*Ra-If.^2*k);
plot(6*If,UUU,'r')
xlabel('I[A]')
ylabel('U[V]')
%--------------------------------------------------------------------------