共轭梯度法,请帮我分析一下那里出错了!

shinemin

clear
x0=[1:1];
A=[2 0;0 200]
[M,x,k]=FR(x0,A)
________________________________________________________________________
function [M,x,k]=FR(x0,A)
x=x0;
g=A*x;
s=g(1)^2+g(2)^2;
k=0;
while s~=0
    if k==0
        p=-g;
    else
        bats=s/t;
        p=-g+bats*p;
    end
    alpha=-(p.')*A*x/((p.')*A*p);
    x=x+alpha*p;
    g=A*x;
    t=s;
    s=g(1)^2+g(2)^2;
    k=k+1;
    M(:,k)=x;
end

cesario

try cgs/pcg
分别为共轭梯度方法以及预处理的CG方法
edit cgs可以看到他的源码
你的算法。。。
有点乱的感觉，推荐你看看数值计算原理里相关的章节
使用CG方法，存储量p_k(搜索方向),A*p_k(A共轭向量),r_k(剩余向量),x_k(对应的结果)。另外，用内积函数dot比较简洁明
虽然CG方法可以是直接方法，但是在一般应用里，大多当作迭代方法，用norm(residual,2)2-范数判断收敛比较常用