《《汽车理论》第五版13题matlabcode》由会员分享，可在线阅读，更多相关《《汽车理论》第五版13题matlabcode（5页珍藏版）》请在金锄头文库上搜索。
1、运行结果 Untitled3Result1 = 99.0589km/h-最高车速Result2 = 8.6862km/h-最大爬坡度时的车速A = 0.3386rad-最大坡角I = 0.3522-最大爬坡度Cf1 = 1.2041-前驱附着率（1，舍弃）Cf2 = 0.4978-后驱克服该坡度的附着率time70 = 26.2088s-2档起步加速至70km/h的时间图1 汽车驱动力-行驶阻力平衡图图2 汽车行驶加速度曲线图3 汽车行驶加速度倒数曲线MATLAB-Code%1.3.1绘制汽车驱动力与行驶阻力平衡图财务盈利能力分析采用的主要评价指标包括静态评价指标和动态评价指标两类。其中静态评
2、价指标主要有投资回收期，投资利润率，投资利税率和资本金利润率；动态评价指标主要有投资回收期，净现值、净现值率，内部收益率。m=3880;r=0.367;nT=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;ig=5.56,2.769,1.644,1.00,0.793;L=3.2;a=1.947;hg=0.9;g=9.8;n=1:1:4000;%为了得到完整的图线，避免矢量长度不符，故从转速1开始Tq=-19.313+295.27*(n./1000)-165.44*(n./1000).2+40.874*(n./1000).3-
3、. 3.8445*(n/1000).4;%发动机转矩拟合公式for k=1:5%求出各档位时发动机600-4000转对应的驱动力、车速和空气阻力 for t=1:4000 Ft(k,t)=Tq(t)*ig(k)*i0*nT/r;%求各档的驱动力 ua(k,t)=0.377*n(t)*r/ig(k)/i0;%求各档的车速范围 Fw(k,t)=CdA*(ua(k,t).2)/21.15;%求不同车速对应的空气阻力 if t=Fx(k)%驱动力小于行驶阻力的数据无效 min1=Ft5(k)-Fx(k);%当前车速下驱动力与行驶阻力差值 end if minmin1 min=min1;%找到最小差值
4、result1=k,Ft5(k),Fx(k);%存储最高车速对应的转速（3200+k）、驱动力和行驶阻力 end endtext(0.377*(3200+result1(1)*r/ig(5)/i0,result1(2),M)%标出最高车速点Result1=0.377*(3200+result1(1)*r/ig(5)/i0%得出最高车速值%1.3.2.2求最大爬坡度%由图得，1档时有最大驱动力，在车速为7-10km/h之间，Ft13100N,查得转速在1850-2280之间max=0;%存储驱动力与行驶阻力最大差值，以找出最大爬坡度点for k=1:430%find the fastest sp
5、eed Ft1(k)=Tq(1850+k)*ig(1)*i0*nT/r;%1档时的驱动力 Fx(k)=CdA*(0.377*(1850+k)*r/ig(1)/i0).2)/21.15+m*g*f;%1档的行驶阻力 max1=Ft1(k)-Fx(k);%当前车速下的驱动力与行驶阻力差值 if maxmax1 max=max1;%找出最大差值 result2=k,Ft1(k),Fx(k);%存储最大爬坡度对应的转速（1850+k）、驱动力和行驶阻力 end endtext(0.377*(1850+result2(1)*r/ig(1)/i0,result2(2),I)%标出最大爬坡度点Result2
6、=0.377*(1850+result2(1)*r/ig(1)/i0%得出爬最大坡时的车速值A=asin(max/m/g)%能爬上的最大坡角I=tan(A)%最大爬坡度I%1.3.2.3求克服该坡度相应的附着率q=I+(1/cos(A)*(1/g)*(max-m*g*I*cos(A)/m);%等效坡度Cf2=q/(a/L+hg/L*q);%后轮附着率Cf1=q/(L-a)/L-hg/L*q);%前轮附着率Cf1Cf2%1.3.3.1绘制汽车行驶加速度倒数曲线for k=1:5 de(k)=1+1/m*(Iw1+Iw2)2/r+1/m*If*ig(k)2*i02*nT/r2;%汽车旋转质量换算系
7、数 for t=1:4000 acc(k,t)=1/(de(k)*m)*(Ft(k,t)-Ff-Fw(k,t);%各档的加速度 if t3310 acc(k,t)=NaN;%最高车速对应的转速为3310，之后就没有加速能力了 end endendfigure(2)%图2plot(ua(1,:),acc(1,:),ua(2,:),acc(2,:),ua(3,:),acc(3,:),ua(4,:),acc(4,:),. ua(5,:),acc(5,:);%输出加速度图形legend(1,2,3,4,5)%图例title(汽车行驶加速度曲线);xlabel(u_a/km*h-1)ylabel(a/m
8、*s-2)for k=1:5 for t=1:4000 if t3310 acc1(k,t)=NaN;%最高车速对应的转速为3310，之后就没有加速能力了 end endend figure(3)%图3plot(ua(1,:),acc1(1,:),ua(2,:),acc1(2,:),ua(3,:),acc1(3,:),. ua(4,:),acc1(4,:),ua(5,:),acc1(5,:);%输出加速度倒数曲线图legend(1,2,3,4,5)%图例title(汽车行驶加速度倒数曲线);axis(0 100 0 10);%设置图形界限，最高车速之后就没有加速能力了xlabel(u_a)yl
9、abel(1/a)%1.3.3.2积分法求2档起步加速至70公里的时间time=0;%二档起步加速至70公里用时；z=1;%循环标记for k=2:4%由图中看出2档起步，4档即可加速至70公里 for t=600:3999 if k=2%二档加速至2档最高车速然后换挡 time1=(ua(k,t+1)-ua(k,t)/3.6*(acc1(k,t+1)+acc1(k,t)/2;%时间积分微元 time(z+1)=time(z)+time1;%时间累加 z=z+1; end if k=3&ua(k-1,4000)ua(k,t)%3档加速至3档最高车速然后换挡，从3档对应的2档最高车速开始 time1=(ua(k,t+1)-ua(k,t)/3.6*(acc1(k,t+1)+acc1(k,t)/2; time(z+1)=time(z)+time1; z=z+1; end if k=4&ua(k-1,4000)ua(k,t)&ua(k,t)70%4档加速至70公里，从4档对应的3档最高车速开始 time1=(ua(k,t+1)-ua(k,t)/3.6*(acc1(k,t+1)+acc1(k,t)/2; time(z+1)=time(z)+time1; z=z+1; end endend time70=time(z-1)
