4.8 实例：矩阵乘法

实例：矩阵乘法

问题分析：

算法描述：

编程：
#include<iostream>

using namespace std;

int main()

{

       const int M=20,N=20,K=20;

       double A[M][N],B[N][K],C[M][K];

       int M1,N1,N2,K1;  //两个矩阵的实际维数

       int i,j,k;

       cout<<"请输入第1个矩阵的维数M N："<<endl;

       cin>>M1>>N1;  //输入第1个矩阵的行数和列数

       cout<<"请按行输入第1个矩阵的元素："<<endl;

       for(i=0;i<M1;i++){

              for(j=0;j<N1;j++){

                     cin>>A[i][j];

              }

       }

       cout<<"请输入第2个矩阵的维数N K："<<endl;

       cin>>N2>>K1;

       while(N2!=N1){

              cout<<"第2个矩阵的行数应等于第1个矩阵的列数，请重输："<<endl;

       }

       cout<<"请按行输入第2个矩阵的元素："<<endl;

       for(i=0;i<N1;i++){  //按行输入第2个矩阵的元素

              for(j=0;j<K1;j++){

                     cin>>B[i][j];

              }

       }

       for(i=0;i<M1;i++){

              for(j=0;j<K1;j++){

                     C[i][j]=0;  //i行j列元素赋初值0

                     for(k=0;k<N1;k++){

                            C[i][j]=C[i][j]+A[i][k]*B[k][j];

                     }

              }

       }

      

       //输出乘积矩阵的元素

       for(i=0;i<M1;i++){

              for(j=0;j<K1;j++){

                     cout<<C[i][j]<<"\t";

              }

              cout<<endl;

       }

       return 0;

}

运行结果：