can someone help me with this code? there is something wrong and I don't know what it's.

조회 수: 1 (최근 30일)
%prameters
L = 160
w = 1000 * 2 * pi/60
r = 60
N = 100
th0 = 0
thl = 180
dth = (thl - th0)/N
q=r/L
l = 1
%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

답변 (1개)

KALYAN ACHARJYA
KALYAN ACHARJYA 2021년 11월 21일
편집: KALYAN ACHARJYA 2021년 11월 21일
Replace all lines
axis ([th0 thl -inf inf]);
More
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 ([th0 thl -inf inf]);
  댓글 수: 2
Omar Ahmad AL Bloushi
Omar Ahmad AL Bloushi 2021년 11월 21일
Error using plot
Vectors must be the same length.
still not working
KALYAN ACHARJYA
KALYAN ACHARJYA 2021년 11월 21일
편집: KALYAN ACHARJYA 2021년 11월 21일
%prameters
L = 160;
w = 1000 * 2 * pi/60;
r = 60;
N = 100;
th0 = 0;
thl = 180;
dth = (thl - th0)/N;
q=r/L;
l = 1;
%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 ([th0 thl -inf inf]);
It's working, please do the memory allocation in arrays.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 2-D and 3-D Plots에 대해 자세히 알아보기

태그

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by