Spectral Substraction for removing noise

Question

Amirhosein Khanlari 2020-1-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/502021-spectral-substraction-for-removing-noise

编辑： Thiago Henrique Gomes Lobato 2020-1-26

i have this code buy i get index out of range.how can i fix it

[x,fs]=audioread('NoisyTest.wav');
y=x(1:350,1);
Y=fft(y);
length(Y)
magY=abs(Y);
b=[];
for i=0:2000;
n=350;
x1=x(1+n*i:n+n*i);
X1=fft(x1);
magX=abs(X1);
S=(magX.^2-magY.^2);
S1=abs(S).^0.5;
s1=ifft(S1);
a=s1';
b=[b a];
end
x2=b';
plot(x2);
sound(x2,fs);
audiowrite('13.wav',x2,fs)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Thiago Henrique Gomes Lobato 2020-1-26

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/502021-spectral-substraction-for-removing-noise#answer_412082

编辑：Thiago Henrique Gomes Lobato 2020-1-26

在 MATLAB Online 中打开

Probably your signal is not long enough for you to use 2000 blocks of 350 samples, you can make your loop signal dependent to avoid this problem:

[x,fs]=audioread('NoisyTest.wav');
y=x(1:350,1);
Y=fft(y);
length(Y)
magY=abs(Y);
b=[];
for i=0:round(length(x)/350)-1;
n=350;
x1=x(1+n*i:n+n*i);
X1=fft(x1);
magX=abs(X1);
S=(magX.^2-magY.^2);
S1=abs(S).^0.5;
s1=ifft(S1);
a=s1';
b=[b a];
end
x2=b';
plot(x2);
sound(x2,fs);
audiowrite('13.wav',x2,fs)