用共轭梯度法解方程,用Jacobi方法求矩阵的全部特征值和特征向量.docVIP

用共轭梯度法解方程，用Jacobi方法求矩阵的全部特征值和特征向量 二、代码 clear %输入矩阵阶数 n=input(矩阵阶数n=); A=zeros(n,n); b=zeros(n,1); for i=1:n/2 b(2*i-1,1)=5; b(2*i,1)=6; end for i=2:n-1 A(i,i)=4; A(i,i-1)=1; A(i,i+1)=1; end A(1,1)=4; A(n,n)=4; A(n,n-1)=1; X=zeros(n,1); for i=1:n X(i,1)=1; end %用共轭梯度法求解方程 fprintf(方程的精确解\n); X fprintf(用共轭法求解方程\n); x=cg(A,b) %用方法求解方程的特征值和特征向量 fprintf(用Jacobi方法求解方程的特征值和特征向量\n); [D,V]=tezhengJaco(A) 三、数值结果 baogaoer 矩阵阶数n=10 方程精确解 X = 1 1 1 1 1 1 1 1 1 1 用共轭梯度法求解方程 k = 100 x = 1.2500 1.0104 0.7085 1.1555 0.6694 1.1668 0.6634 1.1794 0.6188 1.3453 用方法Jacobi求解矩阵的全部特征值及特征向量 D = Columns 1 through 7 4.0000 0 0 0 0 0 0 0.1382 5.9021 0.0000 0.0000 -0.0000 0.0000 0.0000 -0.2629 0.0000 5.6180 0.0000 0.0000 -0.0000 -0.0000 -0.3618 0.0000 0.0000 5.1756 -0.0000 0 -0.0000 -0.4253 -0.0000 -0.0000 -0.0000 4.6180 -0.0000 -0.0000 0.4472 -0.0000 0.0000 0.0000 -0.0000 4.0000 -0.0000 -0.4253 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 3.3820 0.3618 -0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 0.2629 0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 0.1382 0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 Columns 8 through 10 0 0 0 0 -0.0000 -0.0000 0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 2.8244 -0.0000 0.0000 0.0000 2.3820 -0.0000 -0.0000 -0.0000 2.0979 V = Columns 1 through 7 1.0000 0 0 0 0 0 0 0.1382 -0.2629 -0.3618 -0.4253 0.4472 -0.4253 0 0.2629 -0.4253 -0.4253 -0.2629 -0.0000 0.2629 0 0.3618 -0.4253