我正在编写一个可计算某种矩阵的Matlab应用程序。 我试图用基于矢量的计算替换程序中的for循环，但是，我陷入了困境。

到目前为止，我已经弄清楚了这个简单的代码示例部分：

可以用这样的结构代替：

并给出相同的结果。 但是，如果我需要做一些小改动怎么办：

我怎样才能加快施工速度？ 可以很容易地将其向量化吗？ 我是Matlab的新手，所以我将不胜感激。 非常感谢！

一种选择是使用arrayfun。

说明：

arrayfun带有一个作用在标量输入上的函数以产生标量输出，以及要传递给该函数的输入数组。 然后，它遍历输入数组，在每个数组上运行您的函数，并将每个数组的输出放入输出矩阵(-es)的相应条目中。 因此，arrayfun基本上会执行嵌套循环，这会降低您的工作速度。

在此示例中，我还使用了匿名函数(也称为lambda函数)在最内层的循环中进行工作。 由于lambda函数必须是Matlab中的单个表达式，因此我不得不稍微重写一下内部循环(使用简单的代数操作)。 如果需要将arrayfun与不容易表示为lambda的函数一起使用，则始终可以将该函数编写在单独的.m文件中，并将其传递给arrayfun。

编辑：请注意，您不必再预先分配kernel，并且您也不需要offset。