线性方程组的数值求解.docVIP

高斯消去法 #includestdio.h #includemath.h int main() { printf(高斯消去法求根\n); float a[4][4]={{7,-1,2,1},{3,9,-1,8},{10,-1,7,6},{-1,5,-1,6}}; float b[4]={-11,-11,6,25}; //方程的系数// float x[4]; float temp; int k,i,j; for(k=0;k4;k++) //消元，将系数转化为上三角阵// { for(i=k+1;i4;i++) {temp=a[i][k]/a[k][k]; b[i]=b[i]-temp*b[k]; for(j=k;j4;j++) {a[i][j]=a[i][j]-temp*a[k][j]; } } } for(j=3;j=0;j--) //回代依次求出x4,x3,x2,x1// {temp =0; for(i=j+1;i4;i++) {temp =temp+a[j][i]*x[i]; } x[j]=(b[j]-temp)/a[j][j]; } for(i=0;i4;i++) //依次输出x1,x2,x3,x4// { printf(%.6f\n,x[i]); } } 三角分解法 #include stdio.h #include math.h int main() { printf(三角分解法\n); float a[4][4]={{7,-1,2,1},{3,9,-1,8},{10,-1,7,6},{-1,5,-1,6}}; float b[4]={-11,-11,6,25}; //方程的系数// float l[4][4],u[4][4]; float x[4],y[4]; float x1,x2,x3,x4,sum; int i,j; //求矩阵L和U// for(i=0;i4;i++)l[i][i]=1; for(j=0;j4;j++)u[0][j]=a[0][j]; for(i=1;i4;i++)l[i][0]=a[i][0]/u[0][0]; for(j=1;j4;j++)u[1][j]=a[1][j]-l[1][0]*u[0][j]; for(i=2;i4;i++)l[i][1]=(a[i][1]-l[i][0]*u[0][1])/u[1][1]; for(j=2;j4;j++)u[2][j]=a[2][j]-l[2][0]*u[0][j]-l[2][1]*u[1][j]; for(i=3;i4;i++)l[i][2]=(a[i][2]-l[i][0]*u[0][2]-l[i][1]*u[1][2])/u[2][2]; u[3][3]=a[3][3]-l[3][0]*u[0][3]-l[3][1]*u[1][3]-l[3][2]*u[2][3]; for(j=0;j4;j++) //由L*y=b求y// {sum=0; for(i=j-1;i=0;i--) {sum=sum+l[j][i]*y[i]; } y[j]=b[j]-sum; } for(j=3;j=0;j--) //由U*x=y求x// {sum=0; for(i=j+1;i4;i++) { sum=sum+u[j][i]*x[i]; } x[j]=(y[j]-sum)/u[j][j]; } x1=x[0];x2=x[1];x3=x[2];x4=x[3]; printf(x1=%f\nx2=%f\nx3=%f\nx4=%f\n,x1,x2,x3,x4); } 雅克比迭代法 #includestdio.h #include math.h int main() { float a[4][4]={{7,-1,2,1},{3,9,-1,8},{10,-1,7,6},{-1,5,-1,6}}, b[4]={-11,-11,6,25}; //方程的系数// float x0[4]={0,0,0,0},x[4]; float e=1e-5; int K=1