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).
: 请问有经验的大侠是如何处理这个问题的？
|
