C ++程序执行矩阵乘法
矩阵是数字的矩形阵列，以行和列的形式排列。
矩阵的示例如下。
一个3*2矩阵具有3行2列，如下所示-
8 1 4 9 5 6
执行矩阵乘法的程序如下。
示例
#include<iostream>
using namespace std;
int main() {
int product[10][10], r1=3, c1=3, r2=3, c2=3, i, j, k;
int a[3][3] = { {2, 4, 1} , {2, 3, 9} , {3, 1, 8} };
int b[3][3] = { {1, 2, 3} , {3, 6, 1} , {2, 4, 7} };
if (c1 != r2) {
cout<<"Column of first matrix should be equal to row of second matrix";
} else {
cout<<"第一个矩阵是："<<endl;
for(i=0; i<r1; ++i) {
for(j=0; j<c1; ++j)
cout<<a[i][j]<<" ";
cout<<endl;
}
cout<<endl;
cout<<"第二个矩阵是："<<endl;
for(i=0; i<r2; ++i) {
for(j=0; j<c2; ++j)
cout<<b[i][j]<<" ";
cout<<endl;
}
cout<<endl;
for(i=0; i<r1; ++i)
for(j=0; j<c2; ++j) {
product[i][j] = 0;
}
for(i=0; i<r1; ++i)
for(j=0; j<c2; ++j)
for(k=0; k<c1; ++k) {
product[i][j]+=a[i][k]*b[k][j];
}
cout<<"这两个矩阵的乘积为："<<endl;
for(i=0; i<r1; ++i) {
for(j=0; j<c2; ++j)
cout<<product[i][j]<<" ";
cout<<endl;
}
}
return 0;
}输出结果
第一个矩阵是： 2 4 1 2 3 9 3 1 8 第二个矩阵是： 1 2 3 3 6 1 2 4 7 这两个矩阵的乘积为： 16 32 17 29 58 72 22 44 66
在上面的程序中，两个矩阵a和b初始化如下：
int a[3][3] = { {2, 4, 1} , {2, 3, 9} , {3, 1, 8} };
int b[3][3] = { {1, 2, 3} , {3, 6, 1} , {2, 4, 7} };如果第一矩阵中的列数不等于第二矩阵中的行数，则不能执行乘法。在这种情况下，将显示错误消息。给出如下。
if (c1 != r2) {
cout<<"Column of first matrix should be equal to row of second matrix";
}矩阵a和b都使用嵌套的for循环显示。下面的代码段对此进行了演示。
cout<<"第一个矩阵是："<<endl;
for(i=0; i<r1; ++i) {
for(j=0; j<c1; ++j)
cout<<a[i][j]<<" ";
cout<<endl;
}
cout<<endl;
cout<<"第二个矩阵是："<<endl;
for(i=0; i<r2; ++i) {
for(j=0; j<c2; ++j)
cout<<b[i][j]<<" ";
cout<<endl;
}
cout<<endl;此后，将product[][]矩阵初始化为0。然后使用嵌套的for循环查找2个矩阵a和b的乘积。下面的代码段对此进行了演示。
for(i=0; i<r1; ++i)
for(j=0; j<c2; ++j) {
product[i][j] = 0;
}
for(i=0; i<r1; ++i)
for(j=0; j<c2; ++j)
for(k=0; k<c1; ++k) {
product[i][j]+=a[i][k]*b[k][j];
}获得产品后，将其打印。如下所示。
cout<<这两个矩阵的乘积为："<<endl;
for(i=0; i<r1; ++i) {
for(j=0; j<c2; ++j)
cout<<product[i][j]<<" ";
cout<<endl;
}