主要内容

cheby2

切比雪夫II型滤波器设计

描述

例子

[<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-ba" class="intrnllnk">b,一个= cheby2(<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-n" class="intrnllnk">n,<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-Rs" class="intrnllnk">Rs,<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj-Ws" class="intrnllnk">Ws返回an的传递函数系数n边频归一化止带低通数字切比雪夫二型滤波器Ws而且Rs阻带衰减的分贝从峰值通带值下降。

例子

[<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-ba" class="intrnllnk">b,一个= cheby2(<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-n" class="intrnllnk">n,<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-Rs" class="intrnllnk">Rs,<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj-Ws" class="intrnllnk">Ws,<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj-ftype" class="intrnllnk">ftype根据的值,设计低通、高通、带通或带阻切比雪夫II型滤波器ftype元素的个数Ws.所得到的带通和带阻设计是2阶的n

注意:看到<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/cheby2.html" class="intrnllnk">限制有关影响传递函数形成的数值问题的信息。

例子

[<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-zpk" class="intrnllnk">z,磷、钾= cheby2(<年代p一个nclass="argument_placeholder">___设计了一个低通、高通、带通或带阻数字切比雪夫II型滤波器,并返回其零点、极点和增益。此语法可以包括以前语法中的任何输入参数。

例子

[<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-ABCD" class="intrnllnk">A, B, C, D= cheby2(<年代p一个nclass="argument_placeholder">___设计一个低通、高通、带通或带阻数字切比雪夫II型滤波器,并返回指定其状态空间表示的矩阵。

例子

[<年代p一个nclass="argument_placeholder">___= cheby2(<年代p一个nclass="argument_placeholder">___, ' s ')设计了具有阻带边角频率的低通、高通、带通或带阻模拟切比雪夫II型滤波器<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj-Ws" class="intrnllnk">Ws而且<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-Rs" class="intrnllnk">Rs阻带衰减的分贝。

例子

全部折叠

设计一个六阶低通切比雪夫II型滤波器,阻带衰减为50 dB,阻带边频为300 Hz,对于1000hz采样的数据,对应<年代p一个nclass="inlineequation"> 0 6 π rad /样品。画出它的幅值和相位响应。用它来过滤1000个样本的随机信号。

Fc = 300;Fs = 1000;[b,a] = cheby2(6,50,fc/(fs/2));Freqz (b,a,[],fs) subplot(2,1,1) ylim([-100 20])

图中包含2个轴对象。axis对象1,标题为Phase, xlabel Frequency (Hz), ylabel Phase (degrees),包含一个类型为line的对象。标题为Magnitude, xlabel Frequency (Hz), ylabel Magnitude (dB)的坐标轴对象2包含一个类型为line的对象。

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

设计一个六阶切比雪夫II型带阻滤波器,其边频归一化为<年代p一个nclass="inlineequation"> 0 2 π 而且<年代p一个nclass="inlineequation"> 0 6 π rad/sample和50 dB的阻带衰减。画出它的幅值和相位响应。使用它来过滤随机数据。

[b,a] = cheby2(3,50,[0.2 0.6],<年代p一个n年代tyle="color:#A020F0">“停止”);freqz (b)

图中包含2个轴对象。坐标轴对象1,标题为Phase, xlabel归一化频率(\times\pi rad/sample), ylabel Phase (degrees)包含一个类型为line的对象。坐标轴对象2,标题为Magnitude, xlabel归一化频率(\times\pi rad/sample), ylabel Magnitude (dB)包含一个类型为line的对象。

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

设计一个9阶高通切比雪夫II型滤波器,阻带衰减为20db,阻带边频为300hz,对于1000hz采样的数据,对应<年代p一个nclass="inlineequation"> 0 6 π rad /样品。绘制幅值和相位响应。将零点、极点和增益转换为二阶部分,以供使用fvtool

[z,p,k] = cheby2(9,20,300/500,<年代p一个n年代tyle="color:#A020F0">“高”);SOS = zp2sos(z,p,k);fvtool (sos,<年代p一个n年代tyle="color:#A020F0">“分析”,<年代p一个n年代tyle="color:#A020F0">“频率”)

图1:幅度响应(dB)和相位响应包含一个轴对象。axis对象的标题为幅值响应(dB)和相位响应,xlabel归一化频率(乘以pi空白r d / s a m p le), ylabel幅值(dB)包含一个类型为line的对象。

设计一个20阶切比雪夫II型带通滤波器,低阻带频率为500hz,高阻带频率为560hz。指定阻带衰减为40 dB,采样率为1500 Hz。使用状态空间表示。设计一个相同的过滤器使用designfilt

[A,B,C,D] = cheby2(10,40,[500 560]/750);D = designfilt(<年代p一个n年代tyle="color:#A020F0">“bandpassiir”,<年代p一个n年代tyle="color:#A020F0">“FilterOrder”, 20岁,<年代p一个n年代tyle="color:#0000FF">...“StopbandFrequency1”, 500,<年代p一个n年代tyle="color:#A020F0">“StopbandFrequency2”, 560,<年代p一个n年代tyle="color:#0000FF">...“StopbandAttenuation”现年40岁的<年代p一个n年代tyle="color:#A020F0">“SampleRate”, 1500);

将状态空间表示转换为二阶部分。可视化的频率响应使用fvtool

sos = ss2sos(A,B,C,D);FVT = fvtool(sos,d,<年代p一个n年代tyle="color:#A020F0">“Fs”, 1500);传奇(fvt<年代p一个n年代tyle="color:#A020F0">“cheby2”,<年代p一个n年代tyle="color:#A020F0">“designfilt”)

图1:幅度响应(dB)包含一个轴对象。axis对象的标题为Magnitude Response (dB), xlabel Frequency (Hz), ylabel Magnitude (dB),其中包含2个类型为line的对象。这些对象表示cheby2、designfilt。

设计一个截止频率为2 GHz的五阶模拟巴特沃斯低通滤波器。乘以<年代p一个nclass="inlineequation"> 2 π 将频率转换为每秒弧度。计算滤波器在4096点处的频率响应。

N = 5;Fc = 2e9;[zb,pb,kb] =黄油(n,2*pi*fc,<年代p一个n年代tyle="color:#A020F0">“s”);[bb,ab] = zp2tf(zb,pb,kb);[hb,wb] = freqs(bb,ab,4096);

设计一个具有相同边频和3db通带纹波的五阶切比雪夫I型滤波器。计算它的频率响应。

[z1,p1,k1] = cheby1(n,3,2*pi*fc,<年代p一个n年代tyle="color:#A020F0">“s”);[b1,a1] = zp2tf(z1,p1,k1);[h1,w1] = freqs(b1,a1,4096);

设计具有相同边频和30 dB阻带衰减的五阶切比雪夫II型滤波器。计算它的频率响应。

[z2,p2,k2] = cheby2(n,30,2*pi*fc,<年代p一个n年代tyle="color:#A020F0">“s”);[b2,a2] = zp2tf(z2,p2,k2);[h2,w2] = freqs(b2,a2,4096);

设计一个边频相同,通带纹波为3 dB,阻带衰减为30 dB的五阶椭圆滤波器。计算它的频率响应。

[ze,pe,ke] = ellip(n,3,30,2*pi*fc,<年代p一个n年代tyle="color:#A020F0">“s”);[be,ae] = zp2tf(ze,pe,ke);[he,we] = freqs(be,ae,4096);

设计一个具有相同边缘频率的五阶贝塞尔滤波器。计算它的频率响应。

[zf,pf,kf] = besself(n,2*pi*fc);[bf,af] = zp2tf(zf,pf,kf);[hf,wf] = freqs(bf,af,4096);

以分贝为单位绘制衰减图。以千兆赫表示频率。比较过滤器。

Plot ([wb w1 w2 we wf]/(2e9*pi),<年代p一个n年代tyle="color:#0000FF">...Mag2db (abs([hb h1 h2 he hf])))轴([0 5 -45 5])<年代p一个n年代tyle="color:#A020F0">“频率(GHz)”) ylabel (<年代p一个n年代tyle="color:#A020F0">“衰减(dB)”)传说([<年代p一个n年代tyle="color:#A020F0">“黄油”“cheby1”“cheby2”“ellip”“besself”])

图中包含一个轴对象。具有xlabel Frequency (GHz), ylabel Attenuation (dB)的坐标轴对象包含5个类型为line的对象。这些对象表示butter、cheby1、cheby2、ellip、besself。

巴特沃斯和切比雪夫II型滤波器有平坦的通带和宽的过渡带。切比雪夫I型和椭圆滤波器滚出更快,但有通带纹波。切比雪夫II型设计函数的频率输入设置止带的开始,而不是通带的结束。贝塞尔滤波器沿通带具有近似恒定的群延迟。

输入参数

全部折叠

筛选器顺序,指定为整数标量。对于带通和带阻设计,n表示筛选器顺序的一半。

数据类型:

止带衰减从峰值通带值下降,指定为以分贝表示的正标量。

如果您的规格ℓ是线性单位,则可以使用Rs= -20 log<年代ub>10ℓ。

数据类型:

止带边频率,指定为标量或两元矢量。阻带边缘频率是滤波器的幅值响应所在的频率<年代p一个nclass="inlineequation">- - - - - -<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-Rs" class="intrnllnk">Rs分贝。阻带衰减值越大,Rs,导致更宽的跃迁带。

  • 如果Ws是标量吗cheby2设计了带边频的低通或高通滤波器Ws

    如果Ws是二元向量吗(w1 w2),在那里w1<w2,然后cheby2设计了一种带通或带阻滤波器,具有较低的边频w1以及更高的边缘频率w2

  • 对于数字滤波器,阻带边缘频率必须位于0和1之间,其中1对应于奈奎斯特速率-采样速率的一半或πrad /样品。

    对于模拟滤波器,阻带边缘频率必须以弧度/秒表示,并且可以取任何正值。

数据类型:

过滤器类型,指定为以下之一:

  • “低”指定带阻边频率的低通滤波器<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj-Ws" class="intrnllnk">Ws“低”是否默认为标量Ws

  • “高”指定带阻边频率的高通滤波器Ws

  • “带通”指定阶为2的带通滤波器<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-n" class="intrnllnk">n如果Ws是一个二元向量。“带通”是默认值,当Ws有两个元素。

  • “停止”指定阶为2的带阻滤波器n如果Ws是一个二元向量。

输出参数

全部折叠

滤波器的传递函数系数,作为长度的行向量返回<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-n" class="intrnllnk">n+ 1为低通和高通滤波器,2n带通和带阻滤波器为+ 1。

  • 对于数字滤波器,传递函数表示为b而且一个作为

    H z B z 一个 z b (1) + b (2) z 1 + + b (n + 1) z n (1) + (2) z 1 + + (n + 1) z n

  • 对于模拟滤波器,传递函数用b而且一个作为

    H 年代 B 年代 一个 年代 b (1) 年代 n + b (2) 年代 n 1 + + b (n + 1) (1) 年代 n + (2) 年代 n 1 + + (n + 1)

数据类型:

滤波器的零点、极点和增益,作为长度的两个列向量返回<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-n" class="intrnllnk">n(2n对于带通和带阻设计)和标量。

  • 对于数字滤波器,传递函数表示为zp,k作为

    H z k 1 z (1) z 1 1 z (2) z 1 1 z (n) z 1 1 (1页) z 1 1 (2页) z 1 1 p (n) z 1

  • 对于模拟滤波器,传递函数用zp,k作为

    H 年代 k 年代 z (1) 年代 z (2) 年代 z (n) 年代 (1页) 年代 (2页) 年代 p (n)

数据类型:

过滤器的状态空间表示,以矩阵形式返回。如果=<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj_sep_shared-n" class="intrnllnk">n用于低通和高通设计= 2n对于带通和带阻滤波器,然后一个×B×1,C等于1 ×,D等于1 × 1。

  • 对于数字滤波器,状态空间矩阵与状态向量相关联x,输入u,输出y通过

    x k + 1 一个 x k + B u k y k C x k + D u k

  • 对于模拟滤波器,状态空间矩阵与状态向量相关联x,输入u,输出y通过

    x ˙ 一个 x + B u y C x + D u

数据类型:

更多关于

全部折叠

限制

传递函数语法的数值不稳定性

一般情况下,使用(z,磷、钾)语法来设计IIR过滤器。要分析或实现筛选器,可以使用(z,磷、钾)输出与zp2sos.属性来设计过滤器[b]语法,你可能会遇到数值问题。这些问题是由于舍入错误造成的,可能发生在n低至4。下面的例子说明了这种限制。

N = 6;Rs = 80;Wn = [2.5e6 29e6]/500e6;ftype =<年代p一个n年代tyle="color:#A020F0">“带通”;<年代p一个n年代tyle="color:#228B22">传递函数设计[b,a] = cheby2(n,Rs,Wn,ftype);<年代p一个n年代tyle="color:#228B22">此过滤器不稳定%零极增益设计[z,p,k] = cheby2(n,Rs,Wn,ftype);SOS = zp2sos(z,p,k);<年代p一个n年代tyle="color:#228B22">绘制并比较结果HFVT = fvtool(b,a,sos,<年代p一个n年代tyle="color:#A020F0">“FrequencyScale”,<年代p一个n年代tyle="color:#A020F0">“日志”);传奇(hfvt<年代p一个n年代tyle="color:#A020F0">“特遣部队设计”,<年代p一个n年代tyle="color:#A020F0">“ZPK设计”)

图1:幅度响应(dB)包含一个轴对象。axis对象的标题为Magnitude Response (dB), xlabel归一化频率(乘以pi blank r a d / s a m p le), ylabel Magnitude (dB)包含2个类型为line的对象。这些对象分别代表TF Design, ZPK Design。

算法

切比雪夫II型滤波器在通带是单调的,在阻带是等纹波的。II型滤波器不像I型滤波器那样滚转快,但没有通带纹波。

cheby2使用五步算法:

  1. 它使用该函数找到低通模拟原型极点、零点和增益<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/cheb2ap.html">cheb2ap

  2. 它将极点、零点和增益转换为状态空间形式。

  3. 如果需要,它使用状态空间转换将低通滤波器转换为具有所需频率约束的带通、高通或带阻滤波器。

  4. 对于数字滤波器设计,它使用<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/bilinear.html">双线性通过频率预翘曲的双线性变换将模拟滤波器转换为数字滤波器。小心地调整模拟滤波器和数字滤波器的频率,使之具有相同的频率响应幅度<一个href="//www.tianjin-qmedu.com/se/help/signal/ref/#bucr0qj-Ws" class="intrnllnk">Wsw1而且w2

  5. 它根据需要将状态空间滤波器转换回传递函数或零极增益形式。

扩展功能

版本历史

R2006a之前介绍

另请参阅

|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">|<年代p一个nitemscope itemtype="//www.tianjin-qmedu.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">