Every time I run this code, I get "Error using plot Vectors must be the same length. Error in ex (line 22)

Question

Omar Ahmad AL Bloushi 2021-11-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1591999-every-time-i-run-this-code-i-get-error-using-plot-vectors-must-be-the-same-length-error-in-ex-li

回答： Kumar Pallav 2021-11-24

%prameters
L = 160
w = 1000 * 2 * pi/60
r = 60
N = 100
th0 = 0
thl = 180
dth = (thl - th0)/N
q=r/L
%calculation
for i = 0:1:N;
    th = (th0 + i * dth) * pi/180;
    phi = asin (q * sin(th));
    dphidt = w * q * cos(th)/sqrt(1 - (q * sin(th))^2);
    vd (i + l) = w * r * sin(phi + th)/cos(phi);
    wbd(i + l) = -w * r * cos(th)/(L * cos(phi));
    ad(i + l)= w * r * ((cos(phi) * cos(phi + th) * (w + dphidt) - sin(phi) * sin(phi + th) * dphidt)/cos(phi)^2);
    albd (i + l) = -w * q * ((-cos(phi) * sin(th) * w + sin(phi) * cos(th) *dphidt)/cos(phi)^2);
end
figure
plot ((th0:dth:thl) , wbd , (th0:dth:thl) , albd/100 ,'lineWidth', 2);
xlabel ('\theta ({\circ})') ;
legend ('\omega_{BD} (rad/s)','\alpha_{BD} (10^2 rad/s^2');
axis ([th0 thl - inf inf]);
figure
plot ( [th0:dth:thl] , vd/10 , [th0:dth:thl] , ad/1000 , 'LineWidth' , 2 );
xlabel ('\theta ({\circ})');
legend ('v_D (10^1 in/s)','a_D (10^3 in/s^2)');
axis ([tho thl - inf inf]);
N = 100
thl = 180
N = 100
clc