gongexieliangfa.m

来自「解线性方程组的一种非线性迭代法：共轭斜量法。适用于方程组系数矩阵实对称正定的情况」· M 代码 · 共 25 行

function [X]=GongEXieLiangFa(A,b,x0,epsilon)
x(:,1)=x0;
r(:,1)=b-A*x(:,1);p(:,1)=r(:,1);
k=1;
while 0<1
alpha(k)=r(:,k)'*r(:,k)/((A*r(:,k))'*r(:,k));
x(:,k+1)=x(:,k)+alpha(k)*p(:,k);
r(:,k+1)=r(:,k)-alpha(k)*A*p(:,k);
if norm(r(:,k+1))<epsilon
X=x(:,k+1);
return
end
beta(k)=r(:,k+1)'*r(:,k+1)/(r(:,k)'*r(:,k));
p(:,k+1)=r(:,k+1)+beta(k)*p(:,k);
k=k+1;
end
end
% clc
% clear
%conlum 
% A=[10 -1 -2;-1 10 -2;-1 -1 5];
% b=[7.2;8.3;4.2];
% x0=[0;0;0];
% epsilon=0.0001;
% [X]=GongEXianLiangFa(A,b,x0,epsilon)