《汽车动力性matlab仿真源程序.doc》由会员分享，可在线阅读，更多相关《汽车动力性matlab仿真源程序.doc（2页珍藏版）》请在三一文库上搜索。
1、clear allclcn=1500:500:5500;%转速范围T=78.59 83.04 85.01 86.63 87.09 85.87 84.67 82.50 80.54;%对应各转矩dt=polyfit(n,T,3);%对发动机输出转矩特性进行多项式拟合，阶数取4n1=1000:100:5500;%?t=polyval(dt,n1);figure(1)title(发动机外特性)plot(n1,t,n,T,o),grid on%图示发动机输出转矩特性 %汽车驱动力计算G=input(整车重力/N，G=);%输入970*9.8ig=3.416 1.894 1.28 0.914 0.757;
2、%变速器速比k=1:5;%5个前进档r=0.272;i0=4.388;eta=0.9;ngk=800 800 800 800 800;ngm=5500 5500 5500 5500 5500;ugk=0.377.*r.*ngk(k)./(ig(k).*i0);%计算每一档发动机800rpm时的最低行驶速度ugm=0.377.*r.*ngm(k)./(ig(k).*i0);%计算每一档发动机5400rpm最高行驶速度for k=1:5%依次计算5个档的驱动力u=ugk(k):ugm(k);n=ig(k)*i0.*u./r/0.377;t=54.8179+2.2441.*(n./100)-4.80
3、03.*(n./1000).2+2.815e-10.*n.3Ft=t.*ig(k).*i0*eta/r;figure(2)plot(u,Ft)hold on,grid on %保证K的每次循环的图形都保留显示end %行驶阻力计算f0=0.009;f1=0.002;f4=0.0003;%三者都是轿车滚动阻力系数% disp空气阻力系数Cd=0.3-0.41,迎风面积A=1.7-2.1Cd=input(空气阻力系数Cd=);%输入0.3A=input(迎风面积/m2,A=);%输入2.3u=0:10:180;f=f0+f1.*(u./100)+f4.*(u./100).4;Ff=G*f;%计算滚
4、动阻力Fw=Cd*A.*u.2./21.15;%计算空气阻力F=Ff+Fw;%滚动阻力、空气阻力之和title(驱动力-阻力图（五档速比为3.416 1.894 1.28 0.914 0.757）)plot(u,F,mo-);grid onhold on %图解法求最高车速for u=50:180;k=5; n=ig(k)*i0.*u./r/0.377; t=54.8179+2.2441.*(n./100)-4.8003.*(n./1000).2+2.815e-10.*n.3; Ft=t.*ig(k).*i0*eta/r; f=f0+f1.*(u./100)+f4.*(u./100).4; Ff=G*f; Fw=Cd*A.*u.2./21.15; F=Ff+Fw; if abs(Ft-F)1; umax=u;break endenddisp(= = =汽车动力性能仿真计算结果= = =)disp(驱动力-行驶阻力平衡图及最高车速)fprintf(汽车理论最高车速（驱动力与行驶阻力曲线交点）Vmax=%3.3f km/hn,umax)