《matlab函数学习》-精选·课件.ppt

几种常用的窗函数 矩形窗(Boxcar) 三角形窗(Bartlett) 汉宁窗(Hanning) 汉明窗(Hamming) 布拉克曼(Blackman) 凯萨窗(Kaiser) 各种窗的比较 Name 近似过渡带 准确过渡带 最小阻带衰减 Boxcar 4?/N 1.8?/N 13dB Bartlett 8?/N 6.1?/N 27dB Hanning 8?/N 6.2?/N 32dB Hamming 8?/N 6.6?/N 43dB Blackman 12?/N 11?/N 58dB 窗函数法的设计步骤 给定理想频响函数 ； 根据指标选择窗函数。确定窗函数类型的主要依据是过度带宽和阻带最小衰耗的指标，确定N； 由 求 ，加窗得h(n)=hd(n)·w(n) 检验 频率取样法设计 频率抽样法是指定离散的理想频率响应Hd(k)抽样值，通过内插求得H(z) ，因而是插值法逼近 设计步骤： 给定理想频响 确定采样点数，对理想频响采样得 计算 ，即得FIR系统函数 用Matlab设计FIR滤波器 窗函数 bartlett（三角窗） blackman（布莱克曼窗） boxcar（矩形窗） hamming（哈明窗） hann（汉宁窗） triang（三角窗） chebwin（切比雪夫窗） kaiser（凯赛窗） help window FIR设计函数 fir1调用格式（默认为Hamming窗） b = fir1(N,Wn); %带通, b = fir1(N,Wn,’high’); % 高通 b = fir1(N,Wn,’stop’); % 带阻 N：阶次，滤波器长度为N＋1； Wn：通带截止频率，其值在0～1之间，1对应 Fs/2 b: 滤波器系数。 对于格式(1),若Wn为标量，则设计低通滤波器，若Wn是1×2的向量，则用来设计带通滤波器，若Wn是1×L的向量，则可用来设计L带滤波器。这时，格式（1）要改为: b = fir1(N,Wn,’DC-1’); %第一个带为通带 b = fir1(N,Wn,’DC-0’); % 第一个带为阻带 fir2可设计具有任意幅频响应的FIR 数字滤波器（默认为Hamming窗） b = fir2(N,F,M); N：阶次，滤波器长度为N＋1； F：频率向量，其值在0～1之间 M: 是和F相对应的所希望的幅频响应。 相关FIR设计函数 remez.m：设计Chebyshev最佳一致逼近FIR滤波器、Hilbert变换器和差分器 remezord.m：用来确定在用Chebyshev最佳一致逼近设计FIR滤波器时所需要的滤波器阶次 firls.m：用最小平方法设计线性相位FIR滤波器，可设计任意给定的理想幅频响应 fircls.m：用带约束的最小平方法设计线性相位FIR滤波器，可设计任意给定的理想幅频响应 fircls1.m：用带约束的最小平方方法设计线性相位FIR低通和高通滤波器 sgolay.m：用来设计 Savitzky-Golay FIR 平滑滤波器 firrcos.m：用来设计低通线性相位FIR滤波器，其过渡带为余弦函数形状 End of This Lecture Thank! 信息科学与工程学院 宁波大学 Ningbo University 数字信号处理与MATLAB /jiangzhidi/mediasystem /mediasystem 主讲：蒋志迪 宁波大学信息科学与工程学院 办公室：曹光彪信息楼424 电话：665929 email：bighandjzd@126.com 多媒体信息处理系统设计 Matlab中的信号表示 常见连续信号 离散信号 离散卷积 LTI系统 常见连续信号 指数信号Aeat???????????????y = A*exp(a*t); 指数序列ak??????????????????幂运算a.^k实现 正弦型信号??????????????????内部函数cos( ) 和sin( ) 抽样函数Sa(t)????????????? sinc(t) 矩形脉冲信号?????????????? y = rectpuls(t,width) 三角波脉冲信号?????????? y = tripuls(t, width,skew) 周期信号：正弦信号，周期方波 ex1. 产生一个幅度为2，频率为4Hz,相位为 的正弦信号 A=2; f=4; phi=pi/6;