using "solve" and getting rid of unknown parameters

6 次查看（过去 30 天）
Del
Del 2013-1-14
I am trying to solve a system of equations whose answer will depend on arbitrary values.
For example, consider the following system of 2 equations and 3 unknowns:
syms x1 x2 x3
E=[x1+2*x2+x3;x3-x1+x2];
a=solve(E==0,x1,x2,x3)
[a.x1 a.x2 a.x3]
this gives me an answer that depends on an arbitrary value z.
I would like MATLAB to directly give me an answer (with no arbitrary values). So for instance, I would like to get an answer where the x1 value is positive. Is there a way to do that?
I am using "solve" in a subroutine, and after that, I would like to test the answer to see if it is a non negative array. But whenever the code reaches the answer, it tells me there is an error (obviously) because of the unknown parameter.

请先登录，再回答此问题。

回答（2 个）

Azzi Abdelmalek
Azzi Abdelmalek 2013-1-14
编辑：Azzi Abdelmalek 2013-1-14
You have 2 equations and 3 unknowns variables, which means there are infinite solutions. It's like your are solving
E=[x1+2*x2+z;z-x1+x2];
b=solve(E==0,x1,x2)
  6 个评论
显示 4更早的评论
Del
Del 2013-1-14
编辑：Del 2013-1-14
sure, my question, is what is the correct way in MATLAB " of adding , clearly, the conditions that allows to choose some of them."
Azzi Abdelmalek
Azzi Abdelmalek 2013-1-14
编辑：Azzi Abdelmalek 2013-1-14
There are your conditions, if for example you said that the first solution must be positive. The below code will choose one solution.
syms x1 x2 x3
E=[x1+2*x2+x3;x3-x1+x2];
a=solve(E==0,x1,x2,x3)
n= size(fieldnames(a),1)
z=solve(a.x1==1)
for k=1:n
sol(k)=eval(a.(sprintf('x%d',k)));
end

请先登录，再进行评论。

Walter Roberson
Walter Roberson 2013-1-14
You might be able to use numeric::solve
That routine is not exposed at the MATLAB level so you would have to feval() or evalin() the symengine.

类别

在 Help Center 和 File Exchange 中查找有关 Newton-Raphson Method 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by