西电信号大作业歌曲人声消除

1、信号与系统课程实践报告1内容与要求通过信号分析的方法设计一个软件或者一个仿真程序，程序的主要功能是完成对歌曲中演唱者语音的消除。试分析软件的基本设计思路、基本原理，并通过matlab程序设计语言完成设计。更进一步地，从理论和实用的角度改善软件性能的方法和措施。读入音频信号2 思路与方案 分离左右声道 右声道信号左声道信号 相减信号（频率分析）带阻滤波写出音频信号歌曲的伴奏左右声道相同，人声不同。所以通过左右声道不同处理信号，然后通过频率分析做带阻滤波滤除主要人声信号。3 成果及展示代码：clear;clc;x,fs=audioread(d:文本文档林.wav);ts=1/fs;n=length

2、(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

3、(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);so

4、und_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,

5、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(d:文本文档林_去人声.wav,sound_final,fs);1歌曲原始左右声道的幅度和频率曲线2相减得到的信号幅度和频率曲线3进行消除人声处理后信号的幅度和频率曲线4 总结与感想 在本次实践中，熟悉了matlab的操作，了解了很多命令。实现歌曲的人声消除，了解