直流电动机matlab仿真实验,实验四：直流电动机的MATLAB仿真实验

实验四：直流电动机的MATLAB仿真实验结果

clear

UN=220;PN=5.5;IN=30.4;

nN=1000;Ra=0.55;Rf=430;

IaN=IN-UN/Rf;

CePhi=(UN-Ra*IN)/nN;

CtPhi=9.55*CePhi;

Ia=0:IN;

n=UN/CePhi-Ra/CePhi*Ia;

Te=CtPhi*Ia;

T2N=9500*PN/nN;

figure(1);

plot(Te,n,'*');

xlabel('电磁转矩');

ylabel('转速');

ylim([0,1200]);

figure(2);

plot(Te,n,'rs');

xlabel('电磁转矩');

ylabel('转速');

hold on;

for coef=1:-0.25:0.25;

U=UN*coef;

n=U/CePhi-Ra/(CePhi-Ra)/(CePhi*CtPhi)*T e;

plot(Te,n,'k-');

str=strcat('U=',num2str(U),'V');

sy=1000*coef;

text(50,sy,str);

end

ylim([0,1200]);

figure(3);

Rc=0;

n=UN/CePhi-(Ra+Rc)/(CePhi*CtPhi)*Te; plot(Te,n,'rs');

xlabel('电磁转矩');

ylabel('转速');

hold on;

U=UN;

for Rc=0:2.0:8.0;

n=U/CePhi-(Ra+Rc)/(CePhi*CtPhi)*Te; plot(Te,n,'k-');

str=strcat('R=',num2str(Ra+Rc),'\Omega'); sy=260*(4-Rc*0.5);

text(50,sy,str);

end

ylim([0,1200]);

figure(4);

Rc=0;

n=UN/CePhi-(Ra+Rc)/(CePhi*CtPhi)*Te; plot(Te,n,'rs');

xlabel('电磁转矩');

ylabel('转速');

hold on;

U=UN;

for coef=0.5:0.25:1.3;

Ce=CePhi*coef;

Ct=CtPhi*coef;

n=U/Ce-Ra/(Ce*Ct)*Te;

plot(Te,n,'k-');

str=strcat('\phi=',num2str(coef),'\phiN');

sy=600*(4-coef*2.1);

text(50,sy,str);

end

ylim([0,2500])