fft 失败的噪声消除 - 无法分配元素

Question

我有以下降噪代码：

[z,fs] = audioread('noisy_voices.wav'); % Use >>soundsc(z, fs) to hear the unprocessed signal
zproc_vec=zeros(1,length(z));
tail = zeros(1,256);
for k = 0:128:length(z)-256
Z = fft(z(k +1:k + 256).* hann(256));
[zmax, zl] = max(abs(Z(1:128)));
Z(zl-3: zl +3)=0;
Z(256-(zl-3:zl +3)+2)=0;
zproc = ifft(Z);
zproc = zproc+tail;
tail(1:128) = zproc(129:256);
zproc_vec(k+1:k+256)=zproc;
end
soundsc(zproc_vec , fs)

谁能告诉我为什么会出现这个错误？

由于左右两边的元素数量不同，无法进行赋值

task_one 中的错误（第 12 行） zproc_vec(k+1:k+256)=zproc;

Answer 1

我认为您的Z = fft( ___ ) 行的输出将是一个列向量，但您将tail 初始化为带有tail = zeros(1,256); 的行向量

所以在这一行：

zproc = zproc+tail;

Implicit expansion 将使 zproc 成为一个 256*256 的正方形矩阵。

然后您的索引将失败，如错误消息所指定，因为您尝试将此方阵分配给您使用 zproc_vec(k+1:k+256) 进行索引的 256 个元素。

将tail 初始化为列向量应该可以解决问题。

或者，您可以采取“懒惰”的方式，并确保您只对列向量进行操作以创建 zproc

zproc = zproc(:)+tail(:); % Make both terms column vectors for addition

Answer 2

''无法进行赋值，因为左右两边的元素数量不同''

您不理解错误消息的哪一部分？这意味着等号左边和右边的变量大小不同，所以你不能把右边的东西分配给左边的东西。

检查尺寸并确保它们相同。