矩阵乘法
用户输入Am×n,Bn×k两个矩阵的元素,计算它们的乘积
并输出。其中m,n,k也由用户输入,它们均不超过20。
问题分析
输入矩阵的元素只要根据行数、列数,使用循环即可。设C是矩阵的乘积,则矩阵乘法的运算公式是
算法描述
代码实现
#include<iostream>
using namespace std;
int main()
{
const int M = N = 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;
cout<<"请按行输入第1个矩阵的元素"<<endl;
for(i=0; i<M; ++i)
{
for(j=0; j<N; ++j)
{
cin >> A[i][j];
}
}
cout<<"请输入第2个矩阵的维数N K"<<endl;
cin>>N2>>K1; //输入第2个矩阵的
while(N2 != N1)
{
cout<<“第2个矩阵的行数应等于第1个矩阵的列数,请重输"<<endl;
cin >> N2 >> N1;
}
cout << 请按行输入第2个矩阵的元素"<<endl;
for(i=0; i< N1; ++i)
{
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; // init
for(k=0; k<N1; ++k)
C[i][j] = C[i][j] + A[i][j] * B[k][j];
}
}
// 输出乘积矩阵的元素
for(i=0;i<M1;i++)
{ //M1行
for(j=0;j<K1;j++)
{ //K1列
cout<<C[i][j]<<"\t"; //行中中间用Tab键分隔
}
cout<<endl; //行间换行
}
return 0;
}