c语言实现高斯消元debug

Tom学编程

各位大佬，我i写的高斯消元求方程组的程序计算三阶方阵时就没有算不对，计算二阶时能得出正确结果，请问除了什么问题？谢谢大家！
代码如下：

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<math.h>
   
4. #pragma warning(disable:4996)
   
5. 
   
6. void swap(double **row,double *b,int r,int t)//将矩阵A和向量b的第r行和第t行交换
   
7. {   double *m,n;
   
8.         m=row[r];row[r]=row[t];row[t]=m;
   
9.         n=b[r];b[r]=b[t];b[t]=n;
   
10.         return;
    
11. }
    
12. 
    
13. void gauss(double **row,int n,double *b)//高斯消元
    
14. {   int i,j,flag,k;
    
15.         for(j=0;j<n;j++)
    
16.         {   flag=0;
    
17.                 for(i=j;i<n;i++)
    
18.                 {   if(row[i][j]>0) flag=1;
    
19.                         if(flag==1&&row[j][j]==0&&i!=j) swap(row,b,j,i);//如果第j行第j列为0，通过交换使其不为0
    
20.                         if(flag==1&&row[i][j]!=0&&i!=j)
    
21.                         {   for(k=j+1;k<n;k++)//消去i,j元素
    
22.                                 {   row[i][k]=row[i][k]-row[i][j]/row[j][j]*row[j][k];
    
23.                                         b[i]=b[i]-row[i][j]/row[j][j]*b[j];
    
24.                                 }
    
25.                                 row[i][j]=0;
    
26.                         }
    
27.                 }
    
28.         }
    
29. }
    
30. 
    
31. int main()
    
32. {   double **row,*b,*x;
    
33.         int i,j,n;
    
34.         printf("请输入矩阵阶数:");
    
35.         scanf("%d",&n);
    
36.         row=(double**)malloc(sizeof(double*)*n);
    
37.         for(i=0;i<n;i++)
    
38.                 row[i]=(double*)malloc(sizeof(double)*n);
    
39.         b=(double*)malloc(sizeof(double)*n);
    
40.         x=(double*)malloc(sizeof(double)*n);
    
41.         printf("请输入矩阵A:\n");
    
42.         for(i=0;i<n;i++)
    
43.                 for(j=0;j<n;j++)
    
44.                         scanf("%lf",&row[i][j]);
    
45.         printf("请输入向量b:\n");
    
46.         for(i=0;i<n;i++)
    
47.                 scanf("%lf",&b[i]);
    
48.         gauss(row,n,b);
    
49.         if(fabs(row[n-1][n-1])<1e-8) printf("No Solution or Infinite solutions!\n");
    
50.         else
    
51.         {   for(i=n-1;i>=0;i--)
    
52.                 {   for(j=n-1;j>i;j--)
    
53.                                 b[i]=b[i]-row[i][j]*x[j];
    
54.                         x[i]=b[i]/row[i][i];
    
55.                 }
    
56.                 for(i=0;i<n;i++)
    
57.                         printf("x%d=%.2lf\n",i+1,x[i]);
    
58.         }
    
59.         free(b);free(x);
    
60.         for(i=n-1;i>=0;i--)
    
61.                 free(row[i]);
    
62.         free(row);
    
63.         return 0;
    
64. }

复制代码

wp231957

给两个例子吧，关于向量的部分  不懂

Tom学编程

已经解决了

Tom学编程

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<math.h>
   
4. #pragma warning(disable:4996)
   
5. 
   
6. void swap(double **row,double *b,int r,int t)//将矩阵A和向量b的第r行和第t行交换
   
7. {   double *m,n;
   
8.         m=row[r];row[r]=row[t];row[t]=m;
   
9.         n=b[r];b[r]=b[t];b[t]=n;
   
10.         return;
    
11. }
    
12. 
    
13. void gauss(double **row,int n,double *b)//高斯消元
    
14. {   int i,j,flag,k;
    
15.         for(j=0;j<n;j++)
    
16.         {   flag=0;
    
17.                 for(i=j;i<n;i++)
    
18.                 {   if(row[i][j]>0) flag=1;
    
19.                         if(flag==1&&row[j][j]==0&&i!=j) swap(row,b,j,i);//如果第j行第j列为0，通过交换使其不为0
    
20.                         if(flag==1&&row[i][j]!=0&&i!=j)
    
21.                         {   for(k=j+1;k<n;k++)//消去i,j元素
    
22.                                     row[i][k]=row[i][k]-row[i][j]/row[j][j]*row[j][k];
    
23.                                 b[i]=b[i]-row[i][j]/row[j][j]*b[j];
    
24.                                 row[i][j]=0;
    
25.                         }
    
26.                 }
    
27.         }
    
28.         return;
    
29. }
    
30. 
    
31. int main()
    
32. {   double **row,*b,*x;
    
33.         int i,j,n;
    
34.         printf("请输入矩阵阶数:");
    
35.         scanf("%d",&n);
    
36.         row=(double**)malloc(sizeof(double*)*n);
    
37.         for(i=0;i<n;i++)
    
38.                 row[i]=(double*)malloc(sizeof(double)*n);
    
39.         b=(double*)malloc(sizeof(double)*n);
    
40.         x=(double*)malloc(sizeof(double)*n);
    
41.         printf("请输入矩阵A:\n");
    
42.         for(i=0;i<n;i++)
    
43.                 for(j=0;j<n;j++)
    
44.                         scanf("%lf",&row[i][j]);
    
45.         printf("请输入向量b:\n");
    
46.         for(i=0;i<n;i++)
    
47.                 scanf("%lf",&b[i]);
    
48.         gauss(row,n,b);
    
49.         if(fabs(row[n-1][n-1])<1e-8) printf("No Solution or Infinite solutions!\n");
    
50.         else
    
51.         {   for(i=n-1;i>=0;i--)
    
52.                 {   for(j=n-1;j>i;j--)
    
53.                                 b[i]=b[i]-row[i][j]*x[j];
    
54.                         x[i]=b[i]/row[i][i];
    
55.                 }
    
56.                 for(i=0;i<n;i++)
    
57.                         printf("x%d=%.3lf\n",i+1,x[i]);
    
58.         }
    
59.         free(b);free(x);
    
60.         for(i=n-1;i>=0;i--)
    
61.                 free(row[i]);
    
62.         free(row);
    
63.         return 0;
    
64. }

复制代码