索斯菲尔特

二阶（双二次）IIR数字滤波

折叠所有页面

语法

y=sosfilt（sos，x）

y = sosfilt（sos，x，dim）

描述

实例

Y=sosfilt(SOS.,x)应用二阶分段数字滤波器SOS.输入信号x.

Y=sosfilt(SOS.,x,暗淡的)沿维度进行操作暗淡的.

例子

全部崩溃

二阶分段滤波

打开直播脚本

负载chirp.mat.该文件包含一个信号，Y，它的大部分力量都在上面FS / 4.，或奈奎斯特频率的一半。采样率为8192赫兹。

负载啁啾声t =（0：长度（y）-1）/ fs;

设计一个七阶Butterworth高通滤波器来衰减以下信号的分量FS / 4.. 使用标准化截止频率0.48πrad/sample。用二阶截面表示滤波器系数。

[zhi，phi，khi]=黄油（7,0.48，“高”)；soshi=zp2sos（zhi，phi，khi）；freqz（soshi）

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

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

outhi=sosfilt（soshi，y）；图子地块（2,1,1）地块（t，y）标题('原始信号'）ys = ylim;子图（2,1,2）绘图（T，OUTHI）标题（'高通滤波信号'）xlabel（'时间'')叶利姆（ys）

图中包含2个轴对象。标题为原始信号的轴对象1包含line类型的对象。标题为高通滤波信号的轴对象2包含line类型的对象。

设计一个具有相同规格的低通滤波器。对信号进行滤波，并将结果与原始结果进行比较。使用相同的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（'时间'')叶利姆（ys）

图中包含2个轴对象。标题为原始信号的轴对象1包含line类型的对象。标题为低通滤波信号的轴对象2包含line类型的对象。

输入参数

全部崩溃

`SOS.`—二阶截面数字滤波器
L-by-6矩阵

二阶截面数字滤波器，指定为一个L-by-6矩阵，其中L是二阶节数。矩阵

$SOS. = [\begin{matrix} B_{01} & B_{11} & B_{21} & 1. & {A.}_{11} & {A.}_{21} \\ B_{02} & B_{12} & B_{22} & 1. & {A.}_{12} & {A.}_{22} \\ ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋮ \\ B_{0 L} & B_{1. L} & B_{2. L} & 1. & {A.}_{1. L} & {A.}_{2. L} \end{matrix}]$

表示二阶分段数字滤波器

$H (Z) = \prod_{K = 1.}^{L} H_{K} (Z) = \prod_{K = 1.}^{L} \frac{B_{0 K} + B_{1. K} Z^{- 1.} + B_{2. K} Z^{- 2.}}{1. + {A.}_{1. K} Z^{- 1.} + {A.}_{2. K} Z^{- 2.}} .$

例子：[b，a]=黄油（3,1/32）；sos=tf2sos（b，a）指定标准化频率为3 dB的三阶巴特沃斯滤波器π/32拉德/样品。

数据类型：仅有一个的|双重的

`x`—输入信号
矢量|矩阵|N-D阵列

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

例子：x=[21].*sin（2*pi*（0:127）./[16 64]）指定双通道正弦曲线。

数据类型：仅有一个的|双重的
复数支持：金宝app对

`暗淡的`—操作的维度
正整数标量

维度运行，指定为正整数标量。默认情况下，索斯菲尔特沿第一阵列维度操作x其大小大于1。

数据类型：仅有一个的|双重的

输出参数

全部崩溃

`Y`- 过滤信号
矢量|矩阵|N-D阵列

滤波信号，以向量、矩阵或N-D阵列。Y大小与x.

参考

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

[2] 奥法尼迪斯，索福克勒斯J。信号处理简介恩格尔伍德悬崖，新泽西州：普伦蒂斯大厅，1996年。

扩展能力

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

GPU阵列
使用并行计算工具箱™在图形处理单元（GPU）上运行，加速代码。

使用说明和限制：

输入滤波器SOS.必须稳定。用缺席检查过滤器的稳定性。
输入滤波器的所有二阶子部分必须是IIR。
这个gpuArraySOSFILT版本使用并行算法[1]与matlab不同^®版本这些算法对于具有不同参数的复数输入给出了不同的结果楠或者INF.价值观：
- 在Matlab版本中，楠s和INF.s只在真实的部分传播。
- 在里面gpuArray版本楠s和INF.s在实部和虚部都传播。