The MATLAB function to solve implicit equation
215 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Teja Reddy
am 1 Mär. 2023
Beantwortet: Steven Lord
am 1 Mär. 2023
Hi all,
Thanks once again to answer my questions every time.
I would like to request you tell me the best way to solve an Implicit equation in MATLAB.
Thank you very much
6 Kommentare
Torsten
am 1 Mär. 2023
Verschoben: John D'Errico
am 1 Mär. 2023
You open one time [ and close two times ] ] . Still unbalanced parenthesis.
Akzeptierte Antwort
Steven Lord
am 1 Mär. 2023
@John D'Errico posted how to solve this equation symbolically, if you don't have values for the various symbols used in the code. If you do have values for the symbols (except for a0, the symbol for which you're trying to solve) write a function to evaluate f(a0) and use fzero.
m = 1;
a = 2;
k1 = 3;
k2 = 4;
v1 = 5;
v2 = 6;
EQ = @(a0) [m^2 /( 3*(a+a0)^2)] * (2 * k1 + 3 * k2) - (a0 * (2/v1 + 5/v2));
y = fzero(EQ, 7)
check = EQ(y) % should be small
I changed your EQ slightly. Instead of being [one expression] == [a second expression] I made it f(a0) = [one expression] - [a second expression]. fzero solves for where that f(a0) is equal to 0 which is equivalent to your formulation.
0 Kommentare
Weitere Antworten (1)
John D'Errico
am 1 Mär. 2023
Bearbeitet: John D'Errico
am 1 Mär. 2023
Assuming your question is how to solve the equation
[m^2 /( 3(a+a0)^2)] * (2 * k1 + 3 * k2) = a0 * (2/v1 + 5/v2)
I removed the second closing ], since it seems to have been the spurious one. Even if not, the solution is still little different. Just use solve. This is just a simple polynomial equation at heart. And, you NEED to use * to multiply. MATLAB does not allow implicit multiplication. That is, 3(a+a0) is NOT seen as 3*(a + a0) in MATLAB. It is just a syntax error.
syms m a a0 k1 k2 v1 v2
EQ = [m^2 /( 3*(a+a0)^2)] * (2 * k1 + 3 * k2) == a0 * (2/v1 + 5/v2)
a0sol = solve(EQ,a0,'maxdegree',3)
Yes, the rsult is messy. But what do you expect from what is effectively a cubic polynomial? Not all problems have a nice simple solution. In fact, most implicit equations are written in that form because they do not have a simple solution, and worse yet, most implicit equations will have no algebraic solution at all.
Siehe auch
Kategorien
Mehr zu Boundary Conditions 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!