避免for循环

http://www.360doc.com/content/10/0610/22/1217721_32411251.shtml

1. 避免使用for循环：

在Matlab中，for循环运算效率非常低，因为Matlab是一种解释语言。像矩阵乘法的宏操作与诸如增加标号的微操作差不多一样快，因为代码解释的顶层都存在于这两个情形之中。for循环应该只用于做最后的手段，并且一般用于控制运算，而不是由于计算的原因。

一般的Matlab程序中，90%的for循环都可以用等效但更为快速的向量代码代替。（一般为寻找替代函数或使用

技巧）例如，求和运算就可以使用sum函数，也可以通过将该向量乘以所有元素都为1的列向量来求得。

1. 向量化

转换一个for循环矩阵向量的操作过程可以称为向量化。有时向量化看上去似乎十分低效，因为可能它做了比

for循环更多的计算。然而，向量化后的程序运行速度要快得多，因为重复应用到这个向量的是一个简单运算。

1） 重复行与列经常需要从头至尾重复一个或多个数值来构建矩阵。如果矩阵所有的数值都要相同，则可以使用像ones和zeros这样的函数。但是假设有一个行向量x，并且要产生一个10行的矩阵，其每一行都是x的复制。这里应避免使用循环，而是使用外积矩阵乘法操作。

X = ones(10,1)*x

2） 向量逻辑操作

程序运行速率慢的一个地方是条件语句。表面上看，条件检验无法进行向量化，但实际上，在Matlab中，如大于、等于这些比较函数都能够在向量或矩阵上运算。这样，下面的Matlab代码

[1 2 3 4 5 6] <  4

将返回结果[1 1 1 0 0 0]，其中0代表假，1代表真。由后面产生一个冲激信号向量的技巧给出另一个简单例

子：

nn = [-20:80];

impulse = (nn==0); stem(nn,impulse);

3） Clip函数向量化

Clip函数用来以给定的上限和下限剪切一个输入信号。用常用的语句编写的代码如下：

function y = clip(x,lo,hi)

%Clip --- threshold large and small elements in matrix x

%======>slowest possible version <========= [M,N] = size(x);

for m = 1:M for n = 1:N if x(m,n) > hi x(m,n) = hi;

elseif x(m,n) < lo x(m,n) = lo; end,end,end

y = x;

上面代码具有双重嵌套的for循环，用于经过矩阵所有的元素。为得到速度较快的版本，必须完全放弃该循环而利用逻辑操作的向量性质。进一步，我们可以利用：true和false具有1和0的数值作为屏蔽，例如上述代码可以改写如下：

function y = clip(x,lo,hi)

%=========>fast version<=========

%(uses matrix logicals to replace loops) y = (x .* [x<=hi]) + (hi .* [x>hi]);

y = (y .* [x>=10]) + (lo .* [x<lo]);

如果使用上面的方法，算术运算的次数大于第一个版本。（可以使用etime和flops函数对这两个版本的函数进行计时）。即使是产生第二个版本需要进行10倍次数的大量计算的情形，第二个版本仍运行得比第一个版本快上10 倍。

4）“：”算子

“：”符号可以通过给出起始序号、步长以及结尾序号来产生标号范围。因此，规则相间的整数（或实数）向量可经由下式得到

iii = start:skip:end

另外，“：”可以与矩阵结合起来操作，对矩阵A，A(2:5,1:3)分出一个4*3的子矩阵。而A（：）产生一个列向量，该向量正是A的连接在一起的列。矩阵A的更一般的“再成形”可以用reshape(A,M,N)函数实现。

5）经常使用help指令来查询函数的帮助信息，而用type指令查看函数内容