1、用共轭梯度法解方程用Jacobi方法求矩阵的全部特征值和特征向量用共轭梯度法解方程,用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;endfor i=2:n-1 A(i,i)=4; A(i,i-1)=1; A(i,i+1)=1;endA(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(方程的精确解
2、n);Xfprintf(用共轭法求解方程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 = 100x = 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
3、 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.000
4、0 -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 -
5、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.0979V = 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 -0.1382 0.2629 -0.4472 0.2629 0 0.4253 -0.2629 0.
6、2629 0.4253 -0.0000 -0.4253 0 0.4472 0.0000 0.4472 0.0000 0.4472 0.0000 0 0.4253 0.2629 0.2629 -0.4253 0.0000 0.4253 0 0.3618 0.4253 -0.1382 -0.2629 -0.4472 -0.2629 0 0.2629 0.4253 -0.4253 0.2629 -0.0000 -0.2629 0 0.1382 0.2629 -0.3618 0.4253 0.4472 0.4253 Columns 8 through 10 0 0 0 0.3618 0.2629 0.
7、1382 -0.4253 -0.4253 -0.2629 0.1382 0.4253 0.3618 0.2629 -0.2629 -0.4253 -0.4472 0.0000 0.4472 0.2629 0.2629 -0.4253 0.1382 -0.4253 0.3618 -0.4253 0.4253 -0.2629 0.3618 -0.2629 0.1382 baogaoer矩阵阶数n=20方程精确解X = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1用共轭梯度法求解方程k = 76x = 1.2500 1.0104 0.7085 1.1554 0.66
8、97 1.1659 0.6669 1.1666 0.6667 1.1667 0.6667 1.1667 0.6666 1.1667 0.6664 1.1676 0.6632 1.1795 0.6188 1.3453用方法Jacobi求解矩阵的全部特征值及特征向量D = Columns 1 through 7 4.0000 0 0 0 0 0 0 0.0495 5.9754 0.0000 0.0000 0.0000 0.0000 -0.0000 -0.0977 -0.0000 5.9021 -0.0000 0.0000 0.0000 0.0000 -0.1436 -0.0000 0.0000 5
9、.7820 0.0000 -0.0000 -0.0000 -0.1859 -0.0000 0.0000 0.0000 5.6180 -0.0000 -0.0000 -0.2236 0.0000 -0.0000 0.0000 0.0000 5.4142 0.0000 -0.2558 0.0000 -0.0000 0.0000 0.0000 0.0000 5.1756 0.2818 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 0.3008 0.0000 0.0000 0.0000 0.0000 -0.0000 -0.0000 0.3123 -0.000
10、0 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.3162 -0.0000 -0.0000 -0.0000 0.0000 0.0000 -0.0000 -0.3123 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.3008 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.2818 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.2558 -0.0000 -0.0000 -0.0000 -0.0000 0.00
11、00 0.0000 0.2236 0.0000 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.1859 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.1436 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0977 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0495 0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 Columns 8 through 14 0 0 0 0 0
12、 0 0 -0.0000 -0.0000 -0.0000 0.0000 0.0000 0 -0.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 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 -0.0000 0.0000 0.0000 -0.0000 0.0000 0.
13、0000 0.0000 4.9080 0 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 4.6180 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 4.3129 0.0000 -0.0000 -0.0000 0.0000 0.0000 -0.0000 -0.0000 4.0000 0.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 3.6871 -0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
14、 3.3820 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 3.0920 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 -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.000
15、0 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 Columns 15 through 20 0 0 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 0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 -0.0000 -0.0000
16、 -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 -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 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000
17、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 -0.0000 -0.0000 0.0000 2.5858 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 2.3820 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 2.2180 0.0000 0.0000 -0.0000 0.0000 0.0000 -0
18、.0000 2.0979 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 2.0246V = Columns 1 through 7 1.0000 0 0 0 0 0 0 0 0.0495 -0.0977 -0.1436 -0.1859 -0.2236 -0.2558 0 0.0977 -0.1859 -0.2558 -0.3008 -0.3162 -0.3008 0 0.1436 -0.2558 -0.3123 -0.3008 -0.2236 -0.0977 0 0.1859 -0.3008 -0.3008 -0.1859 -0.0000 0.1859
19、 0 0.2236 -0.3162 -0.2236 -0.0000 0.2236 0.3162 0 0.2558 -0.3008 -0.0977 0.1859 0.3162 0.1859 0 0.2818 -0.2558 0.0495 0.3008 0.2236 -0.0977 0 0.3008 -0.1859 0.1859 0.3008 0.0000 -0.3008 0 0.3123 -0.0977 0.2818 0.1859 -0.2236 -0.2558 0 0.3162 -0.0000 0.3162 0.0000 -0.3162 -0.0000 0 0.3123 0.0977 0.28
20、18 -0.1859 -0.2236 0.2558 0 0.3008 0.1859 0.1859 -0.3008 -0.0000 0.3008 0 0.2818 0.2558 0.0495 -0.3008 0.2236 0.0977 0 0.2558 0.3008 -0.0977 -0.1859 0.3162 -0.1859 0 0.2236 0.3162 -0.2236 -0.0000 0.2236 -0.3162 0 0.1859 0.3008 -0.3008 0.1859 0.0000 -0.1859 0 0.1436 0.2558 -0.3123 0.3008 -0.2236 0.09
21、77 0 0.0977 0.1859 -0.2558 0.3008 -0.3162 0.3008 0 0.0495 0.0977 -0.1436 0.1859 -0.2236 0.2558 Columns 8 through 14 0 0 0 0 0 0 0 0.2818 0.3008 0.3123 0.3162 -0.3123 -0.3008 0.2818 0.2558 0.1859 0.0977 0.0000 0.0977 0.1859 -0.2558 -0.0495 -0.1859 -0.2818 -0.3162 0.2818 0.1859 -0.0495 -0.3008 -0.3008
22、 -0.1859 -0.0000 -0.1859 -0.3008 0.3008 -0.2236 -0.0000 0.2236 0.3162 -0.2236 0.0000 -0.2236 0.0977 0.3008 0.2558 0.0000 0.2558 0.3008 -0.0977 0.3123 0.1859 -0.1436 -0.3162 0.1436 -0.1859 0.3123 0.1859 -0.1859 -0.3008 -0.0000 -0.3008 -0.1859 -0.1859 -0.1436 -0.3008 0.0495 0.3162 -0.0495 0.3008 -0.14
23、36 -0.3162 -0.0000 0.3162 0.0000 0.3162 0.0000 0.3162 -0.1436 0.3008 0.0495 -0.3162 -0.0495 -0.3008 -0.1436 0.1859 0.1859 -0.3008 -0.0000 -0.3008 0.1859 -0.1859 0.3123 -0.1859 -0.1436 0.3162 0.1436 0.1859 0.3123 0.0977 -0.3008 0.2558 0.0000 0.2558 -0.3008 -0.0977 -0.2236 -0.0000 0.2236 -0.3162 -0.22
24、36 0.0000 -0.2236 -0.3008 0.3008 -0.1859 -0.0000 -0.1859 0.3008 0.3008 -0.0495 0.1859 -0.2818 0.3162 0.2818 -0.1859 -0.0495 0.2558 -0.1859 0.0977 0.0000 0.0977 -0.1859 -0.2558 0.2818 -0.3008 0.3123 -0.3162 -0.3123 0.3008 0.2818 Columns 15 through 20 0 0 0 0 0 0 -0.2558 0.2236 0.1859 0.1436 -0.0977 0
25、.0495 0.3008 -0.3162 -0.3008 -0.2558 0.1859 -0.0977 -0.0977 0.2236 0.3008 0.3123 -0.2558 0.1436 -0.1859 -0.0000 -0.1859 -0.3008 0.3008 -0.1859 0.3162 -0.2236 0.0000 0.2236 -0.3162 0.2236 -0.1859 0.3162 0.1859 -0.0977 0.3008 -0.2558 -0.0977 -0.2236 -0.3008 -0.0495 -0.2558 0.2818 0.3008 0.0000 0.3008
26、0.1859 0.1859 -0.3008 -0.2558 0.2236 -0.1859 -0.2818 -0.0977 0.3123 0.0000 -0.3162 0.0000 0.3162 0.0000 -0.3162 0.2558 0.2236 0.1859 -0.2818 0.0977 0.3123 -0.3008 -0.0000 -0.3008 0.1859 -0.1859 -0.3008 0.0977 -0.2236 0.3008 -0.0495 0.2558 0.2818 0.1859 0.3162 -0.1859 -0.0977 -0.3008 -0.2558 -0.3162
27、-0.2236 0.0000 0.2236 0.3162 0.2236 0.1859 0.0000 0.1859 -0.3008 -0.3008 -0.1859 0.0977 0.2236 -0.3008 0.3123 0.2558 0.1436 -0.3008 -0.3162 0.3008 -0.2558 -0.1859 -0.0977 0.2558 0.2236 -0.1859 0.1436 0.0977 0.0495 baogaoer矩阵阶数n=30方程精确解X = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1用共
28、轭梯度法求解方程k = 74x = 1.2500 1.0104 0.7085 1.1554 0.6697 1.1659 0.6669 1.1666 0.6667 1.1667 0.6667 1.1667 0.6667 1.1667 0.6667 1.1667 0.6667 1.1667 0.6667 1.1667 0.6667 1.1667 0.6666 1.1667 0.6664 1.1676 0.6632 1.1795 0.6188 1.3453用方法Jacobi求解矩阵的全部特征值及特征向量D = Columns 1 through 7 4.0000 0 0 0 0 0 0 0.0270
29、 5.9890 -0.0000 0.0000 0.0000 0.0000 0.0000 -0.0537 0.0000 5.9563 -0.0000 -0.0000 -0.0000 0 0.0798 -0.0000 -0.0000 5.9021 -0.0000 0.0000 -0.0000 0.1050 0.0000 -0.0000 0.0000 5.8271 -0.0000 0.0000 0.1291 0.0000 -0.0000 0.0000 0.0000 5.7321 0.0000 0.1518 -0.0000 0.0000 -0.0000 0.0000 0.0000 5.6180 0.1
30、728 0.0000 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.1919 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.2089 -0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.2236 0.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000 0.2359 -0.0000 0.0000 -0.0000 0.0000 0.0000 -0.0000 0.2456 0.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000 -0.2526 0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.2568 0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.2582 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 0.2568 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000 0.2526 -0.000