主要内容

西森

二阶(双层)IIR数字滤波

描述

例子

y= sosfilt(SOS.X应用二阶段数字滤波器SOS.输入信号X

y= sosfilt(SOS.X暗淡沿着维度操作暗淡

例子

全部收缩

加载chirp.mat.。该文件包含一个信号,y,它具有上面的大部分电力FS / 4.或者奈奎斯特频率的一半。采样率为8192Hz。

加载唧唧喳喳t =(0:长度(y)-1)/ fs;

设计第七次Butterworth高通滤波器,可衰减以下信号的组件FS / 4.。使用归一化截止频率为0.48πRad /样品。就二阶部分表示滤波器系数。

[zhi,phi,khi] =黄油(7,0.48,'高的');Soshi = ZP2SOS(ZP2SOS(ZP2SOO);Freqz(Soshi)

图包含2个轴。轴1包含类型线的对象。轴2包含类型线的对象。

过滤信号。显示原始和高通滤波的信号。使用相同y- 两个绘图的轴垢。

Outhi = Sosfilt(Sososfilt(Soshi,Y);图形子图(2,1,1)图(T,Y)标题('原始信号')ys = ylim;子图(2,1,2)绘图(T,OUTHI)标题('高通滤波信号')Xlabel('时间'')ylim(ys)

图包含2个轴。具有标题原始信号的轴1包含类型线的对象。具有标题高通滤波信号的轴2包含类型线的物体。

设计具有相同规格的低通滤波器。过滤信号并将结果与​​原稿进行比较。使用相同y- 两个绘图的轴垢。结果主要是噪音。

[zlo,plo,klo] =黄油(7,0.48);Soslo = ZP2SOS(ZLO,PLO,KLO);Outlo = Sosfilt(Soslo,Y);子图(2,1,1)绘图(T,Y)标题('原始信号')ys = ylim;子图(2,1,2)绘图(T,Outlo)标题('低通滤波信号')Xlabel('时间'')ylim(ys)

图包含2个轴。具有标题原始信号的轴1包含类型线的对象。带标题低通滤波信号的轴2包含类型线的物体。

输入参数

全部收缩

二阶截面数字滤波器,指定为一个L.-6矩阵,在哪里L.是二阶部分的数量。矩阵

SOS. = [ B. 01. B. 11. B. 21. 1 一个 11. 一个 21. B. 02. B. 12. B. 22. 1 一个 12. 一个 22. B. 0. L. B. 1 L. B. 2 L. 1 一个 1 L. 一个 2 L. ]

代表二阶部分数字滤波器

H Z. = π K. = 1 L. H K. Z. = π K. = 1 L. B. 0. K. + B. 1 K. Z. - 1 + B. 2 K. Z. - 2 1 + 一个 1 K. Z. - 1 + 一个 2 K. Z. - 2

例子:[B,A] =黄油(3,1 / 32);SOS = TF2SOS(B,A)指定具有归一化3 dB频率的三阶Butterworth滤波器π/ 32 rad /样品。

数据类型:单身的|双倍的

输入信号,指定为向量,矩阵或N-D阵列。

例子:X = [2 1]。* SIN(2 * PI *(0:127)'./ [16 64])指定双通道正弦曲线。

数据类型:单身的|双倍的
复数支持:万博1manbetx是的

维度运行,指定为正整数标量。默认,西森沿第一阵列维度操作X大小大于1。

数据类型:单身的|双倍的

输出参数

全部收缩

过滤信号,返回矢量,矩阵或N-D阵列。y具有相同的尺寸X

参考

[1]银行,巴拉斯。“将无限脉冲响应滤波器转换为平行形式”。IEEE信号处理杂志。卷。35,第3号,2018年5月,第124-130页。

[2] Orfanidis,Sophocles J.信号处理简介。Englewood Cliffs,NJ:Prentice大厅,1996年。

扩展能力

C / C ++代码生成
使用MATLAB®Coder™生成C和C ++代码。

也可以看看

||

在R2006A之前介绍