用Matlab作出汽车驱动力图

根据汽油发动机使用外特性的$T_q-n$曲线的拟合公式：
$$T_q=-19.313+295.27(\frac{n}{1000})-165.44(\frac{n}{1000}{)}^2+40.874(\frac{n}{1000}{)}^3-3.8445(\frac{n}{1000}{)}^4$$再利用：$$u_a=0.377\frac{rn}{i_g{i}_0}和F_t=\frac{T_q{i}_g{i}_0{\eta }_T}{r}$$
可作出汽车的驱动力图。

Matlab代码：

function plot_F_u
  u=zeros(680,1);
  F=zeros(680,1);
  i=1;
  k=1;
for k=1:5
  for n=600:5:4000;
  u(i)=ua(n,k);
  F(i)=Ft(Tq(n),k);
  i=i+1;
  end
  i=1;
  hold on
  plot(u,F);
end
function[Ft]=Ft(T,i)
  r=0.367;  %车轮半径
  ig=[5.56,2.769,1.644,1.0,0.793];%不同档位的变速器传动比
  i0=5.83;%主减速器传动比
  eta=0.85;%传动系机械效率
  Ft=T*ig(i)*i0*eta/(1000*r);
end
function[uan]=ua(n,i)
  r=0.367;%车轮半径
  ig=[5.56,2.769,1.644,1.0,0.793];%不同档位的变速器传动比
  i0=5.83;%主减速器传动比
  uan=0.377*(r*n)/(ig(i)*i0);
end
function[Tqn]=Tq(n)
  Tqn=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;
end
end

汽车驱动力图：