使用FFT实现任意三个同频带信号的频分复用.docx
《使用FFT实现任意三个同频带信号的频分复用.docx》由会员分享,可在线阅读,更多相关《使用FFT实现任意三个同频带信号的频分复用.docx(22页珍藏版)》请在沃文网上搜索。
1、 使用FFT实现任意三个同频带信号的频分复用1. 课程设计目的数字信号处理课程是一门理论和技术都发展十分迅速、应用非常广泛和前沿性学科,它的理论性和实践性都很强,学生在学习这门课程时,普遍感到数字信号处理的概念抽象,对其中的分析方法与基本理论不能很好地理解与掌握。因此,如何帮助学生理解与掌握课程中的基本概念、基本原理、基本分析方法以及综合应用所学知识解决实际问题的能力,是本课程教学中所要解决的关键问题。为了是配合数字信号处理专业基础课的理论教学,安排数字信号处理课程设计,它是针对数字信号处理的基础理论和算法进行实践环节的一个综合训练,以便学习巩固所学的知识,加强理论和实际结合的能力,培养学生的
2、综合设计能力与实际工作能力。 2. 课程设计要求选择三个不同频段的信号对其进行频谱分析,根据信号的频谱特征设计三个不同的数字滤波器,将三路信号合成为一路信号,分析合成信号的时域和频域特点,然后将合成信号分别通过设计好的三个数字滤波器,分离出原来的三路信号,分析得到的三路信号的时域波形和频谱,与原始信号进行比较说明频分复用的特点。 3 课程设计相关知识介绍 3.1 MATLAB简介 Matlab语言是一种广泛应用于工程计算及数值分析 领域的新型高级语言,Matlab功能强大、简单易学、编程效率高,深受广大科技工作者的欢迎。特别是Matlab还具有信号分析工具箱,不需具备很强的编程能力,就可以很方
3、便地进行信号分析、处理和设计。数字信号处理 常用的 Matlab 命令如下: 绘图相关abs、axis、pi、plot、grid 、title、xlabel、ylabel等。滤波器相关butter、buttord、cheb1ord、cheb2ord、cheby1、cheby2 ellip、ellipord等。3.2 信号与系统分析 实际中遇到的信号一般是模拟信号,对它们等间隔采样便可以得到时域离散信号。时域离散信号可以用集合符号、图形或公式表示。常用的典型序列有:单位采样序列、单位阶跃序列、矩形序列、实指数序列、正弦序列、复指数序列、周期序列等。系统的特性有非线性、时不变性、因果性和稳定性。信
4、号的分析主要是针对其的时域波形以及频谱关系。系统则是根据其特性分析该系统输入与输出的关系。 3.2.1信号的分析 时域分析是指控制系统在一定的输入下,根据输出量的时域表达式,分析系统的稳定性、瞬态和稳态性能。由于时域分析是直接在时间域中对系统进行分析的方法,所以时域分析具有直观和准确的优点。时域分析以线性定常微分方程的解来讨论系统的特性和性能指标。可以通过MATLAB中的plot绘出信号的时域波形图,来实现对信号的时域分析。频谱分析是为了我们能够知道一个混合信号中所有的分量信号的频率,从而我们可以得知信号的组成,频谱图表示时域信号与频谱之间的变换关系。傅立叶变换是数字信号处理领域一种很重要的算
5、法。傅立叶原理表明:任何连续测量的时序或信号,都可以表示为不同频率的正弦波信号的无限叠加。根据该原理创立的傅立叶变换算法利用直接测量到的原始信号,以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。因此,可以说,傅立叶变换将原来难以处理的时域信号转换成了易于分析的频域信号(信号的频谱)。由以上的运算流图可以看出:FFT的基本思想在于将原有的N点序列分成两个较短的序列,这些序列的DFT可以很简单地组合起来得到原有序列的DFT。所以一般我们用FFT来分析信号的的频率分量特征分析与提取。可以通过MATLAB中的fft对原始信号进行变换,再由plot绘出信号的频谱图,来实现对信号的频谱分析。3.
6、2.2系统的分析时域离散系统的频域特性:(2.1)其中、分别是系统的输出序列和输入序列的频域特性(或称为频谱特性),是系统的单位取样响应的频谱,又称为系统的频域响应。输入序列的频谱经过系统处理后即为。因此,只要按照输入信号频谱的特点和处理信号的目的, 适当选择,使得滤波后的满足设计的要求,同时也可以求出系统的输出。3.3 滤波器的设计设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统,并使系统函数H(z)具有指定的频率特性。数字滤波器从实现的网络结构或者从单位冲激响应分类,可以分成无限长单位冲激响应(IIR)数字滤波器和有限长单位冲激响应(FIR)数字滤波器。IIR数字滤波器: IIR数
7、字滤波器的系统函数为 的有理分数,即IIR数字滤波器的逼近问题就是求解滤波器的系数 和 ,使得在规定的物理意义上逼近所要求的特性的问题。如果是在s平面上逼近,就得到模拟滤波器,如果是在z平面上逼近,则得到数字滤波器。FIR数字滤波器,设FIR的单位脉冲响应h(n)为实数,长度为N,则其z变换和频率响应分别为按频域采样定理FIR数字滤波器的传输函数H(z)和单位脉冲响应h(n)可由它的N个频域采样值H(k)唯一确定。 3.3.1 设计滤波器的一般步骤1) 把给出的数字滤波器的性能指标转换为模拟滤波器的性能指标2) 根据转换后的性能指标,通过滤波器阶数选择函数,来确定滤波器的最小阶数N和固有频率W
8、n3) 由最小阶数N得到低通滤波器原型4) 由固有频率Wn把模拟低通滤波器原型转换为低通、高通、带通、带阻滤波器5) 运用脉冲响应不变法或双线性变换法把模拟滤波器转换成数字滤波器3.3.2 巴特沃斯滤波器巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦,没有起伏,而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波得图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。巴特沃斯低通滤波器可用如下振幅的平方对频率的公式表示:(2.2)其中,n = 滤波器的阶数,c =截止频率 =振幅下
9、降为 -3分贝时的 频率,p = 通频带边缘频率。n阶巴特沃斯低通滤波器的振幅和频率关系可用如下的公式表示:(2.3)G 表示滤波器的放大率, H 表示 转移函数, j 是 虚数单位, n 表示滤波器的级数, 是信号的 角频率,以弧度/秒 为单位, c 是振幅下降3分贝时的截止频率。 MATLAB信号处理工具箱提供巴特沃斯滤波器设计函数buttord和butter。 n,Wn=butterd(Wp,Ws,Rp,Rs,s)在给定滤波器性能的情况下(通带临界频率Wp、阻带临界频率Ws、通带内最大衰减Rp和阻带内最小衰减Rs),计算ButterWorth滤波器的阶数n和截止频率Wn。 b,a=but
10、ter(n,Wn,s)根据阶数n和截止频率Wn计算ButterWorth滤波器分子分母系数(b为分子系数的矢量形式,a为分母系数的矢量形式)。 3.3.3 切比雪夫I型滤波器切比雪夫滤波器(又译车比雪夫滤波器)是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“I型切比雪夫滤波器”,在阻带波动的为“II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。n阶第一类切比雪夫滤波器的幅度与频率的关系可用下列公式表示: (2.4)其中:| | 1 ,而 (
11、2.5)是滤波器在截止频率0的放大率 ,是 n阶切比雪夫多项式: (2.6) MATLAB信号处理工具箱提供切比雪夫型滤波器设计函数有cheb1ap、cheb1ord和cheby1:数字域:b,a=cheby1(n,Rp,Wn)可设计出n阶chebyshevI滤波器,其截止频率由Wn确定,通带内的波纹由Rp确定 b,a=cheby1(n,Rp,Wn,ftype)当ftype=high时,可设计出截止频率为Wn的高通滤波器;当ftypestop时,可设计出带阻滤波器z,p,k=cheby1(n,Rp,Wn) zp,k= cheby1 (n,Rp,Wn,ftype) A,B,C,D= cheby1
12、 (n,Rp,Wn)A,B,C,D= cheby1 (n,Rp,Wn,ftype)模拟域:b,a= cheby1 (n,Rp,Wn,s)可设计出截止频率为Wn的n阶chebyshevI型模拟滤波器,其余形式类似于数字域的。3.3.4 椭圆滤波器椭圆滤波器特点:幅值响应在通带和阻带内都是等波纹的,对于给定的阶数和给定的波纹要求,椭圆滤波器能获得较其它滤波器为窄的过渡带宽。椭圆滤波器可以获得对理想滤波器幅频响应的最好逼近。低通椭圆滤波器的频率响应的幅度为:(2.7)MATLAB信号处理工具箱提供椭圆滤波器设计函数ellipap 、ellipord和ellip:z,p,k=ellipap(N,Rp,
13、As) 用于计算N阶归一化中滤波器的零点向量z、极点向量p和增益因子k。Rp和As分别为通带最大衰减和阻带最小衰减。 N,wpo=ellipord(wp,ws,Rp,As,s)。用于计算满足指标的椭圆模拟滤波器的最低阶数N和通带边界频率wpo,指标要求由参数(wp,ws,Rp,As)给定。b,a=ellip(n,Rp,Rs,Wn,ftype,s)计算椭圆滤波器系统函数系数向量b和a。当ftype=high时,可设计出截止频率为Wn的高通滤波器;当ftypestop时,可设计出带阻滤波器3.3.5 以上类型滤波器的比较图3.3.5是各滤波器频域响应图:巴特沃斯滤波器(左上)和同阶第一类切比雪夫滤
14、波器(右上)、第二类切比雪夫滤波器(左下)、椭圆函数滤波器(右下)的频率响应图。巴特沃斯滤波器的衰减速度比其他类型滤波器缓慢,但十分平坦,没有幅度变化。两类切比雪夫滤波器比巴特沃斯滤波器陡峭; 但不如椭圆函数滤波器,然而后者幅度波动较大。椭圆滤波器比其他滤波器更陡,因此在选择滤波器的时候,椭圆滤波器能够以较低的阶数获得较窄的过渡带宽,但是它在通带和阻带上都有波动。 图 2.3.5 滤波器的比较 4. 设计思想选择三个不同频段的信号x1x2x3对其进行频谱分析,根据信号x1x2x3的频谱特征设计三个不同的数字滤波器,再将三路信号x1x2x3合成为一路信号x,分析合成信号的时域和频域特点,然后将合
- 1.请仔细阅读文档,确保文档完整性,对于不预览、不比对内容而直接下载带来的问题本站不予受理。
- 2.下载的文档,不会出现我们的网址水印。
- 3、该文档所得收入(下载+内容+预览)归上传者、原创作者;如果您是本文档原作者,请点此认领!既往收益都归您。
下载文档到电脑,查找使用更方便
20 积分
下载 | 加入VIP,下载更划算! |
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 使用 FFT 实现 任意 三个 频带 信号 频分复用