plotting after a for loop doesn't work

6 次查看（过去 30 天）
Olga Rakvag
Olga Rakvag 2022-2-20
编辑： Olga Rakvag 2022-2-20
close
clear all
clc;
CM=2;
CD=1;
D=1;
k=2*pi/153.9;
h=60;
H=(0:25);
for zeta_a=H/2
y = real(asin(CM/CD*((pi*D^2)/zeta_a)*sinh(k*h)^2/(sinh(2*k*h)+2*k*h)))
end
plot(H,y, '-b');

请先登录，再回答此问题。

采纳的回答

VBBV
VBBV 2022-2-20
编辑：VBBV 2022-2-20
close
clear all
clc;
CM=1;
CD=2;
D=1;
k=2*pi/153.9;
h=60;
H=(0:1:25);I = 1;
for zeta_a=(H/2)
y(I) = real(asin(CM/CD*((pi*D^2)/zeta_a)*sinh(k*h)^2/(sinh(2*k*h)+2*k*h)));
I = I+1;
end
plot(H,y, '-b');
  4 个评论
显示 2更早的评论
VBBV
VBBV 2022-2-20
Welcome, It appears this equation used for propellers , wind mills for finding axial thrust load ?

请先登录，再进行评论。

更多回答（1 个）

Olga Rakvag
Olga Rakvag 2022-2-20
编辑：Olga Rakvag 2022-2-20
Here comes the result with your code, thank you! :-)
close
clear all
clc;
CM=2;
CD=1;
D=1;
k=2*pi/153.9;
h=60;
H=(1:1:153.9/7);
I = 1;
for zeta_a=(H/2)
y(I) = real(asin(CM/CD*((pi*D^2)/zeta_a)*sinh(k*h)^2/(sinh(2*k*h)+2*k*h)))
I = I+1;
end
plot(H,y, '-b');
grid on
hold on
h=legend('$\omega t$');
set(h ,'Interpreter','latex','FontSize', 12);
ylabel(' Phase $\omega t$', 'Interpreter','latex')
xlabel('Wave height H(m) ','Interpreter','latex')
title('Phase \omegat for max total force F_T_m_a_x')

类别

在 Help Center 和 File Exchange 中查找有关 Partial Differential Equation Toolbox 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by