用matlab编程设计切比雪夫低通滤波器..

目录

1摘要 (3)

2设计原理 (4)

2.1 切比雪夫滤波器介绍 (4)

2.2滤波器的分类 (5)

2.3 模拟滤波器的设计指标 (6)

3切比雪夫I型滤波器 (7)

3.1 切比雪夫滤波器的设计原理 (7)

3.2切比雪夫滤波器的设计步骤 (10)

3.3 用matlab编程设计切比雪夫低通滤波器 (11)

4 总结 (18)

5 参考文献 (18)

摘要

随着信息和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中，由于信号中经常混有各种复杂成分，因此很多信号的处理都是基于滤波器而进行的。所以，数字滤波器在数字信号处理中起着举足轻重的作用。

而数字滤波器的设计都要以模拟滤波器为基础的，这是因为模拟滤波器的理论和设计方方法都已发展的相当成熟，且有典型的模拟滤波器供我们选择。，如巴特沃思滤波器、切比雪夫滤波器等。

本次课程设计将运用MATLAB设计一个基于切比雪夫低通滤波器，并出所设计滤波器的幅度及幅度衰减特性。

关键词：滤波切比雪夫模拟低通

1切比雪夫滤波器介绍

在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。因此，若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内，或同时在通带和阻带内都均匀分布，这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质：其频率响应的幅度既可以在通带中是等波纹的，而在阻带中是单调的（称为I 型切比雪夫滤波器），也可以在通带中是单调的，而在阻带中是等波纹的（称为II 型切比雪夫滤波器），其中切比雪夫II 型滤波器又称为逆切比雪夫滤波器。

I 型切比雪夫滤波器的幅度平方函数是

2|)(|Ωj H C =

)/(11

2

2c N C ΩΩ+ε (2.1) 其中ε是一个小于1的正数，它与通带波纹有关，ε越大，波纹也越大，式中为N 阶切比雪夫多项式，定义为

)cos cos()(1x N x C N -= (2.2) 当N 大于或等于1时，从定义切比雪夫多项式可以直接得出由

)

(x C N 和

)(1x C N -求)(1x C N +的递推公式。将三角恒等式代入 (2.2)式，得

)(1x C N +=2x -)(x C N )(1x C N - (2.3)

从 (2.2)式我们注意到，当0

x C N 在0和1之间变化；当x>1时，x 1cos -是虚数，所以)

(x C N 像双曲余弦一样单调地增加。参考(2.1)，

2|)(|Ωj H C 对于0≤p ΩΩ/≤1呈现出在1和1/（21ε+）之间的波动；而对于

p

ΩΩ/〉1单调地减小。需要用三个参量来确定该滤波器：ε，p Ω和N 。

在典型的设计中，用容许的通带波纹来确定ε，而用希望的通带截止频率来确定

c

Ω。然后选择合适的阶次N ，以便阻带的技术要求得到满足。

定义允许的通带最大衰减p α用下式表示： 2

2max ()10lg ,min ()

p p Ha j Ha j αΩ=Ω≤ΩΩ

为了求切比雪夫滤波器在椭圆上极点的位置，我们首先要这样确定，在大圆和小圆上以等角度 等间隔排列的那些点：这些点对于虚轴呈对称分布，并且没有一个点落在虚轴上；但当N 为奇数时要有一个点落在实轴上，而当N 为偶数时，就都不会落在实轴上。切比雪夫滤波器的极点落在椭圆上，起纵坐标由相应的大圆上点的纵坐标来表示，起横坐标由相应的小圆上点的横坐标来表示。

3切比雪夫I 型滤波器

3.1切比雪夫低通滤波器的设计原理

切比雪夫滤波器的幅频特性具有等波纹特性。它有两种形式：振幅特性在通带内是等波纹的，在阻带内是单调递减的切比雪夫I 型滤波器，振幅特性在阻带内是等波纹的，在通带内是单调递减的切比雪夫II 型滤波器，如图所示分别画出了滤波器的幅频特性和衰减函数。

以切比雪夫I 型为例介绍其设计原理 幅度平方函数用2

()H j Ω表示

（2-3-1）

式中，ε为小于1的正数，表示通带内幅度波动的程度，ε越大，波动幅度也越大。p Ω称为通带截止频率。令p

λΩ=Ω ，称为对p Ω的归一化频

率。

定义允许的通带内最大衰减p α用下式表示 2

2

max ()10lg ,min ()

p p Ha j Ha j αΩ=Ω≤ΩΩ （2-3-2）

式中

2

max ()1ha j Ω= （2-3-3）

2

2

1

min ()1ha j εΩ=+ （2-3-4） 因此

210lg(1)p αε=+ （2-3-5）

2

21()1(

)N p

H j C εΩ=

Ω

+Ω

0.1210

1p

αε=- （2-3-6）

这样，可以根据通带内最大衰减p α，可求出参数ε。 阶数N 影响过渡带的宽度，同时也影响通带内波动的疏密，因为N 等于通带内的最大值和最小值的总个数。设阻带的起点频率为s Ω，则有

2

2

2

1()1()

s s N p

H j C εΩ=

Ω

+Ω （2-3-7）

令s

s p λΩ=

Ω，由s λ>1，有

()()()

2

111N s s a s C ch Narch H j λλε

==-Ω （2-3-8）

可以解出

()()

2

111a s s arch H j N arch ελ??

?

?-??

Ω??=

（2-3-9）

()

2

1

111s p a s ch arch N H j ε??

?

???

?

?Ω=Ω-????Ω?????

?

（2-3-10）

3dB 截止频率用c Ω表示，

()2

12

a c H j Ω= （2-3-11）

按照（2-3-1）式，有

1

1c p ch arch N

ε??Ω=Ω ??? （2-3-12）

经过一系列推论

得归一化系统函数为

（2-3-13）

()()

1

1

1

2

a N

N i

i G p p p ε-==

-∏

去归一化的系统函数为 ()()

()

11

2p

N

p a s N

p N i p i H s G p p p αε=

-Ω=Ω==

-Ω∏ （2-3-14）

3.2切比雪夫低通滤波器的设计步骤

（1）确定低通滤波器的技术指标：边带频率p Ω，通带最大衰减p α、阻带最大衰减s α、阻带截至频率s Ω，它们满足

（2-3-15）

（2-3-16）

（2）求滤波器阶数N 和参数ε s

s p λΩ=Ω （2-3-17）

0.11

0.1101

1101

s p

K αα--=- （2-3-18） 1

1s archk N arch λ-= （2-3-19）

这样，先由（2--18）式求出11K -，代入（2-3-19），求出阶数N ，最后取大于或等于N 的最小整数。

（3）求归一化系统函数()a G p

()

2

110lg

p p Ha j α=≤Ω()

2

110lg

s s Ha j α=≤Ω

()()

11

1

2a N

N i i G p p p ε-==

-∏ （2-3-20）

（4）将去归一化，得到实际的()a H s ()()p

a s

p H s G p =Ω= （2-3-21）

3.3用MATLAB 设计切比雪夫低通滤波器

（1）Matlab 的信号处理工具箱提供了频谱分析函数:fft 函数、filter 函数和freqz 函数。

fft 函数

filter 函数功能：利用IIR 滤波器和FIR 滤波器对数据进行滤波。 调用格式：

y=filter(b,a,x) [y,zf]=filter(b,a,x) y=filter(b,a,x,zi)

说明：filter 采用数字滤波器对数据进行滤波，其实现采用移位直接Ⅱ型结构，因而适用于IIR 和FIR 滤波器。滤波器的系统函数为

n

n m

m z a z a z a z a z b z b z b z b b Z H --------+???+++++???++++=33221133221101)(

即滤波器系数a=[a0 a1 a2 ...an],b=[b0 b1 ...bm],输入序列矢量为x 。这里，标准形式为a0=1，如果输入矢量a 时，a0≠1，则MATLAB 将自动进行归一化系数的操作；如果a0=0，则给出出错信息。

y=filter(b,a,x)利用给定系数矢量a 和b 对x 中的数据进行滤波，结果放入y 矢量中，y 的长度取max(N,M)。

y=filter(b,a,x,zi)可在zi 中指定x 的初始状态。

[y,zf]=filter(b,a,x)除得到矢量y 外，还得到x 的最终状态矢量zf 。

freqz 函数功能：离散时间系统的频率响应。

格式：[h,w]=freqz(b,a,n)

[h,f]=freqz(b,a,n,Fs)

h=freqz(b,a,w)

h=freqz(b,a,f,Fs)

freqz(b,a,n)

说明: freqz 用于计算数字滤波器H(Z)的频率响应函数H(ejω)。

[h,w]=freqz(b,a,n)可得到数字滤波器的n点复频响应值，这n个点均匀地分布在[0,π]上,并将这n个频点的频率记录在w中，相应的频响值记录在h中。要求n为大于零的整数,最好为2的整数次幂,以便采用FFT计算，提高速度。缺省时n =512。

[h,f]=freqz(b,a,n,Fs)用于对H(ejω)在[0,Fs/2]上等间隔采样n点，采样点频率及相应频响值分别记录在f 和h中。由用户指定FS（以HZ为单位）值。

h=freqz(b,a,w)用于对H(ejω)在[0,2π]上进行采样，采样频率点由矢量w指定。

h=freqz(b,a,f,Fs) 用于对H(ejω)在[0,FS]上采样，采样频率点由矢量f指定。

freqz(b,a,n) 用于在当前图形窗口中绘制幅频和相频特性曲线。

fft函数函数功能：对信号进行傅里叶变换。

格式：fft(X)

fft(X,N)

fft(X,[],DIM)或fft(X,N,DIM)

说明:fft(X)是对输入信号X的离散傅里叶变换。

fft(X,N)是N点的傅里叶变换，如果X少于N点则补0凑齐位数，长于N点则截断。

如果x是个矩阵，列的长度将会以同样的方式调整，fft会对每列进行傅里叶变换，并返回一个相同维数的矩阵。

fft(X,[],DIM)或fft(X,N,DIM)是离散傅里叶变换在DIM尺度上的应用。DIM 可适应于任意维度的fft运算。

（2）确定低通滤波器的技术指标：

N-椭圆滤波器最小阶数；Wn为椭圆滤波器的带宽；

Wp-椭圆滤波器通带截止角频率；

Ws-椭圆滤波器阻带起始角频率；

Rp-通带波纹（dB）；

Rs-阻带最小衰减(dB)；

设计要求是低通滤波器，需要屏蔽的是15Hz和30Hz的波形，所以可令fp=5Hz，设fs=8Hz，Rp<0.1dB,Rs>40dB。

（3）连续信号的产生及采样:

clear

f1=5;

f2=15;

f3=30;

N=100;

fs=100;

n=(0:N-1); %采样点数为N

x1=sin(2*pi*f1*n/fs);

x2=sin(2*pi*f2*n/fs);

x3=sin(2*pi*f3*n/fs); %模拟信号转化为数字信号

x=x1+x2+x3; %信号叠加

plot(n,x); %作出时间-幅值图像

xlabel('时间(s)');

ylabel('幅值');

grid;

连续信号仿真图

(4)低通滤波器的设计

Wp=8;

Ws=12;

Rp=2;

Rs=40;%设置指标参数

[N,wp1]=cheb1ord(Wp/(fs/2),Ws/(fs/2),Rp,Rs);%利用cheb1ord函数求滤波器最小阶数

[b,a]=cheby1(N,Rp,wp1); %计算切比雪夫I型模拟低通滤波器系统函数系数

[H,f]=freqz(b,a,512,100); %求幅值H，频率f

plot(f,abs(H));

切比雪夫I型模拟低通滤波器仿真图

(5)对滤波后的信号进行分析和变换

sf=filter(b,a,x);%对原信号进行滤波

plot(n,sf);

xlabel('时间 (s)');

ylabel('幅值');

axis([0 1 -1 1]);

X=fft(x,512); %滤波前信号的傅里叶变换

SF=fft(sf,512); %滤波后信号的傅里叶变换

f=(0:255)/256*(fs/2);

plot(f,abs([X(1:256)' SF(1:256)'])); %滤波前后信号图像对比

xlabel('频率(Hz)');

ylabel('傅立叶变换图');

grid;

legend({'before','after'})

信号通过椭圆低通滤波器的仿真图

注：图中蓝色曲线代表滤波前的幅频曲线，绿线代表滤波后的幅频曲线。

3.4 设计结果分析

（1）求得阶数N=6，通带边界频率为wp = 0.1600，表明阶数越大，模拟低通滤波器越理想，且接近矩形。

（2）从图上可以看出：切比雪夫I型滤波器是在通带频率响应幅度等波纹的滤波器，通带内具有等波纹起伏特性，阻带内单调下降且衰减更大；幅

度在低频处通过，随着频率的增加，逐渐衰减至零；综上所看：本设计符合要求。

（3）切比雪夫滤波器的特点是在通带内，具有相等的波纹。截频衰减陡度比同阶数巴特沃斯特性更陡度比同阶数程时的衰减就超过6NdB。在阶数N 一定时，波纹越大，截频衰减陡度越陡。相位响应也是非线性，但较之比巴特沃斯为差。

总结

经过了将近一个周的课程设计，我学会的如何使用matlab来进行滤波器的设计，通过这次课程设计我不仅学到了许多关于matlab的知识，认识到了matlab的强大功能，更重要的是增强了我分析问题的能力，使我深刻地认识到仅仅学习课本上的知识是远远不够的，必须多多实践才能真正理解并掌握所学的知识，达到学以致用的目的，为以后的工作积累了宝贵的经验，同时我也深深地感受到严谨的态度对于科学研究的重要性。

这次课设不仅让我增长了知识，更让我对信号处理有了极大的兴趣，对我今后学习相关知识将会有很大的影响。在设计在学习中进步让我收获了知识，收获了满足，当看见一个题目在自己的手中实现让我充满成就感，也坚定了我继续深入学习的决心。

最后，感谢学校和老师给，让我能够锻炼自己，学习知识，丰富了我的学习生活。

切比雪夫1型数字低通滤波器

目录 1. 数字滤波器的设计任务及要求 (2) 2. 数字滤波器的设计及仿真 (2) 2.1数字滤波器的设计 (3) 2.2数字滤波器的性能分析 (3) 3. 数字滤波器的实现结构对其性能影响的分析 (8) 3.1数字滤波器的实现结构一及其幅频响应 (10) 3.2数字滤波器的实现结构二及其幅频响应 (12) 3.3 数字滤波器的实现结构对其性能影响的小结 (12) 4. 数字滤波器的参数字长对其性能影响的分析 (13) 4.1数字滤波器的实现结构一参数字长及幅频响应特性变化 4.2数字滤波器的实现结构二参数字长及幅频响应特性变化 4.3 数字滤波器的参数字长对其性能影响的小结 (16) 5. 结论及体会 (16) 5.1 滤波器设计、分析结论 (16) 5.2 我的体会 (16) 5.3 展望 (16)

1.数字滤波器的设计任务及要求 1. 设计说明 每位同学抽签得到一个四位数，由该四位数索引下表确定待设计数字滤波器的类型及其设计方法，然后用指定的设计方法完成滤波器设计。 要求：滤波器的设计指标： 低通： (1)通带截止频率πrad (id) pc 32 ln = ω (2)过渡带宽度πrad ) (i d 160 10log tz ≤?ω (3)滚降dB αroll 60= 其中，i d — 抽签得到那个四位数（学号的最末四位数），本设计中i d =0201。 2. 滤波器的初始设计通过手工计算完成； 3. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器 结构进行分析）； 4. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 5. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 6. 课程设计结束时提交设计说明书。 2.数字滤波器的设计及仿真 2.1数字滤波器（编号0201）的设计 数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序，而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。 本次课程设计使用MATLAB 信号处理箱和运用切比雪夫法设计数字滤波器，将手工计算一个切比雪夫I 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。

简单低通滤波器设计及matlab仿真

东北大学 研究生考试试卷 考试科目： 课程编号： 阅卷人: 考试日期： 姓名：xl 学号： 注意事项 1.考前研究生将上述项目填写清楚. 2.字迹要清楚,保持卷面清洁. 3.交卷时请将本试卷和题签一起上交. 4.课程考试后二周内授课教师完成评卷工作,公共课成绩单与试卷交研究生院培养办公室, 专业课成绩单与试卷交各学院,各学院把成绩单交研究生院培养办公室. 东北大学研究生院培养办公室

数字滤波器设计 技术指标： 通带最大衰减： =3dB ， 通带边界频率： ＝100Hz 阻带最小衰减： ＝20dB 阻带边界频率： ＝200Hz 采样频率：Fs=200Hz 目标： 1、根据性能指标设计一个巴特沃斯低通模拟滤波器。 2、通过双线性变换将该模拟滤波器转变为数字滤波器。 原理： 一、模拟滤波器设计 每一个滤波器的频率范围将直接取决于应用目的，因此必然是千差万别。为了使设计规范化，需要将滤波器的频率参数作归一化处理。设所给的实际频 率为Ω（或f ），归一化后的频率为λ，对低通模拟滤波器令λ＝p ΩΩ/，则1 =p λ， p s s ΩΩ=/λ。令归一化复数变量为p ，λj p =，则p p s j j p Ω=ΩΩ==//λ。所以巴 特沃思模拟低通滤波器的设计可按以下三个步骤来进行。 (1)将实际频率Ω规一化 (2)求Ωc 和N 11010/2-=P C α s p s N λααlg 1 10 110lg 10 /10/--= 这样Ωc 和N 可求。 p x fp s x s f

根据滤波器设计要求＝3dB ，则C ＝1，这样巴特沃思滤波器的设计就只剩一个参数N ，这时 N p N j G 222 )/(11 11)(ΩΩ+= += λλ (3)确定)(s G 因为λj p =，根据上面公式有 N N N p j p p G p G 22)1(11 )/(11)()(-+= += - 由 0)1(12=-+N N p 解得 )221 2exp(πN N k j p k -+=，k ＝1，2， (2) 这样可得 1 )21 2cos(21 ) )((1 )(21+-+-= --= -+πN N k p p p p p p p G k N k k 求得)(p G 后，用p s Ω/代替变量p ，即得实际需要得)(s G 。 二、双线性变换法 双线性变换法是将s 平面压缩变换到某一中介1s 平面的一条横带里，再通过标准变换关系)*1exp(T s z =将此带变换到整个z 平面上去，这样就使s 平面与z 平面之间建立一一对应的单值关系，消除了多值变换性。 为了将s 平面的Ωj 轴压缩到1s 平面的1Ωj 轴上的pi -到pi 一段上，可以通过以下的正切变换来实现： )21 tan(21T T Ω= Ω 这样当1Ω由T pi -经0变化到T pi 时，Ω由∞-经过0变化到∞+，也映射到了整个Ωj 轴。将这个关系延拓到整个s 平面和1s 平面，则可以得到

切比雪夫Ⅰ型滤波器和切比雪夫Ⅱ型滤波器的对比

设计流程图如下： 设计思想： 首先设计一个源信号和一个混合信号，通过其频谱对比得出最大和最小通带，最大和最小阻带;然后再根据得到的参数来设计切比雪夫滤波器，最后通过切比雪夫Ⅰ型滤波器和切比雪夫Ⅱ型滤波器的对比来得出那种效果好。 切比雪夫滤波器设计原理：

切比雪夫滤波器的振幅具有等波纹特性，它有两种形式：1)振幅特性在通带内是等波纹的、在阻带内是单调的切比雪夫I 型滤波器；2)振幅特性在通带内是单调的、在阻带内是等波纹的切比雪夫II型滤波器，采用何种形式的切比雪夫滤波器取决于实际用途. 切比雪夫滤波器的设计方法就是将逼近精确度均匀分布在整个通带内，或者均匀分布在整个阻带内，或者均匀分布在两者之内，这样就可以使滤波器阶数大大降低。 切比雪夫I型滤波器平方幅度响应函数表示为: 2 ) (Ωj G=[1+2εC2 N (Ω)]2/1- 其中ε<1(正数)，它与通带波纹有关,ε越大，波纹也越大；C N (Ω) 是切比雪夫多项式，它被定义为： C N (Ω)=cos(Narccos(Ω)),Ω≤1, C N (Ω)=cosh(Narcosh(Ω)),Ω>1. 而切比雪夫II型滤波器平方幅度响应函数表示为： ) (Ωj G2={1+2ε{ C2 N (Ω)/[2 N (Ω/c Ω)]2}}1- 其中ε<1(正数)，表示波纹变化情况；c Ω为截止频率；N为滤波器的阶次,也 是C N ( N Ω Ω/)的阶次。 源信号编码及其图形： t=-1:0.01:1 y=(cos(2*pi*10*t)+cos(2*pi*40*t)); N=length(y); fx=fft(y); df=100/N; n=0:N/2; f=n*df; subplot(2,1,1); plot(f,abs(fx(n+1))*2/N); grid; title('源波形频谱')

切比雪夫模拟高通滤波器matlab实验

题目滤波器设计与实现 学生姓名王朝全学号201005020228 专业班级电信1002 设计内容与要求一、设计内容： 设计一个模拟高通切比雪夫Ⅰ滤波器，技术指标：通带截止频率2000Hz，通带最大衰减1dB；过渡带700Hz，阻带最小衰减60dB，画出其幅度谱和相位谱。 二、设计要求： 1 设计报告一律按照规定的格式，使用A4纸，格式、封面统一给出模版。 2 报告内容 （1）设计题目及要求 （2）设计原理(包括滤波器工作原理、涉及到的matlab函数的说明) （3）设计内容（设计思路，设计流程、仿真结果） （4）设计总结（收获和体会） （5）参考文献 （6）程序清单 起止时间2012年12 月3日至2012年12月11 日指导教师签名2012年12月2日 系（教研室）主任签 名 年月日 学生签名年月日

一、设计题目及要求 设计一个模拟高通切比雪夫Ⅰ滤波器，技术指标：通带截止频率 2000Hz ，通带最大衰减1dB ；过渡带700Hz ，阻带最小衰减60dB ，画出其幅度谱和相位谱。 二、设计原理 (包括滤波器工作原理、涉及到的matlab 函数的说明) 2.1切比雪夫滤波器的内容： 2.1.1、切比雪夫滤波器的幅度平方函数： 从公式中我们可以发现它需要计算的 ε和N ， 110 1.0-=p a ε 边带频率p Ω，通带最大衰减p α=0.1dB 、阻带最大衰减s α=50dB 、 阻带截止频率s Ω，它们满足： ) /cosh() /110cosh(11.0p s a ar ar N s ΩΩ-≥ -ε 根据公式可以求出最小的切比雪夫I 型滤波器阶数的N 最小值。 2.1.2、切比雪夫滤波器的三个参量： Ωc 通带截止频率，ε表征通带内波纹大小由： 2 222 1 ()() 1( )a N p A H j C εΩ=Ω= Ω +Ω

切比雪夫I型低通滤波器设计解读

******************* 实践教学 ******************* 兰州理工大学 计算机与通信学院 2013年春季学期 信号处理课程设计 题目：切比雪夫I型低通滤波器设计 专业班级：通信工程三班 姓名： 学号： 指导教师：蔺莹 成绩：

摘要 本次课程设计将完成一个数字切比雪夫低通IIR 滤波器的设计，利用双线性变换和冲激响应不变法完成设计，并利用MATLAB 进行仿真。 已知数字滤波器的性能指标为：通带截止频率为: 0.4,1,0.45,15p p s P R dB R dB ω=π=ω=π=通带波动为阻带波动为，要求设计满足以上技术指标的切比雪夫I 型低通滤波器。绘制出理想冲激响应和实际冲激响应结果图。并且给出幅度响应结果图。 关键字：数字滤波器 切比雪夫 双线性变换 冲激响应不变

前言 (1) 一．数字滤波器 (2) 1．1 数字滤波器的概念 (2) 1．2数字滤波器的分类 (2) 1．3 IIR数字滤波器设计原理 (3) 二．切比雪夫滤波器 (5) 三．双线性变换法 (8) 四．脉冲响应不变法 (12) 五．切比雪夫低通滤波器的设计 (15) 5．1 程序流程图 (15) 5．2 设计步骤 (15) 六．总结 (18) 七．参考文献 (19) 致谢 (20) 附录 (21)

随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。目前数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器（DF，Digital Filter），根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR（Infinite Impulse Response）滤波器和有限冲激响应FIR（Finite Impulse Response）滤波器。与FIR 滤波器相比，IIR的实现采用的是递归结构，极点须在单位圆内，在相同设计指标下，实现IIR滤波器的阶次较低，即所用的存储单元少，从而经济效率高。MATLAB是英文MATrix LABoratory(矩阵实验室)的缩写。它是美国的MathWorks公司推出的一套用于科学计算和图形处理可视化、高性能语言与软件环境。MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常方便。

各类滤波器的MATLAB程序清单

各类滤波器的MATLAB程序 一、理想低通滤波器 IA=imread(''); [f1,f2]=freqspace(size(IA),'meshgrid'); Hd=ones(size(IA)); r=sqrt(f1.^2+f2.^2); Hd(r>=0; Y=fft2(double(IA)); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=ifft2(Ya); figure subplot(2,2,1),imshow(uint8(IA)); subplot(2,2,2),imshow(uint8(Ia)); figure surf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong'); 二、理想高通滤波器 IA=imread(''); [f1,f2]=freqspace(size(IA),'meshgrid'); Hd=ones(size(IA)); r=sqrt(f1.^2+f2.^2); Hd(r<=0; Y=fft2(double(IA));

Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya)); figure subplot(2,2,1),imshow(uint8(IA)); subplot(2,2,2),imshow(uint8(Ia)); figure surf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong'); 三、B utterworth低通滤波器 IA=imread(''); [f1,f2]=freqspace(size(IA),'meshgrid'); D=; r=f1.^2+f2.^2; n=4; for i=1:size(IA,1) for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1/(t^n+1); end end Y=fft2(double(IA)); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya));

matlab程序之——滤波器(带通-带阻)教学内容

m a t l a b程序之——滤波器(带通-带阻)

matlab程序之——滤波器（带通，带阻） 以下两个滤波器都是切比雪夫I型数字滤波器，不是巴特沃尔滤波器，请使用者注意！ 1.带通滤波器 function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs) %带通滤波 %使用注意事项：通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半 %即，f1,f3,fs1,fsh,的值小于 Fs/2 %x:需要带通滤波的序列 % f 1：通带左边界 % f 3：通带右边界 % fs1：衰减截止左边界 % fsh：衰变截止右边界 %rp：边带区衰减DB数设置 %rs：截止区衰减DB数设置 %FS：序列x的采样频率 % f1=300;f3=500;%通带截止频率上下限 % fsl=200;fsh=600;%阻带截止频率上下限 % rp=0.1;rs=30;%通带边衰减DB值和阻带边衰减DB值 % Fs=2000;%采样率 % wp1=2*pi*f1/Fs; wp3=2*pi*f3/Fs; wsl=2*pi*fsl/Fs; wsh=2*pi*fsh/Fs; wp=[wp1 wp3]; ws=[wsl wsh]; % % 设计切比雪夫滤波器； [n,wn]=cheb1ord(ws/pi,wp/pi,rp,rs); [bz1,az1]=cheby1(n,rp,wp/pi); %查看设计滤波器的曲线 [h,w]=freqz(bz1,az1,256,Fs); h=20*log10(abs(h));

figure;plot(w,h);title('所设计滤波器的通带曲线');grid on; y=filter(bz1,az1,x); end 带通滤波器使用例子 %-------------- %带通滤波器测试程序 fs=2000; t=(1:fs)/fs; ff1=100; ff2=400; ff3=700; x=sin(2*pi*ff1*t)+sin(2*pi*ff2*t)+sin(2*pi*ff3*t); figure; subplot(211);plot(t,x); subplot(212);hua_fft(x,fs,1); % y=filter(bz1,az1,x); y=bandp(x,300,500,200,600,0.1,30,fs); figure; subplot(211);plot(t,y); subplot(212);hua_fft(y,fs,1); %调用到的hua_fft()函数代码如下 function hua_fft(y,fs,style,varargin) %当style=1,画幅值谱；当style=2,画功率谱;当style=其他的，那么花幅值谱和功率谱 %当style=1时，还可以多输入2个可选参数 %可选输入参数是用来控制需要查看的频率段的 %第一个是需要查看的频率段起点 %第二个是需要查看的频率段的终点 %其他style不具备可选输入参数，如果输入发生位置错误 nfft= 2^nextpow2(length(y));%找出大于y的个数的最大的2的指数值（自动进算最佳FFT步长nfft） %nfft=1024;%人为设置FFT的步长nfft y=y-mean(y);%去除直流分量 y_ft=fft(y,nfft);%对y信号进行DFT，得到频率的幅值分布 y_p=y_ft.*conj(y_ft)/nfft;%conj()函数是求y函数的共轭复数，实数的共轭复数是他本身。

切比雪夫Ⅱ型低通滤波器

一、设计一个切比雪夫Ⅱ型低通滤波器 wp=0.2*pi; %通带边界频率； ws=0.4*pi; %阻带截止频率； rp=1; %通带最大衰减； rs=80; %阻带最小衰减； Fs=1000 %假设抽样脉冲1000hz [N,Wn]=cheb2ord(wp,ws,rp,rs,'s'); %Chebyshev II型滤波器参数计算（模拟域）； [Z,P,K]=cheby2(N,rs,Wn,'s'); %构造Chebyshev II型滤波器（零极点模型）；[H,W]=zp2tf(Z,P,K); %将零极点模型转化成传递函数的模型； figure(1); freqs(H,W); %在Figure1上显示滤波器的幅频响应及相频响应；[P,Q]=freqs(H,W); %返回滤波器的冲击响应的复数形式； figure(2); plot(Q*Fs/(2*pi),abs(P));grid; %在Figure2上显示幅频特性曲线； xlabel('频率/Hz'); ylabel('幅值');

二、设计一个高通Chebyshow型数字滤波器 wp=100;ws=80;Fs=300;rp=1;rs=45; %数字滤波器的各项指标； WP=100*2*pi; %把数字滤波器的频率特征转换成模拟滤波器的频率特征； WS=300*2*pi; [N,Wn]=cheb2ord(WP,WS,rp,rs,'s'); %Chebyshev II型滤波器参数计算（模拟域）； [Z,P,K]=cheb2ap(N,rs); %创建Chebyshev滤波器原型； [A,B,C,D]=zp2ss(Z,P,K); %表达式从零极点增益形式转换成状态方程形式； [AA,BB,CC,DD]=lp2hp(A,B,C,D,Wn); %实现低通到高通滤波器类型的转换；[a,b,c,d]=bilinear(AA,BB,CC,DD,Fs); %采用双线性变换法，从模拟高通到数字高通； [P,Q]=ss2tf(a,b,c,d); %表达式从状态方程形形式转换成传输函数形式；figure(1); freqz(P,Q); %绘出频率响应； [H,W]=freqz(P,Q); figure(2); plot(W*Fs/(2*pi),abs(H));grid; xlabel('频率/Hz'); ylabel('幅值');

基于matlab-的巴特沃斯低通滤波器的实现

基于matlab 的巴特沃斯低通滤波器的实现 一、课程设计的目的 运用MATLAB实现巴特沃斯低通滤波器的设计以及相应结果的显示，另外还对多种低通滤波窗口进行了比较。 二、课程设计的基本要求 1）熟悉和掌握MATLAB 的基本应用技巧。 2）学习和熟悉MATLAB相关函数的调用和应用。 3）学会运用MATLAB实现低通滤波器的设计并进行结果显示。 三、双线性变换实现巴特沃斯低通滤波器的技术指标： 1.采样频率10Hz。 2.通带截止频率fp=0.2*pi Hz。 3.阻带截止频率fs=0.3*pi Hz。 4.通带衰减小于1dB，阻带衰减大于20dB 四、使用双线性变换法由模拟滤波器原型设计数字滤波器 程序代码： T=0.1; FS=1/T; fp=0.2*pi;fs=0.3*pi; wp=fp/FS*2*pi; ws=fs/FS*2*pi; Rp = 1; % 通带衰减 As = 15; % 阻带衰减 OmegaP = (2/T)*tan(wp/2); % 频率预计 OmegaS = (2/T)*tan(ws/2); % 频率预计 %设计巴特沃斯低通滤波器原型

N = ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS))); OmegaC = OmegaP/((10^(Rp/10)-1)^(1/(2*N))); [z,p,k] = buttap(N); %获取零极点参数 p = p * OmegaC ; k = k*OmegaC^N; B = real(poly(z)); b0 = k; cs = k*B; ds = real(poly(p)); [b,a] = bilinear(cs,ds,FS);% 双线性变换 figure(1);% 绘制结果 freqz(b,a,512,FS);%进行滤波验证 figure(2); % 绘制结果 f1=50; f2=250; n=0:63; x=sin(2*pi*f1*n)+sin(2*pi*f2*n); subplot(2,2,1);stem(x,'.'); title ('输入信号'); y=filter(b,a,x); subplot(2,2,2);stem(y,'.') ; title('滤波之后的信号'); figure(3) ; stem(y,'.') title（'输出的信号')）

基于切比雪夫I型的高通滤波器设计Matlab

设计题目基于切比雪夫I型的数字高通滤波器的设计 设计要求 设计一个9阶切比雪夫I型高通滤波器,通带纹波为10dB,下边界频率为400 / rad s ,并绘出其幅频响应曲线 设计过程1．系统设计方案 1.1 Matlab的简介和主要功能： 简介：MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。使用 MATLAB，您可以较使用传统的编程语言（如 C、C++ 和 Fortran）更快地解决技术计算问题。 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。 MATLAB 提供了很多用于记录和分享工作成果的功能。可以将您的 MATLAB 代码与其他语言和应用程序集成，来分发您的 MATLAB 算法和应用。 主要功能:1.此高级语言可用于技术计算 2.此开发环境可对代码、文件和数据进行管理 3.交互式工具可以按迭代的方式探查、设计及求解问题 4.数学函数可用于线性代数、统计、傅立叶分析、筛选、优化以及数值 积分等 5.二维和三维图形函数可用于可视化数据 6.各种工具可用于构建自定义的图形用户界面 7.各种函数可将基于 MATLAB 的算法与外部应用程序和语言（如 C、 C++、Fortran、Java、COM 以及 Microsoft Excel）集成 1.2 开发算法和应用程序 开发算法和应用程序 MATLAB 提供了一种高级语言和开发工具，使您可以迅速地开发并分析算法和应用程序。

基于MATLAB的巴特沃思和切比雪夫低通滤波器的设计

第一章数字滤波器及MATLAB语言概述 数字滤波是语音和图像处理、模式识别、谱分析等应用中的一个基本处理算法，在数字信号处理中占有极其重要的地位。研究基于Matlab 环境下的IIR数字滤波器的设计与实现，给出了相应的Matlab函数命令，并将滤波器应用于图像噪声的去除，取得了不同的效果，就其结果做出了进一步的解释和说明。数字滤波器是具有一定传输选择性的数字信号处理装置，其输入、输出均为数字信号，实质上是一个由有限精度算法实现的线性不变离散系统。它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号输出分量。 MATLAB是美国MathWorks公司推出的一套用于工程计算的可视化高性能语言与软件环境。MATLAB为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。MATLAB推出的工具箱使各个领域的研究人员可以直观方便地进行科学研究、工程应用，其中的信号处理(signalproeessing)、图像处理(imageproeessing)、小波(wavelet)等工具箱为数字滤波研究的蓬勃发展提供了有力的工具。 数字滤波器与模拟滤波器相比，具有精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配以及能实现模拟滤波器无法进行的特殊滤波等优点。本文主要介绍有限冲激响应数字滤波器（FIR）和无限冲激响应数字滤波器（IIR）的设计原理、方法、步骤以及在MATLAB中的实现，并以实例形式列出设计程序和仿真结果。

第二章 方案提取和技术要求 第一节 方案提取 IIR 数字滤波器是一种离散时间系统，其系统函数为 假设M ≤N ，当M ＞N 时,系统函数可以看作一个IIR 的子系统和一个(M-N)的FIR 子系统的级联。IIR 数字滤波器的设计实际上是求解滤波器的系数和 ，它是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S 平面上去逼近，就得到模拟滤波器；如果在z 平面上去逼近，就得到数字滤波器。 设计高通、带通、带阻等数字滤波器通常可以归纳为如图所示的两种常用方法： 方法1 方法2 图一 数字滤波器设计的两种方法

FIR低通滤波器+matlab编程+滤波前后图形

Matlab实现振动信号低通滤波 附件txt中的数字是一个实测振动信号，采样频率为5000Hz，试设计一个长度为M=32的FIR低通滤波器，截止频率为600Hz，用此滤波器对此信号进行滤波。要求： （1）计算数字截止频率； （2）给出滤波器系数； （3）绘出原信号波形； （4）绘出滤波后的信号波形； 解答过程： 第一部分：数字截止频率的计算 =600/5000/2=0.24 数字截止频率等于截止频率除以采样频率的一半，即 n 第二部分：滤波器系数的确定 在matlab中输入如下程序，即可得到滤波器系数： n=32 Wn=0.24 b=fir1(n,Wn) 得到的滤波器系数b为 Columns 1 through 9 -0.0008 -0.0018 -0.0024 -0.0014 0.0021 0.0075 0.0110 0.0077 -0.0054 Columns 10 through 18 -0.0242 -0.0374 -0.0299 0.0087 0.0756 0.1537 0.2166 0.2407 0.2166 Columns 19 through 27 0.1537 0.0756 0.0087 -0.0299 -0.0374 -0.0242 -0.0054 0.0077 0.0110 Columns 28 through 33 0.0075 0.0021 -0.0014 -0.0024 -0.0018 -0.0008 第三部分：原信号波形 将附件4中的dat文件利用识别软件读取其中的数据，共1024个点，存在TXT 文档中，取名bv.txt，并复制到matlab的work文件夹。 在matlab中编写如下程序： x0=load('zhendong.txt'); %找到信号数据地址并加载数据。 t=0:1/5000:1023/5000; %将数据的1024个点对应时间加载

切比雪夫低通滤波器讲解

课程设计 课程名称：数字信号处理 题目编号： 0202 题目名称：切比雪夫Ⅱ型IIR低通滤波器 专业名称：电子信息工程 班级：电子1204班 学号： 20124470411 学生姓名：刘春阳 任课教师：黄国玉 2015年09月30日

课程设计任务书

目录 1. 数字滤波器的设计任务及要求（编号202） (2) 2. 数字滤波器的设计及仿真 (3) 2.1数字滤波器（编号202）的设计 (3) 2.2数字滤波器（编号202）的性能分析 (6) 3. 数字滤波器的实现结构对其性能影响的分析 (7) 3.1数字滤波器的实现结构一（直接型）及其幅频响应 (8) 3.2数字滤波器的实现结构二（级联型）及其幅频响应 (10) 3.3 数字滤波器的实现结构对其性能影响的小结 (10) 4. 数字滤波器的参数字长对其性能影响的分析 (11) 4.1数字滤波器的实现结构一（直接型）参数字长及幅频响应特性变化 (12) 4.2数字滤波器的实现结构二（级联型）参数字长及幅频响应特性变化 (14) 4.3 数字滤波器的参数字长对其性能影响的小结 (16) 5. 结论及体会 (16) 5.1 滤波器设计、分析结论 (16) 5.2 我的体会 (16) 5.3 展望 (17)

1.数字滤波器的设计任务及要求（0202） 每位同学抽签得到一个四位数，由该四位数索引下表一确定待设计数字滤波器的类型及其设计方法， 然后用指定的设计方法完成滤波器设计。 要求： （1）滤波器设计指标：通带截止频 pc ln ()32 d rad i πω= , 过渡带宽度 10 tz () 160 log d rad i πω?≤ ，滚降roll 60dB α=; 其中， id —抽签得到那个四位数（题目编号) （2）滤波器的初始设计通过手工计算完成； （3）在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种 以上合适的滤波器结构进行分析）； （4）在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； （5） 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； （6）课程设计结束时提交设计说明书。 2.数字滤波器的设计及仿真 2.1 数字滤波器（编号202）的设计 随着信息和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中，由于信号中经常混有各种复杂的成分，因此很多信号的处理都是基于滤波器而进行的。所以数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的，这是因为模拟滤波器的设计方法都已发展的相当成熟，且有典型的模拟滤波器供我们选择。如巴特沃斯滤波器、切比雪夫滤波器等。本次课程设计将手工计算一个切比雪夫II 型的IIR 的低通模拟滤波器的系统函数，并在MATLAB 的FDATool 设计工具分析其性能指标。

关于滤波器设计的matlab函数简表

关于滤波器设计、实现的Matlab函数分类函数名功能说明 滤波器分析 （求幅频、相频响应）abs求模值 angle求相角 freqs模拟滤波器的频率响应freqz数字滤波器的频率响应grpdelay群延迟 impz脉冲响应（离散的）zplane画出零极点图 fvtool滤波器可视化工具 滤波器实现（求输入信号通过滤波器的响应）conv/conv2卷积/二维卷积 filter求信号通过滤波器的响应 IIR滤波器阶数估算buttord巴特沃斯滤波器阶数估算 cheb1ord切比雪夫Ⅰ型滤波器阶数估算 cheb2ord切比雪夫Ⅱ型滤波器阶数估算 ellopord椭圆滤波器阶数估算 IIR数字滤波器设计（求系统函数H(z)）butter cheby1 cheby2 ellip 模拟低通滤波器原型（归一化的）buttap模拟低通巴特沃斯滤波器原型cheb1ap模拟低通切比雪夫Ⅰ型滤波器原型cheb2ap模拟低通切比雪夫Ⅱ型滤波器原型ellipap模拟低通椭圆滤波器原型besselap模拟低通贝塞尔滤波器原型 模拟低通滤波器设计（求系统函数H(s)）butter巴特沃斯滤波器设计cheby1切比雪夫Ⅰ型滤波器设计cheby2切比雪夫Ⅱ型滤波器设计ellip椭圆滤波器设计besself贝塞尔滤波器设计 模拟滤波器频带变换lp2bp低通→带通 lp2bs低通→带阻 lp2hp低通→高通 lp2lp低通→低通 滤波器离散化（由模拟滤波器得到数字滤波器）bilinear脉冲响应不变法impinvar双线性变换法 FIR滤波器设计fir1基于窗函数的FIR滤波器设计 fir2基于窗函数的任意响应FIR滤波器设 计 窗函数boxcar矩形窗 rectwin矩形窗 bartlett三角窗

切比雪夫低通滤波器设计

摘要 随着信息与数字时代的到来,数字信号处理已成为当今一门极其重要的学科与技术领域。在现代通信系统中,由于信号中经常混有各种复杂成分,因此很多信号的处理都就是基于滤波器而进行的。所以,数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的,这就是因为模拟滤波器的理论与设计方方法都已发展的相当成熟, 且有典型的模拟滤波器供我们选择。,如巴特沃思滤波器、切比雪夫滤波器等。 本次课程设计将运用MATLAB设计一个基于切比雪夫低通滤波器,并出所设计滤波器的幅度及幅度衰减特性。 关键词:模拟低通滤波切比雪夫

1课题描述 数字滤波器就是数字信号处理的重要工具之一,它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序,而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。故本课题使用MATLAB 信号处理箱与运用切比雪夫法设计数字低通滤波器。 2设计原理 2、1切比雪夫滤波器介绍 在巴特沃兹滤波器中,幅度响应在通带与阻带内都就是单调的。因此,若滤波器的技术要求就是用最大通带与阻带的逼近误差来给出的话,那么,在靠近通带低频端与阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径就是使逼近误差均匀地分布于通带或阻带内,或同时在通带与阻带内都均匀分布,这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不就是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质:其频率响应的幅度既可以在通带中就是等波纹的,而在阻带中就是单调的(称为I 型切比雪夫滤波器),也可以在通带中就是单调的,而在阻带中就是等波纹的(称为II 型切比雪夫滤波器)。I 型切比雪夫滤波器的幅度平方函数就是 2|)(|Ωj H C =)/(11 22c N C ΩΩ+ε (2、1) 式中为N 阶切比雪夫多项式,定义为 )cos cos()(1x N x C N -= (2、2) 从定义切比雪夫多项式可以直接得出由 )(x C N 与)(1x C N -求)(1x C N +的递推公式。将三角恒等式代入 (2、2)式,得 )(1x C N +=2x -)(x C N ) (1x C N - (2、3) 从 (2、2)式我们注意到,当01

基于MATLAB的巴特沃斯滤波器

数字信号处理课程设计 2015年 6 月25 日

目录 一．设计目的： (3) 二．设计要求： (3) 三．设计内容： (4) 3.1选择巴特涡斯低通数据滤波器及双线性变换法的原因 (4) 3.2巴特沃思低通滤波器的基本原理 (4) 3.3双线性变换法原理 (5) 3.4数字滤波器设计流程图 (7) 3.5数字滤波器的设计步骤 (7) 四．用matlab实现巴特沃斯低通数字滤波器的仿真并分析 (9) 4.1巴特沃斯低通数字滤波器技术指标的设置 (9) 4.2用matlab实现巴特沃斯低通数字滤波器的仿真 (9) 4.3波形图分析： (12) 五.总结与体会 (13) 六．附录参考文献 (14) 2

一．设计目的： 该课程设计是测控技术与仪器专业的必修课，开设课程设计的目的使学生掌握数字信号处理的基本概念和基本理论，能够利用辅助工具进行FIR和IIR数字滤波器的设计，进行一维信号的频谱分析，并进行仿真验证。加强实践教学环节，加强学生独立分析、解决问题的能力，培养学生动手能力和解决实际问题的能力，实现宽口径教育。 (1)理解低通滤波器的过滤方法。 (2)进一步熟悉低通滤波器的基本应用。 (3)用仿真工具matlab软件对设计的滤波器进行软件和硬件仿真。 (6)将对仿真结果进行比较，从而检验滤波器滤波性能的准确性。 二．设计要求： 地震发生时，除了会产生地震波，还会由地层岩石在断裂、碰撞过程中所发生的震动产生次声波。它的频率大约在每秒十赫兹到二十赫兹之间（可以用11Hz和15Hz的两个信号的和进行仿真，幅度可以分别设定为1、2）。大气对次声波的吸收系数很小，因此它可以传播的很远，而且穿透性很强。通过监测次声波信号可以监测地震的发生、强度等信息，因为自然界中广泛存在着各种次声波，这就对地震产生的次声波产生了干扰（可以用白噪声模拟，方差为5），需要采取一定的处理方法，才能检测到该信号，要求设计检测方案；并处理方法给出具体的软件（可以以51系列单片机、STM32F407、TMS320F28335或TMS320F6745为例）。 假设地震次声波信号为x,输入x=sin(2*π*11*t)+2*sin(2*π*15*t)和伴有白噪声的合成信号，经过滤波器后滤除15Hz以上的分量，即只保留x=sin(2*π*11*t)+2*sin(2*π*15*t)的分量信号，来验证设计的滤波器是否达到了设计要求。 3

巴特沃斯、切比雪夫滤波器的仿真与实现

华北科技学院课程设计任务书 2013 — 2014 学年第二学期 电子信息工程学院（系、部）通信工程专业 B111 班级课程名称：移动通信 设计题目：巴特沃斯、切比雪夫滤波器的仿真与实现完成期限：自16 周至 18 周共 3 周

目录 1．前言 (3) 1.1 MATLAB (3) 1.2 滤波器的概念 (5) 1.2.1滤波器的原理 (6) 1.2.2理想滤波器与实际滤波器 (6) 1.2.3 滤波器的分类 (7) 2.设计目的 (9) 3.设计原理 (9) 3.1.模拟低通滤波器的设计指标及逼近方法 (9) 3.2.巴特沃斯低通滤波器的设计方法 (10) 3.3.切比雪夫滤波器的设计方法 (14) 4.详细设计与系统分析 (21) 4.1程序设计 (21) 4.1.1巴特沃斯滤波器 (21) 4.1.2切比雪肤滤波器 (23) 4.2同一滤波器不同参数的比较 (25) 4.2.1巴特沃斯滤波器 (25) 4.2.2切比雪夫滤波器 (27) 4.3不同滤波器同一阶数的比较 (30) 4.3.1低通滤波器 (30) 4.3.2高通滤波器 (30) 4.3.3带通滤波器 (31) 4.3.4带阻滤波器 (31) 5.心得体会 (32) 6.参考文献 (32)

摘要：利用MATLAB设计滤波器,可以按照设计要求非常方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计的最优化。MATLAB因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。本文介绍了在MATLAB R2009a 环境下滤波器设计的方法和步骤。关键词：滤波器，MATLAB 1．前言 1.1 MATLAB MATLAB是美国MathWorks公司开发的一种功能极其强大的高技术计算语言和内容极其丰富的软件库，集数值计算、矩阵运算和信号处理与显示于一身。该软件最初是由美国教授Cleve Moler 创立的。1980年前后，他在教线性代数课程时，发现用其他高级语言编程时极不方便，便构思开发了MATLAB，即矩阵实验室（Matrix Laboratory）。该软件利用了当时代表数值线性代数领域最高水平的EISPACK和LINPACK两大软件包，并且利用Fortran 语言编写了最初的一套交互式软件系统，MATLAB的最初版本便由此产生了。 最初的MATLAB由于语言单一，只能进行矩阵的运算，绘图也只能用原始的描点法，内部函数只有几十个，因此功能十分简单。1984年该公司推出了第一个MATLAB的商业版，并用C语言作出了全部改写。现在的MATLAB程序是MathWorks公司用C语言开发的，第一版由steve Bangert主持开发编译解释程序，Steve Kleiman完

高通滤波器(matlab编写)切比雪夫

设计一个高通滤波器，并检验它的性能 采样率为10kHZ 阻带边缘为1.5Khz,衰减为40bB 通带边缘为2kHz，波纹为3Db 采用切比雪夫2型滤波器 Fs=1e4 ; fs=1.5e3; fp=2e3; As=40; Rp=3 ;wp=2*fp/Fs; ws=2*fs/Fs; 用MA TLAB工具： 设计的滤波器性能 >>Fs=1e4; >>fs=1.5e3; >>fp=2e3; >>As=40; >>Rp=3; >>wp=2*fp/Fs; >>ws=2*fs/Fs; >>[N,wn]=cheb2ord(wp,ws,Rp,As); >>[b,a]=cheby2(N,As,wn,'high'); >>[db,mag,pha,grd,w]=freqz_m(b,a); >> subplot(2,2,1);plot(w/pi,mag); >> axis([0,1,0,1]); >> setX([0 0.3 0.4 1]); >>setY([0.01 0.7279 1]) >> title('Magnitude Response'); >> subplot(2,2,2);plot(w/pi,db); >> axis([0 1 -70 0]) >> setX([0 0.3 0.4 1]) >> setY([-40 -2.7589]) >> title('Magnitude Response in dB'); 输入信号x=cos(0.6*pi*n); 取200个抽样值，为了便于观看，在画图是裁剪掉100个值

不直观，直接做频谱图 为了说明该滤波器的性能，做一组对比。换上另一种频率的信号x=cos(0.2*pi*n); 由上图可以看到该信号完全被滤去了