拟牛顿法求解轴承静力学模型

kerb

你的方程最好写出来，如果是非线性的方程组，你在牛顿法中，雅克比矩阵往往接近奇异，高斯消元法往往不能进行下去，即便能进行下去，求解结果也未必正确，看你的程序里面好象没有雅克比矩阵，你说的拟牛顿法也不像，里面也没看见海塞矩阵。

牛顿法是这样的：假设你的方程组用F表示，S = X_{n+1}-X_n;假设J是雅克比矩阵，迭代是这样：S =- (J^tJ)^{-1}J^tF(X_n);X_{n+1} = X_n +Ｓ