当前位置：文档之家 › MATLAB汽车仿真技术复习答案

MATLAB汽车仿真技术复习答案

  • 格式：doc
  • 大小：270.00 KB
  • 文档页数：9

汽车仿真技术复习题 一、术语解释：

1．计算机仿真:是利用计算机对一个实际系统的结构和行为进行动态演示,以评价或预测该系统的行为效果。是解决复杂问题的一条有效途径。 2. 矩阵的秩：矩阵线性无关的行数与列数称为矩阵的秩。 3．传递函数: 在线性定常系统中，当初始条件为零时，系统输出量的拉普拉斯变换与输入量的拉普拉斯变换之比。 4．伯德图:幅频特性曲线图和相频特性曲线图两张图合起来称为伯德图。 5．根轨迹法: 用作图的方法表示特征方程的根与系统某一参数的全部数值关系，当这一参数取特定值时，对应的特征根可在上述关系图中找到。这种方法叫根轨迹法 二、填空题

1．58.135.0ea=（ 4.29 ）； 程序为： >> format bank; >> a= exp(-0.5)+sqrt(13.58) a = 4.29

2．建立M文件求取表达式?2009998321的值。 sum=0; for i=1:1:200 sum=sum+i; end 在命令平台上键入sum后回车（键入Enter键），记录计算结果。

3. 方程的求解：0312512.5657.1069.3212345xxxxx， >> d=[1 -21 3.69 -10.57 56.12 -3125]; >> roots(d) ans = 20.8577 2.4747 + 2.6023i 2.4747 - 2.6023i -2.4035 + 2.4168i -2.4035 - 2.4168i

4．一圆柱形金属构件，直径为20cm,长为110cm,质量为362kg,求此金属的密度（3g/cm）。试补充划线部分的程序。 >> d=20; >> L=110; >> v=pi*(d/2)*(d/2)*L; >> m=362000; >> rou=m/v rou = 10.4753

5．复数绘图，已知复数g的表达式)6)(5)(2)(1(60jjjjg，的范围200,200，增量取0.01。试补充划线部分的程序。

w=-200:0.01:200;

g=60./((1+j*w).* (2+j*w) .* (5+j*w) .* (6+j*w));

plot(g); grid;

6．数值计算10195.212)3sin(dyyxdx）（ >> dblquad(inline('sin(3*x)+y.^2'),1,10,-1,2.95) ans = 78.51053721163348 7．求1)sin(5.0xeyx在[2,8]内的最小值。

>>f=inline('exp(-0.5*x).*sin(x)+1'); >> [xmin,fval]=fminbnd(f,2,8) xmin = 4.24872619592456 fval = 0.89310865860724

8．利用函数绘图命令绘制)8/2sin(5.0xeyx，其中]5,0[x； >> fplot( 'exp(-0.5*x).*sin(2*x-pi/8)',[0,5,-1.1,1.1] ) 9.请补充完整以下程序： x=0; n=1; while n<6 x= n+x ; n= n+1 ; end x 运行结果 x= 15 .

三、简答题

1.绘制三维曲面图)2.0sin())cos(sin(xyxz，其中 [x,y]=meshgrid(0:0.2:13)。试写出程序。 （已给出结果图，如图1） [x,y]=meshgrid(0:0.2:13); z=sin(x+cos(y))-sin(0.2*x); mesh(x,y,z); title('三维曲面图'); 2．写出以下系统的多项式模型，并将其转换为零极点模型 2153173261552115.35291)(23452341ssssssssssG

>> n1=[91,-52,3.5,-11,52]; d1=[1,15,26,73,31,215]; sys1=tf(n1,d1) [z1,p1,k1]=tf2zp(n1,d1) sys1zp=zpk(z1,p1,k1)

3．.已知单位负反馈系统的开环传递函数为 10.5s10.2s1s1.0ksG）（

要求编程绘制50k时的极坐标图，确定曲线与负实轴的交点坐标及频率值。 （已给出结果图2） n=50; d=conv([0.1,1],conv([0.2,1],[0.5,1])); sys=tf(n,d); nyquist(sys) 曲线与负实轴的交点坐标为-3.76； 曲线与负实轴的交点频率值9.2； 4．单位反馈系统的开环传递函数为 )3)(2()1()(ssssKsGg

试编程绘制闭环系统的根轨迹（根轨迹图已给出，如图3）。并回答 （1）闭环系统稳定的gK的取值范围；

（2）系统的阶跃响应有超调的gK的取值范围； （3）分离点的坐标。 编程: %ggj01.m z=[-1]; p=[0;-2;-3]; k=[1]; sys=zpk(z,p,k); rlocus(sys) (1)系统稳定的gK的取值范围：0gK；

(2)系统的阶跃响应有超调的gK的取值范围：419.0gK,(3)分离点的坐标：47.2d 5.有二个系统如图（a）、（b）所示，编程分别绘制各系统的单位阶跃响应和单位脉冲响应。

5352ss )5)(3()2(20sss

s 1211.0ss )(sR

)(sC

（a） 1523512223ssss

)10)(3()6(2.0ssss

)173)(12()15.1)(16.0(52sssss )(sR )(sC

(b) )(sC )(sR ss2102 15.01s

11.01s ss632 )(2sG

（c） 图5-17 6．增量与绘图，通过仿真结果回答问题 （1）程序1 t1=0:4*pi; y1=sin(2*t1); plot(t1,y1); 程序1可得到图形正弦曲线A（增量为默认值1）。 （2）程序2 t2=0:0.2:4*pi; y2=sin(2*t2); plot(t2,y2); 程序2可得到图形正弦曲线B（增量为0.2）。 （3）程序3 t3=0:0.01:4*pi; y3=sin(2*t3); plot(t3,y3); 程序3可得到图形正弦曲线C（增量为0.01）。 （4）程序4 t4=0:0.00001:4*pi; y4=sin(2*t4); plot(t4,y4); 程序4可得到图形正弦曲线D（增量为0.00001）：

回答问题： ① 通过MATLAB仿真得到的正弦曲线A、B、C有何区别？ 答：曲线A失真严重，曲线B略有失真，曲线C几乎没有失真。 ② 通过MATLAB仿真得到的正弦曲线C、D有何区别？ 答：曲线C和曲线D区别不大。 ③ 你得到的结论是通常情况下增量（步长）取多少为宜？ 答：通常情况下增量（步长）取0.01即可满足一般的仿真要求和仿真精度。 7. 绘制)3cos(tr，3,0t的极坐标（polar）图，试写出程序。 t=0:0.01:3*pi; r=cos(3*t); polar(t,r)

8．采用图形保持，在同一坐标内绘制曲线)sin(1xy ,)sin(25.0xeyx,

)cos(3xy,]10,0[x。

试写出程序。（已给出结果图2） x=0:0.01:10; y1=sin(x); y2=exp(-0.5*x).*sin(x); y3=cos(x); plot(x,y1, 'k'); hold on; plot(x,y2, 'r'); plot(x,y3, 'b'); hold off;

9. 试编写计算程序画出曲线与曲面图形 .22,22,)1((2).);20(cossin(1).223322yxyxexztty

tx

yx

%1 %2 t=0:0.1: 2*pi; [x,y]=meshgrid(-2:0.1:2); x= sin(t).^3; z=(x.^2+1).* exp(-x.^2- y.^2+)+x.^2.*y; y= cos (t).^3; mesh(x,y,z) plot(x,y)

合集下载

matlab与控制系统仿真期末复习

matlab与控制系统仿真期末复习

matlab与控制系统仿真期末复习考试题型一、填空题（20分）矩阵建立A=[1 2 3;4 5 6;7 8 9]A([1 3],:) 表示取第一三行A(1:2,:)表示取第一二行ay=a*sin(a)例如:Y=A.*sin(Y)表示点乘点乘和乘是两种不同的矩阵计算符号。

（1）点乘表示两个矩阵对应位置元素相乘，所以这两个矩阵应该是尺寸等大的（这里不是说元素等大，而是行列数分别相等，都是m 行n列的矩阵）。

例如：>> a = [1 2 3;2 3 4];>> b = [1 2 3;2 3 4];>> a.*b1 4 94 9 16>> c = [1 2 3;1 2 3;1 2 3];>> a.*cError using ==> timesMatrix dimensions must agree.这里a无法点乘c，因为这两个不是相同尺寸的矩阵。

（2）而乘* ，这个和我们在数学中学到的矩阵相乘是一个意义。

例如：>> a*cans =>> a*bError using ==> mtimesInner matrix dimensions must agree.由于a的列不等于b的行，所以两者不可相乘。

Shift+enter:一行多命令输入…：连接符可以把断开的命令连起来，但是…后的不在被识别，是字符要单独用‘隔开;：表示不显示结果A=5;A=5,B=6\矩阵左除/矩阵右除x=a/b即x=a*inv（b）a乘以b的逆矩阵。

x=a\b,即x=inv（a)*b a的逆矩阵乘以b矩阵。

AB是标量A/B等于B\AAB是矩阵A/B不等于A\B标量：单个数据a=2；向量：一行或一列数据a=[3,5,7,8]或a=[3 5 7 8] 行向量，在方括号内的一行数据%：注释close all 关闭图形窗口clear清除工作空间clc清屏who：who命令只列出当前工作空间的变量名whos：以列表的形式，给当前工作空间的变量名以及变量的其它信息(大小，字节，类等）figure（n）：figure是建立图形的意思，系统自动从1，2，3，4...来建立图形，数字代表第几幅图形，figure(1)，figure(2)就是第一第二副图的意思save acb.mat A B –append将A B变量保存成文件名为acb.mat的文件。

汽车理论课后作业答案MATLAB

汽车理论课后作业答案MATLAB

汽车理论作业MA TLAB过程汽车驱动力与阻力平衡图加速度倒数-速度曲线图0 10 20 30 40 50 60 70 80 90 100u汽车功率平衡图u/(km/h)最高档等速百公里油耗曲线Ua/(km/h)燃油积极性-加速时间曲线源程序：《第一章》m=3880;g=9.8;r=0.367;x=0.85;f=0.013;io=5.83;CdA=2.77;lf=0.218;Iw1=1.798;Iw2=3.598;Iw=lw1+lw2;ig=[6.09 3.09 1.71 1.00]; %变速器传动比L=3.2;a=1.947;hg=0.9;n=600:1:4000;T=-19.313+295.27* n/1000-165.44*(门/1000)人2+40.874*(门/1000)人3-3.8445*( n/IOOO).%;Ft1=T*ig(1)*io*x/r; %计算各档对应转速下的驱动力Ft2=T*ig(2)*io*x/r;Ft3=T*ig(3)*io*x/r;Ft4=T*ig(4)*io*x/r;u1=0.377*r*n/(io*ig(1));u2=0.377*r*n/(io*ig(2));u3=0.377*r*n/(io*ig(3));u4=0.377*r*n/(io*ig(4));u=0:130/3400:130;F仁m*g*f+CdA*u”2/21.15;%计算各档对应转速下的驱动阻力F2=m*g*f+CdA*u2.A2/21.15;F3=m*g*f+CdA*u3.A2/21.15;F4=m*g*f+CdA*u4.A2/21.15;figure(1);plot(u1,Ft1, '-r' ,u2,Ft2, '-m' ,u3,Ft3, '-k' ,u4,Ft4, '-b' ,u1,F1, '-r' ,u2,F2, '-m' ,u3,F3, ' k' ,u4,F4, '-b' , 'LineWidth' ,2) title( ' 汽车驱动力与阻力平衡图' );xlabel( 'u_{a}/km.hA{-1}' )ylabel( 'F/N' )gtext( 'F_{t1}' )gtext( 'F_{t2}' )gtext( 'F_{t3}' )gtext( 'F_{t4}' )gtext( 'F_{f}+F_{w}' )%由汽车驱动力与阻力平衡图知，他们无交点，u4在最大转速时达到最大umax=u4(3401)Ft1max=max(Ft1);imax=(Ft1max-m*g*f)/(m*g)disp( ' 假设是后轮驱动' );C=imax/(a/L+hg*imax/L) % 附着率delta1=1+(Iw1+Iw2)/(m*rA2)+If*ig(1)*rA2*ioA2*x/(m*rA2);delta2=1+(Iw1+Iw2)/(m*rA2)+If*ig(2)*rA2*ioA2*x/(m*rA2);delta3=1+(Iw1+Iw2)/(m*rA2)+If*ig(3)*rA2*ioA2*x/(m*rA2);delta4=1+(Iw1+Iw2)/(m*rA2)+If*ig(4)*rA2*ioA2*x/(m*rA2);a1=(Ft1-F1)/(delta1*m); %加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);h1=1./a1; %加速度倒数h2=1./a2;h3=1./a3;h4=1./a4;figure(2);plot(u1,h1,u2,h2,u3,h3,u4,h4, 'LineWidth',2);title( ' 加速度倒数- 速度曲线图xlabel( 'u' )ylabel( '1/a' )gtext( '1/a1' )gtext( '1/a2' )gtext( '1/a3' )gtext( '1/a4' )%由加速度倒数- 速度曲线图可知u1min=min(u1);u1max=max(u1); u2min=u1max; u2min=min(u2);u2max=max(u2); u3min=u2max;u3max=max(u3); u4min=u3max; u4max=70;x1=[]; x2=[]; x3=[]; x4=[];y=3401; for i=1:3401;if u3(i)<=u3min;x1=[i];end end q1=max(x1); ua3=u3(q1:y); a3=h3(q1:y); for i=1:3401;if u4(i)<=u4min;x2=[i];elseif u4(i)<=u4max;x3=[i];end end q2=max(x2); q3=max(x3); ua4=u4(q2:q3); a4=h4(q2:q3); s1=trapz(h2,u2 ); %二挡运行时间s2=trapz(ua3,a3);s3=trapz(ua4,a4);s=[s1 s2 s3];disp( ' 积分得' ) t=sum(s)*1000/3600 %总时间《第二章》Pe1=Ft1.*u1./3600; %计算各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;P1=F1.*u1./(3600*x); %计算各档对应的各个车速下的行驶功率P2=F2.*u2./ (3600*x);P3=F3.*u3./ (3600*x);P4=F4.*u4./ (3600*x);figure(3);plot(u1,Pe1, '-r' ,u2,Pe2, '-m' ,u3,Pe3, '-k' ,u4,Pe4, '-b' ,u1,P1, 'k' ,u2,P2, 'k' ,u3,P3, 'k' u4,P4, 'k' ,'linewidth' ,2);gtext( 'Pe1' )gtext( 'Pe2' )gtext( 'Pe3' )gtext( 'Pe4' )xlabel( 'u/(km/h)' );ylabel( 'P/kW' );title( ' 汽车功率平衡图' );n=[815 1207 1614 2012 2603 3006 3403 3804];Ua=[];Ua=0.377*r*n./(io*ig(4))ft=[];ft=m*g*f+(2.77/21.15)*Ua.A2; 算各转速对应的各个车速下的行驶阻力Pe(1)=ft(1).*Ua(1)./(3600*x); %计算各转速对应的各个车速下的行驶阻力功率Pe(2)=ft(2).*Ua(2)./(3600*x);Pe(3)=ft(3).*Ua(3)./(3600*x);Pe(4)=ft(4).*Ua(4)./(3600*x);Pe(5)=ft(5).*Ua(5)./(3600*x);Pe(6)=ft(6).*Ua(6)./(3600*x);Pe(7)=ft(7).*Ua(7)./(3600*x);Pe(8)=ft(8).*Ua(8)./(3600*x)B0=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7];B1=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291];B2=[72.739 36.657 14.525 7.0035 4.4763 2.8593 2.9788 0.7113];B3=[-5.8629 -2.0533 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215]; B4=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568];b1=(B0⑴)+(B1(1)*Pe(1))+(B2(1)*Pe(1f2)+(B3(1)*Pe⑴T)+(B4(1)*Pe⑴人4); b2=(B0(2))+(B1(2)*Pe (2))+(B2(2)*Pe(2)A2)+(B3(2)*Pe(2)A3)+(B4(2)*Pe(2)A4);b3=(B 0(3) )+(B1(3)*Pe(3))+(B2 (3)*Pe(3)A2)+(B3(3)*Pe(3)A3)+(B4(3)*Pe(3)A4);b4=(B0(4))+(B1(4)*Pe(4))+(B2(4)*Pe(4)A2)+(B3(4)*Pe(4)A3)+(B4(4)*Pe(4)A4);b5=(B0(5))+(B1(5)*Pe(5))+(B2(5)*Pe(5)A2)+(B3(5)*Pe(5)A3)+(B4(5)*Pe(5)A4);b6=(B0(6))+(B1(6)*Pe(6))+(B2(6)*Pe(6)A2)+(B3(6)*Pe(6)A3)+(B4(6)*Pe(6)A4);b7=(B0(7))+(B1(7)*Pe(7))+(B2(7)*Pe(7)A2)+(B3(7)*Pe(7)A3)+(B4(7)*Pe(7)A4);b8=(B0(8))+(B1(8)*Pe(8))+(B2(8)*Pe(8)A2)+(B3(8)*Pe(8)A3)+(B4(8)*Pe(8)A4);p=0.7;Qs=[];Qs(1)=(Pe(1)*b1)/(1.02*Ua(1).*p*g);Qs(2)=(Pe(2)*b2)/(1.02*Ua(2).*p*g);Qs(3)=(Pe(3)*b3)/(1.02*Ua(3).*p*g);Qs(4)=(Pe(4)*b4)/(1.02*Ua(4).*p*g);Qs(5)=(Pe(5)*b5)/(1.02*Ua(5).*p*g);Qs(6)=(Pe(6)*b6)/(1.02*Ua(6).*p*g);Qs(7)=(Pe(7)*b7)/(1.02*Ua(7).*p*g);Qs(8)=(Pe(8)*b8)/(1.02*Ua(8).*p*g);M=polyfit(Ua,Qs,2);UA=0.377*r*600/(ig(4)*io):1:0.377*r*4000/(ig 4)*io); %卩表示车速QS=polyval(M,UA); %0表示油耗figure(4);plot(UA,QS, 'linewidth' ,2);title( '最高档等速百公里油耗曲线' );xlabel( 'Ua/(km/h)' );ylabel( 'Qs/L' );《第三章》io=[5.17 5.43 5.83 6.17 6.33];Va=0.377*r*n(7)./(io.*ig(4));Ps=46.9366;disp('假设以最高档，较高转速(n取3401 )，最经济负荷(即90%负荷大约46.9366KW )行驶时油耗' ) qs=[];qs(1)=(Ps*b7)/(1.02*Va(1).*p*g); qs(2)=(Ps*b7)/(1.02*Va(2).*p*g);qs(3)=(Ps*b7)/(1.02*Va(3).*p*g);qs(4)=(Ps*b7)/(1.02*Va(4).*p*g); qs(5)=(Ps*b7)/(1.02*Va(5).*p*g);st=[ 17.5813 16.2121 14.5126 13.3775 12.9185]; %加速时间：(这里以最高档〈四档〉、速度由0加速到94.93Km/h 的时间) 因与题 1.3 第三问求法相同，这里不在累述，可直接有计算机求得：figure(5);plot(qs,st, '+' ,'lineWidth' ,2)hold on plot(qs,st,); gtext( 517') gtext( 543') gtext( 583') gtext( 617') gtext( 633')title('燃油积极性-加速时间曲线’); xlabel('燃油经济性(qs/L)'); ylabel('动力性--原地起步加速时间(st/s)');《第四章》4-31)时：前轮先抱死 0时：后轮先抱死利用MATLAB 作图得到下图:前轴利用附着系数为:L z b zh g后轴利用附着系数为:L 1 zazh g空载时： 03.95 0.38 1.85所以满载时：ohg0.8450.413空载时后轮总是先抱死。

基于matlab的车辆工程仿真实例 -回复

基于matlab的车辆工程仿真实例 -回复

基于matlab的车辆工程仿真实例-回复一个基于MATLAB的车辆工程仿真实例是车辆悬挂系统的分析与优化。

在汽车设计中，悬挂系统起着至关重要的作用，影响着车辆的舒适性、稳定性和操控性。

使用MATLAB可以通过建立车辆的动力学模型，对悬挂系统进行仿真分析和优化。

以下是一个可能的仿真实例步骤：1. 建立车辆的动力学模型：使用MATLAB 建立车辆的多体动力学模型，包括车体、车轮、悬挂系统等。

模型可以考虑车辆的质量分布、悬挂系统的刚度和阻尼等参数。

2. 仿真悬挂系统的响应：在建立好动力学模型后，可以通过输入不同的激励，如车速、道路条件等，模拟悬挂系统的动态响应。

通过仿真可以了解车辆在不同工况下的悬挂系统行为，如车身姿态、轮胎力等。

3. 分析悬挂系统参数的影响：在模拟悬挂系统的响应之后，可以通过改变悬挂系统的参数，如刚度和阻尼，来分析这些参数对悬挂系统响应的影响。

通过分析可以确定最优的悬挂系统参数，以满足特定的需求，如舒适性、稳定性等。

4. 优化悬挂系统设计：基于分析结果，可以对悬挂系统的设计进行优化。

通过MATLAB的优化工具箱，可以使用不同的优化算法，如遗传算法、粒子群算法等，来搜索最优的悬挂系统参数组合。

5. 评估优化结果：在优化过程中，评估优化结果的有效性是必要的。

可以使用MATLAB的仿真工具再次运行优化后的悬挂系统模型，并进行性能评估。

性能评估可以包括悬挂系统的响应、车辆稳定性、舒适性等方面。

通过以上步骤，可以使用MATLAB对车辆悬挂系统进行仿真分析和优化，以改善车辆的悬挂系统性能。

这个实例展示了MATLAB 在车辆工程中的应用，可以帮助工程师更好地理解和优化车辆的悬挂系统。

汽车理论习题-matlab

汽车理论习题-matlab

《汽车理论》习题一、确定一轻型货车的动力性能。

1）绘制汽车驱动力与行驶阻力平衡图；2）求汽车最高车速与最大爬坡度；3）绘制汽车行驶加速度倒数曲线；用计算机求汽车用Ⅱ档起步加速行驶至70km/h所需的加速时间。

已知数据略。

（参见《汽车理论》习题第一章第3题）解题程序如下:用Matlab语言m1=2000; m2=1800; mz=3880;g=9.81; r=0.367; CdA=2.77; f=0.013; nT=0.85;ig=[5.56 2.769 1.644 1.00 0.793]; i0=5.83;If=0.218; Iw1=1.798; Iw2=3.598; #输入已知参数Iw=2*Iw1+4*Iw2;for i=1:69n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;end #从600~4000rpm油拟合公式计算发动机转距for j=1:5for i=1:69Ft(i,j)=Ttq(i)*ig(j)*i0*nT/r;ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Fz(i,j)=CdA*ua(i,j)^2/21.15+mz*g*f;endend #计算各档位下的驱动力和行驶阻力plot(ua,Ft,ua,Fz)title('汽车驱动力与行驶阻力平衡图');xlabel('ua(km/h)');ylabel('Ft(N)'); #作汽车驱动力与行驶阻力平衡图for k=1:175n1(k)=3300+k*0.1;Ttq(k)=-19.313+295.27*(n1(k)/1000)-165.44*(n1(k)/1000)^2+40.874*(n1(k)/1000)^33.8445*(n1(k)/1000)^4;Ft(k)=Ttq(k)*ig(5)*i0*nT/r;ua(k)=0.377*r*n1(k)/(ig(5)*i0);Fz(k)=CdA*ua(k)^2/21.15+mz*g*f;E(k)=abs((Ft(k)-Fz(k)));end #从3300rpm开始以0.1rpm为步长搜索直到找到for k=1:175 #驱动力与行驶阻力之差绝对值最小点(Ⅴ档处)if(E(k)==min(E)) #通过其他检验程序可估计出该点大致位置(3310rpm附近)disp('汽车最高车速=');disp(ua(k));disp('km/h');end #输出结果endfor p=1:150 #从2000rpm开始以0.5为步长搜索找驱动力n2(p)=2000+p*0.5; #与行驶阻力之差最大点(Ⅰ档处)Ttq(p)=-19.313+295.27*(n2(p)/1000)-165.44*(n2(p)/1000)^2+40.874*(n2(p)/1000) ^3-3.8445*(n2(p)/1000)^4;Ft(p)=Ttq(p)*ig(1)*i0*nT/r;ua(p)=0.377*r*n2(p)/(ig(1)*i0);Fz(p)=CdA*ua(p)^2/21.15+mz*g*f;af(p)=asin((Ft(p)-Fz(p))/(mz*g));end #通过其他检验程序可估计出该点大致位置(2000rpm附近)for p=1:150if(af(p)==max(af))i=tan(af(p));disp('汽车最大爬坡度=');disp(i);endendfor i=1:69 #开始计算绘制加速度倒数曲线所需数据n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4;endfor j=1:5for i=1:69deta=1+Iw/(mz*r^2)+If*ig(j)^2*i0^2*nT/(mz*r^2);ua(i,j)=0.377*r*n(i)/(ig(j)*i0);a(i,j)=(Ttq(i)*ig(j)*i0*nT/r-CdA*ua(i,j)^2/21.15-mz*g*f)/(deta*mz);if(a(i,j)<=0)a(i,j)=a(i-1,j);endif(a(i,j)>0.05) #剔除加速度值小于0.05的数据b1(i,j)=a(i,j);u1(i,j)=ua(i,j);elseb1(i,j)=a(i-1,j);u1(i,j)=ua(i-1,j);endb(i,j)=1/b1(i,j);endendx1=u1(:,1);y1=b(:,1);x2=u1(:,2);y2=b(:,2);x3=u1(:,3);y3=b(:,3);x4=u1(:,4);y4=b(:,4);x5=u1(:,5);y5=b(:,5);plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5);title('加速度倒数曲线');axis([0 120 0 30]);xlabel('ua(km/h)');ylabel('1/aj'); #作加速度倒数曲线for i=1:69 #从作出的加速度倒数曲线知,各档之间无交点A=ua(i,3)-ua(69,2); #因此换档时刻按低档最高转速对应的车速计if (A<1&A>0)j=i;endB=ua(i,4)-ua(69,3);if(B<2&B>0)k=i;endif(ua(i,4)<=70)m=i;endend #以上是计算三个换档时刻t=ua(1,2)*b(1,2); #开始计算加速时间,这里是以二档起步得到的原始时间for p1=2:69t1(p1)=(ua(p1,2)-ua(p1-1,2))*(b(p1,2)+b(p1-1,2))*0.5;t=t+t1(p1);end #为减小计算误差,采用梯形算法累加,以下均是for p2=j:69t2(p2)=(ua(p2,3)-ua(p2-1,3))*(b(p2,3)+b(p2-1,3))*0.5;t=t+t2(p2);endfor p3=k:mt3(p3)=(ua(p3,4)-ua(p3-1,4))*(b(p3,4)+b(p3-1,4))*0.5;t=t+t3(p3);endt=t+(ua(j,3)-ua(69,2))*b(69,2)+(ua(k,4)-ua(69,3))*b(69,3)+(70-ua(m,4))*b(m,4); #加上了换档处未计算的那部分时间tz=t/3.6; #换算单位disp('加速时间=');disp(tz);disp('s'); #输出计算结果计算结果及图形：汽车最高车速=99.0679km/h汽车最大爬坡度=0.3518加速时间=29.0585s二、计算与绘制题1 中货车的1）汽车功率平衡图；2）最高档与次高档的等速百公里油耗曲线。