语音信号的数字滤波处理课设.docVIP

1 绪论 1.1 课题背景 数字滤波是数字信号处理的重要内容，数字滤波器可分为IIR和FIR两大类。对于IIR数字滤波器的设计，需要借助模拟原型滤波器，再将模拟滤波器转化为数字滤波器，文中采用的设计方法是脉冲响应不变法、双向性变换法和完全函数设计法；对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，文中采用的设计方法是窗函数法。根据IIR滤波器和FIR滤波器的特点，在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并对采集的语音信号进行分析，最后给出了IIR和FIR对语音滤波的效果。 1.2 课题要求 掌握数字信号处理的基本概念，基本理论和基本方法。 熟悉离散信号和系统的时域特性。 掌握序列快速傅里叶变换方法。 学会MATLAB的使用，掌握MATLAB的程序设计方法。 掌握利用MATLAB对语音信号进行频谱分析。 掌握滤波器的网络结构。 　　　　　　　　 　2 课程设计预习与原理 2.1 卷积运算的演示 2.1.1 线性卷积 序列x1(n)=[2 0 1 2 5 7 0 5 0 2 0 3]，序列x2(n)=[2 0 0 1 0 1 1 0]。动态演示序列进行线性卷积x1(n)﹡x2(n)翻转、移位、乘积、求和的过程 n=-7:18; M=17; yn=zeros(1,19); figure(1) stem(yn); xlabel(n) ylabel(y(n)) xn1=[2 0 1 2 5 7 0 5 0 2 0 3]; xm1=[zeros(1,7) xn1 zeros(1,7)];%为26个数字的矩阵 figure(2) stem(n,xm1) xlabel(m) ylabel(x1(m)) xn2=[2 0 0 1 0 1 1 0]; xm2=[fliplr(xn2) zeros(1,18)]; %移位，补零为26个数字的矩阵 figure(3) stem(n,xm2) xlabel(m) ylabel(x2(N-m)) title(n=0) yn(1)=sum((xm1.*xm2));%对xm1与xm2进行对应原素乘方之后进行数组转置，求和;即为求卷积 figure(4) stem(yn) xlabel(n) ylabel(y(n)) title(n=N) for N=1:17 xm3=[zeros(1,N) fliplr(xn2) zeros(1,M)]; figure(); stem(n,xm3) xlabel(m) ylabel(x2(N-m)) title(n=N) M=M-1; yn(N+1)=sum((xm1.*xm3)); figure() stem(yn) xlabel(n) ylabel(y(n)) title(n=N) end xm3=[zeros(1,18) fliplr(xn2)] figure() stem(xm3) xlabel(m) ylabel(x2(N-m)) title(n=N); yn(19)=sum((xm1.*xm3)); figure() stem(yn) xlabel(n) ylabel(y(n)) 2.1.2循环卷积 序列x1(n)=[2 0 1 2 5 7 0 5 0 2 0 3]，序列x2(n)=[2 0 0 1 0 1 1 0]，N=12动态演示序列进行圆周卷积x1(n)⊙x2?(n)翻转、移位、乘积、求和的过程 yn=zeros(1,12); figure(1) stem(yn) xlabel(n) ylabel(y(n)) %图1,12个0 xn1=[2 0 1 2 5 7 0 5 0 2 0 3]; figure(2) stem(n,xn1) xlabel(m) ylabel(x1(m)) xn2=[2 0 0 1 0 1 1 0]; xm2=[xn2 zeros(1,length(xn1)-length(xn2))]; figure(3) stem(n,xm2) xlabel(m) ylabel(x2(m)) title(n=0); yn(1)=sum((xn1.*xm2)); figure(4) stem(yn) xlabel(n) ylabel(y(n)) title(n=N) for N=1:11 xm1=[fliplr(xn1(1:N)) fliplr(xn1(N+1:12))]; figure() stem(n,xm1) xlabel(m) ylabel(x1(N-m)) title(n=N) yn(N)=sum((xm1.*xm2)); fig