MATLAB画驱动力-行驶阻力图,两个程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
m1=2000
m2=1800;
m=3880;
r0=0.367
gt=0.85;
f=0.013;
CDA=2.77;
i0=5.83;
If=0.218;
Iw1=1.798;
Iw2=3.598;
Ig5=[5.56 2.769 1.644 1.00 0.793];
for i=1:3401;
n(i)=i+599;
Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).^2+40.874*(n(i)/1000).^3-3.8445*(n(i)/100 0).^4;
end
for i=1:3401;
for j=1:5;
Ft(i,j)=Ttq(i)*i0*Ig5(j)*gt/r0;
ua(i,j)=0.377*r0*n(i)./(Ig5(j)*i0);
F(i,j)=f*m*9.8+CDA*ua(i,j).^2/21.15;
end
end
plot(ua,Ft,ua,F)
xlabel('ua/(km/h)');
ylabel('F/N');
title('汽车驱动力-行驶阻力平衡图');
gtext('Ft1')
gtext('Ft2');
gtext('Ft3');
gtext('Ft4');
gtext('Ft5');
gtext('Ff+Fw');
>> m=3880;
>> g=9.8;
>> G=m*g;
>> ig=[5.56 2.769 1.644 1.00 0.793];
>> nT=0.85;
>> r=0.367;
>> f=0.013;
>> CDA=2.77;
>> i0=5.83;
>> L=3.2;
>> a=1.947;
>> hg=0.9;
>> If=0.218;
Iw1=1.798;
Iw2=3.598;
nmin=600;nmax=4000;
>> n=600:10:4000;
>> Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; >> Ft1=Tq*ig(1)*i0*nT/r;
>> Ft2=Tq*ig(2)*i0*nT/r;
>> Ft3=Tq*ig(3)*i0*nT/r;
>> Ft4=Tq*ig(4)*i0*nT/r;
>> Ft5=Tq*ig(5)*i0*nT/r;
>> ua1=0.377*r*n/ig(1)/i0;
>> ua2=0.377*r*n/ig(2)/i0;
>> ua3=0.377*r*n/ig(3)/i0;
>> ua4=0.377*r*n/ig(4)/i0;
>> ua5=0.377*r*n/ig(5)/i0;
>> ua=0:5:120;
>> Ff=G*f;
>> Fw=CDA*ua.^2/21.15;
>> Fz=Ff+Fw;
>> plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);
>> title('汽车驱动力-行驶阻力平衡图');
>> xlabel('ua(km*h^-1)');
>> ylabel('Ft/N');
>> gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');