在MATLAB中介绍矢量化 - 任何好的教程？

Question

我正在寻找关于MATLAB中矢量化(循环)的任何好教程.

我有非常简单的算法,但它使用两个for循环.我知道向它进行矢量化应该很简单,我想学习如何做而不是要求你提供解决方案.

但是为了让你知道我有什么问题,所以你可以建议最好的教程,展示如何解决类似的问题,这是我的问题的大纲:

B = zeros(size(A));    % //A is a given matrix.
for i=1:size(A,1)
   for j=1:size(A,2)
      H = ... %// take some surrounding elements of the element at position (i,j) (i.e. using mask 3x3 elements)
      B(i,j) = computeSth(H); %// compute something on selected elements and place it in B
   end
end

所以,我不是要求解决方案.我要求一个很好的教程,在MATLAB中矢量化循环的例子.我想学习如何做,并自己做.

Answer 1

以下是我经常链接到的几个MathWorks教程作为主题的参考:

• 代码矢量化指南
• 提高性能的技术:矢量化循环

这是Loren的博客文章之一,它针对特定的示例问题进行了很好的代码矢量化演练:

• 加快MATLAB应用程序的速度

您作为样本提供的特定类型的问题(涉及处理给定矩阵的子矩阵)可以以不同的方式进行矢量化,具体取决于您正在执行的操作类型.您可以使用CONV2或FILTER2代替嵌套for循环.图像处理工具箱中还有许多函数可以处理矩阵的邻域和块处理,例如NLFILTER和BLOCKPROC.这些函数的文档应该可以帮助您弄清楚如何将它们用作向量化代码的方法.