[信息与通信]数字信号处理-FIR系统设计上.pptVIP

FIR系统与理想系统间的误差 总平方误差： FIR系统与理想系统间的误差 总平方误差可以通过系统阶数进行控制，对于给定的阶数，此类滤波器具有最小平方积分误差（LS优化滤波器）； FIR系统与理想系统间的误差 频域误差的来源分析 时间窗口与理想冲激响应乘积 窗口频谱与理想频率响应卷积 误差来自于窗口频谱的卷积效应： 主瓣宽度决定过渡带宽度，与阶数有关； 主瓣与旁瓣的相对面积比决定带内纹波幅度，与阶数无关。 频域误差的来源分析 窗口法设计：频域乘积 窗口频谱自乘，等效于时域窗口卷积，窗口由矩形窗口改变为三角窗口； 卷积会导致窗口宽度增加，滤波器阶数增加； 保持滤波器阶数不变时，窗口频谱的主瓣旁瓣宽度加倍(过渡带加宽），旁瓣相对高度下降。 窗口法设计：频域乘积 N阶三角窗口的函数表达： 三角窗 triang(N) 巴特利特窗 bartlett(N) 窗口法设计：频域乘积 窗口频域乘积法的推广 时域采用m个宽度为N/m的矩形窗口进行卷积； 频域表现为m个宽度m倍的sinc函数乘积； 效果： 主瓣与旁瓣变宽m倍，旁瓣衰减加大为m倍。 利用sinc函数的周期与对称关系，使用平移的频谱叠加来抵消旁瓣的面积！ 窗口法设计：频域调制 频谱叠加使第一旁瓣并入主瓣（主瓣加宽），其余旁瓣宽度不变，幅度相互抵消，面积削减效果明显。 窗口法设计：频域调制 该类窗口的频谱可以看作是矩形窗口频谱与三个冲激频谱卷积的结果。 汉明窗 汉宁窗 hanning(N) hamming(N) 窗口法设计：频域调制 窗口法设计：频域调制 窗口法设计：频域调制 Hamming窗口设计仿真 Hanning窗：主瓣宽度 旁瓣衰减 -31dB Hamming窗：主瓣宽度 旁瓣衰减 -41dB Blackman 窗：主瓣宽度 旁瓣衰减 -57dB 窗口法设计：频域调制 常用固定窗口设计的结果对比 窗口类型 阻带衰减 通带插损 最小阶数 矩形 21 0.811 0.91 /d 汉宁 44 0.055 3.32 /d 汉明 55 0.0155 3.44 /d 布莱克曼 75 0.0015 5.98 /d 上表中阻带衰减和通带插损均以dB为单位，d为相对过渡带宽度。 窗口法设计的MATLAB函数 fir1(n,f, ‘标识’ , window) n阶FIR滤波器，截止频率f； 输入参数： 阶数n，截止频率f， 滤波器类型，窗口类型； 输出：滤波器冲激响应序列。 阶数n：设计低通和带通时可取任意正整数，设计高通和带阻时只能取正偶数。 截止频率f： f的取值区间为[0,1]; 对低通和高通，f为单一值； 对带通和带阻，f为2元素向量f=[f1 f2]; 对多带滤波器，f为多元素向量。 窗口法设计的MATLAB函数 滤波器标识： 缺省时表达低通、带通或最低带为阻带； high表达高通或最高带为通带； stop表达带阻或最低带为通带； 窗口法设计的MATLAB函数 窗口类型：可以直接使用已有的窗口函数，也可以自行编制相关的函数。 矩形窗 boxcar(N) 三角窗 triang(N) 巴特利特窗 bartlett(N) 汉明窗 hamming(N) 汉宁窗hanning(N) 布莱克曼窗blackman(N) 凯泽窗 kaiser(N，b) Dolph-Chebyshev窗 chebwin(N，r) 窗函数中的参数N为窗口序列长度 窗口法设计的MATLAB函数 窗口法设计的MATLAB仿真程序 N=41;n=-5:N+4;x=boxcar(N); x1=[0 0 0 0 0 x 0 0 0 0 0]; subplot(2,2,1),stem(n,x1),title(窗口函数); axis([-5,N+4,-0.2,1.2]);h=fir1(N-1,1/pi ,x); h1=[0 0 0 0 0 h 0 0 0 0 0]; subplot(2,2,3),stem(n,h1),title(冲激响应); axis([-5,N+4,-0.2,0.4]);[H1,w1]=freqz(x,1,40000); subplot(2,2,2),plot(w1,abs(H1)),title(窗口幅频特性); [H2,w2]=freqz(h,1,40000); subplot(2,2,4),plot(w2,abs(H2)),title(滤波器幅频特性); 矩形窗口设计仿真 FIR滤波