汽车理论matlab作业;

1、 一、确定一轻型货车的动力性能。 1） 绘制汽车驱动力与行驶阻力平衡图； 2） 求汽车最高车速与最大爬坡度； 3） 绘制汽车行驶加速度倒数曲线；用计算机求汽车用档起步加速行驶至 70km/h 所需的加速时间。 已知数据略。（参见汽车理论习题第一章第3题） 解题程序如下:用matlab语言 (1)绘制汽车驱动力与行驶阻力平衡图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

2、=3.598; iw=2*iw1+4*iw2; for 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; end for j=1:5 for i=1:69 ft(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; end end plot(ua,ft,ua,ff,ua,ff+

3、fw) title(汽车驱动力与行驶阻力平衡图); xlabel(ua(km/h); ylabel(ft(n); gtext(ft1)gtext(ft2)gtext(ft3)gtext(ft4)gtext(ft5)gtext(ff+fw) (2)求最大速度和最大爬坡度for k=1:175 n1(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.3

4、77*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 for k=1:175 if(e(k)=min(e) disp(汽车最高车速=); disp(ua(k); disp(km/h); end end for p=1:150 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

5、)*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 for p=1:150 if(af(p)=max(af) i=tan(af(p); disp(汽车最大爬坡度=); disp(i); end end 汽车最高车速=99.0679km/h 汽车最大爬坡度=0.3518 （3） 计算2档起步加速到70km/h所需时间 for i=1:69 n(i)=(i+11)*50; ttq(i)=-19.313+295.27*(n(i)/100

6、0)-165.44*(n(i)/1000)2+40.874*(n(i)/1000)3-3.8445*(n(i)/1000)4; end for j=1:5 for i=1:69 deta=1+iw/(mz*r2)+if*ig(j)2*i02*nt/(mz*r2); 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.05) b1(i,j)=a(i,j); u1(i,j)=ua(i,j); else b1(i,j)=a(i-1

7、,j); u1(i,j)=ua(i-1,j); end b(i,j)=1/b1(i,j); end end x1=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); gtext(1/a1)gtext(1/a2)gtext(1/a3

8、)gtext(1/a4)gtext(1/a5) for i=1:69 a=ua(i,3)-ua(69,2); if (a0) j=i; end b=ua(i,4)-ua(69,3); if(b0) k=i; end if(ua(i,4)=70) m=i; end end t=ua(1,2)*b(1,2); for p1=2:69 t1(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:69 t2(p2)=(ua(p2,3)-ua(p2-1,3)*(b(p2,3)+b(p2-1,3)*0.5; t=

9、t+t2(p2); end for p3=k:m t3(p3)=(ua(p3,4)-ua(p3-1,4)*(b(p3,4)+b(p3-1,4)*0.5; t=t+t3(p3); end t=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); 加速时间=29.0585s 二、计算与绘制题1 中货车的1）汽车功率平衡图； 2）最高档与次高档的等速百公里油耗曲线。 已知数据略。（参见汽车理论习题第二章第7题） 解题程序如

10、下:用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; n1=815 1207 1614 2012 2603 3006 3403 3804; iw=2*iw1+4*iw2; nd=400; qid=0.299; for j=1:5 for i=1:69 n(i)=(i+11)*50; ttq(i)=-19.313+295.27*(n(i)/

11、1000)-165.44*(n(i)/1000)2+40.874*(n(i)/1000)3-3.8445*(n(i)/1000)4; pe(i)=n(i)*ttq(i)/9549; ua(i,j)=0.377*r*n(i)/(ig(j)*i0); pz(i,j)=(mz*g*f*ua(i,j)/3600.+cda*ua(i,j)3/76140.)/nt; end end plot(ua,pe,ua,pz); title(汽车功率平衡图); xlabel(ua(km/h); ylabel(pe,pz(kw); gtext(i)gtext(ii)gtext(iii)gtext(iv)gtext(v

12、)gtext(p阻) for j=1:5 for i=1:8 td(i)=-19.313+295.27*(n1(i)/1000.0)-165.44*(n1(i)/1000.0)2+40.874*(n1(i)/1000.0)3-3.8445*(n1(i)/1000.0)4; pd(i)=n1(i)*td(i)/9549; u(i,j)=0.377*n1(i)*r/(ig(j)*i0); end end b(1)=0.17768*pd(1)4-5.8629*pd(1)3+72.379*pd(1)2-416.46*pd(1)+1326.8; b(2)=0.043072*pd(2)4-2.0553*p

13、d(2)3+36.657*pd(2)2-303.98*pd(2)+1354.7; b(3)=0.0068164*pd(3)4-0.51184*pd(3)3+14.524*pd(3)2-189.75*pd(3)+1284.4; b(4)=0.0018555*pd(4)4-0.18517*pd(4)3+7.0035*pd(4)2-121.59*pd(4)+1122.9; b(5)=0.00068906*pd(5)4-0.091077*pd(5)3+4.4763*pd(5)2-98.893*pd(5)+1141.0; b(6)=0.00035032*pd(6)4-0.05138*pd(6)3+2.8

14、593*pd(6)2-73.714*pd(6)+1051.2; b(7)=0.00028230*pd(7)4-0.047449*pd(7)3+2.9788*pd(7)2-84.478*pd(7)+1233.9; b(8)=-0.000038568*pd(8)40.00075215*pd(8)3+0.71113*pd(8)245.291*pd(8) +1129.7; u1=u(:,1); u2=u(:,2); u3=u(:,3); u4=u(:,4); u5=u(:,5); b1=polyfit(u1,b,3); b2=polyfit(u2,b,3); b3=polyfit(u3,b,3); b

15、4=polyfit(u4,b,3); b5=polyfit(u5,b,3); for q=1:69 bh(q,1)=polyval(b1,ua(q,1); bh(q,2)=polyval(b2,ua(q,2); bh(q,3)=polyval(b3,ua(q,3); bh(q,4)=polyval(b4,ua(q,4); bh(q,5)=polyval(b5,ua(q,5); end for i=1:5 for q=1:69 q(q,i)=pz(q,i)*bh(q,i)/(1.02*ua(q,i)*7.05); end end plot(ua(:,4),q(:,4),ua(:,5),q(:,5

16、); title(四档五档等速百公里油耗图); xlabel(ua(km/h); ylabel(qs(l/100km); 三、改变1.3题中轻型货车的主减速器传动比，做出为5.17、5.43、5.83、6.17、6.33时的燃油经济性加速时间曲线，讨论不同值对汽车性能的影响。matlab程序：m1=2000 ; m2=1800 ; m=3880 ; r0=0.367 ; gt=0.85 ; f=0.013 ; cda=2.77 ; i0=5.83 ; if=0.218 ; iw1=1.798 ; iw2=3.598 ; ig5=5.56 2.769 1.644 1.00 0.793; ig0=

17、5.17 5.43 5.83 6.17 6.33; b=1326.8 -416.46 72.379 -5.8629 0.17768; 1354.7 -303.98 36.657 -2.0553 0.043072; 1284.4 -189.75 14.524 -0.51184 0.0068164; 1122.9 -121.59 7.0035 -0.18517 0.0018555; 1141.0 -98.893 4.4763 -0.091077 0.00068906; 1051.2 -73.714 2.8593 -0.05138 0.00035032; 1233.9 -84.478 2.9788

18、-0.047449 0.00028230; 1129.7 -45.291 0.71113 -0.00075215 -0.000038568;n=815 1207 1614 2012 2603 3006 3403 3804;for i=1:5 for k=1:8ua(i,k)=0.377*0.367*n(k)/(ig0(i)*ig5(5); 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; f5(i,k)=0.013*3880*9.8+2.77

19、.*ua(i,k)2/21.15;pe(i,k)=f5(i,k)*ua(i,k)/(3600*0.85);b5(i,k)=b(k,1)+b(k,2)*pe(i ,k)+b(k,3)*pe(i,k)2+b(k,4)*pe(i,k)3+b(k,5)*pe(i,k)4; endend ua1=25; s1=50;fa5=0.013*3880*9.8+2.77.*ua1.2/21.15;pe5=fa5.*ua1/(3600*0.85);d1=polyfit(pe(1,:),b5(1,:),3);ba1=polyval(d1,pe5);d2=polyfit(pe(2,:),b5(2,:),3);ba2=

20、polyval(d2,pe5);d3=polyfit(pe(3,:),b5(3,:),3);ba3=polyval(d3,pe5);d4=polyfit(pe(4,:),b5(4,:),3);ba4=polyval(d4,pe5);d5=polyfit(pe(5,:),b5(5,:),3);ba5=polyval(d5,pe5) ;ba=ba1 ba2 ba3 ba4 ba5;qa1=pe5.*ba*50/(ua1*102*7)ua2=25:40;q2=1+(2*iw1+4*iw2)/(m*r02)+if*ig5(2).2*ig0(2)2*gt/(m*r02)fb5=0.013*3880*9.

21、8+2.77.*ua2.2/21.15+q2*m*0.25;pb5=fb5.*ua2/(3600*0.85);db1=polyfit(pe(1,:),b5(1,:),3);bb1=polyval(db1,pb5);db2=polyfit(pe(2,:),b5(2,:),3);bb2=polyval(db2,pb5);db3=polyfit(pe(3,:),b5(3,:),3);bb3=polyval(db3,pb5);db4=polyfit(pe(4,:),b5(4,:),3);bb4=polyval(db4,pb5);db5=polyfit(pe(5,:),b5(5,:),3);bb5=po

22、lyval(db5,pb5);bb=bb1 bb2 bb3 bb4 bb5;pb=pb5 pb5 pb5 pb5 pb5;qb=pb.*bb/(367.1*7);for i=1:5 for j=1:15 qb(i,j)=qb(i,j)+qb(i,j+1); end end qb2=sum(qb)ua3=40; s2=250;fc5=0.013*3880*9.8+2.77.*ua3.2/21.15;pc5=fc5.*ua3/(3600*0.85);dc1=polyfit(pe(1,:),b5(1,:),3);bc1=polyval(dc1,pc5);dc2=polyfit(pe(2,:),b5(

23、2,:),3);bc2=polyval(dc2,pc5);dc3=polyfit(pe(3,:),b5(3,:),3);bc3=polyval(dc3,pc5);dc4=polyfit(pe(4,:),b5(4,:),3);bc4=polyval(dc4,pc5);dc5=polyfit(pe(5,:),b5(5,:),3);bc5=polyval(dc5,pc5) ;bc=bc1 bc2 bc3 bc4 bc5;qc3=pc5.*bc*250/(ua3*102*7)ua4=40:50;q4=1+(2*iw1+4*iw2)/(m*r02)+if*ig5(4).2*ig0(4)2*gt/(m*r

24、02)fd5=0.013*3880*9.8+2.77.*ua4.2/21.15+q4*m*0.2;pd5=fd5.*ua4/(3600*0.85);dd1=polyfit(pe(1,:),b5(1,:),3);bd1=polyval(dd1,pd5);dd2=polyfit(pe(2,:),b5(2,:),3);bd2=polyval(dd2,pd5);dd3=polyfit(pe(3,:),b5(3,:),3);bd3=polyval(dd3,pd5);dd4=polyfit(pe(4,:),b5(4,:),3);bd4=polyval(dd4,pd5);dd5=polyfit(pe(5,:

25、),b5(5,:),3);bd5=polyval(dd5,pd5);bd=bd1 bd2 bd3 bd4 bd5;pd=pd5 pd5 pd5 pd5 pd5;qd=pd.*bd/(367.1*7);for i=1:5 for j=1:10 qd(i,j)=qd(i,j)+qd(i,j+1); end end qd4=sum(qd)ua5=50; s2=250;ff5=0.013*3880*9.8+2.77.*ua5.2/21.15;pf5=ff5.*ua5/(3600*0.85);df1=polyfit(pe(1,:),b5(1,:),3);bf1=polyval(df1,pf5);df2=

26、polyfit(pe(2,:),b5(2,:),3);bf2=polyval(df2,pf5);df3=polyfit(pe(3,:),b5(3,:),3);bf3=polyval(df3,pf5);df4=polyfit(pe(4,:),b5(4,:),3);bf4=polyval(df4,pf5);df5=polyfit(pe(5,:),b5(5,:),3);bf5=polyval(df5,pf5) ;bf=bf1 bf2 bf3 bf4 bf5;qf5=pf5.*bf*250/(ua3*102*7)qi=0.299;qg=(50-25)/3.6/0.36*0.299;qg6=qg qg

27、qg qg qgq=qa1 qb2 qc3 qd4 qf5 qg6;qz=sum(q)/1075*100for k=1:5for i=1:3401; for j=1:5; n(i)=i+599; ua(i,j)=0.377*r0*n(i)./(ig5(j)*ig0(k); q(j)=1+(2*iw1+4*iw2)/(m*r02)+if*ig5(j).2*ig0(k)2*gt/(m*(r0)2); endenduamax=max(ua);uamin=min(ua);ua2=uamin(2):uamax(2);n=ua2*ig0(k)*ig5(2)/(0.377*0.367);ttq=-19.31

28、3+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4; ft2=ttq*ig0(k)*ig5(2)*0.85/0.367;f2=0.013*3880*9.8+2.77*ua2.2/21.15;a2=(ft2-f2)./(q(2)*m);t2=trapz(ua2,1./a2)/3.6;ua3=uamax(2):uamax(3);n3=ua3*ig0(k)*ig5(3)/(0.377*0.367);ttq3=-19.313+295.27*(n3/1000)-165.44*(n3/1000).2+40.874*

29、(n3/1000).3-3.8445*(n3/1000).4; ft3=ttq3*ig0(k)*ig5(3)*0.85/0.367;f3=0.013*3880*9.8+2.77*ua3.2/21.15;a3=(ft3-f3)./(q(3)*m);t3=trapz(ua3,1./a3)/3.6;ua4=uamax(3):70;n4=ua4*ig0(k)*ig5(4)/(0.377*0.367);ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).2+40.874*(n4/1000).3-3.8445*(n4/1000).4; ft4=ttq4*ig0(k

30、)*ig5(4)*0.85/0.367;f4=0.013*3880*9.8+2.77*ua4.2/21.15;a4=(ft4-f4)./(q(4)*m);t4=trapz(ua4,1./a4)/3.6;t(k)=t2+t3+t4;end f=polyfit(qz,t,2);qq=qz(1):0.001:qz(5);b=polyval(f,qq)plot(qq,b,qz,t,*) gtext(5.17)gtext(5.43)gtext(5.83)gtext(6.17)gtext(6.33)title(燃油经济性-动力性曲线)xlabel(燃油经济性/l(100km)-1)ylabel(t/s)四

31、、一中型货车装有前后制动器分开的双管路制动系，其有关参数如下：载荷质量（kg）质心高hg/m轴距l/m质心至前轴距离a/m制动力分配系数空载40800.8453.9502.1000.38满载92901.1703.9502.9500.381） 计算并绘制利用附着系数曲线和制动效率曲线2） 求行驶车速ua30km/h，在0.80路面上车轮不抱死的制动距离。计算时取制动系反应时间0.02s，制动减速度上升时间0.02s。3） 求制动系前部管路损坏时汽车的制动距离s，制动系后部管路损坏时汽车的制动距离。matlab程序：m1=4080;hg1=0.845;a1=2.100; m2=9290;hg2=1

32、.17;a2=2.95; beta=0.38;l=3.95;z=0:0.05:1gf1=beta.*z*l./(l-a1+z*hg1);gf2=beta.*z*l./(l-a2+z*hg2);gr1=(1-beta).*z*l./(a1-z*hg1);gr2=(1-beta).*z*l./(a2-z*hg2);g=z;for i=1:21 if (z(i)0.15); g3(i)=z(i)+0.08; end if(z(i)=0.3); g3(i)=0.38+(z(i)-0.3)/0.74; end end z1=0.15:0.01:0.3; g4=z1-0.08; plot(z,gf1,-.

33、,z,gf2,z,gr1,-.,z,gr2,z,g,z,g3,xk,z1,g4,x)axis(0 1 0 1.2)title(利用附着系数与制动强度的关系曲线)xlabel(制动强度z/g)ylabel(利用附着系数g)gtext(空车前轴)gtext(空车后轴)gtext(满载前轴)gtext(满载后轴)gtext(ece法规)c=0:0.05:1;er1=(a1/l)./(1-beta)+c*hg1/l)*100;ef=(l-a2)/l./(beta-c*hg2/l)*100;er=(a2/l)./(1-beta)+c*hg2/l)*100;plot(c,er,c,ef,c,er1)axi

34、s(0 1 0 100)title(前后附着效率曲线)xlabel(附着系数c)ylabel(制动效率（%）)gtext(满载)gtext(ef)gtext(er)gtext(空载)gtext(er)c1=0.8e1=(ak1/l)./(1-beta)+c1*hg1/l);e2=(am2/l)/(1-beta)+c1*hg2/l);a1=e1*c1*9.8;a2=e2*c1*9.8;ua=30;i21=0.02;i22=0.02;s1=(i21+i22/2)*ua/3.6+ua2/(25.92*ak1);s2=(i21+i22/2)*ua/3.6+ua2/(25.92*am2);disp(满载

35、时不抱死的制动距离=)disp(s2)disp(空载时不抱死的制动距离=)disp(s1)满载时不抱死的制动距离=5.3319空载时不抱死的制动距离=6.8406beta3=1 beta4=0 ekr=(a1/l)/(1-beta4)+c1*hg1/l);ekf=(l-a1)/(beta3*l-c1*hg1);emf=(l-a2)/l./(beta3-c1*hg2/l);emr=(a2/l)./(1-beta4)+c1*hg2/l);akr=0.8*9.8*ekr;akf=0.8*9.8*ekf;amr=0.8*9.8*emr;amf=0.8*9.8*emf;skr=(i21+i22/2)*u

36、a/3.6+ua2/(25.92*akr);skf=(i21+i22/2)*ua/3.6+ua2/(25.92*akf);smf=(i21+i22/2)*ua/3.6+ua2/(25.92*amf);smr=(i21+i22/2)*ua/3.6+ua2/(25.92*amr);disp(空车后管路失效时制动距离)disp(skf)disp(空车前管路失效时制动距离)disp(skr)disp(满载后管路失效时制动距离)disp(smf)disp(满载前管路失效时制动距离)disp(smr)运行结果为：空车后管路失效时制动距离8.0879空车前管路失效时制动距离10.0061满载后管路失效时制动

37、距离 13.5986满载前管路失效时制动距离7.5854五、二自由度轿车模型的有关参数如下： 总质量 m=1818.2kg绕oz轴转动惯量 轴距 l=3.048m质心至前轴距离 a=1.463m质心至后轴距离 b=1.585m前轮总侧偏刚度 k1=-62618n/rad后轮总侧偏刚度 k2=-110185n/rad转向系总传动比 i=20试求：1) 稳定性因数k、特征车速uch。2) 稳态横摆角速度增益曲线、车速u=22.35m/s时的转向灵敏度。3) 静态储备系数s.m.，侧向加速度为0.4g时的前、后轮侧偏角绝对值之差与转弯半径的比值r/r0(r0=15m)。4) 车速u=30.56m/s

38、时，瞬态响应的横摆角速度波动的固有（圆）频率、阻尼比、反应时间与峰值反应时间matlab程序：m=1818.2;iz=3885;l=3.048;a=1.463;b=1.585;k1=-62618;k2=-110185;i=20;g=9.8;r0=15;u1=30.56;k=m/l2*(a/k2-b/k1);uch=1/k(1/2);u=0:0.05:25;s=u/l./(1+k*u.2);disp(稳定因素k=);disp(k);disp(特征车速uch=);disp(uch);plot(u,s);xlabel(ua/(m/s);ylabel(稳态横摆增益);title(汽车的稳态横态摆角速度

39、增益曲线);disp(ua=22.35m/s时，转向灵敏度为);disp(s(448);sm=k2/(k1+k2)-a/l;a=k*0.4*g*l;q=l/r0;r=l/(q-a);r=r/r0;disp(静态储备系数s.m.=);disp(sm);disp(前后轮侧偏角绝对值之差(a1-a2)=);disp(a);disp(转弯半径的比值=);disp(r);w0=l/u1*(k1*k2*(1+k*u12)/(m*iz)(1/2);n=(-m*(a2*k1+b2*k2)-iz*(k1+k2)/(2*l*(m*iz*k1*k2*(1+k*u12)(1/2);t=atan(1-n2)(1/2)/

40、(-m*u1*a*w0/(l*k2)-n)/(w0*(1-n2)(1/2);e=atan(1-n2)(1/2)/n)/(w0*(1-n2)(1/2)+t;disp(横摆角速度波动时的固有频率为);disp(w0);disp(阻尼比为);disp(n);disp(反应时间为);disp(t);disp(达到第一峰值的时间为);disp(e);运算结果：稳定因素k=0.0024特征车速uch=20.6053ua=22.35m/s时，转向灵敏度为 3.3690静态储备系数s.m.=0.1576前后轮侧偏角绝对值之差(a1-a2)=0.0281转弯半径的比值= 1.1608横摆角速度波动时的固有频率为

41、5.5758阻尼比为0.5892反应时间为0.1811达到第一峰值的时间为0.3899六、车身-车轮双质量系统参数：。“人体-座椅”系统参数：。车速，路面不平度系数，参考空间频率n0=0.1m-1。计算时频率步长，计算频率点数。1) 计算并画出幅频特性、和均方根值谱、谱图。进一步计算值2) 改变“人体-座椅”系统参数：。分析值随的变化。3) 分别改变车身-车轮双质量系统参数：，。绘制三个响应量均方根值随以上四个系统参数变化的曲线。matlab程序：f0=1.5;g=0.25;r=9;u=10;fs=3;gs=0.25;v=20;gqn0=2.56e-008;n0=0.1;df=0.2;n=18

42、0;f=0:0.2:36;ff0=f/f0;d=(1-ff0.2).*(1+r-1/u*ff0.2-1).2+4*g2*ff0.2.*r-(1/u+1)*ff0.2.2;z1q=r*sqrt(1-ff0.2).2+4*g2*ff0.2)./d);z2z1=sqrt(1+(2*g*ff0).2)./(1-ff0.2).2+(2*g*ff0).2);ffs=f/fs;pz2=sqrt(1+(2*gs*ffs).2)./(1-ffs.2).2+(2*gs*ffs).2)loglog(f,z1q,*,f,z2z1,.,f,pz2);axis(0 100 0 10)set(gca,xtick,0 0.1

43、 1 10 100)set(gca,ytick,0 0.1 1 10)xlabel(激振频率f/hz) w=2*pi*f;gz1=w.2.*z1q.*(2*3.14./w).*sqrt(gqn0*n02*v);loglog(f,gz1)hold ongz2=w.2.*z1q.*z2z1.*(2*3.14./w).*sqrt(gqn0*n02*v);title(车轮部分加速度均方根值)xlabel(激振频率f/hz)ylabel(|gz1|)loglog(f,gz2)title(车身部分加速度均方根值)xlabel(激振频率f/hz)ylabel(|gz2|)gza=w.2.*z1q.*z2z1

44、.*pz2.*(2*3.14./w).*sqrt(gqn0*n02*v);loglog(f,gza)title(人体加速度均方根值)xlabel(激振频率f/hz)ylabel(|gza|)format longf0=1.5;g=0.25;r=9;u=10;fs=3;gs=0.25;v=20;gqn0=2.56*10(-8);n0=0.1;df=0.2;n=180;f=0.000001:0.2:36;ff0=f/f0;ffs=f/fs;w=2*pi*f;d=(1-ff0.2).*(1+r-ff0.2/u)-1).2+4.*g2.*ff0.2.*(r-(1/u+1).*ff0.2).2;z1q=

45、r*sqrt(1-ff0.2).2+4*g2*ff0.2)./d);z2z1=sqrt(1+(2*g*ff0).2)./(1-ff0.2).2+(2*g*ff0).2);ffs=f/fs;pz2=sqrt(1+(2*gs*ffs).2)./(1-ffs.2).2+(2*gs*ffs).2);ez=sqrt(sum(16.*(3.14)4.*f.2.*gqn0.*n02.*v*df);disp(路面不平度加速度均方根值=)disp(ez)ez1=sqrt(sum(z1q.2.*16.*pi4.*gqn0.*v.*f.2*df*n02);disp(车轮加速度均方根值=)disp(ez1)ez2=sqrt(sum(z2z1.2.*z1q.2.*16.*pi4.*gqn0.*v.*f.2*n02*df);disp(车身加速度均方根值=)disp(ez2)gza=w.2.*z1q.*z2z1.*pz2.*(2*3.14./w).*sqrt(gqn0*n02*v);eza=sqrt(sum(gz