主要内容

sosfilt

二阶(双二次)IIR数字滤波

在页面中全部折叠

语法

x y = sosfilt (sos)
x, y = sosfilt (sos暗)

描述

例子

y= sosfilt (紧急求救信号，x）应用二阶截面数字滤波器紧急求救信号到输入信号x．

  • 如果x是一个矩阵，则该函数沿第一个维度操作，并为每一列返回经过筛选的数据。

  • 如果x是多维数组，则该函数沿大小大于1的第一个数组维度操作。

y= sosfilt (紧急求救信号，x，昏暗的）沿着维度进行运算昏暗的．

例子

全部折叠

打开生活的脚本

负载chirp.mat．该文件包含一个信号，y，它的大部分能量都在上面Fs / 4或者Nyquist频率的一半。采样速率为8192hz。

负载尖声地说t =(0:长度(y) 1) / Fs;

设计一个七阶巴特沃斯高通滤波器来衰减下面信号的分量Fs / 4．使用归一化截止频率0.48πrad /样品。用二阶截面表示滤波器系数。

(智,φ,川崎)=黄油(7、0.48“高”）;soshi = zp2sos(智,φ,川崎);freqz (soshi)

图中包含2个轴对象。标题为Phase的Axes对象1包含一个类型为line的对象。标题为Magnitude的Axes对象2包含一个类型为line的对象。

过滤信号。显示原始信号和高通滤波信号。使用相同的y-轴尺度为两个地块。

outhi = sosfilt (soshi y);Figure subplot(2,1,1) plot(t,y)原始信号的) ys = ylim;次要情节(2,1,2)情节(t, outhi)标题(“Highpass-Filtered信号”)包含(“时间(s)”) ylim (y)

图中包含2个轴对象。标题为Original Signal的Axes对象1包含一个类型为line的对象。标题为“高通滤波信号”的Axes对象2包含一个类型为line的对象。

设计一个相同规格的低通滤波器。对信号进行滤波，并将结果与原始信号进行比较。使用相同的y-轴尺度为两个地块。结果大多是噪音。

[zlo,巴解组织,klo] =黄油(7,0.48);soslo = zp2sos(巴解组织,zlo klo);outlo = sosfilt (soslo y);次要情节(2,1,1)情节(t, y)标题(原始信号的) ys = ylim;次要情节(2,1,2)情节(t, outlo)标题(“Lowpass-Filtered信号”)包含(“时间(s)”) ylim (y)

图中包含2个轴对象。标题为Original Signal的Axes对象1包含一个类型为line的对象。标题为low - pass- filtered Signal的Axes对象2包含一个类型为line的对象。

输入参数

全部折叠

二阶截面数字滤波器，指定为l6矩阵l为二阶分段数。矩阵

紧急求救信号 ＝ ［ 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] =黄油(1/32);sos = tf2sos (b)指定一个三阶巴特沃斯滤波器，其归一化3db频率为π/ 32 rad /样品。

数据类型:单|双

输入信号，指定为向量、矩阵或N- d数组。

例子:X = [2 1].*sin(2*pi*(0:127)'。/ [64])指定双通道正弦信号。

数据类型:单|双
复数的支持:是的

要操作的维度，指定为正整数标量。默认情况下，该函数沿着的第一个数组维度操作xsize大于1。

数据类型:单|双

输出参数

全部折叠

经过过滤的信号，返回为向量、矩阵或N- d数组。y尺寸和x．

参考文献

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

奥法尼迪斯，索福克勒斯。信号处理概论．恩格尔伍德悬崖，新泽西州:Prentice-Hall, 1996。

扩展功能

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

版本历史

之前介绍过的R2006a

另请参阅

||

Baidu
map