📄 cg.m
字号:
function x=cg(A,b)
x=0*b;
r=b;
p=b;
normr0=norm(r);
while norm(r)>1e-3*normr0
Ap=A*p;
rdotr=r'*r;
alpha=rdotr/(p'*Ap);
x=x+alpha*p;
r=r-alpha*Ap;
beta=(r'*r)/rdotr;
p=r+beta*p;
end
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -