使用FFT实现任意三个同频带信号的频分复用

使用FFT实现任意三个同频带信号的频分复用第1页共30页1引言DSP数字信号处理(DigitalSignalProcessing，简称DSP)是一门涉及许多学科而又广泛应用于许多领域的新兴学科。20世纪60年代以来，随着计算机和信息技术的飞速发展，数字信号处理技术应运而生并得到迅速的发展。数字信号处理是一种通过使用数学技巧执行转换或提取信息，来处理现实信号的方法，这些信号由数字序列表示。在过去的二十多年时间里，数字信号处理已经在通信等领域得到极为广泛的应用。DSP（digitalsignalprocessor）是一种独特的微处理器，是以数字信号来处理大量信息的器件。其工作原理是接收模拟信号，转换为0或1的数字信号，再对数字信号进行修改、删除、强化，并在其他系统芯片中把数字数据解译回模拟数据或实际环境格式。它不仅具有可编程性，而且其实时运行速度可达每秒数以千万条复杂指令程序，远远超过通用微处理器，是数字化电子世界中日益重要的电脑芯片。它的强大数据处理能力和高运行速度，是最值得称道的两大特色。DSP芯片，也称数字信号处理器，是一种特别适合于进行数字信号处理运算的微处理器，其主要应用是实时快速地实现各种数字信号处理算法。根据数字信号处理的要求，DSP芯片一般具有如下主要特点：（1）在一个指令周期内可完成一次乘法和一次加法；（2）程序和数据空间分开，可以同时访问指令和数据；（3）片内具有快速RAM，通常可通过独立的数据总线在两块中同时访问；（4）具有低开销或无开销循环及跳转的硬件支持；（5）快速的中断处理和硬件I/O支持；（6）具有在单周期内操作的多个硬件地址产生器；（7）可以并行执行多个操作；（8）支持流水线操作，使取指、译码和执行等操作可以重叠执行。使用FFT实现任意三个同频带信号的频分复用第2页共30页1.2课程设计目的通过本次课程设计，综合运用《DSP原理及应用》课程和其他有关先修课程的理论和生产实际知识去分析和解决具体问题，并使所学知识得到进一步巩固、深化和发展。初步培养学生对工程设计的独立工作能力，学习设计的一般方法。通过课程设计树立正确的设计思想，提高学生分析问题、解决问题的能力。进行设计基本技能的训练，如查阅设计资料和手册、程序的设计、调试等。是针对《DSP原理及应用》的基础理论和运用编程进行实践环节的一个综合训练，以便学习巩固所学的知识，加强理论和实际结合的能力，培养综合设计能力与实际工作能力。1.3课程设计内容通过TLC320AD50采集音频信号(f.max10kHz)，编写DSP的FFT处理程序（自定频谱分辨力），获得幅频特性之后，在点阵液晶(128*64)中大致显示出幅频图。并在液晶中用文字显示频率幅值前三的频率值。1、DSP与TLC320AD50接口电路的原理图绘制；2、DSP控制TLC320AD50的程序编写与调试；3、TLC320AD50进行语音模拟量到数字信号的转换，实现声音的采集，在CCS软件中分析信号的幅频特性；4、编写DSP的FFT处理程序；5、控制点阵液晶，实现绘图功能，将幅频图显示出来6、按要求编写课程设计报告书，正确、完整的阐述设计和实验结果。7、在报告中绘制程序的流程图，并文字说明。使用FFT实现任意三个同频带信号的频分复用第3页共30页2设计原理2.1频分复用所谓频分复用技术，是指按照频率的不同来复用多路信号的方法。在频分复用中，信道的带宽被分成若干个相互不重叠的频段，每路信号占用其中一个频段，因而在接收端可以采用适当的带通滤波器将多路信号分开，从而恢复出所需要的信号。一个简单的频分复用系统如图2.1所示。图2.1频分复用系统图中，各路基带信号首先通过低通滤波器限制基带信号的带宽，避免它们的频谱出现相互混叠。然后，各路信号分别对各自的载波进行调制、合成后送入信道传输。在接收端，分别采用不同中心频率的带通滤波器分离出各路已调信号，解调后恢复出基带信号。频分复用是利用各路信号在频率域不相互重叠来区分的。若相邻信号之间产生相互干扰，将会使输出信号产生失真。为了防止相邻信号之间产生相互干扰，应合理选择各路信号的载波频率，并使各路已调信号频谱之间留有一定的保护间隔。若基带信号是模拟信号，则调制方式可以是DSB-SC、AM、SSB、VSB或FM等，其中SSB方式频带利用率最高。若基带信号是数字信号，则调制方式可以是ASK、FSK、PSK等各种数字调制方式。使用FFT实现任意三个同频带信号的频分复用第4页共30页频分复用技术的特点是所有子信道传输的信号以并行的方式工作，每一路信号传输时可不考虑传输时延，因而频分复用技术取得了非常广泛的应用。频分复用技术除传统意义上的频分复用(FDM)外，还有一种是正交频分复用(OFDM)。传统的频分复用典型的应用莫过于广电HFC网络电视信号的传输了，不管是模拟电视信号还是数字电视信号都是如此，因为对于数字电视信号而言，尽管在每一个频道(8MHz)以内是时分复用传输的，但各个频道之间仍然是以频分复用的方式传输的。OFDM(OrthogonalFrequencyDivisionMultiplexing)实际是一种多载波数字调制技术。OFDM全部载波频率有相等的频率间隔，它们是一个基本振荡频率的整数倍，正交指各个载波的信号频谱是正交的。OFDM系统比FDM系统要求的带宽要小得多。由于OFDM使用无干扰正交载波技术，单个载波间无需保护频带，这样使得可用频谱的使用效率更高。另外，OFDM技术可动态分配在子信道中的数据，为获得最大的数据吞吐量，多载波调制器可以智能地分配更多的数据到噪声小的子信道上。目前OFDM技术已被广泛应用于广播式的音频和视频领域以及民用通信系统中，主要的应用包括：非对称的数字用户环线(ADSL)、数字视频广播(DVB)、高清晰度电视(HDTV)、无线局域网(WLAN)和第4代(4G)移动通信系统等。2.2MATLAB简介Matlab语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，Matlab功能强大、简单易学、编程效率高，深受广大科技工作者的欢迎。特别是Matlab还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行信号分析、处理和设计。数字信号处理常用的Matlab命令如下：绘图相关——abs、axis、pi、plot、grid、title、xlabel、ylabel等。滤波器相关——butter、buttord、cheb1ord、cheb2ord、cheby1、cheby2ellip、ellipord等。2.3信号与系统分析使用FFT实现任意三个同频带信号的频分复用第5页共30页实际中遇到的信号一般是模拟信号，对它们等间隔采样便可以得到时域离散信号。时域离散信号可以用集合符号、图形或公式表示。常用的典型序列有：单位采样序列、单位阶跃序列、矩形序列、实指数序列、正弦序列、复指数序列、周期序列等。系统的特性有非线性、时不变性、因果性和稳定性。信号的分析主要是针对其的时域波形以及频谱关系。系统则是根据其特性分析该系统输入与输出的关系。2.3.1信号的分析时域分析是指控制系统在一定的输入下，根据输出量的时域表达式，分析系统的稳定性、瞬态和稳态性能。由于时域分析是直接在时间域中对系统进行分析的方法，所以时域分析具有直观和准确的优点。时域分析以线性定常微分方程的解来讨论系统的特性和性能指标。可以通过MATLAB中的plot绘出信号的时域波形图，来实现对信号的时域分析。频谱分析是为了我们能够知道一个混合信号中所有的分量信号的频率，从而我们可以得知信号的组成，频谱图表示时域信号与频谱之间的变换关系。傅立叶变换是数字信号处理领域一种很重要的算法。傅立叶原理表明：任何连续测量的时序或信号，都可以表示为不同频率的正弦波信号的无限叠加。根据该原理创立的傅立叶变换算法利用直接测量到的原始信号，以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。因此，可以说，傅立叶变换将原来难以处理的时域信号转换成了易于分析的频域信号（信号的频谱）。由以上的运算流图可以看出：FFT的基本思想在于将原有的N点序列分成两个较短的序列，这些序列的DFT可以很简单地组合起来得到原有序列的DFT。所以一般我们用FFT来分析信号的的频率分量特征分析与提取。可以通过MATLAB中的fft对原始信号进行变换，再由plot绘出信号的频谱图，来实现对信号的频谱分析。2.3.2系统的分析时域离散系统的频域特性:（2.1）其中、分别是系统的输出序列和输入序列的频域特性（或称为使用FFT实现任意三个同频带信号的频分复用第6页共30页频谱特性）,是系统的单位取样响应的频谱，又称为系统的频域响应。输入序列的频谱经过系统处理后即为。因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择，使得滤波后的满足设计的要求，同时也可以求出系统的输出。2.4滤波器的设计设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统，并使系统函数H(z)具有指定的频率特性。数字滤波器从实现的网络结构或者从单位冲激响应分类，可以分成无限长单位冲激响应（IIR）数字滤波器和有限长单位冲激响应（FIR）数字滤波器。IIR数字滤波器:IIR数字滤波器的系统函数为的有理分数，即IIR数字滤波器的逼近问题就是求解滤波器的系数和，使得在规定的物理意义上逼近所要求的特性的问题。如果是在s平面上逼近，就得到模拟滤波器，如果是在z平面上逼近，则得到数字滤波器。FIR数字滤波器，设FIR的单位脉冲响应h(n)为实数，长度为N，则其z变换和频率响应分别为按频域采样定理FIR数字滤波器的传输函数H(z)和单位脉冲响应h(n)可由它的N个频域采样值H(k)唯一确定。2.4.1设计滤波器的一般步骤1）把给出的数字滤波器的性能指标转换为模拟滤波器的性能指标2）根据转换后的性能指标，通过滤波器阶数选择函数，来确定滤波器的最小阶数N和固有频率Wn3）由最小阶数N得到低通滤波器原型4）由固有频率Wn把模拟低通滤波器原型转换为低通、高通、带通、带阻滤波器5）运用脉冲响应不变法或双线性变换法把模拟滤波器转换成数字滤波器使用FFT实现任意三个同频带信号的频分复用第7页共30页2.4.2巴特沃斯滤波器巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。在振幅的对数对角频率的波得图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。巴特沃斯低通滤波器可用如下振幅的平方对频率的公式表示:（2.2）其中，n=滤波器的阶数，ωc=截止频率=振幅下降为-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]=butter(n,Wn，’s’)——根据阶数n和截止频率Wn计算ButterWorth滤波器分子分母系数（b为分子系数的矢量形式，a为分母系数的矢量形式）。2.4.3切比雪夫I型滤波器切比雪夫滤波器（又译车比雪夫滤波器）是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“I型切比雪夫滤波器”，在阻带波动的为“II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃使用FFT实现任意三个同频带信号的频分复用第8页共30页斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。n阶第一类切比雪夫滤波器的幅度与频率的关系