C语言求矩阵的乘法！！！！！！！！！ - C语言论坛

楼主

已结贴√ 问题点数：5 回复次数：3

C语言求矩阵的乘法！！！！！！！！！

有没有一种我输入N，程序就求出矩阵M的N次方的程序？？？？我自己写了一个，多一次方就多一个for循环，理论上何以求到N次，但是觉得不够简洁。。。。。。。

程序代码：

#include<stdio.h>
#include<conio.h>
void main(void){
int M[5][5],Mk[5][5],MKK[5][5],MKKK[5][5],i,j,p;
M[0][0]=0;M[0][1]=0;M[0][2]=0;M[0][3]=0;M[0][4]=1;
M[1][0]=1;M[1][1]=0;M[1][2]=1;M[1][3]=0;M[1][4]=0;
M[2][0]=0;M[2][1]=0;M[2][2]=0;M[2][3]=0;M[2][4]=1;
M[3][0]=1;M[3][1]=0;M[3][2]=1;M[3][3]=0;M[3][4]=0;
M[4][0]=0;M[4][1]=1;M[4][2]=0;M[4][3]=1;M[4][4]=0;
for(i=0;i<5;i++){
  for(j=0;j<5;j++){Mk[i][j]=0;
    for( p=0;p<5;p++)
        Mk[i][j]=Mk[i][j]+M[i][p]*M[p][j];
  }
}
for(i=0;i<5;i++){
  for(j=0;j<5;j++){MKK[i][j]=0;
    for( p=0;p<5;p++)
        MKK[i][j]=MKK[i][j]+Mk[i][p]*M[p][j];
  }
}
for(i=0;i<5;i++){
  for(j=0;j<5;j++){MKKK[i][j]=0;
    for( p=0;p<5;p++)
        MKKK[i][j]=MKKK[i][j]+MKK[i][p]*M[p][j];
  }
}
for(i=0;i<5;i++){
  for(j=0;j<5;j++)
     printf("%8d",MKKK[i][j]);
     printf("\n");}
getch();
}

搜索更多相关主题的帖子: C语言　 color　

第 2 楼

得分:0

回复 2 楼 wyw19880809

谢谢，我在根据这个思路想想

第 3 楼

得分:0

回复 5 楼魑魅魍魉！

一个个赋值，便于查找修改，等程序写的成功了，就该用定义时赋值，如果写函数的话，矩阵的幂不知道怎么表示？有的人说用循环，有的人说递归，搞不懂了。。。

第 4 楼

得分:0

经过几天思考，现在来结贴，以下是关于矩阵N次幂的求法：

程序代码：

#include<stdio.h>
#include<conio.h>
#define A  5
#define B  5//A,B的值可根据具体的矩阵来确定；
void main(void){
    int M[A][B]={0,0,0,0,1,1,0,1,0,0,0,0,0,0,1,1,0,1,0,0,0,1,0,1,0},MT[A][B],
       ME[A][B]={0,0,0,0,1,1,0,1,0,0,0,0,0,0,1,1,0,1,0,0,0,1,0,1,0},i,j,p,c,k;//M[A][B]为初始矩阵，MT[A][A]为中转矩阵；ME[A][B]为结果矩阵
printf("M的()次幂:");
scanf("%d",&k);
for(c=0;c<k-1;c++){
    for(i=0;i<5;i++){
      for(j=0;j<5;j++){
         MT[i][j]=0;
         for( p=0;p<5;p++)
             MT[i][j]=MT[i][j]+ME[i][p]*M[p][j];//迭代
          }
          } 
for(i=0;i<5;i++){
      for(j=0;j<5;j++){
        ME[i][j]=*(MT[i]+j);}
      
    }
}
for(i=0;i<5;i++){
      for(j=0;j<5;j++){
        printf("%4d",ME[i][j]);}
      printf("\n");}
getch();
}

[ 本帖最后由追寻XT 于 2014-9-17 14:34 编辑 ]