动力性分析源程序

梦遗落在地上

以前做过整车的动力性分析，参考了一些程序，然后自己改进了一下，但是水平有限，程序还不完善。
希望高手能进行指点一下
MATLAB语言，数据可修改。
n0=[9026 8481 8010 7509 7026 6469 6019 5491 5023 4521 4000 3500 3000 2500 2000 1400];
n1=n0/1000;
T1=[8.34 8.79 8.93 9.04 9.2 9.13 8.9 8.46 8.1 7.68 7.45 7.23 7.08 6.64 6.40 6.04];
p=polyfit(n1,T1,4);
n=(1400:50:9200);
n2=n/1000;
Tq=polyval(p,n2);
figure(1)
plot(n0,T1,'o',n,Tq)
axis([1400 9500 5 12]);
title('125ml发动机转矩外特性曲线')
xlabel('n   (r/min)'); ylabel('T_q   (N.m)')
i1=4.055;
ig=[2.846 2.063 1.400 1.130 0.923];
i2=3.455;
ik=i1*ig*i2;
r=0.248;
G=480*10;
xl=0.9;
f=0.0105;
CD=0.7;
A=1.25;
for i=1:5
ua=0.377*n*r/ik(i);
Ft=Tq*ik(i)*xl/r;
Ff=G*f;
Fw=2.15*CD*A.*ua.*ua/21.15;
figure(2)
if i==1
    alpha=max(asin((Ft-(Ff+Fw))/G));
    sprintf('最大爬坡度:')
    I=tan(alpha)
end
if  i==5   
    Ft==Ff+Fw;
      sprintf('最高车速:')
     Umax=max(ua)
  end
plot(ua,Ft,ua,Ff+Fw),hold on,
xlabel('u_a   (km/h)'); ylabel('F   (N)')
axis([0 100 0 1600]);
end
gtext('F_t_1 '),gtext('F_t_2 '),gtext('F_t_3 ') ,gtext('F_t_4 '),gtext('F_t_5 '),gtext('F_f+F_w ')
title('正三轮摩托车车驱动力-行驶阻力平衡图')
for i=1:5
ua=0.377*n*r/ik(i);
Ft=Tq*ik(i)*xl/r;
Ff=G*f;
Fw=2.25*CD*A.*ua.*ua/21.15;
Fi=Ft-Ff-Fw;
a=asin(Fi/G);
i=tan(a)*100;
figure(3)
plot(ua,i),hold on
end
%axis([0 70 0 30]);
xlabel('u_a   (km/h)'); ylabel('i(%)')
title('正三轮摩托车的爬坡度图')
gtext('i_1 '),gtext('i_2 '),gtext('i_3 ') ,gtext('i_4 '),gtext('i_5 '),
for i=1:5
ua=0.377*n*r/ik(i);
s=1.074+0.0000126*(ik(i)/r);
Ft=Tq*ik(i)*xl/r;
Ff=G*f;
Fw=2.25*CD*A.*ua.*ua/21.15;
Fj=Ft-Ff-Fw;
a=(Fj*10/(G*s));
figure(4)
plot(ua,a),hold on,
figure(4)
    plot(ua,a),hold on
   figure(5)
   j0=length(a);
   for j=1:j0
       if a(j)>=0.00001
   plot(ua(j),1/a(j)),hold on
end
end
end
axis([0 55 0 8])
xlabel('u_a   (km/h)'); ylabel('1/j   (s^2/m)')
title('正三轮摩托车的加速度倒数图')

wendao

谢谢楼主，参考学习一下。