mathlab汽车动力匹配程序Word格式.docx

1、i0=4.22,4.68,4.88,5.01,5.28,5.48;for x=1:1:6n=1000:100:2800; %等距取数ua5=（0.377*r.*n）/i0（x）; %求转数对应车速F5=f*G+CDA*（ua5.2）/21.15; %求阻力P_fw5=F5.*ua5./（yita*3.6*1000）; %阻力功率ua0=50:60; %50到60公里加速ua1=50:delta=1+（Iw1+Iw2）/（m*r2）+（If*ig2*i0（x）2*yita）/（m*r2）; %求旋转质量换算系数P0=（G*f.*ua0./3600+CDA.*ua0.3/76140+（delta*

2、m.*ua0/3600）*a）/yita; %初速度为ua0时的阻力功率P=（G*f.*ua1/3600+CDA.*ua1.3/76140+（delta*m.*ua1/3600）*a）/yita; %末速度为ua1时的阻力功率b1=0.0133.*P.2-2.0715.*P+287.47; %燃油消耗率Qt=P.*b1./（367.1.*8）; %燃油消耗i1=size（Qt）; %求Qt向量长度 i=i1（2）;Qt1=Qt（2:i-1）;dt=1/（3.6*a）; %速度每增加1km/h所需要的时间q=（Qt（1）+Qt（i）*dt./2+sum（Qt1）*dt; %加速段燃油消耗量qen

3、d%减速油耗ua1=60;ua2=40;a=0.2584;Qi=191;Q=（8*Qi）/（367.1*8）; %计算怠速燃油消耗Qd=（ua1-ua2）*Q）/（3.6*a）; %减速燃油消耗Qd%六工况油耗Q=237.1422 238.6344 241.21566 241.6732 244.3986 246.7512;Qs=（Q./（125+175+1000）*100 %六工况燃油消耗量n1=linspace（0,50000）; %先求各个档位的驱动力nmax=2800;nmin=1000;CDA=4,64;G=（9160）*9.8;ig=8.69,5.01,2.91,1.7,1.00;i

4、0=5.01;for i=1:5 %i为档数 uamax（i）=chesu（nmax,r,ig（i）,i0）; %计算各个档位的最大速度与最小速度 uamin（i）=chesu（nmin,r,ig（i）,i0）; ua（i,:）=linspace（uamin（i）,uamax（i）,30）; n（i,:）=zhuansu（ua（i,:）,r,ig（i）,i0）; %计算各个档位的转速范围 Ttq（i,:）=zhuanju（n（i,:）; %求出各档位的转矩范围 Ft（i,:）=qudongli（Ttq（i,:）,ig（i）,i0,yita,r）; %求出驱动力 F（i,:）=f*G+CDA*（

5、ua（i,:）.2）/21.15; %求出滚动阻力和空气阻力的和 delta（i,:）=1+（1.798+3.598+0.218*（ig（i）2）*（i02）*yita）/（9160*r2）; %转动质量换算系数 a（i,:）=1./（delta（i,:）.*14000./（Ft（i,:）-F（i,:）; %求出加速度 c=1./a（i,:）; plot（ua,c）; title（加速度倒数曲线 xlabel（ua %x轴名称 ylabel（1/a %y轴名称 F2（i,:）=Ft（i,: temp1（1,:）=ua（2,:）/3.6; %下面分各个档位进行积分，求出加速时间 temp1（2,

6、:）=1./a（2,: n1=1; for j1=1:30 if ua（3,j1）max（ua（2,:）&ua（3,j1）max（ua（3,:ua（4,j1）max（ua（4,:ua（5,j1） temp4（1,n3）=ua（5,j1）/3.6; temp4（2,n3）=1./a（5,j1）; n3=n3+1;y=temp1（1,1）*temp1（2,1）+qiuji（temp1（1,:）,temp1（2,:）+qiuji（temp2（1,:）,temp2（2,:）+qiuji（temp3（1,:）,temp3（2,:）+qiuji（temp4（1,:）,temp4（2,:ygtext（ig1

7、）,gtext（ig2ig3ig4ig5function ua=chesu（n,r,ig,i0）; %由转速计算车速ua=0.377*r.*n/（ig*i0）; function n=zhuansu（ua,r,ig,i0）; %求转速n=ig*i0.*ua./（0.377*r）;function y=zhuanju（n）; %求转矩函数y=-6.4e-5*n.2+0.2453*n+139.69;function y=qudongli（Ttq,ig,i0,yita,r）; %求驱动力函数y=（ig*i0*yita.*Ttq）/r;%C曲线b=18.2417 18.5902 18.7999 18.

8、9809;t=36.1848 33.8979 32.0154 28.4538; %加速时间plot（b,t,+r） %绘制i0散点图用+表示hold on; %保持图像b1=linspace（b（1）,b（4）,100）; %b1到b4划100点t1=spline（b,t,b1）; %三次样条差值plot（b1,t1）; %绘制c曲线title（燃油经济性加速时间曲线 %图表标题xlabel（百公里油耗（L/100km）ylabel（加速时间si0=4.22i0=5.01i0=5.28i0=5.48%驱动力图uamax=90; %最大车速imax=0.3; %最大爬坡度Ttqmax=380;

9、%发动机最大转矩nPemax=2800; %最大功率转数d=0.855; %轮胎直径r=d/2; %轮胎半径Ttq=331,341,355,367,375,386,370,351,340,323; %发动机转矩n=1000,1200,1400,1600,1800,2000,2200,2400,2600,2800; %发动机转数io=（0.377*r*2800）/uamax; %最小传动比afamax=atan（imax）; %最大爬坡度角度aitT=0.89;g=9.8;ig1=8.69;ig5=1;ig4=1.7ig3=2.91;ig2=5.01;ua1=（0.377*r*n）/（ig1*i

10、o）; %一档对应车速ua2=（0.377*r*n）/（ig2*io）;ua3=（0.377*r*n）/（ig3*io）;ua4=（0.377*r*n）/（ig4*io）;ua5=（0.377*r*n）/（ig5*io）;Ft1=（Ttq*ig1*io*aitT）/（1000*r）; %一档对应驱动力Ft2=（Ttq*ig2*io*aitT）/（1000*r）;Ft3=（Ttq*ig3*io*aitT）/（1000*r）;Ft4=（Ttq*ig4*io*aitT）/（1000*r）;Ft5=（Ttq*ig5*io*aitT）/（1000*r）;figure,plot（ua1,Ft1,ua2,F

11、t2,ua3,Ft3,ua4,Ft4,ua5,Ft5）;驱动力图车速/（km/h）,ylabel（驱动力（kN） ;%动力特性图io=5.01;nt=0.89; %机械效率% G=m*g; %发动机转数 ua=0.377*r*n/ig1/io; %1档对应车速D=（Ttq.*ig1.*io.*nt/r-CDA*ua.*ua/21.15）/G; %求1档动力因数plot（ua,D）;ua=0.377*r*n/ig2/io;D=（Ttq.*ig2.*io.*nt/r-CDA*ua.*ua/21.15）/G;ua=0.377*r*n/ig3/io;D=（Ttq.*ig3.*io.*nt/r-CDA*ua.*ua/21.15）/G;ig4=1.7;ua=0.377*r*n/ig4/io;D=（Ttq.*ig4.*io.*nt/r-CDA*ua.*ua/21.15）/G;ho