【信号处理】语音采集及处理含matlab源码.zip

clear clc [pyr,fs]=audioread('a0007.wav');%声音读取 sound(pyr,fs); %声音回放 n=length(pyr); pyr1=fft(pyr,n); %快速傅里叶变换 figure(1) subplot(2,1,1); plot(pyr); xlabel('时间'); ylabel('幅度'); title('初始信号波形'); %绘出时域波 subplot(2,1,2); %绘出频域频谱 plot(abs(fftshift(pyr1))); title('初始信号频谱'); xlabel('频率'); ylabel('幅度'); noise=0.01*randn(n,1);%加噪声 s=pyr+noise; pause; sound(s,fs); %播放加噪的语音 n=length(s); %画出加噪之后，其时域频域 S=fft(s,n); figure(2) subplot(2,1,1) plot(s); title('加噪声后信号波形') xlabel('时间'); ylabel('幅度'); grid; subplot(2,1,2) plot(abs(fftshift(S))); xlabel('频率'); ylabel('幅度'); title('加噪声后信号信号频谱'); grid; pause %设计IIR低通滤波器 rp=2,rs=80; Ft=8000; Fp=1000; Fs=2000; 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域的变换 %低通滤波器特性 figure(3); [h,w]=freqz(bz,az); title('IIR低通滤波器'); plot(w*fs/(2*pi),abs(h)); %滤波 z=filter(bz,az,s); %回放滤波后的信号 pause; sound(z,fs); %滤波后画出其频谱，波形 Z=fft(z); %滤波后的信号频谱 figure(4) subplot(2,2,1); plot(z); title('低通滤波后的信号波形'); xlabel('时间'); ylabel('幅度'); grid; subplot(2,2,2) plot(s) title('加噪声后信号波形'); xlabel('时间'); ylabel('幅度'); grid; subplot(2,2,3); plot(abs(fftshift(S))); title('加噪声后信号频谱'); xlabel('频率'); ylabel('幅度'); grid; subplot(2,2,4); plot(abs(fftshift(Z))); title('低通滤波后信号的频谱'); xlabel('频率'); ylabel('幅度'); grid;