基于MATLAB软件的语音采集与处理.rar

function sheji61=sheji61 fs=8000; %抽样频率 T=1/fs; duration=3; [y,fs]=audioread('设计六.wav'); %原始语音信号 sound(y,fs); N=length(y); %信号长度 t=(1:N)/fs;df=fs/N;n1=1:N/2;f=(n1-1)*df; fs=8000;fp=1000;fst=1200; %通带截止频率 ap=1;ast=100;wp=2*fp/fs;wst=2*fst/fs; [N,wc]=cheb1ord(wp,wst,ap,ast); [num,den]=cheby1(N,ap,wc, 'low' ); [h,w]=freqz(num,den); figure(1) plot(w/pi,abs(h)); %低通滤波器的幅频响应 grid;xlabel( '\omega/\pi' );ylabel( ' 振幅（幅值） ' ); title( ' 低通滤波器的幅频响应' ); y0=filter(num,den,y);fftwav=fft(y0); figure(2) subplot(2,2,1); plot(y); %原始语音信号 title( ' 原始信号的波形' ) xlabel( ' 时间/t' );ylabel( ' 幅值' ) subplot(2,2,3); plot(y0); %经过低通滤波器后的语音信号 title( ' 通过低通滤波器后的信号图形' ) xlabel( ' 时间/t' );ylabel( ' 幅值' ) subplot(2,2,2);q=fft(y); plot(f,20*log10(abs(q(n1)))); %滤波前的信号频谱 title( ' 滤波前的信号频谱图' ); xlabel( ' 频率/Hz' );ylabel( ' 幅值/db' ) subplot(2,2,4); plot(f,20*log10(abs(fftwav(n1)))); %滤波后信号的频谱 title( ' 滤波后信号的频谱图' ); xlabel( ' 频率/Hz' );ylabel( ' 幅值/db' ) audiowrite( '存储.wav',y0,fs ); IIRchoose end