matlab 怎样矩阵行优先,矩阵基础 (1). 行优先和列优先的问题

╰+攻爆jí腚メ 2023-01-13 13:57 282阅读 0赞

摘要

本文主要探讨的是“行优先”原则和“列优先”原则的问题。

1. 背景

首先了解“行优先”和“列优先”的知识，这两种方式在数学上的直观描述如下，给定如下矩阵：

![0818b9ca8b590ca3270a3433284dd417.png][]

根据行优先的原则，其排序方式为

![0818b9ca8b590ca3270a3433284dd417.png][]

根据列优先的原则，其排序方式为

![0818b9ca8b590ca3270a3433284dd417.png][]

2. 计算机领域的应用

行列优先原则在计算机领域的应用主要如下。行优先或者列优先没有好坏，但其直接涉及到对内存中数据的最佳存储访问方式。因为在内存使用上，程序访问的内存地址之间连续性越好，程序的访问效率就越高；相应地，程序访问的内存地址之间连续性越差。所以，我们应该尽量在行优先机制的编译器，比如C/C++，CUDA等等上，采用行优先的数据存储方式；在列优先机制的编译器，比如Fortune, Matlab等等上，采用列优先的数据存储方式。但这种思想渗透到编程中之后，代码的质量就会提高一个档次。

3. 以矩阵计算为例(Matlab编译器下测试)

% data

A = \[ 1 1

2 2

3 3

4 4

5 5

6 6

7 7

8 8

9 9\];

B = \[ 1 2 3 4 5 6 7 8 9

1 2 3 4 5 6 7 8 9\];

C = zeros(9,9);

% The method of matrix multiplication in Matlab

tic

C = A\*B;

toc

%　Our impletation method of matrix multiplication

tic

for ra = 1:9 % raws of the matrix A

for cb = 1:9 % columns of the matrix B

for len = 1:2

C(ra,cb) = A(ra,len)\*B(len,cb)+C(ra,cb);

end

end

end

toc

% Optimal method 1

tic

for cb = 1:9 % columns of the matrix B

for ra = 1:9 % raws of the matrix A

for len = 1:2

C(ra,cb) = A(ra,len)\*B(len,cb)+C(ra,cb);

end

end

end

toc

% Advanced optimal method 2

A = A'; % you can also directly given A = \[ 1 2 3 4 5 6 7 8 9

% 1 2 3 4 5 6 7 8 9\];

B = \[ 1 2 3 4 5 6 7 8 9

1 2 3 4 5 6 7 8 9\];

tic

for i = 1:9 % columns of the matrix A

for j = 1:9 % columns of the matrix B

for len = 1:2

C(i,j) = A(len,i)\*B(len,j)+C(i,j);

end

end

end

toc

4. 测试和分析

![0818b9ca8b590ca3270a3433284dd417.png][]

测试结果如上图所示，第一个时间为Matlab自带的乘法运算，第二个为我们原始实现的乘法计算，第三个为循环中行列变换(适应列优先编译器的处理)。

最重要的是第四个是本人原创的矩阵乘法方法，简单地说就是将A矩阵转置，然后设计相应的算法实现矩阵乘运算。在这个点上，希望在理解原理的基础上能给读者一些启发。在本例中，这样做效率最高，原因其一是本例中原始数据结构上适合我这样处理；原因其二是这样做的目的是使得任何一个子乘法的处理上，两乘数所在的内存空间上都是连续，而不仅仅是一个连续(注意：这是本文的核心，读者理解透了一定会很有收获，认真看我给出的程序实现。这是核心，不懂的可以交流思想)！

另外，本文中我给出的这个方法是矩阵乘法里面最优的方法，至少数学逻辑上是这样。之所以Matlab自带的乘法计算之所以性能还不错，是因为Matlab自带的运算都是经过优化的，包括硬件加速，系统加速等自己设计的应用很能调用加速方法。

[0818b9ca8b590ca3270a3433284dd417.png]: /images/20221022/2887c67e45584c76a4c6e714bde05497.png

发表评论取消回复

表情：

评论列表（有 0 条评论，282人围观）

还没有评论，来说两句吧...

相关阅读

相关 matlab 怎样矩阵行优先,矩阵基础 (1). 行优先和列优先的问题

摘要本文主要探讨的是“行优先”原则和“列优先”原则的问题。 1. 背景首先了解“行优先”和“列优先”的知识，这两种方式在数学上的直观描述如下，给定如下矩阵： ![0

╰+攻爆jí腚メ/ 2023年01月13日 13:57/ 0 赞/ 283 阅读

相关邻接矩阵实现的广度优先搜索

\include <stdio.h> \include <stdlib.h> \include <string.h> //循环队列最多 \define M

朴灿烈づ我的快乐病毒、/ 2022年09月18日 11:59/ 0 赞/ 136 阅读

相关 matlab 矩阵去掉重复行

问题背景在做多目标离散问题优化算法时，利用多个不同的进化算子产生了重复的解，所以有必要记录一下如何去掉重复的解，精简问题规模，尤其时迭代到后期，很多相同的解解决

痛定思痛。/ 2022年08月30日 04:29/ 0 赞/ 247 阅读

相关邻接矩阵存储的无向图深度优先(DFS)广度优先(BFS)遍历

图的两种存储方式：邻接矩阵和邻接表；两种遍历方式：深度优先和广度优先；首先以一个结构体存储一个图： struct MGraph { int

比眉伴天荒/ 2022年08月10日 09:40/ 0 赞/ 206 阅读

相关 JAVA将M行N列的矩阵交换行和列

题目要求：编写一段代码，打印出一个M行N列的二维数组的转置（交换行和列对应位置的元素）代码实现如下所示： package com.wonders.test;

r囧r小猫/ 2022年06月06日 09:42/ 0 赞/ 181 阅读

相关图的邻接矩阵以及深度优先遍历 + 广度优先遍历

图的邻接矩阵表示法非常简单，一个定点数组，一个二维数组搞定，类似与这样 ![这里写图片描述][SouthEast] 下面简单实现一个邻接矩阵表示的方法的图，以及遍历的两

左手的ㄟ右手/ 2022年06月01日 08:40/ 0 赞/ 317 阅读

相关图：图的邻接矩阵创建、深度优先遍历和广度优先遍历详解

邻接矩阵介绍直接说，邻接矩阵是图的一种存储结构。那么图是什么呢？图是一种逻辑结构，和线性结构、树形结构、集合结构一样是一种逻辑结构用来描述数据对象中的数据元素之间的关

亦凉/ 2022年05月31日 10:09/ 0 赞/ 214 阅读

相关深度优先和广度优先搜索

目录前言 1 深度优先搜索 2 广度优先搜索 3 深度优先和广度优先的比较 -------------------- 前言最近面试，被问到了深度优先和广度

冷不防/ 2022年05月22日 13:39/ 0 赞/ 286 阅读

相关 matlab中的矩阵是行优先还是列优先

分享一下我老师大神的人工智能教程！零基础，通俗易懂！[http://blog.csdn.net/jiangjunshow][http_blog.csdn.net_jiangju

爱被打了一巴掌/ 2022年04月16日 00:59/ 0 赞/ 309 阅读

相关树的深度优先和广度优先

1.深度优先算法----采用栈(非递归) ![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aH

叁歲伎倆/ 2021年09月21日 00:56/ 0 赞/ 343 阅读