system of non linear equations
2 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
sepideh
am 16 Okt. 2023
Kommentiert: sepideh
am 24 Okt. 2023
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 Kommentare
Akzeptierte Antwort
Torsten
am 16 Okt. 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 Kommentare
Torsten
am 24 Okt. 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.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Systems of Nonlinear Equations 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!