| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1628 人关注过本帖
标题:想用循环语句和二维数组来编写程序计算行列式,但不知哪错了
只看楼主 加入收藏
c语言求助
Rank: 1
等 级:新手上路
帖 子:8
专家分:0
注 册:2016-9-24
结帖率:66.67%
收藏
已结贴  问题点数:4 回复次数:2 
想用循环语句和二维数组来编写程序计算行列式,但不知哪错了
程序代码:
#include <stdio.h>
int main()
{
   int a[50][50];
   int i,j,n,v,s,p,q;
   printf("阶数为:\t");
   scanf("%d",&n);
   printf("please input %d data:\n",n);
   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("please input a data as a[%d][%d]:\t",i,j);
         scanf("%d",&a[i][j]);
      }
   }//构建行列式
   
   printf("please check\n");
   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("%5.d",a[i][j]);      
      }
      printf("\n");
   }//检查行列式

   for(j=1;j<=(n-1);j++)
   {
      for(i=j+1;i<=(n-1);i++)
      {
          for(v=1;v<=n;v++)
          {
              p=(a[i][j]/a[i+1][j])*a[i+1][v];
              a[i][v]-=p;
          }      
      }
      for(v=1;v<=n;v++)
      {
          q=(a[i][j]/a[j][j])*a[j][v];
         a[i][v]-=q;
      }
   }//转换为三角行列式

   printf("\n");

   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("%5.d",a[i][j]);      
      }
      printf("\n");
   }//检查三角行列式
   
   s=1;
   for(i=1;i<=n;i++)
   {
     s*=a[i][i];  
   }//计算三角行列式

   printf("结果为 %d\n",s);//输出结果
   

   return 0;
}

想把行列式变成三角行列式,然后主对角线上元素的乘积即为结果。
变成三角行列式的方法是一列一列地运用行列式性质“一列上的元素加另一列的对应元素的常数倍,行列式值不变”来弄出左下方的那堆0
程序在转化为三角行列式那里错了,但检查不出,望各大神指正。谢谢!
搜索更多相关主题的帖子: 编写程序 please 行列式 
2016-10-26 00:12
word123
Rank: 10Rank: 10Rank: 10
等 级:贵宾
威 望:13
帖 子:333
专家分:1622
注 册:2014-4-5
收藏
得分:4 
#include <stdio.h>
int main()
{
   int a[50][50];
   int i,j,n,v,s,p,q;
   printf("阶数为:\t");
   scanf("%d",&n);
   printf("please input %d data:\n",n);
   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("please input a data as a[%d][%d]:\t",i,j);
         scanf("%d",&a[i][j]);
      }
   }//构建行列式
   
   printf("please check\n");
   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("%5.d",a[i][j]);      
      }
      printf("\n");
   }//检查行列式
/*
这是可以的,满足两个条件
    1 2     基础行第一个数不为0
    3 4     后面行第一个数3除基础行第一个数1是整数
*/
   for(j=1;j<=(n-1);j++)       //我感觉 你是想行变换,我改成转换为行阶梯型矩阵了
   {
      for(i=j+1;i<=n;i++)      //你这个程序会有几个问题,1、若两行第一个相除不是整数倍,第一个数就不能变为0
      {                                               2、若当前行第一个数a[j][j]为0,就会报错,应该寻找当前行当前第一个数不为0的行来作为基础行
          p=a[i][j]/a[j][j];
          for(v=j;v<=n;v++)
          {
              a[i][v]-=p*a[j][v];
          }      
      }
   }//转换为三角行列式

   printf("\n");

   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("%5.d",a[i][j]);      
      }
      printf("\n");
   }//检查三角行列式
   
   s=1;
   for(i=1;i<=n;i++)
   {
     s*=a[i][i];  
   }//计算三角行列式

   printf("结果为 %d\n",s);//输出结果
   

   return 0;
}
2016-10-26 10:39
c语言求助
Rank: 1
等 级:新手上路
帖 子:8
专家分:0
注 册:2016-9-24
收藏
得分:0 
二楼给了我一些启发,后来我发现我的错误在于把p和q在循环内部赋值。
程序代码:
for(j=1;j<=(n-1);j++)
   {
      for(i=j+1;i<=(n-1);i++)
      {
          for(v=j;v<=n;v++)
          {
              p=(a[i][j]/a[i+1][j])*a[i+1][v];
              /*  错误:我们希望的倍数是一个常数,恒等于第一列的两项的比值。若将a[i][j]/a[i+1][j]放入循环中,
              则会出现以下情况:我们希望第二行每个元素都加上倍的第三行的对应元素,在第一个
              循环后a[2][1]已经变成了0,但由于每次循环都重新给p赋(a[2][1]/a[3][1])*a[3][v],所以之后的
              p皆变为了0(下面q同理)  */
              a[i][v]-=p;
          }      
      }
      printf("a[2][2]=%f\t",a[2][2]);//检错程序

      for(v=1;v<=n;v++)
      {
          q=(a[i][j]/a[j][j])*a[j][v];
         a[i][v]-=q;
      }
   }//转换为三角行列式


把p和q放在循环外赋值就可以了
程序代码:
#include <stdio.h>
int main()
{
   float a[50][50];
   int i,j,n,v;
   float s,p,q;
   printf("阶数为:\t");
   scanf("%d",&n);
   printf("please input %d data:\n",n);
   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("please input a data as a[%d][%d]:\t",i,j);
         scanf("%f",&a[i][j]);
      }
   }//构建行列式
   
   printf("please check\n");
   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("%5.2f",a[i][j]);      
      }
      printf("\n");
   }//检查行列式

   for(j=1;j<=(n-1);j++)
   {
      for(i=j+1;i<=(n-1);i++)
      {
          p=a[i][j]/a[i+1][j];
          for(v=j;v<=n;v++)
          {
              
              a[i][v]-=p*a[i+1][v];
          }      
      }
      

      q=a[i][j]/a[j][j];
      for(v=1;v<=n;v++)
      {
                 a[i][v]-=q*a[j][v];
      }
   }//转换为三角行列式

   printf("\n");

   for(i=1;i<=n;i++)
   {
      for(j=1;j<=n;j++)
      {
         printf("%5.2f",a[i][j]);      
      }
      printf("\n");
   }//检查三角行列式
   
   
   
   s=1;
   for(i=1;i<=n;i++)
   {
     s*=a[i][i];  
   }//计算三角行列式

   printf("结果为 %f\n",s);
         return 0;
}


该程序现阶段还不适用于一些存在0元素的行列式,以及行列间成倍数的行列式。
之后会改进一下。
谢过二楼!
2016-10-27 11:38
快速回复:想用循环语句和二维数组来编写程序计算行列式,但不知哪错了
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.027206 second(s), 10 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved