实验四:直流电动机的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])