Matlab观察预处理共轭梯度法(PCG)效果

题目

题目来自高等数值分析作业题：

程序

随机给b的取值，可以改成指定，由于这几种预处理方法条件数改变不多，b的随机性会导致结果相差较大。

clc
clear
n=100; %生成矩阵阶数
A = hilb(n); %生成hilb矩阵
con_A=cond(A); %计算条件数b=rand(n,1); %方式一：生成随机矩阵作为y% for i=1:n
%     b(i,1)=1; %方式二：指定y
% enddisp(['无预处理的cg结果，条件数为' , num2str(con_A)]);
x = pcg(A,b,1e-6,10000); %cg计算收敛性%提取对角线矩阵
A