汽车运用作业matlab程序
1.3
1)绘制驱动力与行驶阻力平衡图
利用Matlab画图，具体语句如下：
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig1=5.56;
ua=0.377*r*n/(ig1*io);
Ft1=Ttq*ig1*io*t/r;
plot(ua,Ft1)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig2=2.769;
ua=0.377*r*n/(ig2*io);
Ft2=Ttq*ig2*io*t/r;
plot(ua,Ft2)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig3=1.644;
ua=0.377*r*n/(ig3*io);
Ft3=Ttq*ig3*io*t/r;
plot(ua,Ft3)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig4=1.00;
ua=0.377*r*n/(ig4*io);
Ft4=Ttq*ig4*io*t/r;
plot(ua,Ft4)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig5=0.793;
ua=0.377*r*n/(ig5*io);
Ft5=Ttq*ig5*io*t/r;
plot(ua,Ft5)
hold on;
m=3880;f=0.012;g=9.8;CdA=2.77;
u=0:1:120;
Ff=m*g*f;
Fw=CdA*u.^2/21.15;
Fr=Ff+Fw;
Ft5=Ttq*ig5*io*t/r;
Ff=m*g*f;
Fw=CdA*u.^2/21.15;
Fr=Ff+Fw;
plot(u,Fr)
运行生成图像如下:
2）求汽车最高车速,最大爬坡度及克服该坡度时相应的附着率当Ft=Ff+Fw时,汽车处于稳定状态,并达到最高车速
由图可知
Ｕmax=101km/h
最大爬坡度及附着率：
一档时有最大爬坡度,由驱动力与阻力图
当ua=8.037km/h时: Gsin错误！未找到引用源。Ft-(Ff+Fw)最大。
可知：Ft=13058N,Fr=316,F=Ft-Fr=12742N.
由: Gsinα=F可得：
α=19.2°.
i=tanα=0.348
货车为后轮驱动C2=i/(a/L+hg*i/L)
其中L=3.2m, a=1.947m, hg=0.9m.
解得C2=0.4927
3）绘制加速度倒数曲线：
利用Matlab画图，具体语句如下：
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; i0=5.83;t=0.844;r=0.367;ig1=5.56;Iw1=1.798;Iw2=3.598;If=0.218;m=3880;r=0.367;
g=9.8;f=0.012;CdA=2.77;
ua=0.377*r*n/(ig1*io);
Ft1=Ttq*ig1*io*t/r;
dert=1+(Iw1+Iw2)/(m*r.^2)+If*ig1.^2*i0.^2*t/(m*r.^2);
Ff=m*g*f;
Fw1=CdA*ua.^2/21.15;
Fr1=Ff+Fw1;
a1=(dert*m)./(Ft1-Fr1);
plot(ua,a1)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig2=2.769;Iw1=1.798;Iw2=3.598;If=0.218;m=3880;r=0.367;
g=9.8;f=0.012;CdA=2.77;
ua=0.377*r*n/(ig2*io);
Ft2=Ttq*ig2*io*t/r;
dert=1+(Iw1+Iw2)/(m*r.^2)+If*ig2.^2*io.^2*t/(m*r.^2);
Ff=m*g*f;
Fw2=CdA*ua.^2/21.15;
Fr2=Ff+Fw2;
a2=(dert*m)./(Ft2-Fr2);
plot(ua,a2)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig3=1.644;Iw1=1.798;Iw2=3.598;If=0.218;m=3880;r=0.367;
g=9.8;f=0.012;CdA=2.77;
ua=0.377*r*n/(ig3*io);
Ft3=Ttq*ig3*io*t/r;
dert=1+(Iw1+Iw2)/(m*r.^2)+If*ig3.^2*i0.^2*t/(m*r.^2);
Ff=m*g*f;
Fw=CdA*ua.^2/21.15;
Fr=Ff+Fw;
a3=(dert*m)./(Ft3-Fr);
plot(ua,a3)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig4=1.00;Iw1=1.798;Iw2=3.598;If=0.218;m=3880;r=0.367;
g=9.8;f=0.012;CdA=2.77;
ua=0.377*r*n/(ig4*io);
Ft4=Ttq*ig4*io*t/r;
dert=1+(Iw1+Iw2)/(m*r.^2)+If*ig4.^2*io.^2*t/(m*r.^2);
Ff=m*g*f;
Fw=CdA*ua.^2/21.15;
Fr=Ff+Fw;
a4=(dert*m)./(Ft4-Fr);
plot(ua,a4)
hold on;
运行程序生成图形如下：
求加速时间语句如下：
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.87*(n/1000).^3-3.8445*(n/1000).^4; io=5.83;t=0.844;r=0.367;ig2=2.769;ig3=1.644;ig4=1.00;
ua2=0.377*r*n/(ig2*i0);ua3=0.377*r*n/(ig3*i0);ua4=0.377*r*n/(ig4*i0);
Ft2=Ttq*ig2*i0*t/r;Ft3=Ttq*ig3*i0*t/r;Ft4=Ttq*ig4*i0*t/r;
m=3880;f=0.012;g=9.8;CdA=2.77;
Ff2=m*g*f;Fw2=CdA*ua2.^2/21.15;Fr2=Ff2+Fw2;F2=Ft2-Fr2;a2=m./F2;
Ff3=m*g*f;Fw3=CdA*ua3.^2/21.15;Fr3=Ff3+Fw3;F3=Ft3-Fr3;a3=m./F3;
Ff4=m*g*f;Fw4=CdA*ua4.^2/21.15;Fr4=Ff4+Fw4;F4=Ft4-Fr4;a4=m./F4;
m2=5.142415:0.001:34.282777;m3=34.28277:0.001:57.74269;m4=57.7427:0.001:70;
t2=interp1(ua2,a2,m2,'spline');t3=interp1(ua3,a3,m3,'spline');t4=interp1(ua4,a4,m4,'spline'); plot(m2,t2,m3,t3,m4,t4);
hold on;
%面积
A=0.001*sum(t2-0)+0.001*sum(t3-0)+0.001*sum(t4-0)
T=86.8332/3.6
运行结果：
A= 88.5961；
T= 24.1203
即加速时间为24.1203 S
各档位加速时间段图像如下：
2.7
1）汽车功率平衡图：
Matlab语句如下：
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
io=5.83;r=0.367;ig1=5.56;
ua=0.377*r*n/(ig1*io);
Ttqn=-19.313*n+295.27*1000*(n/1000).^2-165.44*1000*(n/1000).^3+40.874*1000*(n/1000).^ 4-3.8445*1000*(n/1000).^5;
P1=0.377*Ttqn/3600;
plot(ua,P1)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
io=5.83;r=0.367;ig2=2.769;
ua=0.377*r*n/(ig2*io);
Ttqn=-19.313*n+295.27*1000*(n/1000).^2-165.44*1000*(n/1000).^3+40.874*1000*(n/1000).^ 4-3.8445*1000*(n/1000).^5;
P2=0.377*Ttqn/3600;
plot(ua,P2)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
io=5.83;r=0.367;ig3=1.644;
ua=0.377*r*n/(ig3*io);
Ttqn=-19.313*n+295.27*1000*(n/1000).^2-165.44*1000*(n/1000).^3+40.874*1000*(n/1000).^ 4-3.8445*1000*(n/1000).^5;
P3=0.377*Ttqn/3600;
plot(ua,P3)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
io=5.83;r=0.367;ig4=1.00;
ua=0.377*r*n/(ig4*io);
Ttqn=-19.313*n+295.27*1000*(n/1000).^2-165.44*1000*(n/1000).^3+40.874*1000*(n/1000).^ 4-3.8445*1000*(n/1000).^5;
P4=0.377*Ttqn/3600;
plot(ua,P4)
hold on;
n=600:1:4000;
Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
io=5.83;r=0.367;ig5=0.793;
ua=0.377*r*n/(ig5*io);
Ttqn=-19.313*n+295.27*1000*(n/1000).^2-165.44*1000*(n/1000).^3+40.874*1000*(n/1000).^ 4-3.8445*1000*(n/1000).^5;
P5=0.377*Ttqn/3600;
plot(ua,P5)
hold on;
ua=0:1:110;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
Pr=(m*g*f*ua+CdA*ua.^3/21.15)/(1000*3.6*t);
plot(ua,Pr)
生成图形如下：
2）最高档与次高档百公里油耗
次高档百公里油耗:
利用Matlab计算出八组数据，Matlab语句如下：
n=815;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1426.8;B1=-416.46;B2=72.379;B3=-5.8629;B4=0.17768;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs1=P*b/(1.02*ua*roug)
n=1207;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1354.7;B1=-303.98;B2=36.657;B3=-2.0553;B4=0.043072; ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs2=P*b/(1.02*ua*roug)
n=1614;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1284.4;B1=-189.75;B2=14.524;B3=-0.51184;B4=0.0068164; ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs3=P*b/(1.02*ua*roug)
n=2012;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1122.9;B1=-121.59;B2=7.0035;B3=-0.18517;B4=0.0018555;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs4=P*b/(1.02*ua*roug)
n=2063;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1141.0;B1=-98.893;B2=4.4763;B3=-0.091077;B4=0.00068906;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs5=P*b/(1.02*ua*roug)
n=3006;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1051.2;B1=-73.714;B2=2.8593;B3=-0.05138;B4=0.00035032;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs6=P*b/(1.02*ua*roug)
n=3403;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77; t=0.844;
B0=1233.9;B1=-84.478;B2=2.9788;B3=-0.047449;B4=0.00028230;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs7=P*b/(1.02*ua*roug)
n=3884;r=0.367;ig4=1.00;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1129.7;B1=-45.291;B2=0.71113;B3=-0.00075215;B4=-0.000038568; ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs8=P*b/(1.02*ua*roug)
舍去第一组数据
利用EXCEL画出图像：
最高档百公里油耗：
利用Matlab计算出八组数据，Matlab语句如下：
n=815;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1426.8;B1=-416.46;B2=72.379;B3=-5.8629;B4=0.17768;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs1=P*b/(1.02*ua*roug)
n=1207;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1354.7;B1=-303.98;B2=36.657;B3=-2.0553;B4=0.043072; ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs2=P*b/(1.02*ua*roug)
n=1614;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1284.4;B1=-189.75;B2=14.524;B3=-0.51184;B4=0.0068164; ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs3=P*b/(1.02*ua*roug)
n=2012;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1122.9;B1=-121.59;B2=7.0035;B3=-0.18517;B4=0.0018555;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs4=P*b/(1.02*ua*roug)
n=2603;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1141.0;B1=-98.893;B2=4.4763;B3=-0.091077;B4=0.00068906;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs5=P*b/(1.02*ua*roug)
n=3006;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1051.2;B1=-73.714;B2=2.8593;B3=-0.05138;B4=0.00035032;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs6=P*b/(1.02*ua*roug)
n=3403;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1233.9;B1=-84.478;B2=2.9788;B3=-0.047449;B4=0.00028230;
ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs7=P*b/(1.02*ua*roug)
n=3884;r=0.367;ig4=0.793;io=5.83;roug=7.00;
m=3880;f=0.012;g=9.8;CdA=2.77;t=0.844;
B0=1129.7;B1=-45.291;B2=0.71113;B3=-0.00075215;B4=-0.000038568; ua=0.377*n*r/(ig4*io)
P=(m*g*f*ua+CdA*ua.^3/21.15)/(3600*t);
b=B0+B1*P+B2*P.^2+B3*P.^3+B4*P.^4;
Qs8=P*b/(1.02*ua*roug)
舍去最后一组数值，利用EXCEL画图：
汽车理论作业五
2、 某轿车总质量为1900kg ，轴距为3.1m ，单个轮胎的侧偏刚度数值为 40000 N/rad 。已知该车稳态转向特性为过多转向，临界车速为180km/h 。 （ 1）计算该车的稳定性因数和静态储备系数； （ 2）计算该车的轴荷分配比例； （ 3）在质心位置、轴距和前、后轮胎的型号已定的情况下，试找出五种 改善其转向特性的方法。 答：1， 026.01041900/1.380000800008000080000../100.4) 6.3/180(11 42112122 2422-=???+?-=-+=-+=?-=-=-=--K m L K K K L a K K K m s m s U k Cr 2,2242122/100.4)40000 21.3)400002/((1.3/1900)(m s a a K b K a L m K -?-=?---?-=-= 解得：a=1.63m 后轴荷分布比：a/L=1.63/3.1=52.6% 前轴荷分布比：47.4% 3，前轮气压减小，后轮气压增大。 前轴加装横向稳定杆。 前悬架蚕蛹双横臂式等类型 悬架，后悬架采用单横臂式或非独立式。 采用前轮驱动。 合理利用变形转向，如后轮随动 合理利用侧倾转向 3、教材课后习题5.2。 答:轿车前悬架加装横向稳定杆后，前悬架侧倾角刚度1r K ?增大，整车侧倾角刚度增大，车厢侧倾角r φ减小;在分析侧倾时垂直载荷在左、右车轮上的重新分配时，可以得到: 当前悬架增加横向稳定杆后汽车前悬架的侧倾角刚度增大，后悬架侧倾角刚度不变， 前悬架作用于车厢的恢复力矩增加（总侧倾力矩不变），而后悬架作用于车厢的恢复力矩减 小，所以汽车前轴左、右车轮载荷变化量较后轴大。
MATLAB编程作业
《Matlab 编程训练》 作业 专 业 学生姓名 班级 学 号 指导教师 完成日期
实训一 MATLAB 语言介绍和数值计算 1．先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存变量。 12 2sin851z e =+ . 2. 已知 1234413134787,2033657327A B --???? ????==???? ????-???? ，求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） A+6*B:
A-B+I: (2)A*B和A.*B A*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] c=A*B 结果： A.*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=A.*B 结果：
(3)A^3和A.^3 A^3程序： A=[12 34 -4;34 7 87;3 65 7] E=A^3 结果： A．^3程序： A=[12 34 -4;34 7 87;3 65 7] C=A.^3 (4)A/B及B\A A/B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] C=A/B 结果：
B\A程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=B\A 结果： (5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并（3）保存变量（mat文件）程序： A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; C=B*inv(A); D=C(2:3,2:3) 结果：
汽车理论课后作业matlab编程详解(带注释)分析解析
1.3matlab程序： （1） %驱动力-行驶阻力平衡图 %货车相关参数。 m=3880; g=9.8; nmin=600; nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];y=0.85;r=0.367;f=0.013;CdA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598; n=600:10:4000; %发动机转数转换成汽车行驶速度。 ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; ua5=0.377*r*n/ig(5)/i0; %计算各档位驱动力。 Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.84 45*(n/1000).^4; Ft1=Tq*ig(1)*i0*y/r; Ft2=Tq*ig(2)*i0*y/r; Ft3=Tq*ig(3)*i0*y/r; Ft4=Tq*ig(4)*i0*y/r; Ft5=Tq*ig(5)*i0*y/r; %计算行驶阻力。 Fz1=m*g*f+2.77*ua1.^2/21.15; Fz2=m*g*f+2.77*ua2.^2/21.15; Fz3=m*g*f+2.77*ua3.^2/21.15; Fz4=m*g*f+2.77*ua4.^2/21.15; Fz5=m*g*f+2.77*ua5.^2/21.15; %驱动力-行驶阻力平衡图。 plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua1,Fz1,'k',ua2,Fz2,'k', ua3,Fz3,'k',ua4,Fz4,'k',ua5,Fz5,'k'); title('驱动力-行驶阻力平衡图'); xlabel('ua(km/s)');
MATLAB程序设计作业
Matlab程序设计 班级 姓名 学号
《MATLAB程序设计》作业 1、考虑如下x-y 一组实验数据： x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot的原始数据、一次拟合曲线和三次拟合曲线，给出MATLAB代码和运行结果。 代码如下： x=[1,2,3,4,5,6,7,8,9,10]; y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]; plot(x,y); title('原始数据'); p=polyfit(x,y,1); q=polyval(p,x); figure,plot(x,q); title('一次拟合'); p=polyfit(x,y,2); q=polyval(p,x); figure,plot(x,q); title('二次拟合'); 运行结果如下：
1 2 3 4 5 6 7 8 9 10 12 3 4 5 6 7 8 原始数据 123 456789 102 2.5 3 3.54 4.5 55.56 6.57一次拟合 123456789 101 2 3 4 5 6 7 二次拟合 2、在[0，3π]区间，绘制y=sin(x)曲线（要求消去负半波，即(π，2π)区间内的函数值置零），求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 代码如下： clear clc x=(0:0.01:3*pi); y=sin(x); plot(x,y); y1=(y>=0).*y; figure,plot(x,y1);
汽车理论大作业
二自由度轿车模型的有关参数如下： 总质量 m ＝1818．2kg 绕z o 轴转动惯量 23885z I kg m =? 轴距 L=3.048m 质心至前轴距离 a=1.463m 质心至后轴距离 b=1.585m 前轮总侧偏刚度 rad N k /626181-= 后轮总侧偏刚度 2110185k =- /N rad 转向系总传动比 i=20 试求： 1)稳定性因数K 、特征车速ch u 。 2)稳态横摆角速度增益曲线r ωδ)s ----a u 车速u=22.35m/s 时的转向灵敏度r sw ωδ。 3)静态储备系数S.M.，侧向加速度为0．4g 时的前、后轮侧偏角绝对值之差12a a -与转弯半径的比值R/R 0(R 0=15m)。 4)车速u=30.56m/s,瞬态响应的横摆角速度波动的固有(圆)频率0ω、阻尼比ζ、反应时间τ与峰值反应时间ε。 提示： 1) 稳定性系数：??? ? ??-=122k b k a L m K =0.002422/m s 特征车速K u ch 1= =20.6s m /=74.18km/h 2) 转向灵敏度21Ku L u s r +=???δω=0.618 3) ()211αα-=L a K y ? 21αα-=0.0281rad δL R = 0 () 21ααδ--=L R ?0R R =1.16 4) 固有圆频率 m c '= 0ω=5.58rad/s.0f =0.8874Hz 阻尼比m h ' =02ωξ=0.5892 反应时间ω τΦ-== 0.1811s 峰值反应时间 ωξωωεΦ-=0 arctg =0.3899s
汽车理论课后习题Matlab程序
1.3确定一轻型货车的动力性能（货车可装用 4挡或5挡变速器，任选 其中的一种进行整车性能计算）： 1） 绘制汽车驱动力与行驶阻力平衡图。 2） 求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3） 绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用 2档起步加速行 驶至70km/h 的车速一时间曲线，或者用计算机求汽车用 2档起步加速行驶至 70km/h 的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n 曲线的拟合公式为 19.313 295.27(孟o )165,44(^)2 40-874(^)3 环45為4 式中，Tq 为发动机转矩（N?m ） ;n 为发动机转速（r/min ）。 发动 机的最低转速n min =600r/min,最高转速n max =4000r/min 。 装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 n =0.85 滚动阻力系数 f=0.013 空气阻力系数泌风面积 C D A=2.77m 2 解：Matlab 程序： (1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000]; Tq=-19.313+295.27*( n/1000)-165.44*(门/1000)八2+40.874*(门/1000)八3-3.8445*( n/10 00).A 4; m=3880;g=9.8; nmi n=600; nm ax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793]; nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;lf=0.218;lw1= 1.798;Iw2=3.598; Ft 仁 Tq*ig(1)*i0* nT/r; Ft2=Tq*ig(2)*i0* nT/r; 主减速器传动比 飞轮转动惯量 二前轮转动惯量 四后轮转动惯量 i 0=5.83 l f =0.218kg?m I w1=1.798kg?nf I w2=3.598kg?m ig(数据如下表) 轴距 质心至前轴距离（满载） 质心高（满载） a=1.974m hg=0.9m
matlab程序设计第三章课后习题答案
1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。 >> t=0:0.01:pi; >> x1=t; >> y1=2*x1-0.5; >> x2=sin(3*t).*cos(t); >> y2=sin(3*t).*sin(t); >> plot(x1,y1,'r-',x2,y2,'g-') >> axis([-1,2,-1.5,1]) >> hold on >> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)'); >> plot(double(s.x),double(s.y),'*'); 截图：
p366 第4题绘制极坐标曲线，并分析对曲线形状的影响。 function [ output_args ] = Untitled2( input_args ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here theta=0:0.01:2*pi; a=input('请输入a的值：'); b=input('请输入b的值：'); n=input('请输入n的值：'); rho=a*sin(b+n*theta); polar(theta,rho,'k'); end 下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。
对a的值进行改变：对比发现a只影响半径值的整倍变化 对b的值进行改变：对比发现b的值使这个圆转换了一定的角度
对n的值进行改变：对比发现当n>=2时有如下规律 1、当n为整数时，图形变为2n个花瓣状的图形 2、当n为奇数时，图形变为n个花瓣状的图形 分别让n为2、3、4、5
汽车理论大作业.
《汽车理论实习》实习报告 别克凯越1.6LE-AT 2011款 综合性能分析 学院： 专业班级： 指导老师： 实习时间： 姓名：学号：成绩： 姓名：学号：成绩： 组员任务分配： 动力性，燃油经济性—— 制动性，操纵稳定性——
目录 一、别克凯越1.6LE-AT 2011款动力性分析 (2) 1.发动机主要参数 (2) 2.参数计算 (3) 3.驱动力和行驶阻力平衡图 (6) 4.动力特性图 (7) 5.功率平衡图 (8) 二、别克凯越1.6LE-AT 2011款燃油经济性分析 (9) 1.百公里油耗估算 (9) 2.等速行驶百公里燃油消耗量计算 (12) 3.等加速行驶工况燃油消耗量的计算 (13) 4.等减速行驶工况燃油消耗量的计算 (15) 5.数据分析 (16) 三、别克凯越1.6LE-AT 2011款制动性分析 (18) 1.结构参数 (18) 2.参数分析 (18) 四、别克凯越1.6LE-AT 2011款操纵稳定性分析 (22) 1.结构参数 (23) 2.参数分析 (23)
一、别克凯越1.6LE-AT 2011款动力性分析 1.发动机主要参数 整车技术参数 动力参数
2.参数计算 （1）转矩和功率计算 根据发动机的最大功率max e P 和最大功率时的发动机转速p n ，则发动机的外特性的功率e P n --曲线可用下式估算： 23 max 12e e p p p n n n P P C C n n n ?? ??????=+- ? ? ? ???????? ? 汽油机中C1=C2=1， n 为发动机转速(r ／min), Pe max =81kw ， p n =6000r/min ； 发动机功率Pe 和转矩tq T 之间有如下关系：9549e tq P T n = 可得发动机外特性中的功率与转矩曲线：
汽车理论课后习题Matlab程序
1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选 其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。 2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n 曲线的拟合公式为 234 19.313295.27()165.44()40.874() 3.8445()1000100010001000 q n n n n T =-+-+- 式中，Tq 为发动机转矩（N?m）;n 为发动机转速（r/min ）。 发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。 装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 ηt =0.85 滚动阻力系数 f =0.013 空气阻力系数×迎风面积 C D A =2.77m 2 主减速器传动比 i 0=5.83 飞轮转动惯量 I f =0.218kg?m 2 二前轮转动惯量 I w1=1.798kg ?m 2 四后轮转动惯量 I w2=3.598kg?m 2
变速器传动比ig(数据如下表) 轴距L=3.2m 质心至前轴距离（满载）a=1.974m 质心高（满载）hg=0.9m 解：Matlab程序： (1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000]; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3 .8445*(n/1000).^4; m=3880;g=9.8;nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r;
汽车理论大作业3
汽车理论大作业3 汽车理论大作业 题目:燃油经济性计算 指导老师:侯永平 作者:徐宁 学号:081828 2011年11月 题目内容: 负荷特性曲线的拟合公式为: 234 b,B，BP，BP，BP，BP eeee01234 式中，b为燃油消耗率[g,(kw. h)]; Pe为发动机净功率(kw) 拟合式中的系数为 Q,0.299mL/s怠速油耗 (怠速转速400r/min)。 id 计算与绘制题1.3中货车的 1)汽车功率平衡图。 2)最高档与次高挡的等速百公里油耗曲线。 3)利用计算机求货车按JB3352,83规定的六工况循环行驶的百公路油耗。计算中确定燃油消耗率值b时，若发动机转速与负荷特性中给定的转速不相等，可由相邻转速的两根曲线用插值法求得。 一、绘制汽车功率平衡图 有效转速n=600—4000(r/min)。
ua=0.377rn/ii(km/h)。不同档位取不同。i g0g根据拟合公式分别求出各转速对应的转矩Tq=-19.313+295.27(n/1000)-165.44(n/1000)2+40.874 (n/1000)3-3.8445(n/1000)4(N/m)。再根据公式Pe=Ttq×n/9550(kw)求出净功率。然后依次描点就得到汽车各档功率曲线。发动机输出功率与阻力功率相平衡。 Pe=1/η(Gfu/3600+CDAua3/76140+Giu/3600+δmua/3600) aaa绘制功率平衡图时只考虑P和P,所以 fw Pe=1/η(Gfu/3600+CDAua3/76140) a 利用公式分别求出各点阻力功率，并描点画图，得到阻力曲 线。 二、求最低档和最高档的等速百公里曲线 由已知条件 234 b,B，BP，BP，BP，BPeeee01234 “计算中确定燃油消耗率值b时，若发动机转速与负荷特性中给定的转速不相等，可由相邻转速的两根曲线用插值法求得。”得: b的线性差值等同于B0，B1,B2,B3的线性差值，并依次为原理，编写了函数:
实验二--MATLAB程序的设计(含实验报告)
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1
汽车理论大作业
汽车理论大作业 计算任务书 题目： 第一题汽车动力性能的计算 确定一轻型货车的动力性能（货车装用5档变速器） 1.绘制汽车驱动力与行驶阻力平衡图。 2.求汽车最高时速，最大爬坡度及克服该爬坡度时对应的附着率。 3.绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速至70km/h的车速- 时间曲线，或者用计算机求汽车用2档起步加速至70km/h的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n曲线的拟合公式位： Tq=-19.313+295.27(n/1000)-165.44(n/1000)2+40.874(n/1000)3-3.8455(n/1000)4 式中，Tq为发动机的转矩(N*m);n位发动机的转速(r/min)。 发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min 装备质量2000kg 整车装备质量1800kg 总质量3880kg 车轮半径0.367m 传动系机械效率nT=0.85 滚动阻力系数f=0.013 空气阻力系数*迎风面积CDA=2.77m2 主减速器传动比i0=5.83 飞轮转动惯量If=0.218kg*m2 2前轮转动惯量Iw1=1.798kg*m2 4后轮转动惯量Iw2=3.598kg*m2 变速器传动比Ig（数据如下表） 1档2档3档4档5档 4档变速器 6.03 3.03 1.67 1 ---- 轴距L=3.2m 质心至前轴距离(满载) a=1.947m 质心高(满载) hg=0.9m 第二题： 第三题： 计算 1.确定一轻型货车的动力性能 1.1绘制汽车驱动力与行驶阻力平衡图 1.1.1发动机转矩 1.汽油发动机使用外特性的Tq-n曲线拟合公式为： Tq=-19.313+295.27(n/1000)-165.44(n/1000)2+40.874(n/1000)3-3.8455(n/1000)4 2.式中，Tq为发动机的转矩(N*m);n位发动机的转速(r/min)。 发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min 2.发动机的转矩：将转速带入汽油发动机使用外特性的Tq-n曲线拟合公式，取转速间隔 为100r/min，得发动机转矩表如下：
汽车理论图形MATLAB程序
功率平衡图 m=1230;g=9.8; ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75; r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2; np=6000;Pemax=83; %绘制汽车驱动力与行驶阻力平衡图 for i=1:56; n=500:100:6000; Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3); Tq(i)=9549*Pe(i)/n(i); end for j=1:5 for i=1:56 Ft(i,j)=Tq(i)*ig(j)*i0*yt/r; ua(i,j)=0.377*r*n(i)/(ig(j)*i0); Fz(i,j)=m*g*f+CD*A*(ua(i,j)^2)/21.15; end end figure plot(ua,Ft,ua,Fz); title('汽车驱动力与行驶阻力平衡图'); xlabel('ua(km/h)'); ylabel('Ft(N)'); text(20,6700,'Ft1'); text(40,4000,'Ft2'); text(50,2800,'Ft3'); text(80,2000,'Ft4'); text(100,1600, 'Ft5'); text(100,800,'Ff+Fw'); for k=1:56; n=500:100:6000; Pe(k)=Pemax*(n(k)/np+(n(k)/np)^2-(n(k)/np)^3); Tq(k)=9549*Pe(k)/n(k); Ft(k)=Tq(k)*ig(4)*i0*yt/r; ua(k)=0.377*r*n(k)/(ig(4)*i0); Fz(k)=m*g*f+CD*A*(ua(k)^2)/21.15; E(k)=abs((Ft(k)-Fz(k))); end [Emin,kmin]=min(E); Umax=ua(kmin)
吉林大学汽车理论第一次作业
汽车理论第一次作业 1-3. 确定该轻型货车的动力性能（货车可装用4挡或 5挡变速器，任选其中的一种进行整车性能计算）。 1) 绘制汽车驱动力—行驶阻力平衡图； 2) 求汽车最高车速、最大爬坡度及克服该坡度时相应的附着率； 3) 绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用 2 挡起步加速行驶至70km/h 的车速—时间曲线，求汽车用2挡起步加速行驶至70km/h 的时间。 轻型货车的有关数据： 汽油发动机使用外特性的 Tq ? n 曲线的拟合公式为 tq T =-19.313+295.27( 1000n )-165.44(1000n )2+40.8747(1000n )3-3.8445(1000 n )4 式中，q T 为发动机转矩（N·m）；n 为发动机转速（r/min ）。 发动机的最低转速600r/min ，最高转速4000r/min 。 该车的其他基本参数如表 1-3 所示。 表 1-3某轻型货车的基本参数
变速器（4挡和 5挡）的传动比g i 如表 1-4所示。 解： 1）绘制汽车驱动力与行驶阻力平衡图 。 驱动力矩 tq T =-19.313+295.27( 1000n )-165.44(1000n )2+40.8747(1000n )3-3.8445(1000 n )4 主减速器传动比 0i =5.83 车速 Ua=0.377n*r/(g i *i0) 驱动力 Ft=tq T *g i *0i *ηt /r
行驶阻力（不考虑爬坡因素，加速阻力） Ff+Fw=Gf+A C D *Ua^2/21.15 整合以上信息可得 当一档行驶时：g i =5.56时， 2.561≦Ua ≦17.074 当二挡行驶时：g i =2.769时，5.142≦ Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ?? ????++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 （mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g) 孙野 20081268 1.3（1）绘制汽车驱动力与行驶阻力平衡图 选用5挡变速器进行整车性能计算 发动机转速与汽车行驶速度之间的关系：0 377.0i i rn u g a = 发动机使用外特性的T q -n 曲线拟合公式： 4 32)1000 (8445.3)1000(874.40)1000(44.165)1000(27.295313.19n n n n T q -+-+-= 汽车驱动力：r i i T F T g q t η0= 行驶阻力：2 215 .2115.21a D a D j w i f u A C mgf dt du m u A C Gi Gf F F F F +=++ +=+++δ 用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: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 %从600~4000rpm 油拟合公式计算发动机转 距 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; 汽车理论m a t l新编作 业 Document number：BGCG-0857-BTDO-0089-2022 一、确定一轻型货车的动力性能。 1）绘制汽车驱动力与行驶阻力平衡图； 2）求汽车最高车速与最大爬坡度； 3）绘制汽车行驶加速度倒数曲线；用计算机求汽车用Ⅱ档起步加速行驶至 70km/h 所需 的加速时间。 已知数据略。（参见《汽车理论》习题第一章第3题） 解题程序如下:用Matlab语言 (1)绘制汽车驱动力与行驶阻力平衡图 m1=2000; m2=1800; mz=3880; g=; r=; CdA=; f=; nT=; ig=[ ]; i0=; If=; Iw1=; Iw2=; Iw=2*Iw1+4*Iw2; for i=1:69 n(i)=(i+11)*50; Ttq(i)=+*(n(i)/1000)*(n(i)/1000)^2+*(n(i)/1000)^*(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)=*r*n(i)/(ig(j)*i0); Fz(i,j)=CdA*ua(i,j)^2/+mz*g*f; end end plot(ua,Ft,ua,Ff,ua,Ff+Fw) title('汽车驱动力与行驶阻力平衡图'); xlabel('ua(km/h)'); ylabel('Ft(N)'); gtext('Ft1') gtext('Ft2') gtext('Ft3') gtext('Ft4') gtext('Ft5') gtext('Ff+Fw') Matlab程序设计作业 姓名: 学号: 专业: ? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10 7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合 x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >> 1.3 n=600:1:4000; r=0.367; i0=5.83; eff=0.85; f=0.013; m=3880; g=9.8; G=m*g; CdA=2.77; a=1.947; hg=0.9; L=3.2; Iw1=1.798; Iw2=3.598; Iw=Iw1+Iw2; If=0.218; Ttq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000) .^4; %驱动力行驶阻力平衡图 for ig=[5.56,2.769,1.644,1.00,0.793] Ua=0.377*r*n/ig/i0; Ft=Ttq*ig*i0*eff/r; plot(Ua,Ft); hold on; end Ff=G*f; ua=0:0.1:max(Ua); Fw=CdA*ua.^2/21.15; plot(ua,(Ff+Fw)); title('驱动力-行驶阻力平衡图'); xlabel('Ua/(km/h)');ylabel('Ft/N'); gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw') [x,y]=ginput(1); disp('汽车的最高车速');disp(x);disp('km/h'); %最大爬坡度及最大爬坡度时的附着率 Ua=0.377*r*n/5.56/i0; Ft=Ttq*5.56*i0*eff/r; Fw=CdA*Ua.^2/21.15; i=tan(asin((Ft-(Ff+Fw))/G)); disp('汽车的最大爬坡度');disp(max(i)); C=max(i)/(a/L+hg/L*max(i)); disp('克服最大爬坡度时的附着率');disp(C); %加速度倒数曲线 figure; for ig=[5.56,2.769,1.644,1.00,0.793] Ua=0.377*r*n/ig/i0; q=1+Iw/(m*r^2)+If*ig^2*i0^2*eff/(m*r^2); Ft=Ttq*ig*i0*eff/r; Fw=CdA*Ua.^2/21.15; as=(Ft-(Ff+Fw))/q/m; plot(Ua,1./as); hold on; end axis([0 98 0 10]); title('行驶加速度倒数曲线');xlabel('Ua/(km/h)');ylabel('1/a'); gtext('1/a1'),gtext('1/a2'),gtext('1/a3'),gtext('1/a4'),gtext('1/a5'); %加速时间曲线 汽车理论作业 汽73 2007010806 许四聪 6.5 解： 车身-车轮双质量系统参数：10925.05.10====μγζ、、、Hz f 。 “人体-座椅”系统参数：25.03==s s Hz f ζ、。 车速u=20m/s ，路面不平度系数3401056.2)(m n G q -?=，参考空间频率101.0-=m n 。 计算时频率步长Hz f 2.0=?，计算频率点数N=180。 1) 计算并画出幅频特性2121///z p z z q z 、、和均方根值谱)(1f G z 、)(2f G z 、 )(f G a 谱图，进一步计算q σ、1z σ、2z σ、a σ、w a 、aw L 值。 2) 改变“人体-座椅”系统参数：5.0~125.0,3~25.0==s s Hz f ζ，分析w a 、aw L 值随s f 、s ζ的变化 3) 分别改变车身-车轮双质量系统参数：Hz f 3~25.00=、5.0~125.0=ζ、 18~5.4=γ、20~5=μ。绘制2z σ、fd σ、G Fd /σ三个响应量均方根值随以上 四个系统参数变化的曲线。 解: 1）、公式 2 12 22214)1(?? ? ????+-=λζλq z ()() ()2 1 2222122121z ??? ?????+-+=ζλλζλz ()() ()2 1222222121??? ????? +-+=s s s s s z p λζλλζ 其中() () ()()()2 20202 2202 0ω/ω1μ1γω/ωζ41ω/ωμ1γ1ω/ω1Δ?? ???????? ??+-+????? ?-???? ??-+-= MATLAB在 汽车理论 中的简单应用 汽车理论 可以使用MATLAB作为科学计算工具的课程三本文档的目的在于向之前没有接触过MATLAB的同学介绍其在 汽车理论 中的简单应用,更多关于MATLAB的知识需要大家结合其他课程,通过进一步的学习和研究来掌握三 需要说明的是: 1)本文档默认读者已掌握C语言和线性代数基本知识三 2)红色字体的语句可以直接复制到MATLAB中进行运算三语句中用到的符号(引号内部分除外)均为英文状态下输入三本文档所有语句均已在MATLAB7.8.0(R2009a)版本下调试过三 3)例程编写样式纯属个人习惯,大家不必拘泥三 4)如对本文档及‘汽车理论“第6版中的附录C一一个学生的 汽车理论 课程MATLAB习题编程思路有任何问题,欢迎发送邮件一起讨论三邮件地址:wang-d07@ https://www.360docs.net/doc/4c18387244.html,三 1.MATLAB使用指南 打开MATLAB后出现的是软件的主界面三几个窗口中最重要的是Command Window,用户可以在里面输入程序,编写的函数的结果也在这里显示三 首先对MATLAB进行配置三如果是在自己的计算机上编程,建议先建一个属于自己的workspace,用来保留用户的工作环境三选择 file Save Workspace As ,在合适的目录下保存用户自己的workspace(.mat文件)三这样,以后用户编程的环境就默认为这个workspace了三 虽然可以直接在Command Window里编写程序,但这样的程序无法被保存和发布三选择 file New Blank M File ,可以新建空白的MATLAB函数文件(.m文件),在这里编写的程序可以被保存和发布三编写好的程序要经过编译之后才能执行三选择 Tools Save File and Run (也可以直接按Matlab编程与应用习题和一些参考答案
汽车理论1.3和2.7matlab编程答案
汽车理论matl新编作业
matlab程序设计作业
汽车理论课后题matlab程序
汽车理论第六章作业3
MATLAB在“汽车理论”中的简单应用