汽车理论课后习题Matlab程序详解

汽车理论课后习题Matlab程序详解
汽车理论课后习题Matlab程序详解

陕西理工学院

杨飞华

1 / 34

1.3 确定一轻型货车的动力性能(货车可装用4挡或5挡变速器,任选

其中的一种进行整车性能计算):

1)绘制汽车驱动力与行驶阻力平衡图。

2)求汽车最高车速,最大爬坡度及克服该坡度时相应的附着率。

3)绘制汽车行驶加速度倒数曲线,用图解积分法求汽车用2档起步加速行

驶至70km/h的车速-时间曲线,或者用计算机求汽车用2档起步加速行驶至

70km/h的加速时间。

轻型货车的有关数据:

汽油发动机使用外特性的Tq-n曲线的拟合公式为23419.313295.27()165.44()40.874()3.8445()

1000100010001000qnnnn

T=?+?+?

式中,

Tq为发动机转矩(N?m);n为发动机转速(r/min)。发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min。装载质量 2000kg 整车整备质量1800kg 总质量3880kg 车轮半径0.367m 传动系机械效率εt=0.85 滚动阻力系数 f=0.013 空气阻力系数×迎风面积CDA=2.77m2 主减速器传动比i0=5.83 飞轮转动惯量 If=0.218kg?m2 二前轮转动惯量Iw1=1.798kg?m2 四后轮转动惯量 Iw2=3.598kg?m2 变速器传动比 ig(数据如下表)

Ⅰ档Ⅱ档Ⅲ档Ⅳ档Ⅴ档四档变速器 6.09 3.09 1.71 1.00 - 五档变速器 5.56 2.769 1.644 1.00 0.793 轴距 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/10

00).^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;

2 / 34

Ft4=Tq*ig(4)*i0*nT/r;

Ft5=Tq*ig(5)*i0*nT/r;

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;

ua=[0:5:120];

Ff=G*f;

Fw=CDA*ua.^2/21.15;

Fz=Ff+Fw;

plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);

title('驱动力-行驶阻力平衡图');

xlabel('ua(km/s)');

ylabel('Ft(N)');

gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw'); zoom on;

[x,y]=ginput(1);

zoom off;

disp('汽车最高车速=');disp(x);disp('km/h');

汽车最高车速=

99.3006

km/h

3 / 34

(2)求汽车最大爬坡度程序:

n=[600:10:4000];

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^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;

ua1=0.377*r*n/ig(1)/i0;

Ff=G*f;

Fw1=CDA*ua1.^2/21.15;

Fz1=Ff+Fw1;

Fi1=Ft1-Fz1;

Zoom on;

imax=100*tan(asin(max(Fi1/G)));

disp('汽车最大爬坡度=');

disp(imax);

disp('%');

汽车最大爬坡度=

35.2197%

(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序:

clear

n=[600:10:4000];

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^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;

Ft5=Tq*ig(5)*i0*nT/r;

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;

Fw1=CDA*ua1.^2/21.15;

Fw2=CDA*ua2.^2/21.15;

Fw3=CDA*ua3.^2/21.15;

4 / 34

Fw4=CDA*ua4.^2/21.15;

Fw5=CDA*ua5.^2/21.15;

Ff=G*f;

deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2);

deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2);

deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2);

deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2);

deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2);

a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;

a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;

a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;

a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;

a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;

plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);

axis([0 99 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'); a=max(a1);

af=asin(max(Ft1-Ff-Fw1)/G);

C=tan(af)/(a/L+hg*tan(af)/L);

disp('假设后轮驱动,最大爬坡度相应的附着率=');

disp(C);

假设后轮驱动,最大爬坡度相应的附着率=

0.4219

5 / 34

(4) >>clear

nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;

Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;

G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];

nmin=600;nmax=4000;

u1=0.377*r*nmin./ig/i0;

u2=0.377*r*nmax./ig/i0;

deta=0*ig;

for i=1:5

deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);

end

ua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;

Ff=G*f;

Fw=CDA*ua.^2/21.15;

for i=1:N

k=i;

if ua(i)<=u2(2)

n=ua(i)*(ig(2)*i0/r)/0.377;

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100

0)^4;

Ft=Tq*ig(2)*i0*nT/r;

inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

elseif ua(i)<=u2(3)

n=ua(i)*(ig(3)*i0/r)/0.377;

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;

Ft=Tq*ig(3)*i0*nT/r;

inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

elseif ua(i)<=u2(4)

n=ua(i)*(ig(4)*i0/r)/0.377;

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;

Ft=Tq*ig(4)*i0*nT/r;

inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

else

n=ua(i)*(ig(5)*i0/r)/0.377;

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100武汉理工大学

6 / 34

0)^4;

Ft=Tq*ig(5)*i0*nT/r;

inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

end

a=delta(1:k);

t(i)=sum(a);

end

plot(t,ua);

axis([0 80 0 100]);

title('汽车2档原地起步换挡加速时间曲线');

xlabel('时间t(s)');

ylabel('速度ua(km/h)');

>> ginput

ans =

25.8223 70.0737

25.7467 70.0737

所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s

2.7已知货车装用汽油发动机的负荷特性与万有特性。负荷特性曲线的拟合公式为:4

4

3

3

2

210eeeeP

BPBPBPBBb++++=

其中,b为燃油消耗率[g/(kW?h)];Pe为发动机净功率(kW);拟合式中的系数随转速n变化。怠速油耗smLQid/

299.0=(怠速转速400r/min)。武汉理工大学

7 / 34

计算与绘制题1.3中货车的

1)汽车功率平衡图。

2)最高档与次高档的等速百公里油耗曲线。或利用计算机求货车按

JB3352-83规定的六工况循环行驶的百公里油耗。计算中确定燃油消耗值b时,若发动机转速与负荷特性中给定的转速不相等,可由相邻转速的两根曲线用插值

法求得。

解:Matlab程序:

(1) 汽车功率平衡图程序:

clear

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;

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;

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;

Pe1=Tq.*ig(1)*i0.*ua1./(3600*r);

Pe2=Tq.*ig(2)*i0.*ua2./(3600*r);

Pe3=Tq.*ig(3)*i0.*ua3./(3600*r);

Pe4=Tq.*ig(4)*i0.*ua4./(3600*r);

Pe5=Tq.*ig(5)*i0.*ua5./(3600*r);

ua=[0:0.35:119];

Ff=G*f;

Fw=CDA*ua.^2/21.15;

Pf=Ff*ua/3600;

Pw=Fw.*ua/3600;

Pe0=(Pf+Pw)./nT;

Pe=max(Pe1);

plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pe0,ua,Pe);

axis([0 119 0 100]);

title('汽车功率平衡图');

xlabel('ua(km/h)');

ylabel('Pe(kw)');

gtext('1'),gtext('2'),gtext('3'),gtext('4'),gtext('5'),gtext('(Pf+Pw)/et'),gtext ('Pe');

8 / 34

(2)最高档与次高档的等速百公里油耗曲线程序:

clear

n=600:1:4000;

m=3880;g=9.8;

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;

n0=[815 1207 1614 2012 2603 3006 3403 3804];

B00=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7];

B10=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291]; B20=[72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113];

B30=[-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449

-0.00075215];

B40=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032

0.00028230 -0.000038568];

B0=spline(n0,B00,n);

B1=spline(n0,B10,n);

B2=spline(n0,B20,n);

B3=spline(n0,B30,n);

B4=spline(n0,B40,n);

Ff=G*f;

ua4=0.377*r*n/ig(4)/i0;

ua5=0.377*r*n/ig(5)/i0;

9 / 34

Fz4=Ff+CDA*(ua4.^2)/21.15;

Fz5=Ff+CDA*(ua5.^2)/21.15;

Pe4=Fz4.*ua4./(nT*3.6*1000);

Pe5=Fz5.*ua5./(nT*3.6*1000);

for i=1:1:3401

b4(i)=B0(i)+B1(i)*Pe4(i)+B2(i)*Pe4(i).^2+B3(i)*Pe4(i).^3+B4(i)*Pe4(i). ^4;

b5(i)=B0(i)+B1(i)*Pe5(i)+B2(i)*Pe5(i).^2+B3(i)*Pe5(i).^3+B4(i)*Pe5(i). ^4;

end

pg=7.0;

Q4=Pe4.*b4./(1.02.*ua4.*pg);

Q5=Pe5.*b5./(1.02.*ua5.*pg);

plot(ua4,Q4,ua5,Q5);

axis([0 100 10 30]);

title('最高档与次高档等速百公里油耗曲线');

xlabel('ua(km/h)');

ylabel('百公里油耗(L/100km)');

gtext('4'),gtext('5');

3.1改变1.3题中轻型货车的主减速器传动比,做出0i

为5.17、5.43、5.83、6.17、6.33

时的燃油经济性—加速时间曲线,讨论不同0i

值对汽车性能的影响。武汉理工大学

10 / 34

解:Matlab程序:

主程序:

i0=[5.17,5.43,5.83,6.17,6.33]; %输入主传动比的数据for i=1:1:5

y(i)=jiasushijian(i0(i)); %求加速时间

end

y;

for i=1:1:5

b(i)=youhao(i0(i)); %求对应i0的六工况百公里油耗end

b;

plot(b,y,'+r')

hold on

b1=linspace(b(1),b(5),100);

y1=spline(b,y,b1); %三次样条插值

plot(b1,y1); %绘制燃油经济性-加速时间曲线

title('燃油经济性—加速时间曲线');

xlabel('百公里油耗(L/100km)');

ylabel('加速时间s');

gtext('i0=5.17'),gtext('i0=5.43'),gtext('i0=5.83'),gtext('i0=6.17'),gtext('i0=6. 33');

子程序:

(1) function y=jiasushijian(i0) %求加速时间的处理函数

n1=linspace(0,5000); %先求各个档位的驱动力

nmax=4000;nmin=600;r=0.367;yita=0.85;CDA=2.77;f=0.013;G=(3880)*9.8;ig=[6.09

,3.09,1.71,1.00];%i0=5.83

for i=1:1:4 %i为档数

uamax(i)=chesu(nmax,r,ig(i),i0); %计算各个档位的最大速度与最

小速度

uamin(i)=chesu(nmin,r,ig(i),i0);

ua(i,:)=linspace(uamin(i),uamax(i),100);

n(i,:)=zhuansu(ua(i,:),r,ig(i),i0); %计算各个档位的转速范围

Ttq(i,:)=zhuanju(n(i,:)); %求出各档位的转矩范围

Ft(i,:)=qudongli(Ttq(i,:),ig(i),i0,yita,r); %求出驱动力

F(i,:)=f*G+CDA*(ua(i,:).^2)/21.15; %求出滚动阻力和空气阻

力的和

delta(i,:)=1+(1.798+3.598+0.218*(ig(i)^2)*(i0^2)*yita)/(3880*r^2); %求转动质量换算系数

a(i,:)=1./(delta(i,:).*3880./(Ft(i,:)-F(i,:))); %求出加速度

F2(i,:)=Ft(i,:)-F(i,:);

end

%下面分各个档位进行积分,求出加速时间

temp1(1,:)=ua(2,:)/3.6;

11 / 34

temp1(2,:)=1./a(2,:);

n1=1;

for j1=1:1:100

if ua(3,j1)>max(ua(2,:))&&ua(3,j1)<=70 temp2(1,n1)=ua(3,j1)/3.6;

temp2(2,n1)=1./a(3,j1);

n1=n1+1;

end

end

n2=1;

for j1=1:1:100

if ua(4,j1)>max(ua(3,:))&&ua(4,j1)<=70; temp3(1,n2)=ua(4,j1)/3.6;

temp3(2,n2)=1./a(4,j1);

n2=n2+1;

end

end

y=temp1(1,1)*temp1(2,1)+qiuji(temp1(1,:),temp1(2,:))+qiuji(temp2(1,:),temp2(2,:) )+

qiuji(temp3(1,:),temp3(2,:));

end

(2) function ua=chesu(n,r,ig,i0); %由转速计算车速

ua=0.377*r.*n/(ig*i0);

(3) function n=zhuansu(ua,r,ig,i0); %求转速

n=ig*i0.*ua./(0.377*r);

end

(4) function y=zhuanju(n); %求转矩函数

y=-19.313+295.27.*(n./1000)-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n .

/1000).^4;

(5) function y=qudongli(Ttq,ig,i0,yita,r); %求驱动力函数

y=(ig*i0*yita.*Ttq)/r;

end

(6) function p=qiuji(x0,y0) %求积分函数

n0=size(x0);

n=n0(2);

x=linspace(x0(1),x0(n),200) ;

y=spline(x0,y0,x); %插值

% figure;plot(x,y);

p=trapz(x,y) ;

end

12 / 34

(7) %求不同i0下的六工况油耗

function b=youhao(i0);

global f G CDA yita m r If Iw1 Iw2 pg B0 B1 B2 B3 B4 n %声明全局变量

ig=[6.09,3.09,1.71,1.00];r=0.367;

yita=0.85;CDA=2.77;f=0.013;%i0=5.83;

G=(3880)*9.8;If=0.218;Iw1=1.798;Iw2=3.598;m=3880; %汽车的基本参数设定n0=[815 1207 1614 2012 2603 3006 3403 3804];

B00=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7];

B10=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291]; B20=[72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113];

B30=[-5.8629 -2.0553 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449

-0.00075215];

B40=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230

-0.000038568];

n=600:1:4000;

B0=spline(n0,B00,n);

B1=spline(n0,B10,n);

B2=spline(n0,B20,n); %使用三次样条插值,保证曲线的光滑连续

B3=spline(n0,B30,n);

B4=spline(n0,B40,n);

ua4=0.377*r.*n./(i0*ig(4)); %求出发动机转速范围内对应的III、IV档车速

F4=f*G+CDA*(ua4.^2)/21.15; %求出滚动阻力和空气阻力的和

P_fw4=F4.*ua4./(yita*3.6*1000); %求出阻力功率

for i=1:1:3401 %用拟合公式求出各个燃油消耗率

b4(i)=B0(i)+B1(i)*P_fw4(i)+B2(i)*(P_fw4(i))^2+B3(i)*(P_fw4(i))^3+B4(i)*(P_fw4( i))^4;

end

pg=7.06; %汽油的重度取7.06N/L

ua4_m=[25,40,50]; %匀速阶段的车速

s_m=[50,250,250]; %每段匀速走过的距离

b4_m=spline(ua4,b4,ua4_m); %插值得出对应速度的燃油消耗率

F4_m=f*G+CDA*(ua4_m.^2)/21.15; %车速对应的阻力

P_fw4_m=F4_m.*ua4_m./(yita*3.6*1000); %发动机功率

Q4_m=P_fw4_m.*b4_m.*s_m./(102.*ua4_m.*pg) ;

汽车理论第五版-课后习题答案

第一章 汽车的动力性 1.1试说明轮胎滚动阻力的定义,产生机理和作用形式。 答:车轮滚动时,由于车轮的弹性变形、路面变形和车辙摩擦等原因所产生的阻碍汽车行驶的力称为轮胎滚动阻力。 产生机理和作用形式: (1)弹性轮胎在硬路面上滚动时,轮胎的变形是主要的,由于轮胎有内部摩擦,产生弹性迟滞损失,使轮胎变形时对它做的功不能全部回收。由于弹性迟滞,地面对车轮的法向作用力并不是前后对称的,这样形成的合力z F 并不沿车轮中心(向车轮前进方向偏移a )。如果将法向反作用力平移至与通过车轮中心的垂线重合,则有一附加的滚动阻力偶矩f z T F a =?。为克服该滚动阻力偶矩,需要在车轮中心加一推力P F 与地面切向反作用力构成一力偶矩。 (2)轮胎在松软路面上滚动时,由于车轮使地面变形下陷,在车轮前方实际形成了具有一定坡度的斜面,对车轮前进产生阻力。 (3)轮胎在松软地面滚动时,轮辙摩擦会引起附加阻力。 (4)车轮行驶在不平路面上时,引起车身振荡、减振器压缩和伸长时做功,也是滚动阻力的作用形式。 1.2滚动阻力系数与哪些因素有关? 答:滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料和气压有关。这些因素对滚动阻力系数的具体影响参考课本P9。 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课后习题解答第二章

第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型,是“双精度” 对象,还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型,而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下,以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1)

ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆,所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8(1)通过符号计算求t t y sin )(=的导数 dt dy 。(2)然后根据此结果,求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t

汽车理论第四章 汽车的制动性课后题答案

第四章 4.1 一轿车驶经有积水层的—良好路面公路,当车速为100km/h 时要进行制动。问此时有无可能出现滑水现象而丧失制动能力?轿车轮胎的胎压为179.27kPa 。 答:假设路面水层深度超过轮胎沟槽深度 估算滑水车速:i h p 34.6=μ i p 为胎压(kPa ) 代入数据得:89.84=h μkm/h 而h μμ > 故有可能出现滑水现象而失去制动能力。 4.2在第四章第三节二中.举出了CA700轿车的制动系由真空助力改为压缩空气助力后的制动试验结果。试由表中所列数据估算''2' 2 2 1ττ+的数值,以说明制动器作用时间的重要性。 提示:由表4-3的数据以及公式max 2 02292.2526.31b a a a u u s +??? ??''+'= ττ 计算' '2'22 1ττ+的数值。 可以认为制动器起作用时间的减少是缩短制动距离的主要原因。 4.3一中型货车装有前、后制动器分开的双管路制功系,其有关参数如下; 1)计算并绘制利用附着系数曲线与制动效率曲线。 2)求行驶车速30km/h ,在.0=? 80 路面上车轮不抱死的制动距离。计算时取制动系反应时间 s 02.0'2=τ,制动减速度上升时间s 02.0''2=τ。 3) 求制功系前部管路损坏时汽车的制功距离,制功系后部管路损坏时汽车的制功距离。 答案:1) 前轴利用附着系数为:g f zh b z L += β? 后轴利用附着系数为: ()g r zh a z L --= β?1 空载时: g h b L -=β?0= 413.0845 .085 .138.095.3-=-? 0??> 故空载时后轮总是先抱死。

Matlab课后习题解答

P16 Q2: 计算表达式()2 tan arccos x x -在0.25x =和0.78x π =时的函数值。 function y=jie(x) y=tan(-x.^2)*acos(x); >> jie(0.25) ans = -0.0825 >> jie(0.78*pi) ans = 0 + 0.4418i Q3:编写M 命令文件,求5010 2 1 1 1k k k k ==+ ∑∑ 的值。 a=0;b=0; for i=1:50 a=a+i*i; end for j=1:10 b=b+1/j; end c=a+b; >> c c = 4.2928e+004 P27 Q2:矩阵1234567 8 9A ????=??????,4 685563 2 2B ?? ?? =? ????? ,计算A B *,.A B *,并比较两者的区别。 >> A=[1 2 3;4 5 6;7 8 9];

>> B=[4 6 8;5 5 6;3 2 2]; >> A*B ans = 23 22 26 59 61 74 95 100 122 >> A.*B ans = 4 12 24 20 25 36 21 16 18 A*B表示A与B两矩阵相乘。 A.*B表示A与B对应元素相乘。 P34 Q2:编写一个转换成绩等级的程序,其中成绩等级转换标准为:考试分数在[] 90,100显示为优秀;分数在[) 0,60的 60,80的显示为及格;分数在[) 80,90的显示为良好;分数在[) 显示为不及格。 if x>=90 disp('优秀'); elseif x>=80 disp('良好'); elseif x>=60 disp('及格'); else disp('不及格'); end >> x=85 x = 85 良好

汽车理论课后作业答案MATLAB

汽车理论作业MA TLAB过程 汽车驱动力与阻力平衡图 加速度倒数-速度曲线图 0 10 20 30 40 50 60 70 80 90 100 u

汽车功率平衡图 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);

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

汽车理论余志生_课后习题答案(正确)剖析

qq 第一章 汽车的动力性278973104 1.1试说明轮胎滚动阻力的定义,产生机理和作用形式。 答:车轮滚动时,由于车轮的弹性变形、路面变形和车辙摩擦等原因所产生的阻碍汽车行驶的力称为轮胎滚动阻力。 产生机理和作用形式: (1)弹性轮胎在硬路面上滚动时,轮胎的变形是主要的,由于轮胎有内部摩擦,产生弹性迟滞损失,使轮胎变形时对它做的功不能全部回收。由于弹性迟滞,地面对车轮的法向作用力并不是前后对称的,这样形成的合力z F 并不沿车轮中心(向车轮前进方向偏移a )。如果将法向反作用力平移至与通过车轮中心的垂线重合,则有一附加的滚动阻力偶矩f z T F a =?。为克服该滚动阻力偶矩,需要在车轮中心加一推力P F 与地面切向反作用力构成一力偶矩。 (2)轮胎在松软路面上滚动时,由于车轮使地面变形下陷,在车轮前方实际形成了具有一定坡度的斜面,对车轮前进产生阻力。 (3)轮胎在松软地面滚动时,轮辙摩擦会引起附加阻力。 (4)车轮行驶在不平路面上时,引起车身振荡、减振器压缩和伸长时做功,也是滚动阻力的作用形式。 1.2滚动阻力系数与哪些因素有关? 答:滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料和气压有关。这些因素对滚动阻力系数的具体影响参考课本P9。 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 质心至前轴距离(满载) a=1.974m 质心高(满载) hg=0.9m 分析:本题主要考察知识点为汽车驱动力-行使阻力平衡图的应用和附着率的计算、等效坡度的概念。只要对汽

matlab课后习题及答案详解

第1章 MATLAB概论 1.1与其他计算机语言相比较,MATLAB语言突出的特点是什么? MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成? MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装? 在安装MATLAB时,安装内容由选择组件窗口中个复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上? 在MATLAB操作桌面上有五个窗口,在每个窗口的右上角有两个小按钮,一个是关闭窗口的Close按钮,一个是可以使窗口成为独立窗口的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面成为独立窗口,在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器? 在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗?如何操作? 存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外,还有什么用途? 命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别? 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径?

汽车理论第五课后习题答案

第一章汽车的动力性 试说明轮胎滚动阻力的定义,产生机理和作用形式。 定义:汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动祖力。 2)产生机理:由于轮胎内部摩擦产生弹性轮胎在硬支律路面上行驶时加载变形曲线和卸载变形曲线不重合会有能全损失,即弹性物质的迟滞损失。这种迟滞损失表现为一种阻力偶。 当车轮不滚动时,地面对丰轮的法向反作用力的分布是前后对称的;当车轮滚动时,由于弹性阻滞现象,处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力,这样,地面法向反作用力的分布前后不对称,而使他们的合力F Z 相对于法线前移一个距离a,它随弹性迟滞损 失的增大变大。即滚动时有滚动阻力偶T f = F Z.? a阻碍车轮滚动。 3]作用形式: T f = Wf,T f = T f /r 滚动阻力系数与哪些因素有关?答:滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料和气压有关。 =+ 由计算机作图有: 空车、满载时汽车动力性有无变化?为什么? 答:动力性会发生变化。因为满载时汽车的质量会增大,重心的位置也会发生改变。质量增大,滚动阻力、坡度阻力和加速阻力都会增大,加速时间会增加,最高车速降低。重心位置的改变会影响车轮附着率,从而影响最大爬坡度。 如何选择汽车发动机功率? 答:发动机功率的选择常先从保证汽车预期的最高车速来初步确定。若给出了期望的最高车速,选择的发动机功率应大体等于,但不小于以最高车速行驶时的行驶阻力功率之和。发动机的最大功率但也不宜过大,否则会因发动机负荷率偏低影响汽车的燃油经济性。

在实际工作中,还利用现有汽车统计数据初步估计汽车比功率来确定发动机应有功率。不少国家还对车辆应有的最小比功率作出规定,以保证路上行驶车辆的动力性不低于一定水平,防止某些性能差的车辆阻碍车流。 超车时该不该换入低一挡的排挡? 答:超车时排挡的选择,应该使车辆在最短的时间内加速到较高的车速,所以是否应该换入低一挡的排挡应该由汽车的加速度倒数曲线决定。如果在该车速时,汽车在此排档的加速度倒数大于低排挡时的加速度倒数,则应该换入低一档,否则不应换入低一挡。 可得、最大爬坡度为: 第二章汽车的燃油经济性 “车开得慢,油门踩得小,就一定省油” ,或者“只要发动机省油,汽车就一定省油”这两种说法对不对? 答:不对。由汽车百公里等速耗油量图,汽车一般在接近低速的中等车速时燃油消耗量最低,并不是在车速越低越省油。由汽车等速百公里油耗算式(2-1)知,汽车油耗量不仅与发动机燃油消耗率有关,而且还与车速、档位选择、汽车的状况、汽车的质量利用系数等使用因素有关,还与汽车的质量和总体尺寸、传动系、轮胎的选择有关,发动机省油时汽车不一定就省油。 试述无级变速器与汽车动力性、燃油经济性的关系。 答:①采用无级变速后,理论上克服了发动机特性曲线的缺陷,使汽车具有与等功率发动机一样的驱动功率,充分发挥了内燃机的功率,大大地改善了汽车动力性。②同时,发动机的负荷率高,用无级变速后,使发动机在最经济工况机会增多,提高了然油经济性。 用发动机的“最小燃油消耗特性”和克服行驶阻力应提供的功率曲线,确定保证发动机在最经济状况下工作的“无级变速器调节特性” 。 答:由发动机在各种转速下的负荷特性曲线的包络线即为发动机提供一定功率时的最低燃油消耗率曲线,如课本图 2-9a。利用此图可以找出发动机提供一定功率时的最经济状况(转速与负荷)。把各功率下最经济状况运转的转速与负荷率表明在外特性曲线上,便得到“最小燃油消耗特性” 。无级变速器的传动比 i'与发动机转速n及汽车行驶速度u a之间关系 如何从改进汽车底盘设计方面来提高燃油经济性? 答:汽车底盘设计应该从合理匹配传动系传动比、缩减尺寸和减轻质量、合理选择轮胎来提高燃油经济性。 为什么汽车发动机与传动系统匹配不好会影响汽车燃油经济性与动力性?试举例说明。答:在一定道路条件下和车速下,虽然发动机发出的功率相同,但传动比大时,后备功率越大,加速和爬坡能力越强,但发动机负荷率越低,燃油消耗率越高,百公里燃油消耗量就越大,传动比小时则相反。所以传动系统的设计应该综合考虑动力性和经济性因素。如最小传动比的选择,根据汽车功率平衡图可得到最高车速 u max(驱动力曲线与行驶阻力曲线的交点处车速),发动机达到最大功率时的车速为 u p。当主传动比较小时,u p>u max,汽车后备功

汽车理论第五版课后习题答案

第一章汽车的动力性 1.1 试说明轮胎滚动阻力的定义,产生机理和作用形式。 定义:汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动祖力。 2)产生机理:由于轮胎内部摩擦产生弹性轮胎在硬支律路面上行驶时加载变形曲线和卸载变形曲线不重合会有能全损失,即弹性物质的迟滞损失。这种迟滞损失表现为一种阻力偶。 当车轮不滚动时,地面对丰轮的法向反作用力的分布是前后对称的;当车轮滚动时,由于弹性阻滞现象,处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力,这样,地面法向反作用力的分布前后不对称,而使他们的合力F Z相对于法线前移一个距离a,它随弹性迟滞损失的增大变大。即滚动时有滚动阻力偶T f = F Z.?a阻碍车轮滚动。 3]作用形式: T f = Wf,T f = T f/r 1.2 滚动阻力系数与哪些因素有关?答:滚动阻力系数与路面的种类、行驶车速以及轮胎的构造、材料和气压有关。 1.3 =494.312+0.13U a2

由计算机作图有:

1.4 空车、满载时汽车动力性有无变化?为什么? 答:动力性会发生变化。因为满载时汽车的质量会增大,重心的位置也会发生改变。质量增大,滚动阻力、坡度阻力和加速阻力都会增大,加速时间会增加,最高车速降低。重心位置的改变会影响车轮附着率,从而影响最大爬坡度。 1.5 如何选择汽车发动机功率? 答:发动机功率的选择常先从保证汽车预期的最高车速来初步确定。若给出了期望的最高车速,选择的发动机功率应大体等于,但不小于以最高车速行驶时的行驶阻力功率之和。发动机的最大功率但也不宜过大,否则会因发动机负荷率偏低影响汽车的燃油经济性。 在实际工作中,还利用现有汽车统计数据初步估计汽车比功率来确定发动机应有功率。不少国家还对车辆应有的最小比功率作出规定,以保证路上行驶车辆的动力性不低于一定水平,防止某些性能差的车辆阻碍车流。 1.6 超车时该不该换入低一挡的排挡? 答:超车时排挡的选择,应该使车辆在最短的时间内加速到较高的车速,所以是否应该换入低一挡的排挡应该由汽车的加速度倒数曲线决定。如果在该车速时,汽车在此排档的加速度倒数大于低排挡时的加速度倒数,则应该换入低一档,否则不应换入低一挡。 1.7

最新MATLAB基础教程-薛山第二版-课后习题答案

《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书 班级:T1243-7 姓名:柏元强 学号:20120430724 总评成绩: 汽车工程学院 电测与汽车数字应用中心

目录 实验04051001 MATLAB语言基础 (1) 实验04051002 MATLAB科学计算及绘图 (18) 实验04051003 MATLAB综合实例编程 (31)

实验04051001 MATLAB语言基础 1实验目的 1)熟悉MATLAB的运行环境 2)掌握MATLAB的矩阵和数组的运算 3)掌握MATLAB符号表达式的创建 4)熟悉符号方程的求解 2实验内容 第二章 1.创建double的变量,并进行计算。 (1)a=87,b=190,计算 a+b、a-b、a*b。 clear,clc a=double(87); b=double(190); a+b,a-b,a*b (2)创建 uint8 类型的变量,数值与(1)中相同,进行相同的计算。 clear,clc a=uint8(87); b=uint8(190); a+b,a-b,a*b 2.计算: (1) () sin 60

(2) e3 (3) 3cos 4??π ??? clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4) 3.设2u =,3v =,计算: (1) 4 log uv v (2) () 2 2 e u v v u +- (3) clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v) 4.计算如下表达式: (1) ()() 3542i i -+ (2) () sin 28i - clear,clc (3-5*i)*(4+2*i) sin(2-8*i) 5.判断下面语句的运算结果。 (1) 4 < 20 (2) 4 <= 20

汽车理论课后习题答案(余志生版)(完全免费版)

D 第一章 1.1、试说明轮胎滚动阻力的定义、产生机理和作用形式? 答:1)定义:汽车在水平道路上等速行驶时受到的道路在行驶方向上的分力称为滚动阻力。 2)产生机理:由于轮胎内部摩擦产生弹性轮胎在硬支撑路面上行驶时加载变形曲线和卸载变形曲线不重合会有能量损失,即弹性物质的迟滞损失。这种迟滞损失表现为一种阻力偶。 当车轮不滚动时,地面对车轮的法向反作用力的分布是前后对称的;当车轮滚动时,由于弹性迟滞现象,处于压缩过程的前部点的地面法向反作用力就会大于处于压缩过程的后部点的地面法向反作用力,这样,地面法向反作用力的分布前后不对称,而使他们的合力Fa 相对于法线前移一个距离a, 它随弹性迟滞损失的增大而变大。即滚动时有滚动阻力偶矩 a F T z f = 阻碍车轮滚动。 3)作用形式:滚动阻力 fw F f = r T F f f = (f 为滚动阻力系数) 1.2、滚动阻力系数与哪些因素有关? 提示:滚动阻力系数与路面种类、行驶车速以及轮胎的构造、材料、气压等有关。 1.3、解答:1)(取四档为例) 由 u F n u n Tq Tq F t t →??? ? ?? →→→ 即 r i i T F T o g q t η= 4 32)1000 (8445.3)1000(874.40)1000(44.165)1000(27.25913.19n n n n Tq -+-+-= o g i i rn u 377.0= 行驶阻力为w f F F +: 2 15 .21a D w f U A C Gf F F +=+ 2 131.0312.494a U += 由计算机作图有

matlab课后习题及答案详解(修改版)

第1章MATLAB概论 1.1与其他计算机语言相比较,MA TLAB语言突出的特点是什么? MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MA TLAB系统由那些部分组成? MATLAB系统主要由开发环境、MA TLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装? 在安装MA TLAB时,安装内容由选择组件窗口中个复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MA TLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上? 在MA TLAB操作桌面上有五个窗口,在每个窗口的右上角有两个小按钮,一个是关闭窗口的Close按钮,一个是可以使窗口成为独立窗口的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面成为独立窗口,在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器? 在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗?如何操作? 存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外,还有什么用途? 命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别? 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成。在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上。 1.9 在MA TLAB中有几种获得帮助的途径?

matlab课后习题答案

第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff 或数值梯度gradient指令计算)(t y'曲线绘制 y',然后把)(t y和)(t 在同一张图上,观察数值求导的后果。(模拟数据从prob_data401.mat获得) 〖目的〗 ●强调:要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 (1)从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat (2)用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') (3)用gradent求导的指令(图形与上相似) dt=t(2)-t(1); yc=gradient(y)/dt;

plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已,不要进行数值求导。 ● 假若一定要计算数值导数,自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法,画出dt t t x y x ? =0 sin )(在]10 ,0[区间曲线,并计算)5.4(y 。 〖提示〗 ● 指定区间内的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间内的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =

matlab课后习题答案第一章

第1章基础准备及入门 习题1及解答 1数字1.5e2,1.5e3中的哪个与1500相同吗? 〖解答〗 1.5e3 2请指出如下5个变量名中,哪些是合法的? abcd-2xyz_33chan a变量ABCDefgh 〖解答〗 2、5是合法的。 3在MATLAB环境中,比1大的最小数是多少? 〖解答〗 1+eps 4设a=-8,运行以下三条指令,问运行结果相同吗?为什么? w1=a^(2/3) w2=(a^2)^(1/3) w3=(a^(1/3))^2 〖解答〗 (1)不同。具体如下 w1=a^(2/3)%仅求出主根 w2=(a^2)^(1/3)%求出(-8)^2的主根 w3=(a^(1/3))^2%求出(-8)主根后再平方 w1= -2.0000+3.4641i w2= 4.0000 w3= -2.0000+3.4641i (2)复数的多方根的,下面是求取全部方根的两种方法: (A)根据复数方根定义 a=-8;n=2;m=3; ma=abs(a);aa=angle(a); for k=1:m%m决定循环次数 sa(k)=(aa+2*pi*(k-1))*n/m;%计算各根的相角 end result=(ma^(2/3)).*exp(j*sa)%计算各根 result= -2.0000+3.4641i 4.0000-0.0000i-2.0000-3.4641i

(B )利用多项式02 3=-a r 求根p=[1,0,0,-a^2];r=roots(p)r =-2.0000+3.4641i -2.0000-3.4641i 4.00005指令clear,clf,clc 各有什么用处? 〖解答〗clear 清除工作空间中所有的变量。clf 清除当前图形。clc 清除命令窗口中所有显示。6以下两种说法对吗?(1)“MATLAB 进行数值的表达精度与其 指令窗中的数据显示精度相同。”(2) MATLAB 指令窗中显示的数值有效位数不超过7位。” 〖解答〗 (1)否;(2)否。 7想要在MATLAB 中产生二维数组???? ??????=987654321S ,下面哪些指令能实现目的? (A )S=[1,2,3;4,5,6;7,8;9] (B )S=[123;456;789] (C )S=[1,2,3;4,5,6;7,8,9] %整个指令在中文状态下输入 〖解答〗 前两种输入方法可以,后一种方法不行。8试为例1.3-5编写一个解题用的M 脚本文件? 〖解答〗 直接点击新文件图标,出现M 文件编辑器窗口;在该M 文件编辑器中,输入例1.3-5中的全部指令;并另存为p109.m ,便得到所需的脚本文件。

汽车理论第五版_课后习题答案

第一章 汽车的动力性 确定一轻型货车的动力性能(货车可装用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 车轮半径 传动系机械效率 ηt = 滚动阻力系数 f = 空气阻力系数×迎风面积 C D A = 主减速器传动比 i 0= 飞轮转动惯量 I f =?m 2 二前轮转动惯量 I w1=?m 2 四后轮转动惯量 I w2=?m 2 变速器传动比 ig(数据如下表) 轴距 L= 质心至前轴距离(满载) a= 质心高(满载) hg= 分析:本题主要考察知识点为汽车驱动力-行使阻力平衡图的应用和附着率的计算、等效坡度的概念。只要对汽车行使方程理解正确,本题的编程和求解都不会有太大困难。常见错误是未将车速的单位进行换算。 2)首先应明确道路的坡度的定义tan i α=。求最大爬坡度时可以对行使方程进行适当简化,可以简化的内容包括两项cos 1α≈和sin tan αα≈,简化的前提是道路坡度角不大,当坡度角较大时简化带来的误差会增大。计算时,要说明做了怎样的简化并对简化的合理性进行评估。 3)已知条件没有说明汽车的驱动情况,可以分开讨论然后判断,也可以根据常识判断轻型货车的驱动情况。 解:1)绘制汽车驱动力与行驶阻力平衡图 汽车驱动力Ft= r i i T t o g tq η 行驶阻力F f +F w +F i +F j =G ?f + 2D 21.12A C a u +G ?i+dt du m δ 发动机转速与汽车行驶速度之间的关系式为:0 g i n r 0.377 ua i ?= 由本题的已知条件,即可求得汽车驱动力和行驶阻力与车速的关系,编程即可得到汽车驱动力与行驶阻力平衡图。

MATLAB 实用教程 课后习题答案

第二章 1.计算复数3+4i与5-6i的乘积。 a=3+4i b=5-6i c=a*b 2.构建结构体Students,属性包含Name、age和Email,数据包括{’Zhang’,18,*‘Zhang@16 https://www.360docs.net/doc/751592759.html,’,’Zhang@https://www.360docs.net/doc/751592759.html,’+}、{’Wang’,21,[]}和{’Li’,[],[]},构建后读取所有Name属性值,并且修改’Zhang’的Age属性值为19。 Students(1).Age=18 Students(1).Email='zhang@https://www.360docs.net/doc/751592759.html,','zhang@https://www.360docs.net/doc/751592759.html,' Students(2).Name='Wang' Students(2).Age=21 Students(2).Email=[] Students(3).Name='Li' Students(3).Age=[] Students(3).Email=[] https://www.360docs.net/doc/751592759.html, Student(1).Age(1)=19 Student.Age 3.用满矩阵和稀疏矩阵存储方式分别构造下属矩阵: A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0] A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0] S=sparse(A) S=sparse([2,1,4],[1,2,4],[1,1,1],4,5) 4.采用向量构造符得到向量[1,5,9....,41]. A=1:4:41 5.按水平和竖直方向分别合并下述两个矩阵:A=[1 0 0;1 1 0;0 0 1],B=[2 3 4;5 6 7;8 9 10] A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B] 6.分别删除第五题两个结果的第2行。 A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B] C(2,:)=[] D(2,:)=[] 7.分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。 A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B]

相关文档
最新文档