黄油

巴特沃斯滤波器设计

折叠所有页面

语法

[b,a] =黄油(n,Wn)

[b,a] =黄油(n,Wn,ftype)

[z,p,k] =黄油(<年代p一个nclass="argument_placeholder">___）

[A,B,C,D] =黄油(<年代p一个nclass="argument_placeholder">___）

［<年代p一个nclass="argument_placeholder">___=黄油<年代p一个nclass="argument_placeholder">___, ' s ')

描述

例子

［<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-ba" class="intrnllnk">b,一个=黄油<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-n" class="intrnllnk">n，<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u-Wn" class="intrnllnk">Wn）返回an的传递函数系数n一阶归一化截止频率低通数字巴特沃斯滤波器Wn．

例子

［<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-ba" class="intrnllnk">b,一个=黄油<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-n" class="intrnllnk">n，<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u-Wn" class="intrnllnk">Wn，<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u-ftype" class="intrnllnk">ftype）的值来设计低通、高通、带通或带阻巴特沃斯滤波器ftype元素的个数Wn．所得到的带通和带阻设计是2阶的n．

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

例子

［<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-zpk" class="intrnllnk">z,磷、钾=黄油<年代p一个nclass="argument_placeholder">___）设计一个低通、高通、带通或带阻数字巴特沃斯滤波器，并返回其零点、极点和增益。此语法可以包括以前语法中的任何输入参数。

例子

［<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-ABCD" class="intrnllnk">A, B, C, D=黄油<年代p一个nclass="argument_placeholder">___）设计低通、高通、带通或带阻数字巴特沃斯滤波器，并返回指定其状态空间表示的矩阵。

例子

［<年代p一个nclass="argument_placeholder">___=黄油<年代p一个nclass="argument_placeholder">___, ' s ')设计一个具有截止角频率的低通、高通、带通或带阻模拟巴特沃斯滤波器<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u-Wn" class="intrnllnk">Wn．

例子

全部折叠

低通巴特沃斯传递函数

打开实时脚本

设计一个6阶低通巴特沃斯滤波器，其截止频率为300hz，对于1000hz采样的数据，对应于<年代p一个nclass="inlineequation"> $0 ． 6 π$ rad /样品。画出它的幅值和相位响应。用它来过滤1000个样本的随机信号。

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

图中包含2个轴对象。Axes对象1包含一个line类型的对象。坐标轴对象2包含一个line类型的对象。

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

带阻巴特沃斯滤波器

打开实时脚本

设计一个边频归一化的六阶巴特沃斯带阻滤波器<年代p一个nclass="inlineequation"> $0 ． 2 π$ 而且<年代p一个nclass="inlineequation"> $0 ． 6 π$ rad /样品。画出它的幅值和相位响应。使用它来过滤随机数据。

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

图中包含2个轴对象。Axes对象1包含一个line类型的对象。坐标轴对象2包含一个line类型的对象。

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

高通巴特沃斯滤波器

打开实时脚本

设计一个9阶高通巴特沃斯滤波器。指定300hz的截止频率，对于1000hz采样的数据，对应于<年代p一个nclass="inlineequation"> $0 ． 6 π$ rad /样品。绘制幅值和相位响应。将零点、极点和增益转换为二阶部分，以供使用fvtool．

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

图幅度响应(dB)和相位响应包含一个轴对象。标题为Magnitude Response (dB)和Phase Response的axis对象包含一个类型为line的对象。

带通巴特沃斯滤波器

打开实时脚本

设计一个20阶巴特沃斯带通滤波器，其较低的截止频率为500hz，较高的截止频率为560hz。指定采样率为1500hz。使用状态空间表示。设计一个相同的过滤器使用designfilt．

[A,B,C,D] =黄油(10，[500 560]/750);D = designfilt(<年代p一个n年代tyle="color:#A020F0">“bandpassiir”，<年代p一个n年代tyle="color:#A020F0">“FilterOrder”, 20岁,<年代p一个n年代tyle="color:#0000FF">.．.“HalfPowerFrequency1”, 500,<年代p一个n年代tyle="color:#A020F0">“HalfPowerFrequency2”, 560,<年代p一个n年代tyle="color:#0000FF">.．.“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">“黄油”，<年代p一个n年代tyle="color:#A020F0">“designfilt”）

图幅度响应(dB)包含一个轴对象。标题为Magnitude Response (dB)的axis对象包含2个类型为line的对象。这些对象代表黄油，设计。

模拟IIR低通滤波器的比较

打开实时脚本

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

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

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

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

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

[z2,p2,k2] = cheby2(n,30,2*pi*f，<年代p一个n年代tyle="color:#A020F0">“年代”）;[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*f，<年代p一个n年代tyle="color:#A020F0">“年代”）;[be,ae] = zp2tf(ze,pe,ke);[he,we] = freqs(be,ae,4096);

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

情节(wb / (2 e9 *π)mag2db (abs (hb)))<年代p一个n年代tyle="color:#A020F0">在图(w1/(2e9*pi)，mag2db(abs(h1)))图(w2/(2e9*pi)，mag2db(abs(h2)))图(we/(2e9*pi)，mag2db(abs(he)))轴([0 4 -40 5])<年代p一个n年代tyle="color:#A020F0">“频率(GHz)”) ylabel (<年代p一个n年代tyle="color:#A020F0">“衰减(dB)”)传说(<年代p一个n年代tyle="color:#A020F0">“黄油”，<年代p一个n年代tyle="color:#A020F0">“cheby1”，<年代p一个n年代tyle="color:#A020F0">“cheby2”，<年代p一个n年代tyle="color:#A020F0">“ellip”）

图中包含一个轴对象。axis对象包含4个line类型的对象。这些对象表示butter、cheby1、cheby2、ellipp。

巴特沃斯和切比雪夫II型滤波器有平坦的通带和宽的过渡带。切比雪夫I型和椭圆滤波器滚出更快，但有通带纹波。切比雪夫II型设计函数的频率输入设置止带的开始，而不是通带的结束。

输入参数

全部折叠

`n`- - - - - -<年代p一个nitemprop="purpose">过滤器订单
整数标量

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

数据类型:双

`Wn`- - - - - -<年代p一个nitemprop="purpose">截止频率
标量|<年代p一个nitemprop="inputvalue">双元素向量

截止频率，指定为标量或两元向量。截止频率是滤波器的幅值响应所在的频率<年代p一个nclass="inlineequation">1 /√2．

如果Wn是标量，那么黄油设计具有截止频率的低通或高通滤波器Wn．
如果Wn是二元向量吗(w1 w2),在那里w1<w2,然后黄油设计一种具有较低截止频率的带通或带阻滤波器w1以及更高的截止频率w2．
对于数字滤波器，截止频率必须介于0和1之间，其中1对应于奈奎斯特速率-采样速率的一半或πrad /样品。
对于模拟滤波器，截止频率必须以弧度/秒表示，并且可以取任何正值。

数据类型:双

`ftype`- - - - - -<年代p一个nitemprop="purpose">过滤器类型
`“低”`|<年代p一个nitemprop="inputvalue">`“带通”`|<年代p一个nitemprop="inputvalue">`“高”`|<年代p一个nitemprop="inputvalue">`“停止”`

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

“低”指定具有截止频率的低通滤波器<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u-Wn" class="intrnllnk">Wn．“低”是否默认为标量Wn．
“高”指定具有截止频率的高通滤波器Wn．
“带通”指定阶为2的带通滤波器<一个href="//www.tatmou.com/nl/help/signal/ref/#bucse3u_sep_shared-n" class="intrnllnk">n如果Wn是一个二元向量。“带通”是默认值，当Wn有两个元素。
“停止”指定阶为2的带阻滤波器n如果Wn是一个二元向量。

输出参数

全部折叠

`b,一个`-传递函数系数
行向量

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

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

$H （ z ）＝ \frac{B （ z ）}{一个（ z ）} ＝ \frac{b (1) + b (2) z^{- 1} + \dots + b (n + 1) z^{- n}}{(1) + (2) z^{- 1} + \dots + (n + 1) z^{- n}} ．$
对于模拟滤波器，传递函数用b而且一个作为

$H （年代）＝ \frac{B （年代）}{一个（年代）} ＝ \frac{b (1) {年代}^{n} + b (2) {年代}^{n - 1} + \dots + b (n + 1)}{(1) {年代}^{n} + (2) {年代}^{n - 1} + \dots + (n + 1)} ．$

数据类型:双

`z,磷、钾`-零，极点和增益
列向量，标量

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

对于数字滤波器，传递函数表示为z，p,k作为

$H （ z ）＝ k \frac{（ 1 - z (1) z^{- 1} ）（ 1 - z (2) z^{- 1} ） \dots （ 1 - z (n) z^{- 1} ）}{（ 1 - (1页) z^{- 1} ）（ 1 - (2页) z^{- 1} ） \dots （ 1 - p (n) z^{- 1} ）} ．$
对于模拟滤波器，传递函数用z，p,k作为

$H （年代）＝ k \frac{（年代 - z (1) ）（年代 - z (2) ） \dots （年代 - z (n) ）}{（年代 - (1页) ）（年代 - (2页) ） \dots （年代 - p (n) ）} ．$

数据类型:双

`A, B, C, D`-状态空间矩阵
矩阵

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

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

$\begin{matrix} x （ k + 1 ）＝一个 x （ k ） + B u （ k ） \\ y （ k ）＝ C x （ k ） + D u （ k ）． \end{matrix}$
对于模拟滤波器，状态空间矩阵与状态向量相关联x，输入u，输出y通过

$\begin{array}{l} \dot{x} ＝一个 x + B u \\ y ＝ C x + D u ． \end{array}$

数据类型:双

算法

巴特沃斯滤波器的幅值响应在通带内最大程度地平坦，总体上是单调的。这种平滑度是以降低滚压陡度为代价的。椭圆滤波器和切比雪夫滤波器通常为给定的滤波器顺序提供更陡的滚转。

黄油使用五步算法:

它使用该函数找到低通模拟原型极点、零点和增益<一个href="//www.tatmou.com/nl/nl/help/signal/ref/buttap.html">buttap．
它将极点、零点和增益转换为状态空间形式。
如果需要，它使用状态空间转换将低通滤波器转换为具有所需频率约束的带通、高通或带阻滤波器。
对于数字滤波器设计，它使用<一个href="//www.tatmou.com/nl/nl/help/signal/ref/bilinear.html">双线性通过频率预翘曲的双线性变换将模拟滤波器转换为数字滤波器。仔细的频率调整使模拟滤波器和数字滤波器具有相同的频率响应幅度Wn或在w1而且w2．
它根据需要将状态空间滤波器转换回其传递函数或零极增益形式。

扩展功能

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

版本历史

R2006a之前介绍

另请参阅

besself|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">buttap|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">buttord|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">cheby1|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">cheby2|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">designfilt|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">ellip|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">过滤器|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">maxflat|<年代p一个nitemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">sosfilt

黄油

语法

描述

例子

低通巴特沃斯传递函数

带阻巴特沃斯滤波器

高通巴特沃斯滤波器

带通巴特沃斯滤波器

模拟IIR低通滤波器的比较

输入参数

n- - - - - -<年代p一个nitemprop="purpose">过滤器订单整数标量

Wn- - - - - -<年代p一个nitemprop="purpose">截止频率标量|<年代p一个nitemprop="inputvalue">双元素向量

ftype- - - - - -<年代p一个nitemprop="purpose">过滤器类型“低”|<年代p一个nitemprop="inputvalue">“带通”|<年代p一个nitemprop="inputvalue">“高”|<年代p一个nitemprop="inputvalue">“停止”

输出参数

b,一个-传递函数系数行向量

z,磷、钾-零，极点和增益列向量，标量

A, B, C, D-状态空间矩阵矩阵

更多关于

限制

算法

扩展功能

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

版本历史

另请参阅

`n`- - - - - -<年代p一个nitemprop="purpose">过滤器订单
整数标量

`Wn`- - - - - -<年代p一个nitemprop="purpose">截止频率
标量|<年代p一个nitemprop="inputvalue">双元素向量

`ftype`- - - - - -<年代p一个nitemprop="purpose">过滤器类型
`“低”`|<年代p一个nitemprop="inputvalue">`“带通”`|<年代p一个nitemprop="inputvalue">`“高”`|<年代p一个nitemprop="inputvalue">`“停止”`

`b,一个`-传递函数系数
行向量

`z,磷、钾`-零，极点和增益
列向量，标量

`A, B, C, D`-状态空间矩阵
矩阵

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