为什么matlab不输出了,求大神帮我看一下这个代码为什么输出不了

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

%% 驱动力和行驶阻力平衡图 (滚阻+风阻状态下)

n_test=[500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 7000 8000 8500];

T_test=[279.6 280 279.8 280 279.9 235.2 200.4 171.7 146.1 124.7 106.8 90.6 65.9 46.2 40.5];

p=polyfit(n_test,T_test,4);

n=[500:8500];

Ttq=polyval(p,n);

r=0.32; M=1970; g=9.8; CD=0.38; A=2.02;i0=9.28;

eta=0.9; %Ⅰ挡传动系统效率，Ⅱ挡传动系统效率

f=0.0138;

for i=1:length(i0);

ua(i,:)=0.377*r*n/i0;

end

%计算滚动阻力系数 f

%良好道路上货车轮胎的滚阻系数,是随车速变化的函数,ua(i,:)中的 i,;表示矩阵 ua 的第 i 行所有元素

%计算滚动阻力 Ff

for i=1:length(i0);

Ff(i,:)=M*g*f/1000;

end

%计算空气阻力 Fw

for i=1:length(i0);

Fw(i,:)=CD*A*(ua(i,:).^2)/21.15/1000;

end

%计算行驶阻力 F

for i=1:length(i0);

F(i,:)=Ff(i,:)+Fw(i,:);

end

%计算汽车驱动力 Ft

for i=1:length(i0);

Ft(i,:)=Ttq.*i0*eta/r/1000; %Ⅰ挡传动系统效率下的驱动力 Ft，Ⅱ挡传动系统效率下的驱动力 Ft，当驱动力 Ft 与行驶阻力 F 平衡，车辆达到最高时速 98.531km/h(高挡)，此时对应电机转速 n=4418rpm

end

figure(4)

hold on

for i=1:length(i0);

plot(ua(i,:), Ft(i,:),color(i)); % 与车辆挡位有关，分为Ⅰ挡驱动力和Ⅱ挡驱动力

plot(ua(i,:), F(i,:),'k'); % 车辆行驶阻力，与车速对应

plot(ua(i,:), Ff(i,:),'b'); % 车辆滚动阻力，与车速有关

hold on

grid on

box on %坐标添加边框

xlabel('车速(km/h)');

ylabel('驱动力、行驶阻力(kN)');

legend('驱动力 Ft','行驶阻力 Ff+Fw','滚动阻力 Ff',1);

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