matlab歌曲人声消除

[X,fs]=audioread('F:\ÉñÌﴨƬ¶Î.wav');ts=1/fs;N=length(X)-1;t=0:1/fs:N/fs;Nfft=N;df=fs/Nfft;fk=(-Nfft/2:Nfft/2-1)*df;a1=1;a2=-1;b1=1;b2=-1;%·ÖÀë×óÉùµÀºÍÓÒÉùµÀSoundLeft=X(:,1);SoundRight=X(:,2);%¶Ô×óÉùµÀºÍÓÒÉùµÀ½øÐпìËÙ¸µÀïÒ¶±ä»»SoundLeft_f=ts*fftshift(fft(SoundLeft,N));SoundRight_f=ts*fftshift(fft(SoundRight,N));%ÏÔʾ×óÓÒÉùµÀ·ù¶È±ä»¯figure(1)subplot(411)plot(t,SoundLeft);subplot(412)plot(t,SoundRight);%ÏÔʾ×óÓÒÉùµÀƵÂʱ仯subplot(413)f_range=[-5000,5000,0,0.1];plot(fk,SoundLeft_f);axis(f_range);subplot(414)plot(fk,SoundRight_f);axis(f_range);NewLeft=a1*SoundLeft+a2*SoundRight;NewRight=b1*SoundLeft+b2*SoundRight;Sound(:,1)=NewLeft;Sound(:,2)=NewRight;Sound_Left_f=ts*fftshift(fft(NewLeft,N));Sound_Right_f=ts*fftshift(fft(NewRight,N));figure(2)subplot(411)plot(t,NewLeft);subplot(412)plot(t,NewRight);f_range=[-5000,5000,0,0.1];subplot(413)plot(fk,Sound_Left_f);axis(f_range);subplot(414)plot(fk,Sound_Right_f);axis(f_range);BP=fir1(300,[800,2200]/(fs/2));%¸ù¾Ý×óÓÒÉùµÀ²îÒì½øÐÐÂ˲¨¡¾800,2200¡¿HzCutDown=filter(BP,1,Sound);Sound_Final=Sound-0.6*abs(CutDown);Sound_Final_f=ts*fftshift(fft(Sound_Final,N));figure(3)subplot(211)plot(t,Sound_Final);subplot(212)f_range=[-5000,5000,0,0.1];plot(fk,Sound_Final_f);axis(f_range);audiowrite('F:\¸èÇúÏûÒô.wav',Sound_Final,fs);