设计数字低通滤波器(用matlab实现)资料
更新时间:2023-12-04 02:10:01 阅读量: 教育文库 文档下载
DSP设计滤波器报告
姓名:张胜男 班级:07级电信(1)班 学号:078319120
一·低通滤波器的设计
(一)实验目的:掌握IIR数字低通滤波器的设计方法。
(二)实验原理: 1、滤波器的分类
滤波器分两大类:经典滤波器和现代滤波器。
经典滤波器是假定输入信号x(n)中的有用成分和希望取出的成分各自占有不同的频带。这样,当x(n)通过一个线性系统(即滤波器)后可讲欲去除的成分有效的去除。
现代滤波器理论研究的主要内容是从含有噪声的数据记录(又称时间序列)中估计出信号的某些特征或信号本身。
经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器(AF)和数字滤波器(DF)。对数字滤波器,又有IIR滤波器和FIR滤波器。 IIR DF 的转移函数是:
Y(z)?X(z)H(z)??bzrr?0Nk?1M?r1??akz?k
FIR DF的转移函数是:
H(z)??h(n)z?nn?0N?1
FIR滤波器可以对给定的频率特性直接进行设计,而IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。 2、滤波器的技术要求 低通滤波器:
?p:通带截止频率(又称通带上限频率) ?s:阻带下限截止频率
?p?:通带允许的最大衰减 ?s:阻带允许的最小衰减 (p,?s:通带上限角频率 ?s:阻带下限角频率
的单位dB)
?p(
?p??pTs,
?s??sTs)即
?p?2??pFC
?s?2??sFC
3、IIR 数字滤波器的设计步骤:
1)按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技
术指标。 2)根据转换后的技术指标设计模拟低通滤波器G(s); 3)再按一定的规则将G(s)转换成H(z)。
4)若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模
拟滤波器的技术指标,然后按上述步骤2)设计出低通G(s),再将G(s)转换为所需的H(z)。
4.几种不同类型的滤波器的介绍:
因为我们设计的滤波器的冲击响应一般都为实数,所以有
G(s)G*(s)?G(s)G(?s)s?j??G(j?s)2
2这样,如果我们能由的G(s)。 不同类型的
G(j?)2?p,
?pG(j?),?s,?s求出,那么就容易得到所需要
的表达式,代表了几种不同类型的滤波器。
(1)巴特沃思(Butterworth)滤波器:
G(j?)?2
11?C2(?2)n
C为待定常数,N为待定的滤波器阶次。 (2)切比雪夫I型(Chebyshev – I )滤波器:
G(j?)?2121??2Cn(?)
5.巴特沃思模拟低通滤波器的设计
由于每一个滤波器的频率范围将直接取决于设计者的应用目的,因此必然是千差万别。为了使设计规范化,我们需要将滤波器的频率参数作归一化处理。设所给的实际频率为?(或f),归一化后的频率为?,对低通模拟滤波器,令
,
?=?/?p
。又令归一化复数变量为p,p?j?,显然
显然,
?p?1
?s??s/?p
p?j??j?/?p?s/?p所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。 (1)将实际频率?规一化
(2)求C和N
C2?10?P/10?1
N?lg10?s/10?110?p/10?1lg?s
这样C和N可求。 若令
?p=3dB,则C=1,这样巴特沃思滤波器的设计就只剩一个参数N,这时
G(j?)?211??2N?11?(?/?p)2N
(3)确定G(s)
因为p?j?,根据上面公式有
G(p)G(?p)?11?1?(p/j)2N1?(?1)Np2N
由
1?(?1)Np2N?0解得
2k?N?12?)2N,k=1,2,···,2N
12k?N?1p2?2pcos(?)?12N
pk?exp(j这样可得
Gk(p)?1?(p?pk)(p?pN?1?k)求得G(p)后,用
s/?p代替变量p,即得实际需要得G(s)。
6.用双线性Z变换法设计IIR数字低通滤波器 s平面到z平面的映射关系
s?2z?1Tsz?1
称为双线性Z变换,由此关系求出
z?1?(Ts/2)s1?(Ts/2)s
j??j2sin(?/2)Tscos(?/2)及
??即
2tan(?/2)Ts
??2arctan(?Ts/2)
2因为设计滤波器时系数Ts会被约掉,所以又有
z?1?s1?s
??tan(?/2)
??2arctan?
(三)实验内容
题一:试用双线性Z变换法设计一低通数字滤波器,给定技术指标是
fp?100Hz,
fs?300Hz,?p?3dB,?s?20 dB,抽样频率Fs?1000Hz。 提示:首先应该得到角频率?,然后再
(1) 将数字滤波器的技术要求转换为模拟滤波器的技术要求。 (2) 设计低通滤波器G(s) 由?=
?/?p依次求出
?p,?s,再求出N,可得G(p)
G(s)?G(p)然后由
p?s?p转换成G(s)
(3) 由G(s)求H(z)
源程序:
fp=100; fs=300;
ap=3; %通带最大衰减 as=20; %阻带最小衰减 Fs1=1000; %抽样频率 wp=2*pi*fp/Fs1;
ws=2*pi*fs/Fs1; %数字角频率 Fs=Fs1/Fs1; T=1/Fs;
Op=2/T*tan(wp/2);
Os=2/T*tan(ws/2); %角频率畸变
[N Wn]=buttord(Op,Os,ap,as,'s'); %N代表滤波器阶数,Wn代表滤波器的
截止频率。简单来说就是在wp处,通带内波
纹系数或者说是通带内达到最大衰减为Rp,如
(3db),而在Ws处,阻带达到最小衰减为Rs(如40db),而我们默认求得的Wn是在(-3db)时的频率
[z p k]=buttap(N); %G(p)的零点,极点,增益用于计算N阶归一化
(3dB截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子。
[b a]=zp2tf(z,p,k); %G(p)的分子、分数系数从零、极点模型得到
系统函数的分子、分母多项式系数向量ba、aa
[B A]=lp2lp(b,a,Op); %H(s)的分子、分数系数改变低通模拟滤波器的
截止频率,原滤波器是以多项式系数Bap,Aap给出的,改后的滤波器是带截止频率Wn的
[Bz Az]=bilinear(B,A,Fs*T/2); %H(z)的分子、分数系数 [H w]=freqz(Bz,Az,256,Fs*Fs1); %H频率响应 plot(w,abs(H),'r'); title('低通滤波器'); grid on;
实验结果:
低通滤波器1.210.80.60.40.20050100150200250300350400450500
同理可以设计出低通滤波器和带通及带阻滤波器。 MATLAB参考程序和仿真内容
%*******************************************************************
正在阅读:
2014.5.28 讲课 护理质量分析12-19
生活大爆炸 经典台词04-22
一、基本公共卫生服务健康教育服务规范解读08-16
解放军陆军兵种师旅编制汇总04-06
财务管理题库05-03
实验8(窗体控件的使用)10-13
美国在辽宁省投资企业名录B05-12
我迷上了玩电脑作文400字07-08
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 滤波器
- 实现
- 数字
- matlab
- 设计
- 资料
- 模拟测试
- 江苏省公共建筑施工图绿色设计专篇(电气)
- ZHFL-008工程工作手册 - 图文
- 2018年1月北京石景山区初三地理期末试卷 - 图文
- 线性代数第五版第一章常见试题及解答
- 机关食堂管理(四办法)
- 财务管理计算题试题和答案
- 我国宽带钢热连轧机的最新发展及其评析 - 图文
- 2012届中考数学四边形专项复习训练题
- 英语词汇学 张维友版经典笔记
- 激光扩散片超声波清洗机
- 2006年高考理综试题(全国卷1)(湖北卷)
- 水湘圩地区成矿规律研究及成矿预测
- 电力系统继电保护与自动化专业毕业实习报告范文
- 国际税收教案
- 衰退与危机下丨被栽赃的劳动者、逃避责任的权贵
- 黑山羊 - 常见疾病的防治
- 北师大版七年级数学下第一章整式的乘除复习课教案
- 尿疗手册
- 电子商务案例分析答案