主要研究希尔伯特矩阵的病态性.docVIP

题目：主要研究希尔伯特矩阵的病态性 Hilbert矩阵的定义为： H 该矩阵是一个对称正定矩阵，其条件数随n的增加迅速增加。 画出ln（cond（H））~n之间的曲线，观察他们的关系。其中cond（H）表示矩阵在2范数下的条件数，它可按如下公式计算： cond（H） ，和分别是H的最大特征值和最小特征值。 设D是H对角元素开方构成的对角矩阵，令 D H D，不难看出依然还是对称矩阵，，并且对角元素都是1，画出ln（cond（H））~之间的曲线，观察他们之间的关系。 对于4n12，给定不同的右端项b，分别用x H；x D（ Db）；x H\b，求解Hx b，比较结果x 、x、 x。 程序如下： 首先，建立eiggg-M文件，求解矩阵的最大特征值和最小特征值。 function [max,min] eiggg s ss size s if ss 1; max s 1 ; min s 1 ; else max s 1 ; min s 1 ; for i 1:1:ss; if s i max max s i ; end if s i min min s i ; end end end 再建立Hilbert.m文件，给出Hilbert矩阵的定义： function H Hilbert n H zeros n,n ; for i 1:1:n; for j 1:1:n; H i,j 1/ i+j-1 ; end end 然后求出D的逆矩阵 function H Hilbert2 n H zeros n,n ; for j 1:1:n; H j,j 1/ 2*j-1 ; end D sqrt H ; D1 inv D ; H D1*H*H; 编译主程序： clear; clc; N 4; mv2 zeros N,1 ; for i 1:1:N; H Hilbert2 i ; s eig H ; [max,min] eiggg s ; mv1 max/min; mv2 i log mv1 ; end n 1:1:N; figure 1 ,subplot 221 ,plot n,mv2 n ; title ln cond Hn 2 与n之间的关系 ; H1 Hilbert N ; H2 Hilbert2 N ; b [1;2;3;4;];; x1 inv H1 *b; x2 H2*b; x3 H1\b; sss size x1 i 1:1:sss 1 ; subplot 222 ,plot i,x1 ; subplot 223 ,plot i,x2 ; subplot 224 ,plot i,x3 ; 运行之后的结果： ss 1 1 ss 2 1 ss 3 1 ss 4 1 sss 4 1 得出图像： 由此可见，其条件数随n的增加而增加，则该矩阵不收敛，故则证得希尔伯特矩阵式的病态性。