?? conjgrad.m
字號:
function [x,n]= conjgrad (A,b,x0)
r1 = b-A*x0;
p = r1;
n = 0;
for i=1:rank(A) %以下過程可參考算法流程
if(dot(p,A*p) < 1.0e-50) %循環結束條件
break;
end
alpha = dot(r1,r1)/dot(p,A*p);
x = x0+ alpha*p;
r2 = r1- alpha*A*p;
if(r2 < 1.0e-50) %循環結束條件
break;
end
belta = dot(r2,r2)/dot(r1,r1);
p = r2+belta*p;
n = n + 1;
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -