当前位置：文档之家 › 用MATLAB设计巴特沃斯低通滤波器

用MATLAB设计巴特沃斯低通滤波器

  • 格式：pdf
  • 大小：169.14 KB
  • 文档页数：4

⽤MATLAB设计巴特沃斯低通滤波器

⽤MATLAB 设计巴特沃斯低通滤波器1 巴特沃斯低通滤波器的特性

⼀个理想低通滤波器的幅频特性如图3-80的阴影部分所⽰。为了实现这个理想低通特性，需要在从0~ωC 的整个频带内增强增益，在ω>ωC 增益要降到0。实际上，理想滤波器是不可能实现的。图3-78是实际滤波器的幅频特性。但是实际滤波器的特性愈接近理想特性愈好，巴特沃斯（Butterworth ）滤波器就是解决这个问题的⽅法之⼀。

巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数，巴特沃斯的低通模平⽅函数为：221

|()|

1,2,,1(/)

N

C H j N j j ωωω==+ (3-138)

式中以C ω是滤波器的电压-3dB 点或半功率点。

不同阶次的巴特沃斯滤波器特性如图3-79(a)所⽰。4阶巴特沃斯滤波器的极点分布如图3-79(b)所⽰。

巴特沃斯滤波器幅频响应有以下特点： 最⼤平坦性：在0=ω附近⼀段范围内是⾮常平直的，它以原点的最⼤平坦性来逼近理想低通滤波器。 通带、阻带下降的单调性。这种滤波器具有良好的相频特性。

3dB 的不变性：随着N 的增加，频带边缘下降越陡峭，越接近理想特性。但不管N 是多少，幅频特性都通过-3dB 点。

极点配置在半径为ωC 的圆上，并且均匀分布。左半平⾯上的N 个极点是)(s H 的极点，右半平⾯上的N 个极点是)(s H -的极点。2 巴特沃斯低通滤波器的实现

为使巴特沃斯滤波器实⽤，我们必须能够实现它。⼀个较好的⽅法是将巴特沃斯滤波器函数化成若⼲⼆阶节级联，其中每⼀节实现⼀对共轭复极点。通过将极点以共轭复数的形式配对，对所有的每⼀个⼆阶节都具有实系数。1图3-78 低通滤波器的幅频特性

图3-80所⽰运算放⼤器电路为实现⼀对共轭极点提供了很好的⽅法。电路的系统函数为

2

02202121121122121)(1)11(1

)(ωωω++=

+++=s Q

s C C R R s C R C R s C C R R s H (3-139)

式中，ω0是S 平⾯原点与极点之间的距离，Q 被称为电路的“品质因数”，它提供了对响应

峰值尖锐程度的⼀种度量。⾼Q 值滤波器的极点距离ω轴较近，它抬⾼了所在频率处的幅度响应值。

电路中元件值的确定，⼀种较简单的⽅法是令R R R ==21，则有 Q RC 012ω=，22

121ω=C C R (3-140) 可得

012ωR Q C =

，0221

ωQR C = (3-141) 巴特沃斯滤波器的极点与原点距离为C ω，所以，0ωω=C 。

巴特沃斯滤波器设计步骤为1．根据需要确定截⽌频率C f 或C ω，巴特沃斯滤波器的阶数。 2．求出巴特沃斯滤波器的系统函数。

3．求出极点，将极点以共轭复数的形式配对。

4．求出⼆阶的系统函数的各分⼦分母系数。并计算Q 值。 5．给R R R ==21分配⼀个实际电阻值。 6．按（3-141）式计算C1、C 2的值。

3 设计举例

例 设计4阶巴特沃斯滤波器。截⽌频率3000=C f Hz ，⽤图3-80的运算放⼤器电路实现。设Ω===5600021R R R 。

解 MATLAB 可以⾃动完成这个设计计算。程序如下% 设计4阶巴特沃斯滤波器 wc=3000*2*pi;R=56000;[b,a]=butter(4,wc,'s'); % 求系统函数的系数

图3-80 ⼆阶节的电路-+i

u 0

sys=tf(b,a) % 系统函数的多项式形式

figure(1)

pzmap(sys), % 画零极点图

hold on;u=0:pi/200:2*pi;

r=wc*exp(i*u); % 画半径为wc的圆

plot(r,':'),axis('equal') % 使x,y轴的刻度相等

title('4阶巴特沃斯滤波器极点分布图')

disp('零点和极点模型')

[z p k]=tf2zp(b,a),format short g % 多项式模型转换成零点和极点模型disp('第1个⼆阶节参数')

a1=poly([p(1),p(2)]),Q1=wc/a1(2),C1=2*Q1/R/wc,C2=1/(2*Q1*R*wc) disp('第2个⼆阶节参数')

a2=poly([p(3),p(4)]),Q2=wc/a2(2),C1=2*Q2/R/wc,C2=1/(2*Q2*R*wc) figure(2)

f=linspace(0,6000,300);

H1=freqs(a1(3),a1,2*pi*f); % 求第1个⼆阶节频率响应值H2=freqs(a2(3),a2,2*pi*f); % 求第2个⼆阶节频率响应值plot(f,abs(H1),'r:',f,abs(H2),'r:',f,abs(H1).*abs(H2))

set(gca,'xtick',[0 3000 6000]);

set(gca,'ytick',[0 0.707 1 1.5]);grid

xlabel('频率f(Hz)');

ylabel('幅度');

title('幅频响应曲线')

运⾏程序可得到如下输出Transfer function:

1.262e017

------------------------------------------------------------

s^4 + 4.926e004 s^3 + 1.213e009 s^2 + 1.75e013 s + 1.262e017

零点和极点模型z =

Empty matrix: 0-by-1

p =

-7213.4 + 17415i

-7213.4 - 17415i

-17415 + 7213.4i

-17415 - 7213.4ik =

1.2624e+017

第1个⼆阶节参数a1 = 1 14427 3.5531e+008

Q1 = 1.3066 C1 = 2.4755e-009 C2 = 3.6254e-010 第2个⼆阶节参数

a2 = 1 34829 3.5531e+008 Q2 = 0.5412 C1 = 1.0254e-009 C2 = 8.7524e-010

运⾏程序可得到如下输出图形如图3-81所⽰。

以上计算的元件值都是实际的，所以这个滤波器是能够实现的。图3-83给出了两级幅度响应曲线（虚线），⽽总的响应曲线（实线）证实了4阶巴特沃斯滤波器响应。第1级的Q 值为1.31较⾼，极点距离ω轴较近，它的幅度响应尖锐。⽽第2级的Q 值为0.54较⼩，极点距离ω轴较远，其幅度响应平缓。在实际中，最好将⾼Q 值⼀级排在后，这样可以降低⾼增益导致滤波器硬件饱和的危险。

除了巴特沃斯滤波器来逼迫理想滤波器之外，常⽤的滤波器还有切⽐雪夫、贝塞尔以及椭圆滤波器。它们的设计⽅法与巴特沃斯滤波器的设计⽅法类似，更详细的介绍请参考有关模拟滤波器设计的资料。

图3-81 4阶巴特沃斯滤波器零极点图和幅频响应 (b) 两个2阶节及乘积的幅频响应(a) 4阶巴特沃斯滤波器极点图

合集下载

基于MATLAB的切比雪夫II型数字低通滤波器设计

基于MATLAB的切比雪夫II型数字低通滤波器设计

科技视界Science&TechnologyVision

Science&TechnologyVision科技视

界渊上接第怨源页冤响袁其中学校和家庭是重要因素袁所以学校与家庭要建立共育机制遥学校是学生们学习尧生活的主要场所袁在积极创造学习生活条件袁开展相应工作的同时袁学校方面还应及时和家长沟通袁了解学生生长环境尧性格爱好袁并反馈学生在校学习尧生活和心理状况袁与家长共同教育管理学生袁必要时要共同商讨学生的成长计划遥对于有心理问题的学生袁学校要加强关注的力度,及时与家长取得联系,共同采取有效的干预措施,将各种心理问题扼杀在萌芽中遥总之袁做好大学生的心理健康教育工作袁高校辅导员应及时了解学生的心理活动袁学校要健全一系列心理健康教育和问题解决机制袁并及时与家长沟通袁针对学生的不同心理问题给予有效的指导袁确保大学生以积极的心态面对学习尧生活袁为今后更顺利地步入社会奠定良好的基础遥揖参考文献铱咱员暂李逸龙,姚海田,等.大学生教育管理与发展指导案例[M].东营:中国石油大学出版社,2012:2-3.咱圆暂陈小菊,丁留贯.高校辅导员参与大学生心理健康教育工作探析[J].文教资料,2009(7):199-199.咱猿暂张东伟.高校辅导员在大学生心理健康教育中的作用[J].教育理论研究,2011年(1):111-112.咱源暂谭平.论高校心理健康教育课程的构建[J].理论探讨,2008(12).咱缘暂张猛,杨琳.新时期高校辅导员工作的创新研究[J].中国科技信息,2007(10):190-192.咱责任编辑院杨扬暂干扰抑制常见的模拟滤波器是巴特沃斯(Butterworth)滤波器和切比雪夫(Chebyshev)滤波器遥巴特沃斯滤波器的特点是具有通带内最大平坦的振幅特性袁且随频率袁升高袁幅频特性单调递减遥切比雪夫滤波器在通带范围内是等幅起伏的袁所以同样的通带衰减袁其阶数较巴特沃斯滤波器要小遥可根据需要对通带内允许的衰减量(波动范围)提出要求袁如要求波动范围小于1dB[1,2]遥MATLAB是美国MathWorks公司推出的一套用于工程计算的可视化高性能语言与软件环境袁是数字信号处理技术实现的重要手段[3]遥本文采用脉冲响应不变法实现Chebyshev数字滤波器的设计遥1程序设计及运行结果根据研究任务袁利用切比雪夫模拟滤波器袁通过脉冲响应不变法设计切比雪夫II型数字低通滤波器遥数字滤波器的技术指标为袁0.90臆H(ejw)臆1袁0臆w臆0.25仔袁H(ejw)臆0.18袁园.35仔臆w臆仔袁采样周期为T=2遥设计程序如下院T=2;%设置采样周期为2fs=1/T;%采样频率为周期倒数Wp=0.25*pi/T;Ws=0.35*pi/T;%设置归一化通带和阻带截止频率Ap=20*log10(1/0.9);As=20*log10(1/0.18);%设置通带最大和最小衰减[N,Wc]=cheb2ord(Wp,Ws,Ap,As);%调用cheb2ord函数确定切比雪夫II型滤波器阶数[B,A]=cheby2(N,Ap,Wc);%调用cheby2函数设计切比雪夫II型滤波器W=linspace(0,pi,400*pi);%指定一段频率值hf=freqs(B,A,W);%计算模拟滤波器的幅频响应subplot(2,1,1);plot(W/pi,abs(hf)/abs(hf(1)));%绘出切比雪夫II型模拟滤波器的幅频特性曲线gridon;title('切比雪夫II型模拟滤波器');xlabel('Frequency/Hz');ylabel('Magnitude');[D,C]=impinvar(B,A,fs);%调用脉冲响应不变法Hz=freqz(D,C,W);%返回频率响应subplot(2,1,2);plot(W/pi,abs(Hz)/abs(Hz(1)));%绘出切比雪夫II型数字低通滤波器的幅频特性曲线gridon;title('切比雪夫II型数字滤波器');xlabel('Frequency/Hz');ylabel('Magnitude');2结果分析运行程序袁可得滤波器阶数为N=3袁Wc=0.5498遥符合设计要求遥切比雪夫II型数字低通滤波器的幅频特性曲线如图1所示遥图1切比雪夫II型数字低通滤波器的幅频特性曲线3结语利用MATLAB设计滤波器方法简单尧快捷直观遥本文运用脉冲响应不变法袁利用切比雪夫模拟滤波器设计了切比雪夫II型数字低通滤波器袁程序运行结果符合设计技术指标要求袁取得了较理想的实验效果遥揖参考文献铱咱员暂李勇.MATLAB辅助现代工程数字信号处理[M].西安:西安电子科技大学出版社袁2002:83-95.咱2暂董霖.MATLAB使用详解[M].北京:科学出版社袁2008:507-513.咱3暂石云霞.MATLAB在滤波器设计中的应用[J].青岛建筑工程学院学报袁2004袁25(2)院93-96.咱责任编辑院汤静暂基于MATLAB的切比雪夫II型数字低通滤波器设计王艳文1史先红2(1.新乡医学院生命科学技术学院袁河南新乡453003曰2.许昌学院公共实验中心,河南许昌461000)揖摘要铱本文利用脉冲响应不变法实现了切比雪夫II型数字滤波器的设计袁设计结果符合数字滤波器技术指标要求遥揖关键词铱MATLAB曰切比雪夫II型曰数字滤波器

matlab数字信号处理85个实用案例精讲

matlab数字信号处理85个实用案例精讲

matlab数字信号处理85个实用案例精讲

MATLAB数字信号处理85个实用案例精讲

MATLAB是一种强大的数学软件，广泛应用于数字信号处理领域。本文将介绍85个实用案例，涵盖了数字信号处理的各个方面，包括信号生成、滤波、频谱分析、时频分析、数字滤波器设计等。

1. 信号生成

案例：生成正弦信号

在MATLAB中，可以使用sin函数生成正弦信号。例如，生成频率为100Hz，幅度为1的正弦信号，代码如下：

t = 0:0.001:1;

f = 100;

x = sin(2*pi*f*t);

2. 滤波

案例：低通滤波

低通滤波器可以滤除高频信号，保留低频信号。在MATLAB中，可以使用fir1函数设计低通滤波器。例如，设计截止频率为100Hz的低通滤波器，代码如下：

fs = 1000;

fc = 100;

N = 100;

b = fir1(N, fc/(fs/2), 'low');

3. 频谱分析

案例：计算功率谱密度

功率谱密度是信号在频域上的能量分布。在MATLAB中，可以使用pwelch函数计算功率谱密度。例如，计算频率为100Hz的正弦信号的功率谱密度，代码如下：

t = 0:0.001:1;

f = 100;

x = sin(2*pi*f*t);

[Pxx, f] = pwelch(x, [], [], [], 1000);

4. 时频分析

案例：计算短时傅里叶变换

短时傅里叶变换可以分析信号在时间和频率上的变化。在MATLAB中，可以使用spectrogram函数计算短时傅里叶变换。例如，计算频率为100Hz的正弦信号的短时傅里叶变换，代码如下：

t = 0:0.001:1;

f = 100;

x = sin(2*pi*f*t);

spectrogram(x, [], [], [], 1000, 'yaxis');

5. 数字滤波器设计

案例：设计巴特沃斯滤波器

巴特沃斯滤波器是一种常用的数字滤波器，可以实现平滑滤波和带通滤波。在MATLAB中，可以使用butter函数设计巴特沃斯滤波器。例如，设计截止频率为100Hz的巴特沃斯低通滤波器，代码如下：

巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计

巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计

巴特沃斯Ⅱ型低通滤波器和切比雪夫Ⅱ型低通滤波器IIR低通数字滤波器设计(总24页)

--本页仅作为文档封面，使用时请直接删除即可--

--内页可以根据需求调整合适字体及大小--

南华大学

课 程 设 计 报 告

课程名称： 数字通信课程设计

设计名称： IIR低通数字滤波器设计

姓 名： XXXX

学 号: xxxx

班 级： xxxx

指导教师： XXXX

起止日期：

南华大学电气工程学院制课 程 设 计 任 务 书

学生班级： xxxx 学生姓名： phatonic 学号： XXXXXXXX

设计名称： IIR低通数字滤波器设计

起止日期： 指导教师： XX

设计要求：

用冲激响应不变法和双线性变换法分别以巴特沃斯滤波器、切比雪夫Ⅰ型滤波器、切比雪夫Ⅱ型滤波器为原型设计IIR低通数字滤波器。

1、IIR低通数字滤波器指标为：

0.20.31dB15dBpspsRdBAdB数字通带截止频率（弧度）数字阻带截止频率（弧度）通带衰减（）阻带衰减（）；

2、采用冲激响应不变法分别设计上述三种数字滤波器；

3、分别绘制这三种数字滤波器的幅度响应曲线和相位响应曲线；

模拟有源滤波器设计的MATLAB实现课程设计（已处理）

模拟有源滤波器设计的MATLAB实现课程设计（已处理）

模拟有源滤波器设计的MATLAB实现课程设计

摘要 本文阐述了滤波器的基本概念,介绍了模拟有源滤波器的设计原理和逼

近理论。其中包括巴特沃思逼近、切比雪夫I型逼近、切比雪夫Ⅱ型逼近、椭圆

函数逼近和贝塞尔逼近。并研究了模拟有源滤波器的设计流程及性能测试。 综合了传统的硬件设计方法与软件编程技术,由MATLAB仿真出了各种滤

波器逼近技术的幅频特性曲线并进行了实例分析。对巴特沃思滤波器实例的研究仿真,由程序快速的得到了最小阶数和截止频率,取代了传统繁复的计算;方便的

实现了由模拟低通滤波器向高通、低通和带阻滤波器的转换;对四运放复杂电路

进行了设计仿真,通过求取其不同点的输出传递函数,模拟了二阶低通、高通、带通和带阻滤波器的幅频特性曲线并得到了较好的仿真结果。

关键词:模拟有源滤波器;逼近理论;幅频特性;MATLAB程序设计

ABSTRACT

This paper describes the basic concepts of filters introduced

analog active filter design principles and approximation theory. Including Butterworth approximation, Chebyshev type I approximation,

Chebyshev Ⅱ type approximation, elliptic function approximation and Bezier approximation. The first author studied the analog active filter

design process and performance testing.

Combines the traditional hardware design methods and software programming techniques, the MATLAB simulation of a variety of filter