第五章上机 文档.docVIP

两个矩阵的相乘操作 【实训内容】 利用数组实现两个矩阵的相乘。 【思路】 前面曾经提到过矩阵的乘法，它的步骤是就是两个矩阵对应位置上的元素相乘，然后再存储到结果矩阵的相同位置上。这种乘法叫做点乘运算。本实训所讲的是另一种矩阵乘法。例如，两个矩阵 ，，则两个矩阵相乘的结果矩阵r等于 （5.14） 在上式中，r00=a00b00+a01b10+ a02b20+…+a0nbn0，r01=a00b01+a01b11+a02b21+…+a0nbn1，…，rmm=am0b0m+am1b1m+am2b2m+…+amnbnm。由此可知，r矩阵的每一个元素都是a矩阵的一行和b矩阵的一列作内积的结果，也就是。因此对于这种乘法，就需要前一个矩阵的列数等于后一个矩阵的行数，否则无法运算。而且两个矩阵的顺序不能颠倒。 【程序源代码】 /* ==================== Program Description ===================*/ /* 程序名称： MatrixMutiply */ /* 程序目的：两个矩阵相乘 */ /* =============================================================*/ #include stdio.h ——————————— M 4 ——————————— N 5 /* ==================== Program Description ===================*/ /*相乘运算 */ /* =============================================================*/ void MatrixMultiply2(int A[M][N], int B[N][M],int result[M][M]) /* 两个二维数组相乘 */ { int i,j,k; int temp; for(i=0;iM;i++) for(j=0;jM;j++) { temp=0; for(k=0;kN;k++) /* 内积运算 */ —————————————————————————————————; —————————————————————————————————; } } /* ==================== Program Description ===================*/ /*主函数 */ /* =============================================================*/ void main() { int i,j; int a[M][N]={{1,3,5,7,9},{0,2,4,6,6},{2,3,7,5,8},{1,7,2,9,0}}; /* 4*5矩阵 */ int b[N][M]={{3,1,4,5},{6,7,5,9},{4,4,3,1},{7,6,8,9},{1,2,3,4}}; /* 5*4矩阵 */ int result[M][M]; printf(矩阵a：\n); for(i=0;i—————;i++) /* 显示矩阵a */ for(j=0;j————————;j++) { printf(%d ,a[i][j]); if(———————————) printf(\n); } printf(矩阵b：\n); for(i=0;iN;i++) /* 显示矩阵b */ for(j=0;jM;j++) { printf(%d ,b[i][j]); if(———————————) printf(\n); } printf(两个数组相乘等于：\n); MatrixMultiply2(a,b,r