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循环,用于经过矩阵所有的元素。为