Quick vectorization question (for loop)

1 次查看（过去 30 天）
Tanmay
Tanmay 2012-9-27
Hi
This loop takes about 0.085 seconds to execute and needs to be executed several thousand times. Is there a way you think I can vectorize this? I really appreciate any help!
wsize = 20;
vin = rand(1200,1);
for idx = wsize:length(vin)
mstd(idx-wsize+1, :) = std(vin(idx-wsize+1:idx));
end
PS: I don't have Parallel Processing toolbox.

请先登录，再回答此问题。

采纳的回答

Matt Fig
Matt Fig 2012-9-27
编辑：Matt Fig 2012-9-27
I hope you are pre-allocating your mstd vector. Other than that, this is about 6 times faster. I didn't compare to stdfilt.
for idx = wsize:length(vin)
x = vin(idx-wsize+1:idx);
x = x-sum(x)./ wsize;
mstd(idx-wsize+1, :) = sqrt(sum(x.^2)/(wsize-1));
end
.
.
.
EDIT add comparison based on IA's recommendation
I find the above to not be faster than using stdfilt when repeatedly run. Here is what I used for comparison:
mstd2 = stdfilt(vin,[0;ones(wsize,1)]);
mstd2 = mstd2(wsize/2:end-wsize/2);
There might be a better way to do this with stdfilt, I have no real experience with the function and had to figure it out based on trial and error and doc reading... Perhaps IA can comment.

更多回答（1 个）

Image Analyst
Image Analyst 2012-9-27
Do you have the Image Processing Toolbox? it looks like your code is a sliding window of a standard deviation filter. This is done by the function stdfilt() in the Image Processing Toolbox. Alternatively (if you don't have the Image Processing Toolbox) you could use an anonymous function and blockproc (I have a demo for that if you need it).

类别

在 Help Center 和 File Exchange 中查找有关 Poisson Distribution 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by