已阅读5页,还剩2页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
DSP设计滤波器报告姓名:张胜男 班级:07级电信(1)班 学号:一低通滤波器的设计(一)实验目的:掌握IIR数字低通滤波器的设计方法。(二)实验原理:1、滤波器的分类滤波器分两大类:经典滤波器和现代滤波器。经典滤波器是假定输入信号中的有用成分和希望取出的成分各自占有不同的频带。这样,当通过一个线性系统(即滤波器)后可讲欲去除的成分有效的去除。现代滤波器理论研究的主要内容是从含有噪声的数据记录(又称时间序列)中估计出信号的某些特征或信号本身。经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器(AF)和数字滤波器(DF)。对数字滤波器,又有IIR滤波器和FIR滤波器。IIR DF 的转移函数是:FIR DF的转移函数是:FIR滤波器可以对给定的频率特性直接进行设计,而IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。2、滤波器的技术要求低通滤波器:通带截止频率(又称通带上限频率):阻带下限截止频率:通带允许的最大衰减 :阻带允许的最小衰减(,的单位dB):通带上限角频率:阻带下限角频率(,)即 3、IIR 数字滤波器的设计步骤:1) 按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。2) 根据转换后的技术指标设计模拟低通滤波器;3) 再按一定的规则将转换成。4) 若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标,然后按上述步骤2)设计出低通,再将转换为所需的。4几种不同类型的滤波器的介绍:因为我们设计的滤波器的冲击响应一般都为实数,所以有这样,如果我们能由,求出,那么就容易得到所需要的。不同类型的的表达式,代表了几种不同类型的滤波器。(1)巴特沃思(Butterworth)滤波器:C为待定常数,N为待定的滤波器阶次。(2)切比雪夫I型(Chebyshev I )滤波器:5巴特沃思模拟低通滤波器的设计由于每一个滤波器的频率范围将直接取决于设计者的应用目的,因此必然是千差万别。为了使设计规范化,我们需要将滤波器的频率参数作归一化处理。设所给的实际频率为(或),归一化后的频率为,对低通模拟滤波器,令显然,。又令归一化复数变量为,显然所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。(1)将实际频率规一化(2)求C和N这样C和N可求。若令3dB,则C1,这样巴特沃思滤波器的设计就只剩一个参数N,这时(3)确定因为,根据上面公式有由解得,1,2,2N这样可得求得后,用代替变量,即得实际需要得。6用双线性Z变换法设计IIR数字低通滤波器s平面到z平面的映射关系称为双线性Z变换,由此关系求出及即因为设计滤波器时系数会被约掉,所以又有(三)实验内容题一:试用双线性Z变换法设计一低通数字滤波器,给定技术指标是Hz,Hz,dB, dB,抽样频率Hz。提示:首先应该得到角频率,然后再(1) 将数字滤波器的技术要求转换为模拟滤波器的技术要求。(2) 设计低通滤波器由依次求出,再求出N,可得然后由转换成(3) 由求源程序: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、aaB 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;实验结果:同理可以设计出低通滤波器和带通及带阻滤波器。MATLAB参考程序和仿真内容%*%functionx=ffts(mode,M)Nfft=2M;x=zeros(1,Nfft); %定义一个长度为Nfft的一维全0数组if mode= =1 for n=0:Nfft-1 x(n+1)=sin(2*pi*n/Nfft); end end %定义一个长度为Nfft的单周期正弦序列if mode= =2 for n=0:Nfft-1 x(n+1)=sin(4*pi*n/Nfft); end end %定义一个长度为Nfft的双周期正弦序列if mode= =3 for n=0:Nfft/2-1 x(n+1)=sin(4*pi*n/Nfft); end end %定义一个长度为Nfft/2的正弦序列,后面一半为0序列。if mode= =4 for n=0:Nfft-1 x(n+1)=square(2*pi*n/Nfft); end endif mode= =5 for n=0:Nfft-1 x(n+1)=square(2*pi*n/Nfft); end endif mode= =6 for n=0:Nfft/2-1 x(n+1)=square(4*pi*n/Nfft); end endn=0:Nfft-1;subplot(2,1,1);stem(n,x);axis(0 Nfft-1 1.1*min(x) 1.1*max(x);xlabel(Point
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 化学人教版九上课后提分练:4.1-4.2 水资源的保护与净化
- 2024年河北省秦皇岛市青龙县金声木铎学校中考数学一模前模拟试卷
- 磨料磨具制造项目可行性研究报告
- 2024-2030全球与中国埋地接线盒市场现状及未来发展趋势
- 山西省2024九年级物理全册第十八章电功率学情评估新版新人教版
- 2024年东营市《高等数学(一)》(专升本)深度自测卷含解析
- 2022年房产估价师考试《相关知识》第七章金融知识习题
- 爬上豆蔓看自己-《爬上豆蔓看自己》的读后感范文
- 士兵突击观后感5篇
- 2024年微波暗室设备项目合作计划书
- 2023学年完整公开课版将军饮马
- 2023年广西壮族崇左市宁明县六年级数学第二学期期末考试模拟试题含解析
- 2023年上海杨浦区初三二模语文试卷(附答案)
- 2022第三届“红旗杯”班组长大赛考试题库大全-上(单选题库2-2)
- 电动游览车设计-毕业论文
- 广州星级酒店名录
- 皮带输送机日常点检表
- 失智老人活动策划与组织-蒙氏四大领域活动
- 花样跳绳校本课程教材内容
- 能源管理体系内审检查表及能源管理体系内部审核检查表
- 非遗文创项目设计与推广合作协议合同模板
评论
0/150
提交评论