Set "not equal to" constraint for a linear system of equations solver:

4 次查看（过去 30 天）
Ash Ash
Ash Ash 2018-7-9
评论： Walter Roberson 2018-7-17
Hi, I'm looking to set a constraint for the solution to my linear system of equations
Cx=D
such that
x~=0
or
|x|>eps
I'm currently looking at lsqlin() but it seems that I am not able to do this. Would you please let me know how can I do this with lsqlin() or another function, if it is possible?
Thank you.
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2018-7-17
What is your situation such that it would not be permitted for any component to be 0, but it would be acceptable for the component to be 4.94065645841247e-324 ?

请先登录，再进行评论。

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

采纳的回答

Walter Roberson
Walter Roberson 2018-7-10
There is no direct way to do that.
You need to run your solver multiple 2^N times where N is the length of your x vector. Each position in the vector needs to be considered with a lower bound of -inf and an upper bound of -eps(realmin), and then again with a lower bound of +eps(realmin) and an upper bound of +inf . Then out of all of those 2^N solutions, you have to somehow figure out which one is "best".

更多回答（1 个）

Matt J
Matt J 2018-7-14
编辑：Matt J 2018-7-14
You may as well just find the solution x without this constraint and then add random noise.
x=x + eps.*(x==0).*rand(size(x));
You are not going to get a better solution by attempting rigorous numerical optimization, subject to |x(i)|>=eps.
Remember that these solvers are iterative approximators - you would have to iterate for ages to get within eps of the actual solution you are looking for, even assuming that floating point arithmetic errors would allow such precise convergence at all. And since you can't get within eps of the ideal solution, you may as well accept any solution that is eps or more away from the ideal solution, e.g., the solution without the |x(i)|>=eps constraint or some random perturbation of it.

类别

在 Help Center 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by