matlab汽车动力仿真源程序代码,汽车动力性matlab仿真源程序

汽车动力性matlab仿真源程序

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=( 整车重力/N，G= );%输入 970*9.8ig=[3.416 1.894 1.28 0.914 0.757];%变速器速比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.8003.*(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.1 Cd=( 空气阻力系数 Cd= );%输入 0.3A=( 迎风面积 /m2,A= );%输入 2.3u=0:10:180;f=f0+f1.*(u./100)+f4.*(u./100).^4;Ff=G*f;%计算滚动阻力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;breakendenddisp( == == =汽车动力性能仿真计算结果== == = )disp( 驱动力 -行驶阻力平衡图及最高车速 )fprintf( 汽车理论最高车速(驱动力与行驶阻力曲线交点)Vmax=%3.3f km/h\n ,umax)