AnyMath 文档

切比2

切比雪夫滤波器类型II的计算。

库::`工程师`

语法

函数调用

  • [参数:ba]=cheby2(<参数:n>>,<参数:Rs>>,<参数:Ws>>) -设计切比雪夫II型数字低通滤波器 [参数:n]`延迟带的归一化边界频率的第一阶 `[参数:Ws] 和延迟带中的衰减 [参数:Rs] dB从峰值带宽值。 功能 cheby2 返回滤波器传递函数的分子和分母的系数。

  • [参数:ba]=cheby2(<参数:n>>,<参数:Rs>>,<参数:Ws>>,<参数:ftype>>) -设计切比雪夫II型数字滤波器:低通,高通,带通或陷波滤波器,具体取决于参数的值 [参数:ftype] 和元素的数量 [参数:Ws]. 带通和陷波滤波器的设计顺序如下 2<参数:n>>.

  • [参数:zpk]=cheby2(_;<参数:nout>>=3) -设计一个切比雪夫II型数字滤波器,并返回其零,极点和增益。 此语法可以包括前面选项中的任何输入参数。

  • A, B, C, D=cheby2(_;nout=4) -设计切比雪夫II型数字滤波器并返回定义其在状态空间中的表示的矩阵。

  • =cheby2(,"s";[参数:nout]) -使用前面语法中的任何输入或输出参数设计模拟切比雪夫II型滤波器。 论点 [参数:nout] 显示函数将返回多少个输出参数。 如果论点是 [参数:nout] 省略,该函数将仅对两个输出参数进行计算。

争论

输入参数

# n — 过滤顺序

+ 通行证:[标量]

Details

指定为整数标量的筛选器顺序小于或等于 500. 用于带通和陷波滤波器 n 表示滤波器阶数的一半。

数据类型

漂浮64</无翻译>

# Rs — 延迟带中的衰减,单位为dB

+ 传:[正标量]

Details

延迟带中的衰减相对于通带中的峰值,以dB为单位设置为正标量。

如果值为 以线性单位表示,您可以使用公式将其转换为dB Rs的 .

数据类型

漂浮64</无翻译>

# Ws — 极限频率 的延迟带+ 通行证:[标量] | 传递:[双元向量]

Details

延迟带的边界频率,指定为标量或双元矢量。 延迟带的边界频率是滤波器的幅度-频率响应(频率响应)等于的频率 [参数:Rs] 在dB中。 大延迟带衰减值 [参数:Rs] 它们导致带宽增加。

如果 Ws -一个标量,然后 *cheby2 设计具有截止频率的低通或高通滤波器 Ws.

+ 如果 Ws -二元素向量 [w1w2],在哪里 w1<w2 然后 cheby2 设计具有较低截止频率的带通或陷波滤波器 w1 和上边界频率 w2.

*对于数字滤波器,延迟带的频率限制应在范围内 0 以前 1,在哪里 1 对应于奈奎斯特频率-采样频率的一半或 rad/倒计时。

+ 对于模拟滤波器,延迟带的边界频率必须以rad/s表示,并且可以取任意正值。

数据类型

漂浮64</无翻译>

# ftype — 过滤器类型

+ 通过:["低"] | 通行证:["带通"] | 通行证:["高"] | 通行证:["停止"]

Details

指定为的筛选器类型:

  • "低" -具有延迟带的边界频率的低通滤波器 [参数:Ws]. 此值默认用于标量 [参数:Ws];

  • "高" -具有延迟带的边界频率的高通滤波器 [参数:Ws];

  • "带通" -带通滤波器 2<参数:n>> 如果 [参数:Ws] -一个双元向量。 默认情况下,此值在以下情况下使用 [参数:Ws] 设置为双元素向量;

  • "停止" -陷波(阻塞)滤波器 2<参数:n>> 如果 [参数:Ws] -一个双元向量。

数据类型

字符串</无翻译>

名称-值输入参数

# nout — 输出参数数

+ 2 (默认情况下)| 3 | 4

Details

输出参数的数量,设置为 2, 34. 如果论点是 nout 未指定,函数默认返回两个输出参数。

输出参数

# *b,a*是 传递函数的系数

+ 传递:[字符串向量]

Details

作为行向量返回的滤波器传递函数的系数。 具有指定的过滤顺序 [参数:n] 函数返回 bar 通过计算在哪里 R=[参数:n]+1 对于低和高通滤波器和 r=2*[参数:n]+1 用于带通和陷波滤波器。

传递函数表示为 :

*用于数字滤波器

+

*模拟滤波器

+

数据类型

漂浮64</无翻译>

# z,p,k — 零、极点和增益

+ 传递:[列向量和标量]

Details

零,极点和滤波器增益,返回为两个向量-列和标量。 具有指定的过滤器顺序 [参数:n] 函数返回 zpr 通过计算在哪里 R=<参数:n>> 对于低通和高通滤波器和 r=2*<参数:n>> 用于带通和陷波滤波器。

传递函数表示为 , :

*用于数字滤波器

+

*模拟滤波器

+

数据类型

漂浮64</无翻译>

# A,B,C,D — 状态空间中滤波器的表示

+ 传递:[矩阵]

Details

滤波器在状态空间中的表示形式,作为矩阵返回。 如果 R=<参数:n>> 对于低和高通滤波器和 R=2<参数:n>> 对于带通和陷波滤波器,则 A 这是矩阵 rr, B 矩阵 r1, C 矩阵 1r,而 D11.

状态空间矩阵与状态向量相关 、入口 和出口 借助方程组:

*用于数字滤波器

+

*模拟滤波器

+

数据类型

漂浮64</无翻译>

例子:

切比雪夫II型低通滤波器

Details

让我们设计一个切比雪夫II型低通滤波器 6-屏障带衰减的阶数 50 dB和带宽截止频率 300 Hz,其对应于 0.6π 频率采样数据的rad/计数 1000 赫兹。 让我们绘制幅频和相位特性。 我们用它来过滤一个频率随机信号 1000 数数。

import EngeeDSP.Functions: cheby2, freqz, randn, filter

fc = 300
fs = 1000

b, a = cheby2(6,50,fc/(fs/2))

freqz(b,a,out=:plot)

cheby2 1

dataIn = randn(1000,1)
dataOut = filter(b,a,dataIn)

切比雪夫II型屏障过滤器

Details

让我们设计一个Ii型切比雪夫屏障滤波器 6-具有归一化边界频率的阶数 0.2π0.6π 吊杆车道的rad/倒计时和衰减 50 分贝。 让我们绘制幅频和相位特性。 我们用它来过滤随机信号。

import EngeeDSP.Functions: cheby2, freqz, randn, filter

b,a = cheby2(6,50,[0.2 0.6],"stop")

freqz(b,a,out=:plot)

cheby2 2

dataIn = randn(1000,1)
dataOut = filter(b,a,dataIn)

切比雪夫II型高通滤波器

Details

让我们设计一个切比雪夫II型高通滤波器 9-带带宽截止频率的阶数 300 Hz,用于以频率采样的数据 1000 Hz,对应于 0.6π rad/倒计时。 设置屏障条中的衰减 20 分贝。 将零、极点和增益转换为二阶部分。 让我们绘制幅频和相位特性。

import EngeeDSP.Functions: cheby2, zp2sos

z,p,k = cheby2(9,20,300/500,"high"; nout=3)
sos = zp2sos(z,p,k)

freqz(sos,out=:plot)

cheby2 3

切比雪夫II型带通滤波器

Details

让我们设计一个切比雪夫II型带通滤波器 10-带宽较低的顺序 500 Hz和上带宽 560 赫兹。 设置屏障带中的衰减 40 dB和采样率 1500 赫兹。 我们在状态空间中使用表示。 我们将状态空间中的表示转换为二阶段的形式。 让我们绘制幅频和相位特性。

import EngeeDSP.Functions: cheby2, ss2sos
fs = 1500

A,B,C,D = cheby2(10,40,[500 560]/(fs/2); nout=4)
sos = ss2sos(A,B,C,D)[1]

freqz(sos,out=:plot)

cheby2 4

算法

切比雪夫II型滤波器在通带中是单调的,在延迟带中具有均匀的脉动。 II型滤波器的性能下降速度不如i型滤波器快,但在带宽上没有波纹。

切比雪夫过滤器II型 cheby2 使用五步算法:

*查找模拟低通原型的极点、零点和增益。

*将极点、零点和增益转换为状态空间。

*如有必要,使用状态空间变换将低通滤波器转换为具有所需频率约束的高通滤波器、带通滤波器或陷波滤波器。

*对于数字滤波器的设计,通过双线性频率预失真转换将模拟滤波器转换为数字滤波器。 精细频率调谐允许模拟和数字滤波器具有相同的频率响应幅度 Wsw1w2.

*必要时将状态空间滤波器转换回传递函数或零极点增益形式。

文学作品

  1. 里昂,理查德G._Understanding数字信号处理._上鞍河,NJ:普伦蒂斯霍尔,2004.