实验九-音频信号采集及处理.pdf

音频信号采集及处理 程序代码及实验结果图：voice,fs=audioread(notify.wav);%声音读取sound(voice,fs);%声音回放n=length(voice);%计算长度voice1=fft(voice,n);%快速傅里叶变换figure(1);subplot(2,1,1);plot(voice);%绘出时域波xlabel(t);ylabel(amp);%坐标名称title(初始音频信号时域波形);grid on;subplot(2,1,2);plot(abs(fftshift(voice1);%绘出原始音频信号频谱title(初始音频信号频域波形);xlabel(f);ylabel(amp);grid on;t=0:1/fs:(n-1)/fs;noise=0.05*sin(2*pi*100000*t);%100kHz 正弦波噪声s=voice+noise;%加噪后的音频信号pause;sound(s,fs);%播放加噪的语音n=length(s);S=fft(s,n);%计算频谱figure(2);subplot(2,1,1)；plot(s);%画出加噪之后的音频信号时域波形title(加噪声后的音频信号时域波形)；xlabel(t);ylabel(amp);grid on;subplot(2,1,2)；plot(abs(fftshift(S);%零频移到频谱中心后，绘制加噪之后的音频信号频谱xlabel(f);ylabel(amp);title(加噪声后的音频信号频域波形);grid on;pause;rp=2;rs=80;Ft=8000;Fp=1000;Fs=1300;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;%求出待设计的模拟滤波器的边界频率n,wn=buttord(wp,ws,rp,rs,s);%低通滤波器的阶数和截止频率b,a=butter(n,wn,s);%S 域频率响应的参数即：滤波器的传输函数bz,az=bilinear(b,a,0.5);%利用双线性变换实现频率响应 S 域到 Z 域的变换h,w=freqz(bz,az);figure(3);plot(w*fs/(2*pi),abs(h);%绘制 IIR 低通滤波器特性曲线title(IIR 低通滤波器特性曲线);grid on;z=filter(bz,az,s);%滤波pause;sound(z,fs);%回放滤波后的信号 Z=fft(z);%滤波后的信号频谱figure(4);subplot(2,2,2);plot(z);%绘制低通滤波后的音频信号时域波形title(低通滤波后的音频信号时域波形);xlabel(t);ylabel(amp);grid on;subplot(2,2,1);plot(s);%绘制加噪声后的音频信号时域波形title(加噪声后的音频信号时域波形);xlabel(t);ylabel(amp);grid on;subplot(2,2,3);plot(abs(fftshift(S);%绘制加噪声后的音频信号频域波形title(加噪声后的音频信号频域波形);xlabel(f);ylabel(amp);grid on;subplot(2,2,4);plot(abs(fftshift(Z);%绘制低通滤波后的音频信号频域波形title(低通滤波后的音频信号频域波形);xlabel(f);ylabel(amp);grid on;