problems using bisection method to find a maximum
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
while abs(a-b)>.001
c = (a + b)/2;
if (abs(y_int(1,find(abs(((c./x_int)-1))<.0001))-maxP)<.001)
break;
elseif (y_int(1,(find(abs(((a./x_int)-1))))<.001) > y_int(1,(find(abs(((b./x_int)-1))))<.001))
b=c;
else
a=c;
end
end
So, I want a and b to converge upon a max, and yet they always seem to go to the minimum in the domain i specify. I'm using a spline, not a specific equation. Just to clarify the if and elseif lines. The if line says, if the power that I get at point c is really close to the maximum Power then stop. The elseif line says, if the power at A is greater than the power at B, then make b equal to c, otherwise make a=c. Can anybody see whats wrong?
1 commentaire
Walter Roberson
le 28 Sep 2011
I recommend you calculate
t = abs(((c./x_int)-1));
right after you assign c in the loop, and that you then replace the occurrences of that expression with the variable.
I also recommend you switch to logical indexing instead of using find().
I think the code will come out cleaner, more readable, and faster.
Réponses (1)
Voir également
Catégories
En savoir plus sur Logical dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!