主要内容gydF4y2Ba

tf2sosgydF4y2Ba

将数字滤波器传递函数数据转换为二阶分段形式gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

[gydF4y2Ba紧急求救信号gydF4y2Ba,gydF4y2BaggydF4y2Ba) = tf2sos (gydF4y2BabgydF4y2Ba,gydF4y2Ba一个gydF4y2Ba)gydF4y2Ba发现一个矩阵gydF4y2Ba紧急求救信号gydF4y2Ba带增益的二阶截面形式gydF4y2BaggydF4y2Ba这相当于用传递函数系数向量表示的数字滤波器gydF4y2BabgydF4y2Ba和gydF4y2Ba一个gydF4y2Ba.gydF4y2Ba

[gydF4y2Ba紧急求救信号gydF4y2Ba,gydF4y2BaggydF4y2Ba) = tf2sos (gydF4y2BabgydF4y2Ba,gydF4y2Ba一个gydF4y2Ba,gydF4y2Ba订单gydF4y2Ba)gydF4y2Ba中行的顺序gydF4y2Ba紧急求救信号gydF4y2Ba.gydF4y2Ba

[gydF4y2Ba紧急求救信号gydF4y2Ba,gydF4y2BaggydF4y2Ba) = tf2sos (gydF4y2BabgydF4y2Ba,gydF4y2Ba一个gydF4y2Ba,gydF4y2Ba订单gydF4y2Ba,gydF4y2Ba规模gydF4y2Ba)gydF4y2Ba指定所有二阶分段的增益和分子系数的比例。gydF4y2Ba

紧急求救信号gydF4y2Ba= tf2sos (gydF4y2Ba___gydF4y2Ba)gydF4y2Ba将整个系统增益嵌入到第一部分。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

利用该函数设计了一个巴特沃斯四阶低通滤波器gydF4y2Ba黄油gydF4y2Ba.指定截止频率为奈奎斯特频率的一半。将过滤器实现为二阶分段。通过比较它们的分子和分母来验证这两种表示是相同的。gydF4y2Ba

(nm, dn) =黄油(0.5 4);(ss, gn) = tf2sos (nm, dn);数= [conv (ss(1:3),学生(1:3))* gn;纳米)gydF4y2Ba
数=gydF4y2Ba2×5gydF4y2Ba0.0940 0.3759 0.5639 0.3759 0.0940 0.0940 0.3759 0.5639 0.3759 0.0940gydF4y2Ba
分母项= [conv (ss(6节),学生(2、6));dn)gydF4y2Ba
分母项=gydF4y2Ba2×5gydF4y2Ba1.0000 0.0000 0.4860 0.0000 0.0177 1.0000 0.0000 0.4860 0.0000 0.0177gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

传递函数系数,指定为向量。将传递函数表示为gydF4y2BabgydF4y2Ba和gydF4y2Ba一个gydF4y2Ba作为gydF4y2Ba

HgydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba =gydF4y2Ba BgydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba 一个gydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba =gydF4y2Ba bgydF4y2Ba 1gydF4y2Ba +gydF4y2Ba bgydF4y2Ba 2gydF4y2Ba zgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba ⋯gydF4y2Ba +gydF4y2Ba bgydF4y2Ba ngydF4y2Ba +gydF4y2Ba 1gydF4y2Ba zgydF4y2Ba −gydF4y2Ba ngydF4y2Ba 一个gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba 2gydF4y2Ba zgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba ⋯gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba 米gydF4y2Ba +gydF4y2Ba 1gydF4y2Ba zgydF4y2Ba −gydF4y2Ba 米gydF4y2Ba .gydF4y2Ba

例子:gydF4y2BaB = [1 3 3 1]/6gydF4y2Ba和gydF4y2BaA = [3 0 1 0]/3gydF4y2Ba指定一个三阶巴特沃斯滤波器,归一化频率为3db, 0.5π rad/sample。gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

行顺序,指定为下列之一:gydF4y2Ba

  • “了”gydF4y2Ba-按顺序排列,第一行gydF4y2Ba紧急求救信号gydF4y2Ba包含离单位圆最远的极点。gydF4y2Ba

  • “下来”gydF4y2Ba-按顺序排列,第一行gydF4y2Ba紧急求救信号gydF4y2Ba包含最接近单位圆的极点。gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba

增益系数和分子系数的比例,指定为以下其中之一:gydF4y2Ba

  • “没有”gydF4y2Ba-不应用缩放。gydF4y2Ba

  • “正”gydF4y2Ba-应用无穷范数缩放。gydF4y2Ba

  • “两个”gydF4y2Ba-应用2-范数缩放。gydF4y2Ba

用无穷范数缩放gydF4y2Ba“了”gydF4y2Ba- ordered最小化了实现中溢出的概率。用2范数缩放gydF4y2Ba“下来”gydF4y2Ba-排序最小化峰值舍入噪声。gydF4y2Ba

请注意gydF4y2Ba

无限范数和二范数缩放仅适用于直接形式II实现。gydF4y2Ba

数据类型:gydF4y2Ba字符gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

二阶截面表示,以矩阵形式返回。gydF4y2Ba紧急求救信号gydF4y2Ba是一个gydF4y2BalgydF4y2Ba6矩阵gydF4y2Ba

紧急求救信号gydF4y2Ba =gydF4y2Ba [gydF4y2Ba bgydF4y2Ba 01gydF4y2Ba bgydF4y2Ba 11gydF4y2Ba bgydF4y2Ba 21gydF4y2Ba 1gydF4y2Ba 一个gydF4y2Ba 11gydF4y2Ba 一个gydF4y2Ba 21gydF4y2Ba bgydF4y2Ba 02gydF4y2Ba bgydF4y2Ba 12gydF4y2Ba bgydF4y2Ba 22gydF4y2Ba 1gydF4y2Ba 一个gydF4y2Ba 12gydF4y2Ba 一个gydF4y2Ba 22gydF4y2Ba ⋮gydF4y2Ba ⋮gydF4y2Ba ⋮gydF4y2Ba ⋮gydF4y2Ba ⋮gydF4y2Ba ⋮gydF4y2Ba bgydF4y2Ba 0gydF4y2Ba lgydF4y2Ba bgydF4y2Ba 1gydF4y2Ba lgydF4y2Ba bgydF4y2Ba 2gydF4y2Ba lgydF4y2Ba 1gydF4y2Ba 一个gydF4y2Ba 1gydF4y2Ba lgydF4y2Ba 一个gydF4y2Ba 2gydF4y2Ba lgydF4y2Ba ]gydF4y2Ba

哪一行包含分子系数和分母系数gydF4y2BabgydF4y2Ba本土知识gydF4y2Ba和gydF4y2Ba一个gydF4y2Ba本土知识gydF4y2Ba的二阶部分gydF4y2BaHgydF4y2Ba(gydF4y2BazgydF4y2Ba):gydF4y2Ba

HgydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ggydF4y2Ba ∏gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba lgydF4y2Ba HgydF4y2Ba kgydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ggydF4y2Ba ∏gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba lgydF4y2Ba bgydF4y2Ba 0gydF4y2Ba kgydF4y2Ba +gydF4y2Ba bgydF4y2Ba 1gydF4y2Ba kgydF4y2Ba zgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba bgydF4y2Ba 2gydF4y2Ba kgydF4y2Ba zgydF4y2Ba −gydF4y2Ba 2gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba 1gydF4y2Ba kgydF4y2Ba zgydF4y2Ba −gydF4y2Ba 1gydF4y2Ba +gydF4y2Ba 一个gydF4y2Ba 2gydF4y2Ba kgydF4y2Ba zgydF4y2Ba −gydF4y2Ba 2gydF4y2Ba .gydF4y2Ba

整体系统增益,返回为实标量。gydF4y2Ba

如果你叫gydF4y2Batf2sosgydF4y2Ba使用一个输出参数,函数将整个系统增益嵌入到第一部分,gydF4y2BaHgydF4y2Ba1gydF4y2Ba(gydF4y2BazgydF4y2Ba),所以,gydF4y2Ba

HgydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba =gydF4y2Ba ∏gydF4y2Ba kgydF4y2Ba =gydF4y2Ba 1gydF4y2Ba lgydF4y2Ba HgydF4y2Ba kgydF4y2Ba (gydF4y2Ba zgydF4y2Ba )gydF4y2Ba .gydF4y2Ba

请注意gydF4y2Ba

当缩放直接形式II结构时,不建议将增益嵌入第一部分,并可能导致不规则缩放。为了避免嵌入增益,使用gydF4y2Batf2sosgydF4y2Ba有两个输出。gydF4y2Ba

算法gydF4y2Ba

tf2sosgydF4y2Ba使用四步算法确定输入传递函数系统的二阶截面表示:gydF4y2Ba

  1. 它找到了系统的极点和零点gydF4y2BabgydF4y2Ba和gydF4y2Ba一个gydF4y2Ba.gydF4y2Ba

  2. 它使用函数gydF4y2Bazp2sosgydF4y2Ba,它首先将零和极点分组成复共轭对,使用gydF4y2BacplxpairgydF4y2Ba函数。gydF4y2Bazp2sosgydF4y2Ba然后根据以下规则将极点与零点对匹配形成二阶分段:gydF4y2Ba

    1. 将最接近单位圆的极点与最接近这些极点的零点匹配起来。gydF4y2Ba

    2. 把离单位圆最近的极点和离这些极点最近的零匹配起来。gydF4y2Ba

    3. 继续,直到所有的极点和零都匹配。gydF4y2Ba

    tf2sosgydF4y2Ba将实极按绝对值最接近的实极分组。同样的规则也适用于真实的零。gydF4y2Ba

  3. 它根据极对与单位圆的接近程度来排列截面。gydF4y2Batf2sosgydF4y2Ba通常情况下,在级联中,最接近单位圆的区段排在最后。你可以告诉gydF4y2Batf2sosgydF4y2Ba通过指定以相反的顺序排列各节gydF4y2Ba订单gydF4y2Ba作为gydF4y2Ba“下来”gydF4y2Ba.gydF4y2Ba

  4. tf2sosgydF4y2Ba按中规定的标准缩放各部分gydF4y2Ba规模gydF4y2Ba.对于任意的gydF4y2BaHgydF4y2Ba(gydF4y2BaωgydF4y2Ba)gydF4y2Ba,比例定义为gydF4y2Ba

    为gydF4y2Ba HgydF4y2Ba 为gydF4y2Ba pgydF4y2Ba =gydF4y2Ba [gydF4y2Ba 1gydF4y2Ba 2gydF4y2Ba πgydF4y2Ba ∫gydF4y2Ba 0gydF4y2Ba 2gydF4y2Ba πgydF4y2Ba |gydF4y2Ba HgydF4y2Ba (gydF4y2Ba ωgydF4y2Ba )gydF4y2Ba |gydF4y2Ba pgydF4y2Ba dgydF4y2Ba ωgydF4y2Ba ]gydF4y2Ba 1gydF4y2Ba /gydF4y2Ba pgydF4y2Ba

    在哪里gydF4y2BapgydF4y2Ba可以是∞或2。有关缩放的详细信息,请参阅参考资料。该算法遵循这种缩放,试图在定点滤波器实现中最小化溢出或峰值舍入噪声。gydF4y2Ba

参考文献gydF4y2Ba

L. B.杰克逊gydF4y2Ba数字滤波器和信号处理gydF4y2Ba.3版。波士顿:克鲁尔学术出版社,1996。gydF4y2Ba

Mitra, s.k。gydF4y2Ba数字信号处理:基于计算机的方法gydF4y2Ba.纽约:麦格劳-希尔出版社,1998年。gydF4y2Ba

[3] Vaidyanathan, P. P.《鲁棒数字滤波器结构》。gydF4y2Ba数字信号处理手册gydF4y2Ba(s.k. Mitra和j.f. Kaiser编)。纽约:约翰·威利父子公司,1993年。gydF4y2Ba

扩展功能gydF4y2Ba

另请参阅gydF4y2Ba

|gydF4y2Ba|gydF4y2Ba|gydF4y2Ba|gydF4y2Ba|gydF4y2Ba|gydF4y2Ba

之前介绍过的R2006agydF4y2Ba