汽车动力性计算matlab程序.docx

%%汽车动力性计算（自己编的动力性计算程序供大家计算动力性时参考具体参数大家根据所给程序对应输入并对坐标轴数值按需要进行修改）clc;clear;closeall;%%根据所给发动机数据拟合外特性曲线（发动机数据按照你所得到的数据进行输入）n_test=[5006007008009001000110012001300140015001600170018001900200021002200];T_test=[97511081298149615461620167017851974197419701889182917481669170015241105];figure(1)plot(n_testT_test'g');holdongridon%p=polyfit(n_testT_test7);p=polyfit(n_testT_test2);n=[450:1:2200];Ttq=polyval(pn);plot(nTtq'k');xlabel('发动机转速n(r/min)');ylabel('发动机转矩Ttq(N*m)');title('发动机转矩曲线');legend('测试曲线''拟合曲线');%%所给车型动力总成相关参数ig=[3.072.161.481.00.82];i0=4.0;eta=0.78;r=0.57;M=25000;g=9.8;c=1.5;f0=0.01;f1=0.0002;f4=0.0005;CD=1;A=8;Iw=3.6;If=0.04;%%发动机外特性曲线图figure(2)holdongridonfori=length(n);Pe=Ttq.*n/9550;end[AXH1H2]=plotyy(nTtqnPe);xlabel('发动机转速n(r/min)');ylabel('发动机转矩Ttq(N*m)');ylabel(AX(2)'发动机功率Pe(Kw)');title('发动机外特性曲线');%%各挡位速度曲线%计算各挡位车速fori=1:length(ig);ua(i:)=0.377*r*n/ig(i)/i0;end%计算各档位最高车速uamax=ua(:length(ua(1:)));figure(3)holdonfori=1:length(ig);plot(nua(i:)'k');endholdongridonxlabel('转速n(r/min)');ylabel('各挡位车速(km/h)');title('各挡位车速-转速表');legend('1挡车速''2挡车速''3挡车速''4挡车速''5挡车速');%%驱动力和行驶阻力平衡图%计算滚动阻力系数fori=1:length(ig);f(i:)=f0+f1*(ua(i:)/100)+f4*(ua(i:)/100).^4;end%计算滚动阻力fori=1:length(ig);Ff(i:)=c*M*g*f(i:);end%计算空气阻力fori=1:length(ig);Fw(i:)=CD*A*(ua(i:).^2)/21.15;end%计算行驶阻力fori=1:length(ig);F(i:)=Ff(i:)+Fw(i:);end%计算汽车驱动力fori=1:length(ig);Ft(i:)=Ttq*ig(i)*i0*eta/r;endfigure(4)holdonfori=1:length(ig);plot(ua(i:)Ft(i:)'k');plot(ua(i:)F(i:)'r');plot(ua(i:)Ff(i:)'b');endholdongridonxlabel('车速(km/h)');ylabel('驱动力、行驶阻力(N)');legend('驱动力Ft''行驶阻力Ff+Fw''滚动阻力Ff');title('驱动力-行驶阻力平衡图');%%汽车功率平衡图%计算各档位功率fori=1:length(ig);P(i:)=Ft(i:).*ua(i:)/(3600*eta);end%计算风阻阻力功率fori=1:length(ig);Pw(i:)=CD*A*ua(i:).^3/(76140*eta);end%计算滚动阻力功率fori=1:length(ig);Pf(i:)=M*g*f(i:).*ua(i:)/(3600*eta);end%计算总阻力功率fori=1:length(ig)