对称矩阵的平方根法.docVIP

数值计算方法 09医软（1）班 本组实验同学：刘康康 秦强 梅世友 马蕾 乔琼 任务分配：刘康康：分析平方根法解对称矩阵的思想，并写出推导过程 秦强：利用c++实现算法，并写出程序 梅世友：调试并运行程序 马蕾：根据解方程组的方法写出流程图 乔琼：负责word的排版（解题流程，程序截图等） 一．实验名称 平方根法解对称方程组 二．实验目的 理解平方根法解对称方程组的基本思想，原理以及公式的推导过程，和用c++实现算法 三．理论分析推导和程序实现 （一）、理论分析推导 设A为对称矩阵，且A的所有顺序主子式均不为0，由定理4（P178）可知，A可以唯一地分解为A=LU，其中L为单位下三角矩阵，U为上三角矩阵，为了利用A的对称性，再将U进行分解，即 = 或写成 U=D 其中D为对角矩阵，U1为单位上三角矩阵。于是 A=LU=LD （2.12） 又 A==D(分别表示A.L的转置矩阵) 由分解的唯一性，得=L。由（2.12）知 A=LD 且这种分解是唯一的。 如果A是对称正定的，利用代数知识可以证明（i=1，2，…，n）皆大于零。事实上，将A=LU按分块形式写出，有 （k=1,2,…，n)， 其中 =，=，= 由于A的各阶主子式皆大于零，所以 det=detdet=0. 上式对k=1,2，…，n成立，故有0（i=1,2，…，n）。 因为A为对称正定，所以0（i=1,2，…，n），于是可将D进一步分解为 D= = = 于是 A=LD=L =（L）=L。 其中L1=L为下三角矩阵。当限定L1的对角元素为正时，这种分解也是唯一的。 由 A=L1=， 用比较法可导出计算L1元素的计算公式，对i=1,2，…，n， =， = 这一方法称为平方根法，又称Cholesky分解。 （二）程序的实现 1.输入方程组的系数和b的值并判断矩阵是否为对称矩阵 #include stdio.h #include math.h #define N 20 int main() { int i,j,k; int size; int n=0; float a[N][N],l[N][N]; float b[N],x[N],y[N]; //float S,H,M,L; printf(\t\t\t\t*Cholesky分解法解方程*\n); printf(请输入方阵A的n：); scanf(%d,size); printf(\n); printf(请输入方程组的系数:\n); //数组a[][]的输入 for(i=0;isize;i++) { for(j=0;jsize;j++) { scanf(%f,a[i][j]); } } printf(\n请输入方程组的y:\n); //数组b[]的输入 for(i=0;isize;i++) scanf(%f,b[i]); printf(\n方阵A[][]为：\n); //2.数组A[][]的输出 for(i=0;isize;i++) { for(j=0;jsize;j++) { printf(%f ,a[i][j]); } printf(\n); } printf(\n方程组y为:\n); //数组b[]的输出 for(i=0;isize;i++) printf(%f ,b[i]); printf(\n\n); //3.判断是否可能用平方根法 for(i=0;isize;i++) { for(j=0;jsize;j++) { if(a[i][j]==a[j][i]) n=n+1; } } if(n!=size*size) { printf(次线性方程组的系数不是对称矩阵，不可以用平方根法，请按键退出\n); return 0; } else { printf(次线性方程组的系数是对称矩阵，可以用平方根法，算法如下：\n\n); //4.平方根法的计算过程 计算Lii l[0][0]=sqrt(a[0][0]); for(i=1;isize;i++) { l[i][0]=a[i][0]/l[0][0]; printf(l[%d][0]=%f\n,i,l[i][0]); } f