[remnant]

你的Fr大小是1x121，而Ft1,Ft2,Ft3,Ft4是1x3410的。

D1=f+1/(g*m)*(Ft1-Fr); %将du/dt=1/(L*m)[Ft-(Ff+Fw)]代入D=Y+L/g*(du/dt)
D2=f+1/(g*m)*(Ft2-Fr); %坡度i忽略，则道路阻力系数Y=f
D3=f+1/(g*m)*(Ft3-Fr);
D4=f+1/(g*m)*(Ft4-Fr);

Ft1,2,3,4 - Fr就是错误来源，数组长度不一致。