matlab滤除音乐中的背景杂音

[yFsbits]=wavread（‘music.wav‘）;%读出信号，采样率和采样位数。 
sigLength=length（y）;
Y = fft（ysigLength）;
%Pyy = Y.* conj（Y） / sigLength;% conj函数，用于求一个复数的复共轭
halflength=floor（sigLength/2）; % 根据模板大小确定要扩展的边界
f=Fs*（1:sigLength）/sigLength;
stor=Y（1）;
Y（1）=0;
absY=abs（Y）;
figure（1）
subplot（221）;plot（fabsY）;xlabel（‘Frequency（Hz）‘）;grid on
axis（[0f（end）/20500]）;
t=（0:sigLength-1）/Fs;
subplot（2