Octaveでの共役勾配法の記述例
出典: フリー百科事典『ウィキペディア(Wikipedia)』 (2021/11/29 01:56 UTC 版)
「共役勾配法」の記事における「Octaveでの共役勾配法の記述例」の解説
Gnu Octaveで書くと以下のようになる。 function [x] = conjgrad(A,b,x0) r = b - A*x0; w = -r; z = A*w; a = (r'*w)/(w'*z); x = x0 + a*w; B = 0; for i = 1:size(A)(1); r = r - a*z; if( norm(r) < 1e-10 ) break; endif B = (r'*z)/(w'*z); w = -r + B*w; z = A*w; a = (r'*w)/(w'*z); x = x + a*w; end end
※この「Octaveでの共役勾配法の記述例」の解説は、「共役勾配法」の解説の一部です。
「Octaveでの共役勾配法の記述例」を含む「共役勾配法」の記事については、「共役勾配法」の概要を参照ください。
- Octaveでの共役勾配法の記述例のページへのリンク