system of non linear equations
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello
I am dealing with a system of two non-linear equations and I tried to solve it by 'fsolve', but the solver stopped and could not fine the answers.
Would you please help me?
I attached the files
0 commentaires
Réponse acceptée
Torsten
le 16 Oct 2023
%solving stsyem of nonlinear equation 3.31 and 3.32 Woods paper with fsolve
% eta_l= x1 , eta_u= x2
r=2;
v=20;
F=13/9;
fun= @(x)rooted(x,r,v,F);
x0 = [0 0];
x = fsolve(fun,x0)
function f = rooted (x,r,v,F)
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
% F<Fd Fd= 2.6 I put F=1.5
f(1)= (r*v-1)/(v*(r-1)) * ((r*v*x(1)^2-x(2)^2)/(r*v-1))^(3/2) - ...
(v-1)/(v*(r-1)) * ((v*x(1)^2-x(2)^2)/(v-1))^(3/2) - 9 ;
f(2)= x(2)^3/(r*v) + (v-1)/(v*(r-1)) * ((v*x(1)^2-x(2)^2)/(v-1))^(3/2) - ...
(r*v-1)/(r*v*(r-1)) * ((r*v*x(1)^2-x(2)^2)/(r*v-1))^(3/2) - 9*F;
end
7 commentaires
Torsten
le 24 Oct 2023
I don't know what the reason for non-convergence is. The equations seem to be correct - and concerning parameters r, v and F and initial values for x - I don't know if they make sense.
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Systems of Nonlinear Equations 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!