数字信号处理第二次讨论课.pptx

数字信号处理第二次讨论课14-通信工程4班目录Part FivePart OnePart TwoPart ThreePart FourPart Six0105Part Seven如何利用数字信号处理知识实现模拟信号采集线性相位滤波器的特点（时域、频域）如何利用窗函数设计FIR滤波器理想低通滤波器和实际低通滤波器的特点如何利用利用数字信号处理知识实现信号的滤波如何利用数字信号处理知识实现信号的时域、频域分析吉布斯效应及改进措施0206070304CONTENTS01Part One线性相位滤波器的特点第一类线性相位特性：满足群延时为常数即可称为线性相位。第二类线性相位特性：线性FIR滤波器的时域约束条件即幅度特性的特点如果要求单位脉冲响应h(n)、长度为N的FIR数字滤波器具有第一类线性相位特性，则h(n)应当关于n=(N-1)/2点偶对称。h(n)=h(N-1-n) n从0到N-1线性FIR滤波器的时域约束条件如果要求单位脉冲响应h(n)、长度为N的FIR数字滤波器具有第二类线性相位特性，则h(n)应当关于n=(N-1)/2点奇对称。h(n)=-h(N-1-n) n从0到N-1线性FIR滤波器的时域约束条件如果要求单位脉冲响应h(n)、长度为N的FIR数字滤波器具有第一类线性相位特性，则h(n)应当关于n=(N-1)/2点偶对称。h(n)=h(N-1-n),N为奇数代入时域约束条件h(n)=h(N-1-n)和和得到：代入因为余弦函数关于w=0、π、2π三点偶对称，所以可以看出幅度特性关于这三点偶对称；因此情况1可以实现低通、高通、带通、带阻滤波器h1=[1,2,3,4,5,4,3,2,1]function [A,w,type,tao]=amplres(h)N=length(h);tao=(N-1)/2; L=floor(tao); %求滤波器阶次及符幅特性的阶次n=1:L+1;w=[0:511]*2*pi/512; % 取滤波器频率向量if all(abs(h(n)-h(N-n+1))<1e-8) % 判断滤波器系数若为对称 A=2*h(n)*cos(((N+1)/2-n)'*w)-mod(N,2)*h(L+1); % 对称条件下计算A（两种类型） % 在N=奇数时h(L+1)项存在，N为偶数时，要取消这项，故乘以mod(N,2) type =2-mod(N,2); % 判断并给出类型elseif all(abs(h(n)+h(N-n+1))<1e-8)&(h(L+1)*mod(N,2)<=1e-8) % 系数若为反对称 % 在N=奇数时要求h(L+1)项必须为零，在N=偶数时不要这条件，故乘以mod(N,2)。 A=2*h(n)*sin(((N+1)/2-n)'*w); % 反对称条件下计算A的公式（两种类型相同） type =4-mod(N,2); % 判断并给出类型else error('错误：这不是线性相位滤波器!') % 滤波器系数非对称，报告错误endh(n)=h(N-1-n),N为偶数同理，可得其中余弦函数关于过零点奇对称，关于w=0和2π偶对称。所以幅频特性关于w=0和2π偶对称。因此，情况2不能实现高通和带阻滤波器。h1=[1,2,3,4,4,3,2,1]h(n)=-h(N-1-n),N为奇数代入时域约束条件h(n)=h(N-1-n)和和代入考虑h[(N-1)/2]=0得到：因为幅度特性函数关于w=0、π、2π三点奇对称，所以可以看出幅度特性关于这三点偶对称；因此情况3只能实现带通滤波器。h1=[1,2,3,4,0,-4,-3,-2,-1]h(n)=-(N-1-n),N为偶数同理，可得幅度特性关于w=0和2π两点奇对称，关于w=π偶对称。所以幅频特性关于w=0和2π偶对称。因此，情况4不能实现低通和带阻滤波器。h1=[1,2,3,4,-4,-3,-2,-1]CONTENTS02Part Four理性低通滤波器与实际低通滤波器理想低通滤波器特性相频特性的频率斜率为常值能让零频到截止频率fc之间的所有信号都能完全通过。而让高于截止频率fc的所有信号都消失。但理想低滤波器只能是理想。相频特性幅频特性CONTENTS03Part One吉布斯效应及改进措施吉布斯效应 吉布斯效应：将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯效应。是什么？ 在用窗函数法设计FIR滤波器时，实际遇到的是与上面相反