我在Matlab中有一个大小为65536 x 192的数据集。如果我想沿着第二维进行一维fft,我可以执行一个for循环:
%pre-allocate ect..
for i=1:65536
F(i,:) = fft(Sig(i,:));
end或者我可以指定维度并在不使用for循环的情况下执行该操作:
F = fft(Sig,[],2);这对我的数据集来说要快20倍。
我一直在寻找类似的离散小波变换(dwt),但没有找到它。所以我想知道是否有人知道在Matlab中有什么方法可以在指定的维度上实现dwt?还是我必须使用循环?
发布于 2011-10-12 14:12:38
我猜想您使用的是小波工具箱中的函数:http://www.mathworks.co.uk/help/toolbox/wavelet/ref/dwt.html
文档似乎并不描述对数组的操作,因此可能不支持它。如果它确实允许您输入一个数组,那么它将对第一个非单例维度进行操作,或者它将忽略该形状并将其视为一个向量。
发布于 2011-10-12 13:26:47
在您的循环FFT例子中,您似乎操作在线上。Matlab使用列-主要顺序。它可以解释性能的差异。如果对列进行操作,性能是否相同?如果这是正确的解释,您可以在循环中使用dwt。
发布于 2011-10-12 14:17:39
如果你真的需要性能,一个解决方案是做你自己的MEX,用你想要的方式调用一个C离散小波变换库。
https://stackoverflow.com/questions/7740294
复制相似问题