楚列斯分解.doc

2013-2014（1）专业课程实践论文 题目：楚列斯基分解 一、算法理论 Cholesky分解可用不选主元的Gauss消去法来实现。然而，更简单而实用的方法是通过直接比较两边的对应元素来计算的。设 比较两边对应的元素，得关系 ??????????????? （1） 首先，由得 再由 这样便得到了矩阵第一列元素。假定已经算出的前列元素，由 得 ?????????????????? （2） 再由 得 ?????????????? (3) 这样便又求出了的第列元素。这种方法称为平方根法。当然，亦可按行来逐次计算由于的元素被用来计算出以后不再使用，所以可的元素存储在的对应位置上，这样我们就得到如下算法。 二、算法框图 三、算法程序 A=input(请输入矩阵A：\nA=); disp(矩阵的阶数为：); n=length(A); lambda=eig(A); flagdc=0; flagzd=0; %判断矩阵A是不是正定矩阵 while (flagdc~=n*(n-1)/2)||(flagzd~=n) for i=1:n; for j=i+1:n; if A(i,j)==A(j,i); flagdc=flagdc+1; end; end; end; if(flagdc~=n*(n-1)/2) disp(矩阵A不是对称矩阵。请重新输入矩阵A:) A=input(A=\n); disp(矩阵的阶数为:) n=length(A); flagdc=0; end %判断矩阵A是不是正定矩阵 if flagdc==n*(n-1)/2 disp(矩阵A的特征值是:) lambda=eig(A); for i=1:n if(lambda(i)0) flagzd=flagzd+1; end end if flagzd~=n disp(矩阵A不是正定矩阵。请重新输入矩阵A:) A=input(A=\n); disp(矩阵的阶数为:) n=length(A) flagzd=0; flagdc=0; end end end %对矩阵A进行分解 l(1,1)=sqrt(A(1,1)); for j=2:n l(j,1)=A(j,1)/l(1,1); end for i=2:n for j=i:n p(j,i)=0; for k=1:i-1 p(j,i)=p(j,i)+l(i,k)*l(j,k); t(j,i)=A(j,i)-p(j,i); if(j==i) l(j,i)=sqrt(t(j,i)); else l(j,i)=t(j,i)/l(i,i); end end end end disp(对正定矩阵A的Cholesky分解为：) l LT=l 四、算法实现 例1．将实对称正交矩阵A=[1,1,1;1,2,3;1,3,6]用楚列斯基分解。 解： 例2．将矩阵A=[1,1,1;1,4,2;1,2,8]进行cholesky分解。 解：