语音信号处理及数字滤波器设计报告.docVIP

信号课程设计报告书 ---语音信号处理 学院：信息工程学院 专业班级：通信工程 081 姓名：XX 学号：XXXXXX 1、信号的抽样和恢复 正弦波2KHz抽样输出 正弦波2KHz抽样恢复 正弦波8KHz抽样输出 正弦波8KHz抽样恢复 正弦波16KHz抽样输出 正弦波16KHz抽样恢复 三角波2KHz抽样输出 三角波2KHz抽样恢复 三角波8KHz抽样输出 三角波8KHz抽样恢复 三角波16KHz抽样输出 三角波16KHz抽样恢复 巴特沃斯低通滤波器图形 低通李沙育图形 巴特沃斯高通滤波器图形 高通李沙育图形 巴特沃斯带通滤波器图形 带通李沙育图形 2、对连续信号进行抽样恢复 t=0:0.001:6; x=cos(5*pi*t); subplot(3,1,1); plot(t,x); axis([0,1,-1.1,1.1]); title(原有x); xlabel(t);ylabel(x(t)) fs=2; %采样频率 ts=1/fs; tp=1; n=0:tp/ts; xn=cos(n*ts*0.5*pi); subplot(3,1,2) stem(n*ts,xn); title(fs=2HZ取样后的x);xlabel(t);ylabel(x(nT)) dt=ts/50; %每个采样间隔上取点数 t1=0:dt:6; %生成序列t1 tp=1; %重构时间区间为1s n=0:tp/ts; tmn=ones(length(n),1)*t1-n*ts*ones(1,length(t1)); xr1=sinc(fs*tmn); x2=xn*xr1; subplot(3,1,3) plot(t1,x2); axis([0,1,-1.1,1.1]); x0=cos(0.5*pi*t1); emax=max(abs(x0-x2)) emax = 1.9492 title(用抽样内插函数恢复的结果); 3、语音信号抽样恢复比较 [x,fs]=wavread(xiang1.wav); subplot(2,2,1) plot(x); title(语音信号原有16k采样波形); x1=fft(x); subplot(2,2,2);plot(abs(x1)); title(语音信号原有16k采样波形频谱); n=1:length(x); n2=1:0.25:length(x); yi=interp1(n,x,n2); sound(yi) subplot(2,2,3) plot(n2,yi); title(内差恢复波形); y2=fft(yi); subplot(2,2,4);plot(abs(y2)); title(语音信号内差恢复波形频谱); figure(2) y2=decimate(yi,2); subplot(2,2,1) plot(y2); title(8k采样波形) y21=fft(y2); subplot(2,2,2);plot(abs(y21)); title(8k采样波形频谱); sound(y21) y3=decimate(yi,4); sound(y3); subplot(2,2,3) plot(y3); title(4k采样波形); y31=fft(y3); subplot(2,2,4);plot(abs(y31)); title(4k采样波形频谱); sound(y31) 4、模拟及数字滤波器设计 N=4;fc=22000;fs=160000;T=1/fs;wc=2*pi*fc; [B,A]=butter(N,wc,s); [mag0,w0]=freqs(B,A); subplot(2,1,1);plot(w0/1000/2/pi,abs(mag0)); title(模拟巴特沃思低通滤波器幅频特性’); xl