调幅发射机的课程设计(共24页).doc





《调幅发射机的课程设计(共24页).doc》由会员分享,可在线阅读,更多相关《调幅发射机的课程设计(共24页).doc(24页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上目录专心-专注-专业第一章 数字滤波器设计的具体要求和目的第一节 摘要数字滤波技术是数字信号处理的一个重要组成部分,滤波器的设计是信号处理的核心问题之一。介绍了用双线性不变法设计IIR数字滤波器的基本流程,以双线性不变法设计了一个高通巴特沃斯IIR数字滤波器。根据FIR滤波器的原理,提出了FIR滤波器的窗函数设计法,并对常用的几种窗函数进行了比较。给出了在MATLAB环境下,用窗函数法设计FIR滤波器的过程和设计实例。然后在MATLAB环境下进行了仿真与调试,实现了设计目标。第二节 关键字双线性变换法、巴特沃斯、凯泽窗、数字高通滤波器、MATLAB第三节 设计目的通过
2、本课程设计,巩固数字信号处理的有关知识,熟悉各种数字滤波器的设计原理,设计方法。学会分析有限字长效应带来的误差,合理设计数字滤波器的结构。学会利用MATLAB来进行滤波器的设计。第四节 设计要求给定技术指标:通带截止频率,通带最大衰减阻带截止频率 ,阻带最小衰减,抽样频率用双线性变换法设计一个巴特沃斯数字高通滤波器。技术指标同上,用凯泽窗法设计一个FIR数字高通滤波器。设计时间为一周。第二章 对数字信号处理的认识和了解 随着信息时代的到来,数字信号处理已经成为当今一门极其重要的学科和技术,并且在通信、语音、图像、自动控制等众多领域得到了广泛的应用。在数字信号处理中,数字滤波器占有极其重要的地位
3、,它具有精度高、可靠性好、灵活性大等特点。现代数字滤波器可以用软件或硬件两种方式来实现。软件方式实现的优点是可以通过滤波器参数的改变去调整滤波器的性能。MATLAB是一种面向科学和工程计算的语言,它集数值分析、矩阵运算、信号处理和图形显示于一体,具有编程效率高、调试手段丰富、扩充能力强等特点。MATLAB的信号处理工具箱具有强大的函数功能,它不仅可以用来设计数字滤波器,还可以使设计达到最忧化,是数字滤波器设计的强有力工具。第三章 数字滤波器的设计原理第一节 IIR数字滤波器的设计一 IIR数字滤波器简介一个数字滤波器可以用系统函数表示为: 由这样的系统函数可以得到表示系统输入与输出关系的常系数
4、线形差分方程为: 可见数字滤波器的功能就是把输入序列x(n)通过一定的运算变换成输出序列y(n)。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h(n)是无限长的,其差分方程如上式所示,是递归式的,即结构上存在着输出信号到输入信号的反馈,其系统函数具有上式的形式,因此在z平面的有限区间(0z)有极点存在。前面已经说明,对于一个给定的线形时不变系统的系统函数,有着各种不同的等效差分方程或网络结构。由于乘法是一种耗时运算,而每个延迟单元都要有一个存储寄存器,因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择,以便提高运算速度和减少存储器。然而,当需要考虑有限
5、寄存器长度的影响时,往往也采用并非最少乘法器和延迟单元的结构。IIR滤波器实现的基本结构有:(一)IIR滤波器的直接型结构;优点:延迟线减少一半,变为N 个,可节省寄存器或存储单元;缺点:其它缺点同直接I型。通常在实际中很少采用上述两种结构实现高阶系统,而是把高阶变成一系列不同组合的低阶系统(一、二阶)来实现。图一 直接型(二)IIR滤波器的级联型结构:特点:系统实现简单,只需一个二阶节系统通过改变输入系数即可完成;极点位置可单独调整;运算速度快(可并行进行);各二阶网络的误差互不影响,总的误差小,对字长要求低。 缺点:不能直接调整零点,因多个二阶节的零点并不是整个系统函数的零点,当需要准确的
6、传输零点时,级联型最合适。图二 串联型(三)IIR滤波器的并联型结构。优点:简化实现,用一个二阶节,通过变换系数就可实现整个系统;极、零点可单独控制、调整,调整1i、2i只单独调整了第i对零点,调整1i、2i则单独调整了第i对极点;各二阶节零、极点的搭配可互换位置,优化组合以减小运算误差;可流水线操作。缺点:二阶阶电平难控制,电平大易导致溢出,电平小则使信噪比减小。 图三 并联型 二 双线性变换法设计IIR滤波器原理将S平面j轴压缩变换到s1平面j轴上的-/T到/T一段,可以采用以下变换关系: 这样变换到,变换到,可将上式写成令,解析延拓到整个s平面和s1平面,可得再将平面通过以下变换关系映射
7、到z平面,即从而得到s平面和z平面的单值映射关系为一般来说,为了使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系,可引入待定常数c,将代入到上式,可得在MATLAB中,双线性Z变换可以通过bilinear函数实现,其调用格式为:Bz,Azbilinear(B,A,Fs);其中B,A为模拟滤波器传递函数G(s)的分子分母多项式的系数向量,而Bz,Az为数字滤波器的传递函数H(z)的分子分母多项式的系数向量。三 方案的设计论证与选择本次滤波器设计的主要方法是先设计低通模拟滤波器,然后转换为高通数字滤波器。对于其他如高通,则通过频率变换转换为设计相应的高通等。在设计的全过程的各个步骤,matl
8、ab都提供相应的工具箱函数,使得IIR数字滤波器设计变得非常简单。总的来说,我的设计思路主要有以下两种:方案一:从归一化模拟低通原型出发,先在模拟域内经频率变换成为所需类型的模拟滤波器;然后进行双线性变换,由S域变换到Z域,而得到所需类型的数字滤波器。数字高通模拟高通归一化模拟低通原型 模拟域 双线性变换法频率变换 图四 先频率变换再离散方案二:先进行双线性变换,将模拟低通原型滤波器变换成数字低通滤波器;然后在Z域内经数字频率变换为所需类型的数字滤波器。 归一化模拟低通原型数字原型低通数字高通 双线性 数字域 变换法 频率变换 图五 先离散再频率变换以上两种方案都可以,我最后选择了第一种方案进
9、行设计,即先在模拟域内经频率变换成为所需类型的模拟滤波器;然后进行双线性变换,由S域变换到Z域,而得到所需类型的数字滤波器。因为高通滤波器的设计方法是先将要设计的滤波器的技术特性指标通过频率转换成模拟低通滤波器的技术指标,再根据这些性能指标设计出低通滤波器的传递函数,为了设计程序的简洁,故选择了方案一。第二节 FIR滤波器的设计一 FIR滤波器的简介FIR滤波器的单位抽样响应为有限长度,一般采用非递归形式实现。通常的FIR数字滤波器有横截性和级联型两种。FIR滤波器实现的基本结构有:(一)FIR滤波器的横截型结构表示系统输入输出关系的差分方程可写作: 直接由差分方程得出的实现结构如图六所示:图
10、六 横截型(直接型卷积型)若h(n)呈现对称特性,即此FIR滤波器具有线性相位,则可以简化加横截型结构,下面分情况讨论:图七 N分别为奇偶数时线形相位FIR滤波器实现结构(二)FIR滤波器的级联型结构将H(z)分解成实系数二阶因子的乘积形式: 这时FIR滤波器可用二阶节的级联结构来实现,每个二阶节用横截型结构实现。如图所示:图八 FIR滤波器的级联结构这种结构的每一节控制一对零点,因而在需要控制传输零点时可以采用这种结构。二 窗函数法设计FIR滤波器如果所希望的滤波器的理想的频率响应函数为,则其对应的单位脉冲响应为 窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近。由于往往是无限长序列,而
11、且是非因果的,所以用窗函数将截断,并进行加权处理,得到:就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数为式中,N为所选窗函数的长度。我们知道,用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见附二这样选定窗函数类型和长度N之后,求出单位脉冲响应,并按照上式求出。是否满足要求,要进行演算。一般在尾部加零使长度满足2的整数次幂,以便用FFT计算。如果要观察细节,补零点数增多即可。如果不满足要求,则要重新选择窗函数类型和长度N ,再
12、次验算,直至满足要求。如果要求线性相位特性,则还必须满足 根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性正确选择其中一类,例如,要设计线性相位低通特性,可以选择这一类,而不能选择这一类。第四章 程序设计及调试第一节 设计要求给定技术指标如下: 通带截止频率Fp=3kHZ,通带最大衰减Rp=3db阻带截止频率Fs=2kHZ,阻带最小衰减Rs=14db抽样频率用双线性变换法设计一个巴特沃斯数字高通滤波器。技术指标同上,用凯泽窗法设计一个FIR数字高通滤波第二节 双线性变换法设计一个巴特沃斯数字高通滤波器一 模拟低通滤波首先应用双线性变换作为数字低通滤波
13、器到模拟高通滤波器的频率预畸变处理: 利用上式确定的技术指标和已知的参数dp、ds,就可以进行模拟原型低通滤波器设计了。该部分的程序设计与注释如下: Fs=10000; Ts=1/Fs; %采样频率和采样周期wp=3000*2*pi;ws=2000*2*pi; %低通、高通数字滤波器截止频率Rp=3;Rs=14; %通带波纹和阻带衰减wpl=Fs*2*tan(ws*Ts/2); %按频率转换公式进行转换wsl=Fs*2*tan(wp*Ts/2);N,Wn=buttord(wpl,wsl,Rp,Rs,s); %计算模拟滤波器的最小阶数z,p,k=buttap(N); %设计模拟原型低通滤波器b,
14、a=zp2tf(z,p,k); %变零极点增益形式为系统传递函数形式H,w=freqs(b,a); %求模拟低通滤波器的频率响应figure;plot(w,abs(H);grid; %绘制低通频率特性title(模拟低通原型);xlabel(w);ylabel(幅度);二 模拟到数字的转换(一)从低通向低通的转换b,a=lp2lp(Bap,Aap,Wn);功能是把模拟滤波器原型转换成截至频率为 Wn 的低通滤波器。再由得到的模拟低通滤波器通过函数Bz,Azbilinear(B,A,Fs)得到数字低通滤波器。该部分的程序设计如下:e1,f1=lp2lp(b,a,2000*2*pi) ; %原型低
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 调幅发射机 课程设计 24

限制150内