主要内容

fircls

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

折叠所有页面

语法

B = fircls(n,f,amp,up,lo)
fircls (n、f、amp,瞧,“design_flag”)

描述

例子

b= fircls (n，f，amp，向上，罗）生成一个长度n+ 1线性相位FIR滤波器。该滤波器的频幅特性与矢量给出的频幅特性相匹配f而且amp．

fircls (n，f，amp，向上，罗，“design_flag”）使您能够为筛选器设计指定可视化显示选项。

例子

全部折叠

打开实时脚本

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

N = 150;F = [0 0.4 1];A = [10];Up = [1.02 0.01];Lo = [0.98 -0.01];B = fircls(n,f,a,up,lo，“两个”）;
界限违规= 0.0788344298966界限违规= 0.0096137744998界限违规= 0.0005681345753界限违规= 0.0000051519942界限违规= 0.0000000348656

图中包含3个轴对象。Axes对象1包含3个line类型的对象。坐标轴对象2包含3个line类型的对象。坐标轴对象3包含3个line类型的对象。

绑定违规= 0.0000000006231

边界违反表示设计收敛时过程的迭代。显示滤波器的幅值响应。

fvtool (b)

{

输入参数

全部折叠

过滤器顺序，指定为实正标量。

请注意

的fircls函数对于具有奈奎斯特频率通带的配置(即高通和带阻滤波器)总是使用偶数滤波器顺序。这是因为对于奇阶，奈奎斯特频率处的频率响应必然为0。如果你指定了一个奇值n，fircls加1。

归一化的频率点，指定为实值向量。转换频率在[0,1]范围内，其中1对应奈奎斯特频率。第一点f必须0最后一点必须是1．频率必须按递增顺序排列。

分段常数频率响应的期望幅值，指定为实值向量。的长度amp等于响应的频带数，长度(f) 1．

每个频带中频率响应的上界，指定为具有相同长度的实值向量amp．

每个频带中频率响应的下界，指定为具有相同长度的实值向量amp．

请注意

通常，阻带中的低值被指定为负值。通过设置罗等于0在停止带中，得到了一个非负的频响振幅。对这种滤波器进行频谱分解以获得最小相位滤波器。

过滤器设计显示，指定为以下之一:

  • “跟踪”-在每个迭代步骤中查看设计错误的文本显示。

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

  • “两个”-查看文本显示和情节。

输出参数

全部折叠

过滤器系数，作为长度的行向量返回n+ 1。

算法

的fircls函数使用迭代最小二乘算法来获得等纹波响应。该算法是一种多次交换算法，每次迭代都使用拉格朗日乘法器和库恩-塔克条件。

参考文献

[1]塞莱斯尼克，i.w.， M.朗，C. S.伯勒斯。无指定过渡带的FIR滤波器的约束最小二乘设计。1995年声学、语音和信号处理国际会议论文集。第2卷，1995年，第1260-1263页。

[2]塞莱斯尼克，i.w.， M.朗，C. S.伯勒斯。无指定过渡带的FIR滤波器的约束最小二乘设计。IEEE®信号处理汇刊．第44卷，第8期，1996年，第1879-1892页。

扩展功能

版本历史

R2006a之前介绍

另请参阅

||