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<<"The first matrix is:"<<endl; for(i=0; i<r1; ++i) { for(j=0; j<c1; ++j) cout<<a[i][j]<<" "; cout<<endl; } cout<<endl; cout<<"The second matrix is:"<<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<<"Product of the two matrices is:"<<endl; for(i=0; i<r1; ++i) { for(j=0; j<c2; ++j) cout<<product[i][j]<<" "; cout<<endl; } } return 0; }
输出
The first matrix is: 2 4 1 2 3 9 3 1 8 The second matrix is: 1 2 3 3 6 1 2 4 7 Product of the two matrices is: 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"; }
使用嵌套 for 循环显示矩阵 a 和 b。以下代码片段对此进行了演示。
cout<<"The first matrix is:"<<endl; for(i=0; i<r1; ++i) { for(j=0; j<c1; ++j) cout<<a[i][j]<<" "; cout<<endl; } cout<<endl; cout<<"The second matrix is:"<<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<<Product of the two matrices is:"<<endl; for(i=0; i<r1; ++i) { for(j=0; j<c2; ++j) cout<<product[i][j]<<" "; cout<<endl; }
广告