mathlab汽车动力匹配程序

附：计算使用的MATLAB程序代码及其含义%不同i0对应车速为40、50、60时的发动机转数clearclccloseallr=0.4275;%车轮半径u0=[405060];%速度i0=[4.224.684.885.015.285.48];%主减速比n=(i0'*u0)/(0.377*r);%对应发动机转数n%加速燃油消耗率clearclccloseallyita=0.89;G=9160*9.8;f=0.02;CDA=4.64;a=0.1704;Iw1=1.798;Iw2=3.598;r=0.4275;m=9160;If=0.218;ig=1;i0=[4.22,4.68,4.88,5.01,5.28,5.48];forx=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:1:60;%50到60公里加速ua1=50:1:60;delta=1+(Iw1+Iw2)/(m*r^2)+(If*ig^2*i0(x)^2*yita)/(m*r^2);%求旋转质量换算系数P0=(G*f.*ua0./3600+CDA.*ua0.^3/76140+(delta*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;%加速段燃油消耗量qend%减速油耗clearclccloseallua1=60;ua2=40;a=0.2584;Qi=191;Q=(8*Qi)/(367.1*8);%计算怠速燃油消耗Qd=((ua1-ua2)*Q)/(3.6*a);%减速燃油消耗Qd%六工况油耗clearclccloseallQ=[237.1422238.6344241.21566241.6732244.3986246.7512];Qs=(Q./(125+175+1000))*100%六工况燃油消耗量clearclcclosealln1=linspace(0,50000);%先求各个档位的驱动力nmax=2800;nmin=1000;r=0.4275;yita=0.89;CDA=4,64;f=0.02;G=(9160)*9.8;ig=[8.69,5.01,2.91,1.7,1.00];i0=5.01;fori=1: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*(ua(i,:).^2)/21.15;%求出滚动阻力和空气阻力的和delta(i,:)=1+(1.798+3.598+0.218*(ig(i)^2)*(i0^2)*yita)/(9160*r^2);%转动质量换算系数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,:)-F(i,:);endtemp1(1,:)=ua(2,:)/3.6;%下面分各个档位进行积分，求出加速时间temp1(2,:)=1./a(2,:);n1=1;forj1=1:1:30ifua(3,j1)max(ua(2,:))&&ua(3,j1)=90temp2(1,n1)=ua(3,j1)/3.6;temp2(2,n1)=1./a(3,j1);n1=n1+1;endendn2=1;forj1=1:1:30ifua(4,j1)max(ua(3,:))&&ua(4,j1)=90;temp3(1,n2)=ua(4,j1)/3.6;temp3(2,n2)=1./a(4,j1);n2=n2+1;endendn3=1;forj1=1:1:30ifua(5,j1)max(ua(4,:))&&ua(5,j1)=90;temp4(1,n3)=ua(5,j1)/3.6;temp4(2,n3)=1./a(5,j1);n3=n3+1;endendy=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'),gtext('ig2'),gtext('ig3');gtext('ig4');gtext('ig5');functionua=chesu(n,r,ig,i0);%由转速计算车速ua=0.377*r.*n/(ig*i0);functionn=zhuansu(ua,r,ig,i0);%求转速n=ig*i0.*ua./(0.377*r);endfunctiony=zhuanju(n);%求转矩函数y=-6.4e-5*n.^2+0.2453*n+139.69;functiony=qudongli(Ttq,ig,i0,yita,r);%求驱动力函数y=(ig*i0*yita.*Ttq)/r;end%C曲线clearclccloseallb=[18.241718.590218.799918.9809];%燃油消耗t=[36.184833.897932.015428.4538];%加速时间plot(b,t,'+r')%绘制i0散点图用+表示holdon;%保持图像b1=linspace(b(1),b(4),100);%b1到b4划100点t1=spline(b,t,b1);%三次样条差值plot(b1,t1);%绘制c曲线title('燃油经济性—加速时间曲线');%图表标题xlabel('百公里油耗(L/100km)');%x轴名称ylabel('加速时间s');%y轴名称gtext('i0=4.22'),gtext('i0=5.01'),gtext('i0=5.28');gtext('i0=5.48');%驱动力图clearclccloseallm=9160;uamax=90;%最大车速imax=0.3;%最大爬坡度Ttqmax=380;%发动机最大转矩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;f=0.02;ig1=8.69;ig5=1;ig4=1.7ig3=2.91;ig2=5.01;ua1=(0.377*r*n)/(ig1*io);%一档对应车速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,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5);title('驱动力图');xlabel('车速/(km/h)'),ylabel('驱动力(kN)');gtext('ig1'),gtext('ig2'),gtext('ig3');gtext('ig4');gtext('ig5');%动力特性图clearclccloseallTtq=[331,341,355,367,375,386,370,351,340,323];%发动机转矩io=5.01;nt=0.89;%机械效率%r=0.4275;g=9.8;m=9160;G=m*g;CDA=4.64;n=[1000,1200,1400,1600,1800,2000,2200,2400,2600,2800];%发动机转数ig1=8.69;ua=0.377*r*n/ig1/io;%1档对应车速D=(Ttq.*ig1.*io.*nt/r-CDA*ua.*ua/21.15)/G;%求1档动力因数plot(ua,D);holdon;ig2=5.01;ua=0.377*r*n/ig2/io;D=(Ttq.*ig2.*io.*nt/r-CDA*ua.*ua/21.15)/G;plot(ua,D);holdon;ig3=2.91;ua=0.377*r*n/ig3/io;D=(Ttq.*ig3.*io.*nt/r-CDA*ua.*ua/21.15)/G;plot(ua,D);holdon;ig4=1.7;ua=0.377*r*n/ig4/io;D=(Ttq.*ig4.*io.*nt/r-CDA*ua.*ua/21.15)/G;plot(ua,D);holdon;ig5=1;ua=0.377*r*n/ig5/io;D=(Ttq.*ig5.*io.*nt/r-CDA*ua.*ua/21.15)/G;plot(ua,D);title('汽车动力特性图D-ua');xlabel('ua(km/h)');ylabel('D');gtext('ig1'),gtext('ig2'),gtext('ig3');gtext('ig4');gtext('ig5');%爬坡度图clearclccloseallTtq=[331,341,355,367,375,386,370,351,340,323];io=5.01;nt=0.89;%机械效率r=0.4275;g=9.8;ma=9160;CDA=4.64;n=[1000,1200,1400,1600,1800,2000,2200,2400,2600,2800];%发动机转数G=ma.*g;%汽车重力f=0.02;ig1=8.69;%1档传动比ua=0.377*r.*n/ig1/io;%相应发动机转数对应的车速a=namelengthmax;%标识符的最长长度a=asin((Ttq.*ig1.*io.*nt/r-(G.*cos(a).*f+CDA*ua.*ua/21.15))/G);%计算1档爬角度i=100*tan(a);%计算爬坡度plot(ua,i);holdon;ig2=5.01;ua=0.3