阵列信号处理--ppt版本

2007-7-131阵列信号处理授课教师：廖桂生[e-mail]2007-7-132课程目的掌握空间传播波携带信号的获取与处理的基本理论和方法，特别是空间多维信号算法，熟悉参数估计和自适应波束形成的常用算法。课程要求™期间：含上机实践™期末：论文、考试2007-7-133参考文献1.PrabhakarS.naidu,SensorArraySignalProcessing2.王永良.空间谱估计理论与算法，清华大学出版社3.Monzingo.R.andMillerT.Introductiontoadaptivearray.WileyInterscience.NewYork,1980.(有中译本)4.HudsonJ.AdaptiveArrayPrinciplesPeterPeregrinusLondon,1981.(有中译本)5.HaykinS.(deitor)AdvancesinSpectrumanalysisandarrayProcessing.VolІП.PrenticeHall.NJ.19916.孙超.加权子空间拟合算法理论与应用，西北工业大学出版社7.刘德树等，空间谱估计及其应用，中国科技大学出版社8.张贤达,保铮.通信信号处理，国防工业出版社2007-7-134课程安排第一章：绪论第二章：数学基础第三章：空域滤波原理及算法第四章：部分自适应处理技术第五章：阵列信号的高分辨处理第六章：相干信源的高分辩处理第七章：最大似然与加权子空间拟和方法估计信号源方向第八章：基于高阶统计量和循环非平稳阵列信号处理简介2007-7-135第一章绪论§1.1引言一、阵列信号处理的简介阵列信号处理是信号处理领域的重要分支，它是将多个传感器设置在空间的不同位置组成的传感器阵列，并利用这一阵列对空间信号场进行接收（多点并行采样）和处理。它与一般的信号处理方式不同，因为其阵列是按一定方式布置在空间不同位置上的传感器组，主要利用信号空域特征来增强信号及有效提取信号空域信息，因此阵列信号处理也称为空域信号处理。2007-7-1361.阵列信号处理的目的：提取阵列所接收的信号及其特征信息（参数）,同时抑制干扰和噪声或不感兴趣的信息。2.阵列信号处理的系统分为两类：有源系统（具有发射传感器阵的系统）无源系统3.阵列信号处理最主要的两个研究方向：自适应空域滤波（自适应阵列处理）空间谱估计（估计信号的空域参数或信源置）2007-7-1374.阵列信号处理的研究对象：空间传播波携带信号（空域滤波）5.阵列信号处理的主要研究内容：•信（号）源定位－－确定阵列到信源的仰角和方位角•信源分离－－确定各个信源发射的信号波形，根据各个信源不同的波达方向加以区分•信道估计－－确定信源与阵列之间的传输信道的参数（多径参数）2007-7-1386.阵列信号处理的主要问题（技术）：•波束形成技术（DBF）－－使阵列方向图的主瓣指向所需的方向•零点形成技术－－使天线的零点对准干扰方向•空间谱估计－－对空间信号波达方向的分布进行超分辨估计基本概念：™传感器(sensor)－－能够感应空间传播信并且能以某种形式传输的功能装置™传感器阵列(sensorarray)－－由一组传感器分布于空间的位置构成2007-7-139由于空间传播波携带信号是空间位置和时间的四维函数，所以传播波的接收可分为：空间采集：连续－－面天线离散－－传感器阵列时间采集：所有传感器同步采样(snapshot)实际阵列空间采样的方式虚拟阵列(合成阵列，如SAR)2007-7-1310传播波的类型与媒质有关，采用的传感器也随之不同:检波器(geophone)换能器(transducer)天线(antenna)传感器大气、大地大气、水中大气（自由空间）媒质地震冲击波声波电磁波传输波传感器的空间检测能力即通常所说的方向性，是由其几何结构的形状和物理特性决定的。2007-7-1311空时采样示意图：空时处理NN元传感器阵列元传感器阵列1N2M次同步采样获取信息：波的到达方(DOA)、波形参数、极化参数估计、空间滤波与检测等2007-7-1312空间阵列传感器实际上就是空域滤波器，类同于时域滤波器。空间角方向可视为空间角频率，信号在各个角方向的功率分布可为空间功率谱或所谓的角功率谱。对比FIR滤波器的时域处理，阵列的空域处理有类似的对偶关系：FIR是在时域对时间信号作离散采样，而阵列则相当于在空域对空间信号作离散采样。因此，和FIR滤波器一样，阵列处理可对信号作一系列的运算，如滤波、分离和参数估计等，与FIR滤波器不同的是它研究的是空域信号。2007-7-1313空时等效性()xnnX对某些方向的信号加强或抑制对某些频率的信号加强或抑制滤波处理方向图传递函数系统函数空间谱（角谱）频谱谱空间采样（快拍）时间采样变元采样空域信号时序信号信号类型名称2007-7-1314时间采样与空间采样时域采样定理：Nyquist理论指出，如果一个信号在频率之外无其它频率分量，那么该信号由其整个持续期内的时间间隔为的信号采样值完全确定，从而使模拟信号可以由无限个离散的点信号来表示（拟合）。空间采样：与时间采样类似，采样频率必须足够高才不会引起空间模糊（即空间混叠），但由于受到实际条件的限制，空间采样的点数不可能无限，这相当于时域加窗，所以会出现旁瓣泄漏。f1/2f2007-7-1315时间谱与空间谱时间谱：表示信号在各个频率上的能量分布空间谱：表示信号在各个方向上的能量分布空间谱实际上就是信号的波达方向(DOA)，故空间谱估计又称为DOA估计，或者方向估计，或角度估计或测向。因为空间谱估计技术具有超高的空间信号的分辨能力，能突破并进一步改善一个波束宽度内的空间不同来向信号的分辨力，所以DOA估计是一种超高分辨的谱估计。2007-7-1316空间谱估计的系统结构目标空间观察空间估计空间信源通道1通道M通道2处理器注意：在观察空间中，通道与阵元并非一一对应，通道可由空间的一个、几个或所有阵元合成的。2007-7-1317二、阵列信号的应用•雷达：相控阵天线系统、波束灵活控制、高分辨测向、干扰置零、成像（SAR/ISAR）•移动通信：波束形成、抗多址干扰、空分多址（SDMA）•声纳：水声工程、宽带阵列处理•地震勘探：爆破、地震检测、地质层结构特征分析、探石油•射电天文：定位、测向•电子医疗工程：层析成像、医学成像2007-7-1318三、阵列信号处理的发展史雷达1936年空域信号处理只有三十多年的历史基本理论：Wiener滤波多维信号处理自60年代以来，经历了三大阶段：¾自适应波束控制IEEETransAP1964.3¾自适应零点控制IEEETransAP1976.9¾空间谱估计IEEETransAP1986.32007-7-1319wiener滤波理论应用于阵列处理（60年代）两个方向滤波方向估计自适应波束控制（指向）近代谱估计（80年代以前）自适应零点控制（70代）参数化模型（基于子空间技术）性能代价，快速算法（80年代以后）稳健算法，盲信号处理（90年代）稳健计算（90年代）2007-7-1320§1.2传播波与阵列信号处理1.传播波信号传播波信号为空时信号，是时间和空间的四维函数，服从物理规律——波动方程Maxwell波动方程：其中：A.直角坐标系中的解：a)一个特解：22221EEct∂∇=∂2222222xyz∂∂∂∇=++∂∂∂(,,,)exp[()]xyzsxyztAjtkxkykzω=−−−(,)exp[()]TsrtAjtkrω=−（（**））2007-7-1321代入波动方程：22222(,)(,)(,)(,)xyzksrtksrtksrtsrtcω++=若约束条件：22222xyzkkkcω++=即222xyzkkkkkcω==++=则：（*）式表示的信号是波动方程的解，称为“单色”或“单频”解。c为传播速度，2TπωΔ=（周期）22Tckkπωπλω=⋅=⋅=2kπλ∴=称为波数矢量，其大小表示单位波长的周期数，单位为弧度/米，其方向为波的传播方向。k2007-7-13222Tπω=2kπω=时间频率空间频率对比：某一时刻（t固定）的恒等相位面，即=常数的平面，该平面与垂直。Tkrk波动方程的任意解可以分解为无穷多个“单频”解的迭加（传播方向和频率分量均任意）。b)任意解：由四维Fourier变换表示：()()()()41,,2Tjtkrsrtsuedkdωωωπ−∞=∫∫()()(),,Tjtkrskstredrdtωω−−∞=∫∫其中2007-7-1323波动方程的单频解可以写成单变量的函数：()(),exp[()exp[]TTsrtAjtkrAjtrωωα=−=−式中，其大小等于传播速度的倒数，其方向与传播方向相同，常称为慢速矢量(slownessvector)。kαω=1cα=∵所以表示从原点传播到位置所需时间。Trαroc)波动方程另一个较复杂的解：()()()0,exp[]TTnnsrtstrsjntrαωα+∞=−∞=−=−∑由Fourier理论可知，任意周期函数，()su周期02Tπω=波形具有基本频率的调和级数形式：2007-7-1324都可以用上述级数表示，其中。()001TjnunssueduTω−=∫()(),Tsrtstrα=−α1cα=0nωω=k0kαω=量有不同的频率和波数矢量，但是各频率这时表示了具有任意波形的传播周期波，波传播方向为，速度为。波的各种分与波数矢量必须满足约束条件，可见，不同频率分量传播速度相同，但是波长不同。利用Fourier理论，波动方程更一般的解，可以表示任意波形（非周期）：2007-7-1325()()()()1,exp[]2TTsrtstrsjtrdαωωαωπ+∞−∞=−=−∫这里函数是任意的，只要其Fourier变换存在即可。该式表达了沿同一方向传播的任意波形（信号），其频率分量任意。()s⋅αB.波动方程球坐标系中的解球坐标系，但是，当波动方程的解具有球形对称时,函数并不依赖于和，使解简化，这时波动方程可简化为：(),,rφθ(),,,srtφθφθ2007-7-1326()()()()22222,,1rsrtrsrtrct∂∂=∂∂单频解为：()(),exp[]Asrtjtkrrω=−jteω近场图1.2该解可以解释为自原点向外传播的球面波，任何时刻恒等相位平面为r=常数的球面上。直角坐标系中的解为平面波，对应远场情况；球坐标系中的解为球面波，对应近场情况，如上图。2007-7-1327几点重要说明：a)对于沿一个特定方向传播的空时信号都可以表示为一元函数的形式。如果是带限信号，则由某一位置上的时间采样信号或某时刻的空间采样信号可重构全部空时信号。b)根据已知传播波的波形以及比较一些位置点上的测量信号，波的传播方向可求得。在某一瞬间空间采样提供了一组数据，用此数据有可能决定波的传播方向（如果空间采样无模糊），这是本课程的一个重要研究内容。c)应用迭加原理，允许多个传播波（不同方向、不同频率）同时出现而无交互作用。d)非理想介质对传播波有影响。（略）()Tstrα−(){}0Tmstrα−α2007-7-13282.阵列的基本原理包含在空间传播信号中的信息可以是信号源的位置，也可以是信号本身的内容。与用频率选择滤波器加强某个频率的信号一样，可以选择集中考虑从一个特定方向来的信号。由于传感器依赖机械的方向性指示，所以在同一时间一个传感器只能提取和跟踪一个方向上的信号，不能同时识别几个方向上的信号且一个传感器不能改变它的响应。但是，可以利用传感器阵列的各种加权来克服这些缺点。×××＋飞机×＋∑*1w*2w*Nw2007-7-13293.阵列信号模型考虑沿某一方向传播的窄带信号:1)窄带信号的定义与时域表示窄带信号：信号的带宽小于其中心频率的信号。在阵列信号处理里，窄带意味着信号在阵列上的延迟比信号的带宽倒数小得多，信号包络沿阵列的延迟可以忽略不计，故阵列孔径内的各阵元处复包络不变。反之，若复包络有变化，则认为是宽带信号。()(