求助一个MATLAB中牛顿法的一个问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我需要做一个牛顿法迭代，即
[at1,bt1]=[at;bt]-([laa,lab;lab,lbb]^(-1))*[da;db];

laa，lab，lbb是二阶导数
但是运行时会提示出现奇异阵：Matrix is singular, close to singular or badly scaled. Results
may be inaccurate. RCOND = NaN.

laa,lab,lbb,da,db都是固定的方程，无法修改，所以不知道怎么办了，求高手帮助下解决。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB matla atlab Mat Atl Matrix close

[at1,bt1]=[at;bt]-inv([laa,lab;lab,lbb])*[da;db];

jyliao 发表于 2016-8-20 07:54
[at1,bt1]=[at;bt]-inv([laa,lab;lab,lbb])*[da;db];

万分感谢。但是Matrix is singular, close to singular or badly scaled.
Results may be inaccurate. RCOND = NaN.
运行后仍会发出警告。不知道什么原因。我可以保证函数本身没有问题。是不是牛顿法迭代有时就是不收敛呢？？

[at1,bt1]=[at;bt]-pinv([laa,lab;lab,lbb])*[da;db];

[at1,bt1]=[at;bt]-pinv([laa,lab;lab,lbb])*[da;db];

pseudo inverse

jyliao 发表于 2016-8-21 07:54
[at1,bt1]=[at;bt]-pinv([laa,lab;lab,lbb])*[da;db];

pseudo inverse

还是有警告，555555555555555555555555555555555