fs=20000; y1=wavread('cctv.wav'); y2=fft(y1,1024); [n m]=size(y1); a=(linspace(0,2,n))'; b=zeros(n,1); t=[a b]'; noise=(0.5*sin(10000*pi*t))'; y3=y1+noise; y4=fft(y3,1024); [b,a]=butter(5,0.3,'high'); y5=filter(b,a,y3); y6=fft(y5); sound(y1,20000) figure(1) subplot(2,1,1),plot(y1) title('原始语音信号波形'); subplot(2,1,2),plot(abs(y2)) title('原始音频信号的频谱') figure(2) subplot(2,1,1),plot(y3) title('加噪音后的波形') subplot(2,1,2),plot(abs(y4)) title('加噪音后的频谱') figure(3) subplot(2,1,1),plot(y5) title('加噪音过滤后的波形') subplot(2,1,2),plot(abs(y6)) title('加噪音过滤后的频谱') %语音信号采样频率为 20000 %读取语音信号的数据，赋给变量 y1 %对原始语音信号进行 FFT 变换 %生成 5000HZ 的正弦噪音 %对原始信号加上噪音 %对加上噪声后的信号进行 FFT 变换 %对加噪信号进行 IIR 滤波 %对加噪信号进行 FFT 变换 %播放音乐 %绘制原始语音信号的波形 %绘制原始语音信号的频谱 %绘制加噪音后的波形 %绘制加噪音后的频谱 %绘制加噪音过滤后的波形 %绘制加噪音过滤后的频谱