汽车理论课后习题1.3

for i=1:5

ual(i,1)=(600*0.377*r)/(ig(i)*i0); if (4000*0.377*r)/(ig(i)*i0)

ual(i,2)=umax; end end

%从二档到五档按照A方法选取合适加速度 for u=1:98 if(u

else if(u<=ual(2,2)) %当在二档ua内的加速度倒数 n=u*ig(2)*i0/(0.377*r);

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4; Ft=Tq*ig(2)*i0*eta/r; Fz=Cd*u^2/21.15+m*g*f; a(u)=(th(2)*m)/(Ft-Fz);

else if(u<=ual(3,2)) %在二档最大速度到三档最大速度内加速度倒数 n=u*ig(3)*i0/(0.377*r);

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4; Ft=Tq*ig(3)*i0*eta/r; Fz=Cd*u^2/21.15+m*g*f; a(u)=(th(3)*m)/(Ft-Fz);

else if(u<=ual(4,2)) %在三档最大速度到四档最大速度内加速度倒数 n=u*ig(4)*i0/(0.377*r);

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4; Ft=Tq*ig(4)*i0*eta/r; Fz=Cd*u^2/21.15+m*g*f; a(u)=(th(4)*m)/(Ft-Fz);

else if(u<=ual(5,2)) %在四档最大速度到五档最大速度内加速度倒数 n=u*ig(5)*i0/(0.377*r);

Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4; Ft=Tq*ig(5)*i0*eta/r; Fz=Cd*u^2/21.15+m*g*f; a(u)=(th(5)*m)/(Ft-Fz); end end end end end end

ux=1:ual(2,1);

%给低速无意义阶段的速度赋予5km/h时同样的加速度

[x,y]=size(ux); for i=1:y a(i)=a(y+1) end t=0; u=1:98; for i=1:98

t=t+1/3.6*a(i);%图解积分法求解时间 ts(i)=t; end

plot(ts,u)

title('速度时间图') ylabel('ua(km/h)') xlabel('t(s)')

text(33,70,'\\leftarrow 70km/h加速时间')

>> ts(70) ans =27.4830