Help with Colebrook equation and fzero command.
11 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
I am trying to do a problem with the colebrook equation in which I have to find the roots of the equation to find "f". I am trying to use the fzero command along with a function file, but it is not working. How do I do this?
Re = linspace(10e4,10e7,NRe); NRe = 31; %Values for Reynold's Numbers
F = 2*log10(Re*sqrt(f))-0.8-(1/sqrt(f)); %Colebrook Equation
fcw = fzero(@fricfactor,1,2); %Attempt at fzero command
2 Kommentare
John D'Errico
am 7 Nov. 2014
Bearbeitet: John D'Errico
am 7 Nov. 2014
This is your second question I've seen about fzero from you. In the last question, you used a loop around fzero (although you do not yet seem to accept that no solution exists for that particular problem.) Is there some reason why a loop did not seem right here?
fzero does not solve multiple problems at once. Yes, it is true that careful use of arrayfun would succeed here, there is no need to use a complicated function syntax when a trivial loop would suffice. Learn to solve problems using the basic tools, and only when you clearly understand the language does it make sense to proceed to the complex.
Antworten (1)
Dong Young Kim
am 10 Mai 2019
Bearbeitet: Dong Young Kim
am 10 Mai 2019
It can be ran with following
% parameterized function
Ret=3000:1000:1000000;
NRe=length(Ret);
for i=1:NRe
F = @(f,Re) 2*log10(Re*sqrt(f))-0.8-(1/sqrt(f));
Re = Ret(i); % parameter
fun = @(f) F(f,Re); % function of x alone
ft(i) = fzero(fun,0.1);
end
0 Kommentare
Siehe auch
Kategorien
Mehr zu Optimization finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!