实验报告。
2016-2017学年第一学期。
课程名称专业综合实验。
--iir滤波器设计。
院（系） 电子信息学院。
专业电子信息工程班级电子1313
学生姓名周立学号 131003430337
实验时间 2017.1.2- 2017.1.13
指导老师胡秀娟。
提交时间2017.1.13
目录。一课题名称 1
二设计任务和要求 1
三实验内容 1
3.1数字滤波器的简介 1
3.2 椭圆数字带通滤波器设计原理 1
3.3实验步骤 3
3.4 程序** 3
四方案的原理框图 4
五电路调试或算法调试 5
5.1实验**结果图 5
5.2实验结果分析 6
六收获与体会 6
参考文献 6
iir滤波器设计。
用双线性变换法设计一个数字滤波器，使其指标接近如下技术指标的模拟带通椭圆滤波器：
wp1=100hz,wp2=200hz,ws1=50hz,ws2=250hz,通带衰减系数为rp=0.5db,阻尼衰减系数rs=60db,采样频率为fs=2000hz。
数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字**进行运算处理，以达到改变信号频谱的目的。
数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（iir）滤波器和有限长冲激响应（fir）滤波器。
iir数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以iir滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。fir数字滤波器的单位脉冲响应是有限长序列。
它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。
3.2 椭圆数字带通滤波器设计原理。
双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，可以适用于高通、带通和带阻滤波器的设计。iir数字滤波器的设计通常要借助于模拟低通滤波器的设计，由原型低通滤波器到其他形式（高通、带通、带阻）iir数字滤波器的频带变换有模拟频带变换法和数字频带变换法。本实验中用的是模拟频带变换法。
首先将给定的对数字滤波器(df)的技术要求转换为一个低通模拟滤波器(af)的技术要求，根据这种要求用某种逼近设计出原型的低通模拟滤波器(lp af)，计算出模拟滤波器的阶数n、极点si和传递函数，再按照双线性变换的变换关系，将模拟滤波器的传递函数转换为数字滤波器的传递函数。表1给出了对数字滤波器(df)的技术要求直接转换为对一个低通模拟滤波器(af)的技术要求的频率预畸变校正关系和转换公式。
图表 1双线性变换公式。
低通到带通的频率变换公式如下：
在p平面与s平面虚轴上的频率关系为：
式中，，表示带通滤波器的带通宽度，分别为带通滤波器的带通下截止频率和通带截止频率，根据公式（2）的映射关系，频率λ=0映射为频率，频率λ=映射为频率。也就是说将低通滤波器g(p)的通带映射为带通带通滤波器的通带，同理，得出将q（p）转换为带通滤波器的系统函数，即：
可以证明 (4)
所以，带通滤波器的通带（阻带）边界频率关于中心频率几何对称，如果原指标给定的边界频率不满足公式（4），就要改变其中一个边界频率，以便满足公式（4），但要保证改变后的指标高于原指标，具体方法是，如果，则减小以使（4）得到满足，具体计算公式为：
或5)减小使带通宽度大于原指标要求的通带宽度，增大或减少都使左边的过渡带宽度小于原指标要求的过渡带宽度，反之，如果，，则减少使式（4）得到满足，而且在关于中心频率几何对称的两个正频率点上，带通滤波器的幅度值相等。
（1）：根据题目要求计算出椭圆模拟带通滤波器的通带频率，阻带频率。
2）：根据参数设计出椭圆低通滤波器。
3）：将椭圆低通滤波器转化为带通滤波器。
4）：归一化模拟带通滤波器的参数，运用双线性变换法计算出数字带通滤波器的参数。
5）：根据参数设计出数字带通滤波器。
6）：使用matlab画出滤波器的频谱分析图和零极点分布图。
设计完成后对已设计的滤波器的频率响应要进行校核，要得到幅频相频响应特性，运算量也是很大的。我们平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候都是要根据设计要求和滤波效果不断的调整，以达到设计的最优化。在这种情况下，滤波器的设计就要进行大量复杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成设计。
利用matlab强大的计算功能进行计算机辅助设计，可以快速有效的设计数字滤波器，大大的简化了计算量，直观简便。
design a ellip digtal bandpass filter
wp1=100;wp2=200;ws1=50;ws2=250;fs=1000*2*pi;rp=0.5;rs=50;
wp1=2*pi*wp1;wp2=2*pi*wp2;ws1=2*pi*ws1;ws2=2*pi*ws2;
wp=[wp1 wp2];ws=[ws1 ws2];nn=128;
n wn]=ellipord(wp/(fs/2),ws/(fs/2),rp,rs,'s');估计滤波器最小阶数。
bw=wp2-wp1;
wo=sqrt(wp2*wp1);
z,p,k]=ellipap(n,rp,rs);
a,b,c,d]=zp2ss(z,p,k);
at,bt,ct,dt]=lp2bp(a,b,c,d,wo,bw);
at1,bt1,ct1,dt1]=bilinear(at,bt,ct,dt,fs);
bz,az]=ss2tf(at1,bt1,ct1,dt1);
h,w]=freqz(bz,az);
semilogy(w*fs/2/pi,abs(h)),grid %画对数幅频谱。nwn
稳定的模拟滤波器经双线性变换后得到的数字滤波器也一定是稳定的。但是它的频率变换是非线性的。这种非线性关系使得通带截止频率，过渡带的边缘频率的相对位置发生了非线性畸变。
这种非线性畸变可以通过预畸变来校正。
用双线性变换来设计数字滤波器时，一般总是先将数字滤波器的各临界频率来预畸变，求得相应的参考模拟滤波器的各临界频率，然后设计参考模拟滤波器的传递函数，最后通过双线性变换公式求得数字滤波器的传递函数。这样通过双线性变换，正好将这些频率截止点映射到我们所需的位置上。完成用双线性变换法滤波器的设计。
采用matlab设计滤波器，使原来非常繁琐复杂的程序设计变成了简单的函数调用，为滤波器的设和实现开辟了广阔的天地，尤其是matlab工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。其中的信号处理工具箱、图像处理工具箱、小波工具箱等更是为数字滤波研究的蓬勃发展提供了可能。
在设计数字滤波器时， 善于应用matlab进行辅助设计， 能够大大提高设计效率。数字滤波器的应用十分广泛，利用matlab语言，很容易设计出数字椭圆滤波器，数字椭圆滤波器可以用较少的阶数获得很高的选择特性，在相位要求不明感的场合，如语音信号处理等，适合用数字椭圆数字滤波器，在设计过程中可以对比滤波器特性，随时更改参数，以达到滤波器设计的最优化。
参考文献。1] 高西全。丁玉美。 数字信号处理。 西安：西安电子科技大学出版社，2008.
2] 管爱红，张红梅，杨铁军等。 matlab基础及其应用教程。 北京：电子工业出版社，2009.
3] 董长虹等。 matlab信号处理与应用。 北京：国防工业出版社， 2005.
4] 楼顺天，李博菡。 基于matlab的系统分析与设计－信号处理。 西安：西安电子科技大学出版社，2005.
5] 蒋志凯。 数字滤波与卡尔曼滤波。 北京：中国科学技术出版社，1993.
6] [美] 海因斯著，张建华等译。 数字信号处理。 北京：科学出版社，2002.
评语及成绩评定记录。
《专业 综合 实训》实训大纲
专业 综合 实训 实训教学大纲。一 实训教学课程的性质与任务。本课程是艺术设计专业的综合实训课程。通过实训，学生完成一份自己设计作品与专业技能展示的个人专业技能宣传册，设计一个让学生宣传与推广自我专业技能的环境，让学生通过动手整理 策划 包装 设计 制作等手段，接受广告设计及策划推广职业技能的训练，...
专业综合实训
实训总结性说明。本系统需要实时采集温度数据，然后经过a d 转换为数字信号，存入at89s52 的内部数据存储器，送液晶显示，并与设定值进行比较，经过pid 算法得到控制量并由单片机输出去控制加热器。进行温度控制程序的设计应考虑如下问题 温度采样，数字滤波 pid 控制程序 按键处理，液晶显示 上位...
《专业CAD》实训周实训计划
专业cad 实训周实训计划。适用班级 建筑工程学院10光伏材料班。一 实训的性质和目的。专业cad 课程实训是学生学习课程 专业cad 之后重要的实践环节，是学生在校学习期间理论联系实际 增强实际操作技能的重要手段和方法之一。通过实训，进一步巩固课堂所学的基础理论知识，提高学生的读图和绘图能力。二 ...