-
大小: 304B文件類(lèi)型: .m金幣: 1下載: 0 次發(fā)布日期: 2021-06-09
- 語(yǔ)言: Matlab
- 標(biāo)簽:
資源簡(jiǎn)介
共軛梯度法matlab程序
function x = cg(A,b)
tol=1e-10;
r = b + A*b;
w = -r;
z = A*w;
s = w'*z;
t = (r'*w)/s;
x = -b + t*w;
for k = 1:numel(b);
r = r - t*z;
if( norm(r) < tol )
return;
end
B = (r'*z)/s;
w = -r + B*w;
z = A*w;
s = w'*z;
t = (r'*w)/s;
x = x + t*w;
end
代碼片段和文件信息
function?x?=?cg(Ab)
tol=1e-10;
r?=?b?+?A*b;
w?=?-r;
z?=?A*w;
s?=?w‘*z;
t?=?(r‘*w)/s;
x?=?-b?+?t*w;
for?k?=?1:numel(b);
????r?=?r?-?t*z;
????i
評(píng)論
共有 條評(píng)論