How to specify tolerance using linsolve ?

2014 11 16
1 个回答
回答已采纳
更新时间：2014 11 17
4 次查看（30 天）

0 个投票

Hi,
I have a linear system that I want to solve with Matlab using linsolve. My problem is that the coefficients' range of the system is huge. The minimum coefficient start at e-20 while the max reach 10e25. When I use linsolve, I get a warning telling me "Warning: Rank deficient, rank = 15, tol = 5.2711e+005.". When I use rank on my system, I effectively observe that my system is not full rank. But when I use rank(sys, tol) with a small tol value (~e-15), I obtain a full rank system.
So I suppose that linsolve choose a tolerance that is to high for my system and my question is if there exists a way to specify a tolerance value to linsolve in ordre to solve system with a high range of coefficients.
Thanks

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

 采纳的回答

pietro
pietro 2014-11-16

0 个投票

It is strongly suggested to normalize your design variables, contraints and so on in order to have a similar variation of the gradient and lagragian functions. Please post your problem and I will help you

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Hugo
Hugo 2014-11-16
I try to solve a system in ordrer to compute the coefficients of the solution of a Laplace problem in a domain that looks like a superposition of annulus. The linear system to solve is described here (from page 8) : http://onlinewww.jpier.org/PIER/pier92/01.09032301.pdf
In summary, I put the coefficients of all equations given in this paper in a matrix representing my system. Having infinite sum and indices this system is theoretically of infinite size but I truncate it to a certain order. My code is pretty long but is nothing else than the coefficients expressions listed in this paper (eq (38)–(49), (52)–(54), (56)–(60), (62) and (64)). I can put it here if needed.
pietro
pietro 2014-11-16
well the design variables may be scaled in this way:
(x-lb)./(ub-lb)
where x is the array with the design variables, ub the array with upper bounds and lb is the array with lower bounds
Hugo
Hugo 2014-11-17
Thanks for your help, normalization indeed solve my problem.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Systems Of Linear Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by