7

（零）填充matlab数组的最简单方法是什么？
即给定[1,2,3,4]和长度6返回[1,2,3,4,0,0]

背景

我有一个数据数组，我想在对数据运行 fft 之前应用一个窗口函数。

我曾经将数据直接传递给 fft，这会将零填充到 2 的下一个幂，但现在我需要在 fft 之前将它的零填充，这样我就可以乘以窗口函数。

fs = 100;                          % Sample frequency (Hz)
t = 0:1/fs:10-1/fs;                % 10 sec sample
x = (1.3)*sin(2*pi*15*t) ...       % 15 Hz component
  + (1.7)*sin(2*pi*40*(t-2)) ...   % 40 Hz component
  + (2.5)*randn(size(t));          % Gaussian noise;

m = length(x);          % Window length
n = pow2(nextpow2(m));  % Transform length
w = barthannwin( n );   % FFT Window

y = fft(data, n);       % DFT

windowed_data = x*w ; % Dimensions do not match as x not padded
y = fft(windowed_data, n); % DFT

我知道padarray作为图像处理工具箱的一部分，我没有。

4

4 回答 4

9

编辑

这种方法可能对向量更好，因为它在转置时不会中断，请注意它会改变原始向量，这可能是不可取的：

myVec = 1:7;
myVec(end+3)=0

或者，您可以连接零和您拥有的向量并用它创建一个新变量。

myVec = 1:7;
requiredpadding = 10-7;
myVecPadded=[myVec zeros(1,requiredpadding)]
于 2013-02-15T16:28:12.023 回答
3

没有内置函数来进行填充，但这里有一个小函数可以在x给定最小长度的情况下填充向量n。

function y = pad(x, n)
y = x;
if length(x) < n
    y(n) = 0;
end
于 2013-02-15T17:50:27.847 回答
1

对于数组，这应该用零填充到最接近的 2 次幂a：

a(2^ceil(log2(length(a))))=0;

于 2013-02-15T17:45:21.823 回答
0

Matlab的图像处理工具箱有一个内置函数来填充数组padarray(A,dim,value

对于您的示例：

A = [1, 2, 3, 4];
dimension = [0 1]; % pad with extra columns
size = 2; % how much to pad
B = padarray(A,size*dimension,'post') % 'post' says to pad at the end
% just for demonstration, let's pre-pad the first dimension (rows)
dimension = [1 0];
C = padarray(A,dimension,'pre') % just as an example
% or pad in both directions
dimension = [1 2];
D = padarray(A,dimension) % by default, it will pad both pre and post

返回

B =

 1     2     3     4     0     0

C =

 0     0     0     0
 1     2     3     4
D =

 0     0     0     0     0     0     0     0
 0     0     1     2     3     4     0     0
 0     0     0     0     0     0     0     0

也可以用于多维数组，只需要对维向量进行相应的扩展，即dimension = [0 0 1]在第三维上扩展。

于 2014-03-16T19:01:01.867 回答