共轭梯度法解线性方程组(Matlab程序)

%-------共轭梯度法解线性方程组-----------
%---Conjugate Gradient method-------
%参考教材《数值分析》李乃成&梅立泉,科学出版社2011
clear;clc;
% A=[10,-1,-2;-1,10,-2;-1,-1,5];
% b=[72,83,42]';
A=[2 0 1;0 1 0;1 0 2];
b=[3 1 3]';
N=length(b);    %解向量的维数
fprintf('库函数计算结果:');
x=inv(A)*b      %库函数计算结果

x=zeros(N,1);%迭代近似向量
eps=0.0000001;%精度
r=b-A*x;d=r;
for k=0:N-1
    fprintf('第%d次迭代:',k+1);
    a=(norm(r)^2)/(d'*A*d)
    x=x+a*d
    rr=b-A*x;    %rr=r(k+1)
    if (norm(rr)<=eps)||(k==N-1)
        break;
    end
    B=(norm(rr)^2)/(norm(r)^2);
    d=rr+B*d;
    r=rr;
end
   
   
   


 

阅读更多
文章标签: matlab 出版
个人分类: Matlab 数值分析
上一篇逐次超松弛迭代法解线性方程组(Matlab程序)
下一篇平方根法和改进的平方根法解线性方程组(Matlab程序)
想对作者说点什么? 我来说一句

共轭梯度法matlab程序

2009年04月14日 304B 下载

没有更多推荐了,返回首页

关闭
关闭