S*******t 发帖数: 97 | 1 在做newton-raphon迭代的时候经常要inverse matrix,
经常是进行到4,5个循环的时候matlab就报warning,
因为matrix基本上已经是singular了(det = 0).
请问有经验的大侠是如何处理这个问题的? |
c*******v 发帖数: 2599 | 2 你是要求解非线性ODE么?
牛顿法主要靠瞎猜初始值和步长。
没有普遍的理论指导。
【在 S*******t 的大作中提到】 : 在做newton-raphon迭代的时候经常要inverse matrix, : 经常是进行到4,5个循环的时候matlab就报warning, : 因为matrix基本上已经是singular了(det = 0). : 请问有经验的大侠是如何处理这个问题的?
|
S*******t 发帖数: 97 | 3 不是,就是普通的matrix reverse.如何才能
避免多步迭代以后变成singular
【在 c*******v 的大作中提到】 : 你是要求解非线性ODE么? : 牛顿法主要靠瞎猜初始值和步长。 : 没有普遍的理论指导。
|
A*******r 发帖数: 768 | |
s*****r 发帖数: 183 | 5 有很多方法处理这类问题:
一种方法是在对角线上加上一个常数,使得矩阵正定。
另外也要注意,即使没有不满秩的情况,也可能出现Hessian不正定的情况,也应该处
理。
【在 S*******t 的大作中提到】 : 在做newton-raphon迭代的时候经常要inverse matrix, : 经常是进行到4,5个循环的时候matlab就报warning, : 因为matrix基本上已经是singular了(det = 0). : 请问有经验的大侠是如何处理这个问题的?
|