数字信号处理在航空航天中的运用

1数字信号处理在航空航天中的运用绪论引言数字信号处理在航空航天、遥测遥感、生物医学、自动控制、振动工程、通信雷达、水文科学等许多领域有着十分广泛的应用。特别是随着计算机及其应用的不断深入发展，数字信号处理计算机软件具有十分广阔的应用前景。如语音识别、数据压缩、医疗信号仪器的核心部分都是信号处理软件。然而信号处理算法及其软件实现包含比较复杂的理论基础知识，很难为从事其他专业领域工作的软件设计者所掌握。在实际应用中，大家面对的是数字信号处理芯片，所要完成的任务是如何对特定的芯片进行编程，但编程所要依据的知识是我们这门课要学的内容。只要明白了数字信号是如何进行处理的，数字系统是如何工作的，我们就很好对芯片进行编程。概述数字信号处理器是在模拟信号变换成数字信号以后进行高速实时处理的专用处理器，其处理速度比最快的CPU还快10-50倍。在当今的数字化时代背景下，DSP已成为通信、计算机、消费类电子产品等领域的基础器件，被誉为信息社会革命的旗手。业内人士预言，DSP将是未来集成电路中发展最快的电子产品，并成为电子产品更新换代的决定因素，它将彻底变革人们的工作、学习和生活方式。1.数字化产品DSP应用广泛，其主要应用市场为3C(Communication、Computer、Consumer-通信、计算机、消费类)领域，合占整个市场需求的90%。数字蜂窝电话、数字电视、数码相机等等，都是采用数字的方式对信号进行处理。数字信号处理就是指用数字方法处理各种信号的技术。例如，分析一下2信号的特性如频谱，就需用软件编程，求出频谱图，然后分析。再如，已知一组数字信号，我们想提取其中满足一定条件的信号，去掉不需要的成分，就可以通过数字滤波器来处理，数字滤波器可以用软件来实现，也可以用专用硬件设备来完成。总之，凡是用数字方式对信号进行滤波、变换、增强、压缩、估计、识别等都是数字信号处理的研究对象。2.数字信号处理技术的发展DSP发展历程大致分为三个阶段：70年代理论先行，80年代产品普及，90年代突飞猛进。1982年世界上诞生了首枚DSP芯片。这种DSP器件功耗和尺寸稍大，但运算速度却比CPU快了几十倍，尤其在语音合成和编码解码器中得到了广泛应用。至80年代中期，随着CMOS技术的进步与发展，第二代基于CMOS工艺的DSP芯片应运而生，其存储容量和运算速度都得到成倍提高，成为语音处理、图像硬件处理技术的基础。80年代后期，第三代DSP芯片问世，运算速度进一步提高，其应用于范围逐步扩大到通信、计算机领域。90年代DSP发展最快，相继出现了第四代和第五代DSP器件。现在的DSP属于第五代产品，它与第四代相比，系统集成度更高，将DSP芯核及外围元件综合集成在单一芯片上。这种集成度极高的DSP芯片不仅在通信、计算机领域大显身手，而且逐渐渗透到人们日常消费领域。经过20多年的发展，DSP产品的应用已扩大到人们的学习、工作和生活的各个方面，并逐渐成为电子产品更新换代的决定因素。目前，对DSP爆炸性需求的时代已经来临，前景十分可观。把本课程的主要内容给梳理一遍。《信号与系统》是这门课的先导课，因为它处理的是连续时间系统，而《数字信号处理》处理的是离散时间系统，有很多相通之处。下面我们先来回忆一下《信号与系统》这门课的主要内容。3第一章离散时间信号与系统1.1数字信号处理系统的基本组成我们来讨论模拟信号的数字化处理系统，此系统先把模拟信号变化为数字信号，然后用数字技术进行处理，最后再还原成模拟信号。这一系统的方框图见图所示。数字信号处理系统的简单方框图当然实际的系统并不一定要包括它的所有框图，例如有些系统只需数字输出，可直接以数字形式显示或打印，那么就不需要D/A变换器。另一些系统，其输入就是数字量，因而就不需要A/D变换器。对于纯数字系统，则只需要数字信号处理器这一核心部分就行了。从图中来区别几种信号：1）连续时间信号，也常称为模拟信号：)(),(tYtXzz2）离散时间信号：在一些离散时刻点有定义的信号。是数值的序列。离散时间信号可以由一个连续时间信号的采样来表示，如)(),(nTYnTXaa，也可以直接由一个离散时间过程产生。3）数字信号：时间和幅度上都离散的信号。)(),(nYnX。和离散时间信号进行区别。1.2本章重点内容本章先认识一些常用离散时间序列，重点是线性时不变离散时间系统，掌握如何判断某系统是否为LTI离散时间系统，判断系统的稳定性、因果性。)(tXa前置滤波器A/D变换器数字信号处理器D/A变换器模拟滤波器)(nX)(nY)(tYa41.3离散时间信号：序列离散时间信号在数学上表示成数值的序列。用)(nx来表示序列的第n个数，其中n为整数。这里不涉及时间，只涉及次序。连续信号经采样后变成离散时间信号，存储在存储器中。这样序列表示为{)(nx}，为方便，用)(nx表示序列。（)(),(),(kpnykx都可以用来表示序列，只要变量是离散即可）这里注意：)(nx仅仅在n为整数时才有定义，认为)(nx在n不为整数时就是零是不正确的。见书上图1-1。1．序列的基本运算方式和、积、移位（或延迟）、翻折、累加、卷积和移位：)(nx----)(mnx，m正则右移。若是)(mx----)(mnx，n正则右移。累加：nkkxny)()(，表示)(0ny在0n上的值等于)(nx在0n上的值及0n以前的所有n上的)(nx值之和。卷积和：求线性时不变离散时间系统的输出响应（零状态响应）的主要方法。定义两个序列)(nx、)(nh的卷积和为mmnhmxny)()()()()(nhnx讲解书上的例题。两种方法：图解，开车法。因为这是离散信号，和连续信号不同，后者可以积分，这里不行。由定义求解，定出求和范围，分段讨论。2．几种常见的序列单位抽样序列)(n5定义式、图形；与)(t进行比较，后者是非现实信号，实现不了；)(n的延迟形式；重要作用：任何序列都可以用一组幅度加权和延迟的冲激序列的和来表示。举例说明。31)()()3()3()2()2()1()1()()0()1()1()(kknkxnxnxnxnxnxnx推出一般式：kknkxnx)()()(单位阶跃序列)(nu定义式、图形；与)(n的关系；延迟的表达式)(mnu，)(0nnu，注意n的取值范围。矩形序列)(nRN-----有限长序列定义式、图形、N的含义；与)(n，)(nu的关系指数序列----实指数序列)(nuan和复指数序列)sin()cos(000nwjnwenjw，0w是数字域频率定义式、图形；讨论a取不同值的图形，1,01,10aaa正弦序列)cos(0nwA63．复指数序列和正弦序列的特点njwjnjwnwjAeeAeAenx0002)2()(，这样频率为)2(0rw的复指数序列相互间无法区别，这一点对正弦序列也成立：)cos())2cos(()(00nwAnrwAnx。这一切和连续时间的复指数和正弦信号不同。结论：对于A为实数的复指数序列或正弦序列，只需考虑长度为2的一段频率区间就可以。如：2000ww或。周期性问题：连续时间情况下，一个复指数信号或正弦信号都是周期的，周期等于2除以频率；离散时间情况下，情况如何？定义：如果对所有n存在一个最小的正整数N，满足)()(Nnxnx则称序列)(nx是周期性序列，周期为N。频率问题：对于连续时间正弦信号)cos()(0tAtx，随着0的增加，)(tx振荡的越来越快；对于离散时间正弦信号)cos()(0nwAnx，当0w从0增加到时，)(nx振荡的越来越快，而当0w从增加到2时，)(nx振荡反而变慢。解释：因为由推导可知0/2wrN，且N为最小正整数，0w从0增加到和0w从增加到2，如0w=（1/3）时，N=6；而0w=（2/3）时，N=3。所以可以认为0w是表示序列复现的速度。因为正弦和复指数序列中0w的周期性，在0w=2周围的频率与0w=0周围的频率区分不开。所以对于正弦和复指数序列，位于kw20邻近的0w值就属于低频范围（相对慢的振荡），而0w在0w=)2(k附近就是高频区域（相对快的振荡）。疑难问题：1）njwjnjwnwjAeeAeAenx0002)2()(，说明频率为)2(0rw的复指数序列相互间无法区别，这是由0w的周期性引起的，而不是由于)(nx的周期性引起的。2）0w为什么说是用来表示序列复现的速度。要和序列的周期公式相关联。71.4离散时间系统离散时间系统在数学上可以理解为将输入序列变换成输出序列的运算或变换。图示：表示为：)]([)(nxTny例1理想延迟系统nnnxnyd),()(其中dn为一个固定的正整数，称为系统的延迟。若dn为一个固定的负整数，则对应于时间超前。例2滑动平均)]()1()()1()([11)(11)(211212121MnxnxnxMnxMnxMMknxMMnyMMk该系统输出序列的第n个样本值等于第n个样本前后的)1(21MM个样本的平均。例2无记忆系统2))(()(nxny)(ny只决定于同一n值的输入)(nx。1．线性系统定义：满足可加性和齐次性则为线性。即：已知)]([)()],([)(2211nxTnynxTny，满足下式是线性：)]([)]([)]()([2121nxbTnxaTnbxnaxT=)()(21nbynay，ba,为任意数。注意点：ba,为任意数，证明时不能特殊。)(nx][T)(ny8例题：证明例1是线性系统，例2是非线性系统。零输入产生零输出：讨论一个输入为x(n)和输出为y(n)的任意线性系统。证明如果对于所有n，x(n)=0，则对于所有n，y(n)必然为零。证明：设y(n)=T[x(n)]。因为对所有n，x(n)=0，所以x(n)=x(n)-x(n)=0，由于线性系统满足叠加原理，因此y(n)=T[x(n)]=T[x(n)-x(n)]=T[x(n)]-T[x(n)]=0例3书上[例1-9]：)](Im[)(nxny非线性；[例1-10]：6)(4)(nxny非线性。（输入输出曲线呈线性关系，并不表示系统是线性系统）2．时不变/移不变系统定义：输入序列的错误！链接无效。将引起输出序列相应的移位和延迟。若)]([)(nxTny，对于移不变系统，对于输入为)()(1mnxnx的序列将产生)()(1mnyny。即)]([)(mnxTmny，m为任意数。例1：6)(4)(nxny是移不变的。令)()(1mnxnx时的输出为)(1ny，有)(6)(4)(1mnymnxny例2：)792sin()()9nnxny不是移不变的。例3下式定义的系统称为压缩器：nMnxny),()(，M为正整数。从M个样本中抛弃（M-1）个，也就是说，输出序列是由输入序列中每隔M个样本选出一个来构成的。该系统是移不变的。证明：当输入为)()(1mnxnx时的输出为)(1ny，有)()()(11mMnxMnxny但是)()]([)(1nymnMxmny所以不是移不变的。例4书上[例1-13]：)()(nnxny不是移不变的。当输入为)()(1mnxnx时的输出为)(1ny，有9)()()()()(1mnxmnmnymnnxny3．线性移不变系统一种特别重要的系统是由线性和移不变性组成的系统LTI，并且将线性和移不变性这两种性质结合起来可以得到LTI系统的特别方便的表示方法。用单位冲激响应表示LTI系统在前面我们得到kknkxnx)()()(，其中)(kx表示在k点上的值，设)(nx作为输入序列，通过LTI系统得到