% 直流電動機轉矩特性分析
% 將該函數定義為dc_mo_tor(dc_motor_torque)
%--------------------------------------------------------------------------
% 下面輸入電機基本數據：
Cm=10;Ra=1.8;k=.1;k1=.2;
% 下面輸入750轉／分鐘時的空載特性試驗數據(Ifdata是勵磁電流,Eadata是感應電勢)：
 Ia=0:.01:15; 
 %-------------------------------------------------------------------------
% 計算并勵電動機外特性：
Temb=Cm*k1*Ia;
 plot(Ia,Temb,'k')
 hold on
 axis([0,20,0,60])
 %-------------------------------------------------------------------------
% 計算串勵電動機外特性
Temc=Cm*k*Ia.^2;
plot(Ia,Temc,'b')
hold on
 %-------------------------------------------------------------------------
% 計算他勵電動機外特性
Temt=Cm*k*Ia;
plot(Ia,Temt,'r')
xlabel('Ia[安培]')
ylabel('Tem[牛頓.米]')
%--------------------------------------------------------------------------