Nonlinear equation numerical solution

16 visualizaciones (últimos 30 días)
Tadej Cigler
Tadej Cigler el 21 de Jul. de 2020
Editada: Matt J el 21 de Jul. de 2020
Hi,
I have a following relation
C = 1.2;
D = 9420;
E = -2.0;
B = 0.1973;
f = D*sin(C*atan(B*x-E*(B*x-atan(B*x))));
and I would like to X if Y = 9250.
When I use below code:
syms x0 y0
y0 = 9250;
eqn = y0 == StiffnessMagicFormula(D,C,B,E,x0);
slip_angle = vpasolve(eqn, x0)
I obtained following result:
slip_angle =
-4.1564971021120118282594483668007e191020734
But actually result should be around 6.8 .
I tryed also with functions: fsolve, roots and find but were all unsuccessful.
BR, Tadej

Respuestas (2)

Matt J
Matt J el 21 de Jul. de 2020
Editada: Matt J el 21 de Jul. de 2020
I tryed also with functions: fsolve, roots and find but were all unsuccessful.
fsolve is overkill for such a simple 1D problem. Use fzero instead:
C = 1.2;
D = 9420;
E = -2.0;
B = 0.1973;
y0 = 9250;
f = @(x) D*sin(C*atan(B*x-E*(B*x-atan(B*x))))-y0;
x0=fzero(f,6.8)
For me, this results in
x0 =
6.9606

Fabio Freschi
Fabio Freschi el 21 de Jul. de 2020
% params
C = 1.2;
D = 9420;
E = -2.0;
B = 0.1973;
y0 = 9250;
% function
fun = @(x)D*sin(C*atan(B*x-E*(B*x-atan(B*x))))-y0;
% solution
x0 = fsolve(fun,0)

Etiquetas

Productos


Versión

R2014a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by