plotting after a for loop doesn't work
显示 更早的评论
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');
采纳的回答
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 个评论
Olga Rakvag
2022-2-20
编辑：Olga Rakvag
2022-2-20
Super! It works, but I want to get a curve like
So I get max x =1.57. But I get max 0.8 don know why?
close
clear all
clc;
CM=1;
CD=2;
D=1;
k=2*pi/153.9;
h=60;
hmax=153.9/7;
H_step =[1:1:hmax];
for zeta_a=1:length(H_step/2)
y(zeta_a) = real(asin(CM/CD*((pi*D^2)/zeta_a)*sinh(k*h)^2/(sinh(2*k*h)+2*k*h)));
end
Olga Rakvag
2022-2-20
VBBV
2022-2-20
Welcome, It appears this equation used for propellers , wind mills for finding axial thrust load ?
更多回答（1 个）
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')
类别
在 帮助中心 和 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 
