改进平方根法.docVIP

2012-2013（1）专业课程实践论文 改进的平方根法 王 震 0818180210 R数学08-2班 张福贵 0818180227 R数学08-2班 一、算法理论 在平方根的基础上，为了避免开方运算,所以用计算；其中，； 得 按行计算的元素及对元素公式 对于 步1 . 步2 . 步3 计算出的第行元素后，存放在的第行相置，然后再计算的第行元素，存放在的第行.的对角元素存放在的相应位置. 例如 对称正定矩阵按分解和按分解计算量差不多，但分解不需要开放计算。求解, 的计算公式分别如步4、步5所述公式。 步4 步5 。 . 二、算法框图 三、算法程序 #includemath.h #includestdio.h main() {int i,j,k; int n=4; float m,a[4][4],t[4][4],l[4][4],b[4],d[4],x[4],y[4]; printf(Please input matrix a:\n); for(i=0;in;i++) for(j=0;jn;j++) {l[j][j]=1.0; scanf(%f,a[i][j]); } printf(\n); printf(Please input matrix b:\n); for(i=0;in;i++) scanf(%f,b[i]); printf(\n); d[0]=a[0][0]; for(i=1;in;i++) for(j=0;ji;j++) {m=0.0; for(k=0;kj;k++) m+=t[i][k]*l[j][k]; t[i][j]=a[i][j]-m; l[i][j]=t[i][j]/d[j]; m=0.0; for(k=0;ki;k++) {m+=t[i][k]*l[i][k]; d[i]=a[i][i]-m; } l[j][i]=l[i][j]; } y[0]=b[0]; for(i=1;in;i++) {m=0.0; for(k=0;ki;k++) {m+=l[i][k]*y[k]; y[i]=b[i]-m; } } x[n-1]=y[n-1]/d[n-1]; for(i=n-1;i=0;i--) {m=0.0; for(k=i+1;kn;k++) {m+=l[k][i]*x[k]; x[i]=y[i]/d[i]-m; } } for(i=0;in;i++) printf(x(%d)=%f \n,i,x[i]); } 四、算法实现 例1.，，求X。 解：运行程序 输入矩阵 6 2 1 -1 2 4 1 0 1 1 4 -1 -1 0 - 1 3 （2）输入矩阵B 得出结果 通过程序解得 例2. ，，，求X。 得出结果