The MATLAB function to solve implicit equation
174 次查看(过去 30 天)
显示 更早的评论
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 个评论
Torsten
2023-3-1
移动:John D'Errico
2023-3-1
You open one time [ and close two times ] ] . Still unbalanced parenthesis.
采纳的回答
Steven Lord
2023-3-1
@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 个评论
更多回答(1 个)
John D'Errico
2023-3-1
编辑:John D'Errico
2023-3-1
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.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Boundary Conditions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!