主要内容

fircls

约束最小二乘FIR多频带滤波器设计

语法

b = fircls (n、f、amp, lo)
fircls (n、f, amp,,,”design_flag”)

描述

b = fircls (n、f、amp, lo)生成一个长度n + 1线性相位FIR滤波器b.该滤波器的频率-幅度特性与向量给出的特性相匹配famp

  • f为转换频率的向量,取值范围为0 ~ 1,其中1为奈奎斯特频率。第一点f必须0最后一点1.频率点必须按递增顺序排列。

  • amp是描述频率响应的分段常数期望振幅的向量。的长度amp等于响应中的频带数应该等于长度(f) 1

  • 向上向量的长度和amp.它们定义了每个频带的频率响应的上下限。

fircls对于通频带在奈奎斯特频率(即高通滤波器和带阻滤波器)的配置,总是使用一个均匀的滤波器顺序。这是因为对于奇阶,奈奎斯特频率处的频率响应必然为0。如果你指定一个奇值nfircls把它加1。

fircls (n、f, amp,,,”design_flag”)使您可以监视过滤器的设计,其中design_flag可以

  • “跟踪”,以文本形式显示每个迭代步骤的设计错误。

  • “阴谋”,用于显示滤波器的全频带幅度响应和每个子频带幅度响应的放大视图的图集合。在每个迭代步骤更新所有的图。图上的O是新迭代的估计极值,X是前一次迭代的估计极值,其中的极值是滤波器波纹的峰值(最大值和最小值)。只有具有相应的O和X的波纹才相等。

  • “两个”,用于文本显示和情节。

请注意

通常,阻带中较低的值将被指定为负值。通过设置等于0在阻带中,可以得到一个非负的频率响应幅值。这样的滤波器可以被频谱分解得到最小相位滤波器。

例子

全部折叠

设计一个150阶低通滤波器,归一化截止频率为 0 4 π rad /样品。指定通带的最大绝对误差为0.02,阻带的最大绝对误差为0.01。显示带图。

n = 150;F = [0 0.4 1];A = [1 0];涨= [1.02 0.01];Lo = [0.98 -0.01];b = fircls (n, f,,,“两个”);
被绑定的违规= 0.0788344298966被绑定的违规= 0.0096137744998被绑定的违规= 0.0005681345753被绑定的违规= 0.0000051519942被绑定的违规= 0.0000000348656

图中包含3个轴。axis 1包含3个类型为line的对象。axis 2包含3个类型为line的对象。axis 3包含3个类型为line的对象。

绑定违例= 0.0000000006231

界限违背表示设计收敛时过程的迭代。显示滤波器的幅度响应。

fvtool (b)

图过滤器可视化工具-幅度响应(dB)包含一个轴和其他类型的uitoolbar, uimenu对象。标题为“大小响应(dB)”的轴包含一个类型为line的对象。

算法

fircls采用迭代最小二乘算法得到等纹波响应。该算法是一种多重交换算法,每次迭代使用拉格朗日乘子和库恩-塔克条件。

参考文献

Selesnick, i.w., M. Lang和c.s. Burrus。无指定过渡带的FIR滤波器的约束最小二乘设计。1995年声学、语音和信号处理国际会议论文集。1995年第2卷,第1260-1263页。

Selesnick, i.w., M. Lang和c.s. Burrus。无指定过渡带的FIR滤波器的约束最小二乘设计。IEEE®信号处理汇刊.1996年第44卷第8期,1879-1892页。

扩展功能

另请参阅

||

之前介绍过的R2006a