连续换挡加速时间计算程序

% clear all
% clc
umax=124;
u_limited=60;
s_limited=8000;
% a=1.444;
% b=1.552;
i1=a*b*2;
i2=2*a;
i3=a*b;
i4=（a*b+a）/2;
i5=a;
ig=[i1 i2 i3 i4 i5];
k=1:5;
r=0.447;i0=5.857;eta=0.9;i0=5.857;
delta=1.045+0.04*ig;
ngk=[500 500 500 500 500];      
ngm=[9000 9000 9000 9000 9000]; 
ugk=0.377.*r.*ngk（k）./（ig（k）.*i0）; %每个档位的最低车速
ugm=0.377.*r.*ngm（k）./（ig（k）.*i0）; %每个档位的最高车速
f=0.01;
Cd=0.65;
A=3.9;
m=5000;
G=m*9.8;
Ff=G*f;
beta=0.9;%换挡转速系数
du=0.1;
u1=0;u2=0;u3=0;u4=0;u5=0;
t1=0;t2=0;t3=0;t4=0;t5=0;
tt1=0;tt2=0;tt3=0;tt4=0;tt5=0;
s1=0;s2=0;s3=0;s4=0;s5=0;
ss1=0;ss2=0;ss3=0;ss4=0;ss5=0;
for k=1:5;
    if k==1
       lu=1;
        for u=ugk（k）:du:ugm（k）;
            n=ig（k）*i0*u./r/0.377;
            for（i=1:length（n））
                if （n（i）<=2864）