我正在嘗試分析 different iterative sub-solvers of the Newton equation system, using a Newton-Fischer reformulation of the LCP(linear complementarity problem) 的結果 . 到目前為止，我已經實現了精確求解器 - Gauss-Siedel和牛頓修改方法，該方法使用 bicg matlab 作為J * h = -p等式的子求解器(其中 J 是雅可比， p 是Fischer函數的值)和 h 是我的實現步驟) .

我實現bicg子解析器和精確求解器的部分代碼：

if itt

% use iterative solver for newton eq

while ~all(fischer(x, A*x+b) == 0) & its < max_it

% compute the Jacobian

J = eval_jacobian(A, b, x);

% compute the value of the Fischer function

p = fischer(x, A*x + b);

% the natural merit function for convergence measure

residual(its) = .5*(p'*p);

% the newton eq, solve J*h = -p

h = bicg(J, -p, eps, s_its);

% update the solution vector

x = x + h;

% increment the iteration counter

its = its + 1;

end

else

% the exact solver for newton equation

while ~all(fischer(x, A*x+b) == 0) & its < max_it

% compute the Jacobian

J = eval_jacobian(A, b, x);

% compute the value of the Fischer function

p = fischer(x, A*x + b);

% the natural merit function for convergence measure

residual(its) = .5*(p'*p);

% the newton eq, solve J*h = -p

h = - J/p;

% update the solution vector

x = x + h;

% increment the iteration counter

its = its + 1;

end

那么，我的問題是你會使用其他迭代子解算器嗎？如果在matlab中沒有它們的功能，我不介意為它們實現代碼 . 我意識到這更像是一個理論問題 . 謝謝 .

《新程序員》：云原生和全面數字化實踐50位技術專家共同創作，文字、視頻、音頻交互閱讀

總結

以上是生活随笔為你收集整理的bicg matlab,除了bicg之外，还可以使用哪些迭代子解算器来解决MATLAB中的牛顿方程？...的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。