天天看点

实现两个N*N矩阵的乘法,矩阵由一维数组表示

实现两个N*N矩阵的乘法,矩阵由一维数组表示。

先介绍一下矩阵的加法:

实现两个N*N矩阵的乘法,矩阵由一维数组表示
实现两个N*N矩阵的乘法,矩阵由一维数组表示
实现两个N*N矩阵的乘法,矩阵由一维数组表示
实现两个N*N矩阵的乘法,矩阵由一维数组表示

若两个矩阵要做乘法运:只有在一个矩阵的行数与另一个矩阵的列数相同时,才能做两个矩阵的乘法。

如何得到矩阵的转置:

矩阵的转置也是一个矩阵,原始矩阵中的行转变为转置矩阵的列。例如,有下述一个3×3矩阵:

1 2 3

6 7 8

4 5 9

那么它的转置矩阵为:

1 6 4

2 7 5

3 8 9

让我们从程序员的角度仔细地考察一下这一现象。假设原始数组为M,转置矩阵为MT。那么M[1][0]=6,在转置矩阵中我们发现MT [0][1]=6。因此,我们能够得到程序化的结论:转置一个矩阵实际上就是对换下标变量。用技术术语讲:

MT[Row][Column] = M[Column][Row]; 

下面是得到转置矩阵的C代码:

void show_transpose(float mat[][10],int row,int col)    

{    

   int i,j;    

   for(i=0;i<row;i++)    

   {    

      for(j=0;j<col;j++)    

         printf("%f\t",mat[j][i]);    

      printf("\n");    

   }    

}   

以上方法显示了矩阵的转置。

#include<iostream>  

using namespace std;  

#define size 2  

int multi(int *a , int *b , int N)  

{  

    int i , j , k , temp;  

    int *c = (int*)malloc(N * N * sizeof(int));  

    for(i = 0 ; i < N ; i++)  

    {  

        for(j = 0 ; j < N ; j++)  

        {  

            temp = i * N + j;  

            *(c + temp) = 0;  

            for(k = 0 ; k < N ; k++)  

            {  

                *(c + temp) += a[i * N + k] * b[k * N + j];  

            }  

            cout<<*(c + temp)<<" ";  

        }  

    }  

    return *c;  

}  

int main()  

    int a[size * size] = {2 , 1 , 4 , 3};  

    int b[size * size] = {1 , -1 , 3 , 2};  

    multi(a , b , size);  

    return 0;  

本文转自夏雪冬日博客园博客,原文链接:http://www.cnblogs.com/heyonggang/p/3185648.html,如需转载请自行联系原作者

继续阅读