河北工业大学计算方法实验6..docVIP

实验六 线性方程组数值解法 实验题目： 1.用列主元消去法解方程组 2.用LU分解法求解方程组Ax=b A= b= ‘ 3.用雅可比迭代法解方程组 4.用高斯-塞尔德迭代法解方程组 1）、 源代码 #includestdio.h #includemath.h void main(void) { double a[3][4]={0.101,2.304,3.555,1.183,-1.347,3.712,4.623,2.137,-2.835,1.072,5.643,3.035},x1,x2,x3; int i,j; printf(原矩阵为：\n); for(i=0;i3;i++) { for(j=0;j4;j++) printf(%12lf,a[i][j]); printf(\n); } for(i=1;i3;i++) for(j=0;j4;j++) a[i][j]=-a[i][0]/a[0][0]*a[0][j]+a[i][j]; for(i=2,j=1;j4;j++) a[i][j]=-a[2][1]/a[1][1]*a[1][j]+a[i][j]; printf(消元后矩阵为：\n); for(i=0;i3;i++) { for(j=0;j4;j++) printf(%12lf,a[i][j]); printf(\n\n); } x3=a[2][3]/a[2][2]; x2=(a[1][3]-a[1][2]*x3)/a[1][1]; x1=(a[0][3]-a[0][2]*x3-a[0][1]*x2)/a[0][0]; printf(方程组的解如下：\n); printf(x1=%lf\n,x1); printf(x2=%lf\n,x2); printf(x3=%lf\n,x3); } 3、实验结果 、 源代码 #includeiostream.h #includestdio.h void main() { float x[4]; int i; float a[4][5]={1,2,-12,8,27, 5,4,7,-2,4, -3,7,9,5,11, -6,-12,-8,3,49}; void LU(float*,int,float[]); LU(a[0],4,x); for(i=0;i=3;i++)printf(x[%d]=%f\n,i,x[i]); } void LU(float*u,int n,float x[]) { int i,r,k; for(r=0;r=n-1;r++){ for(i=r;i=n;i++) for(k=0;k=r-1;k++) *(u+r*(n+1)+i)-=*(u+r*(n+1)+k)*(*(u+k*(n+1)+i)); for(i=r+1;i=n-1;i++){ for(k=0;k=r-1;k++) *(u+i*(n+1)+r)-=*(u+i*(n+1)+k)*(*(u+k*(n+1)+r)); *(u+i*(n+1)+r)/=*(u+r*(n+1)+r); } } for(i=n-1;i=0;i--){ for(r=n-1;r=i+1;r--) *(u+i*(n+1)+n)-=*(u+i*(n+1)+r)*x[r]; x[i]=*(u+i*(n+1)+n)/(*(u+i*(n+1)+i)); } } 实验结果 、 源代码 #includestdio.h void main() { double x[3]={0,0,0}; double y[3]; int i=0,j,n,k; scanf(%d,n); for(k=1;kn;k++) { y[0]=(8-2*x[1]-x[2])/5; y[1]=(21-2*x[0]+3*x[2])/8; y[2]=(1-x[0]+3*x[1])/6; for(j=0;j3;j++) { x[i]=y[j]; printf( %f ,x[i]); i++; } printf(\n); } } 实验结果 、