% 直流電動機機械特性分析
% 將該函數定義為dc_mo_mec(dc_motor_mech)
%--------------------------------------------------------------------------
% 下面輸入電機基本數據：
U=220;Ra=0.17;p=2;N=398;a=1;psi=0.0103;Cpsi=0.0013;
% 下面輸入電磁轉矩的變化范圍：
 Te=0:.01:5; 
 %-------------------------------------------------------------------------
% 計算并勵電動機機械特性：
Ce=p*N/60/a;
Cm=p*N/2/pi/a;
n=U/Ce/psi-Ra*Te/Ce/Cm/psi^2;
subplot(2,1,1) 
plot(Te,n,'k')
 hold on
 xlabel('Te')
 ylabel('n')
 %-------------------------------------------------------------------------
% 計算串勵電動機機械特性
C1=1/Ce*(Cm/Cpsi)^.5;
C2=1/Ce/Cpsi;
n=C1*U*(Te+.001).^(-.5)-C2*Ra;
subplot(2,1,2)
 plot(Te,n,'b')
hold on 
axis([0,5,0,60000])
 xlabel('Te')
 ylabel('n')

%-------------------------------------------------------------------------