why is my code not working?
3 ビュー (過去 30 日間)
古いコメントを表示
Franco Belletti
2019 年 5 月 18 日
コメント済み: Walter Roberson
2019 年 5 月 19 日
I'd like to get the index (i) for wich the condition I set in the foor loop appears on the results
Fm0= linspace(0.01,0.02,100);
%ciclo iterativo
for i=1:100
for Fm= (0.25).*(((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1))))).^-1).*((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1)))))).^-1)
if Fm(i)-Fm0<=10^-6
break
end
end
end
0 件のコメント
採用された回答
Walter Roberson
2019 年 5 月 18 日
Fm0 = linspace(0.01,0.02,100);
Fm= (0.25).*(((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1))))).^-1).*((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1)))))).^-1)
i = find( Fm-Fm0 <= 10^-6, 1);
You will find that the answer is 1 -- Fm0(1) > Fm(1) so the subtraction is negative and so satisfies the condition.
You might perhaps have been thinking,
i = find( abs(Fm-Fm0) <= 10^-6, 1);
but you will find that this is empty.
2 件のコメント
Walter Roberson
2019 年 5 月 19 日
If you are looking for the value of Fm such that
Fm = (0.25).*(((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm).^-1))))).^-1).*((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm).^-1)))))).^-1)
then the value turns out to be complex, at least for positive Rey (I did not experiment with negative Rey.)
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Thermal Analysis についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!