1、LANZHOUUNIVERSITYOFTECHNOLOGY毕业设计题目非对称分布式目标的声成像方法研究学生姓名学号11250101专业班级通信工程（1）班指导教师学院计算机与通信学院答辩日期2015年6月15日非对称分布式目标的声成像方法研究STUDYONTHEMETHODOFASYMMETRICDISTRIBUTEDTARGETACOUSTICIMAGINGI摘要本设计介绍的是一种相干分布式目标模型的声成像理论和方法。针对分布式目标声成像的需要，建立相干的回波模型，目标的本质特征被分析出，然后利用分布式目标成像的DSPE算法，把该算法运用在对目标回波的处理上。通过MATLAB仿真分析，仿真结
2、果与仿真前进行的理论分析一致，进而分布式目标声成像的有效性被有效的验证。为实现精确定位需要构建目标的声图像大力的技术提供，研究成果对于基于分布式目标模型的声成像理论和方法的研究有着特别重要的作用。关键词分布式目标；DSPE；MATLAB仿真IIABSTRACTTHISPAPERDISCUSSESTHETHEORYANDMETHODOFTARGETACOUSTICIMAGINGBASEDONTHECOHERENTDISTRIBUTEDSOURCESMODELCONSIDERINGTHENEEDOFUNDERWATERDISTRIBUTEDSOURCEIMAGING，ACOHERENTECHOMO
3、DELISCREATEDBYANALYZINGTHEUNIQUEFEATURESOFTARGETSANDTHENTHEDSPEMETHODOFDISTRIBUTEDSOURCEIMAGINGISPROVIDEDTOPROCESSTHETARGETECHOOFTARGETIMAGINGTHEANALYSISFORMATLABSIMULATIONISCOMPLETED,ITSRESULTCOMPLYINGWITHTHETHEORETICALANALYSIS，TOPROVETHEEFFECTIVENESSOFTHEDISTRIBUTEDSOURCESIMAGINGMETHODTHETARGETIMA
4、GINGPROVIDESTECHNICALSUPPORTFORUNDERWATERDETECTIONTOLOCATEOBJECTSACCURATELYANDQUICKLYTHERESEARCHRESULTSSHOWTHATTHEYAREVALUABLETOTHETHEORETICALANDMETHODOLOGICALRESEARCHOFTARGETIMAGINGBASEDONTLLEDISTRIBUTEDSOURCESMODELKEYWORDDISTRIBUTEDTARGET；DSPE；MATLABSIMULATION兰州理工大学毕业设计III目录第1章绪论111研究背景112目标声成像研究的
5、历史和现状213本文的内容安排2第2章分布式目标声成像系统的概述321分布式目标声成像特性简介3211目标散射原理与截面积3212目标多散射中心的距离成像模型422分布式目标回波技术4第3章目标声成像回波模型原理631线列阵的概述6311线列阵的历史及发展6312线列阵的原理732亮点模型的简介833单亮点模型934多亮点回波模型9第4章分布式目标估计方法1141MUSIC类算法1142基于子空间的DSPE参数估计方法1143相干分布源的方位估计方法14第5章MATLAB程序仿真1751MATLAB概述1752MATLAB程序仿真18521不同阵元数仿真模型的比较20总结24致谢25参考文献2
6、6附录1英文文献27兰州理工大学毕业设计IV附录2英文翻译34附录3程序41兰州理工大学毕业设计1第1章绪论11研究背景随着成像声呐技术发展，声成像在各个探测领域得到了广泛的应用。利用声图像技术对目标进行识别已经成为了目标探测的一种主要手段，并有效的应用在水下作战领域，为武器智能化作战和精确打击提供充分的依据。声成像通过某些处理方法对目标的声图像处理，进而形成由主要亮点构成的目标伪图像，获取了较为完整的目标空间信息和一些细节的描述。较高质量的目标声像，能对目标准确的定位。已经在实际声纳系统中应用的方法主要包括声透镜、全息成像、波束形成。现在又出现了一种崭新的思维模式，我们可以把目标看成是分布式
7、的目标，这样一来的话就有机会从距离域上充分的反映出目标的几何形状及结构特征。基于分布式目标参数估计理论，利用均匀线阵列接收目标的反射回波，然后再运用一些相对的方法去解决，提取目标特性，进行计算并且进行仿真，对目标进行声成像，为探测体实施精确探测及定位展现了一种全新的方法。目标回波技术作为目标探测的关键技术之一其回波特性特性除与物体固有的物理属性（如外在形状、尺寸大小、固有材质）等因素有关，还与入射信号的时频特性、介质等综合因素有着紧密的联系。随着信号处理技术的不断发展，对研究目标回波的办法也是层出不穷。在其他一些国家有一部分新发明出来的水下航行体声纳但利用目标亮点回波的时域、频域特性对目标进行
8、检测和识别，更加突出的是可以利用目标在空间方位中的回波分布特性完成尺度目标的识别，那就可以有效地防止那些能够有效的模拟回波时频域特性的点源诱饵对水下航行体导引系统的引诱。因此，研究分布式目标回波特性的仿真与识别的办法，有利于提升目标识别与电子反对抗的性能。日常情况下，模拟尺度目标的亮点经常会使用不同空间方位的点源靶标延时重发脉冲的方法。与之相似的在国内有基于延时应答和多普勒频移处理的声纳目标模拟与应答及应用于雷达场合的技术研究，在此模拟条件下声纳检测出来的目标时方位有着很大的离散性，不方便对目标尺度的检测精度。目标探测和武器制导主要依赖于声探测设备来实现。然而随着声技术的发展，迅速发展迅速的还
9、有对抗和饭对抗技术对抗和饭对抗技术，还有很多人工对抗干扰器材在模拟潜艇或潜艇的辐射噪声、反射回波、运动状态等一系列目标特征方面取得了相当明显的进步。现在软杀伤和硬杀伤是主要的对抗手段，软杀伤主要是指依靠施放假目标如潜艇模拟器、声学靶、干扰器、声诱饵、气幕弹等干扰或者诱骗鱼雷，使鱼雷偏离所要攻击的目标，损失兰州理工大学毕业设计2鱼雷航程，最后能源耗尽最终无法运行。硬杀伤主要是指采用一定的装备来将鱼雷摧毁，或者使其丧失战斗能力。由于鱼雷武器自动检测、跟踪、向目标导引的能力很大一部分是依靠声学的自导，所以对抗模拟潜艇特性的声诱饵才是鱼雷目标识别真正的目标。12目标声成像研究的历史和现状目标识别时在第
10、一次世界大战末期随着声呐技术、计算机技术、信息检测等技术的发展而发展起来的。但是由于涉及到军事所以各国相对保守保密度也比较高，再加上环境比较复杂，这些都影响了该技术的发展多以目标识别技术一直发展比较缓慢。目前世界各国军队都特别重视水下目标识别技术的研究也有很多先进的技术已经用于实践，美国将目标自动识别技术列为十大关键技术之一并且为其投入了大量的财力和人力，足见该技术的重要性和意义。我国是一个资源大国，国土面积广阔。当然我国也是一个海洋大国，我国领海面积300万平方公里海岸线长18万公里长水产资源十分丰富。为了保护我国的领土完整和海上资源我国对水下目标识别技术的发展也十分重视。伴随着水下信号处理
11、技术的不断发展，声诱饵对鱼雷声自导的发射信号回波的模拟越来越真实，所以目前为止鱼雷目标识别方法的研究十分的重要，目标识别和分类的基本则是目标特征提取。虽然声诱饵等一系列的对抗器材对鱼雷声自导的发射信号回波等部分特征的模拟越来越真实，但毕竟是纸上谈兵，它的尺度信息却是模拟不了的，所以空间特征是根本模拟不了的，于是研究分布式目标回波声成像的空间特征提取对分布式目标的识别技术有着相当重要的意义。13本文的内容安排本文的主要内容和章节安排如下第一章绪论简单介绍了分布式目标声成像的研究背景。第二章目标回波系统的概述，本章主要介绍了目标回波系统的原理以及研究目标回波的关键技术。第三章目标回波模型原理简介，
12、本章核心介绍了本文的关键技术即线列阵和亮点模型以及其各自的原理发展状和况应用前景等。第四章介绍对分布是目标的估计方法，并利用该方法进行MATLAB仿真。本文最后总结了全文的内容，分析了论文中的不足之处，并对非对称分布式目标声成像方法研究做出一系列的展望。兰州理工大学毕业设计3第2章分布式目标声成像系统的概述21分布式目标声成像特性简介分布式目标回波信号是分布式目标和入射声波经过相互作用之后产生的，它的影响因素有很多，不但和入射声波的方向、分布式目标的运动姿态还与声源与分布式目标的相对位置等一系列因素有关，因此分布式目标回波与入射声波在一般情况下存在着相当大的差异，即使是均匀规则形状的物体，回波
13、与入射波也不是相同的。在主动声纳体制下，声纳与分布式目标距离较近时，在声纳发射高频声脉冲下，分布式目标的尺度特性主要包含在分布式目标反射回波中。在分布式目标回波中有分布式目标的三个重要尺度特征分布式目标回波的时间展宽；分布式目标回波的亮点起伏；分布式目标回波空间方位分布。分布式目标的特征信息包括分布式目标的位置情况，例如分布式目标的距离、方位角、府仰角和径向速度等，以及分布式目标自身特征的状况。分布式目标位置状况的参量通常用于测量目标位置，而分布式目标自己的特性参量通常用在目标识别的途径。研究分布式目标的回波就必须研究目标的特性，镜面反射特性对于分布式目标结构的不同，当声波从正横方向射入的时候
14、因为它的曲率半径比声波的波长大，回波基大概都是由镜反射产生的，这种回波的强度相当大，与目标的尺寸有关，所以能用大球体产生的回波来替代它。而这部分在很多时候都存在，而且很强，对回波的贡献也是很大的；对于因为表面不规则而引起的散射波由于表面不连续的结构，曲率半径小于波长，当入射波照到它上面时，会产生不规则散射。多数分布式目标表面都有着这种微小的不规则散射，从而形成的散射波是经过大量散射中心发出的散射波的集合。当分布式目标表面上只存在少数具有作用的散射体时，在回波的包络上能方便的看出他们形成的亮点。211目标散射原理与截面积经过理论上的计算和实验的测量都表示，在高频区，运动目标的电磁散射可以当做是某
15、些局部位置的电磁散射的组合，这些局部性的散射源一般情况下叫做多散射中心，或者称为多散射中心。分布式目标散射中心是分布式目标在高频区散射的特性之一。伴随着雷达技术的发展，人们可以利用宽带信号技术从而获取到分布式目标散射中心在径向距离上的高分辨率，然后得到赖以对多散射中心目标进行分类的大量信息，按照一定的识别算法能对这些目标进行分类。目标回波的强弱主要取决于目标散射截面积RADARCROSSSECTION（亦称为雷达截面积）的大小。对相对单一的目标而言，能看成当发生各向同性的散射时和接受到的实际分布兰州理工大学毕业设计4式目标产生相同能量回波的幻想中拦截面积。雷达截面积是下列因素的函数（1）目标结
16、构；（2）频率；（3）入射场极化形式；（4）接受天线极化形式；（5）目标对于来波方向的角向位置。因此，一般来说，可以表示为IJQ,F,式中I和J表示入射场和接受天线的极化方向，而Q,F则表示球坐标下的视角。通常，运动目标被雷达探测的时候，飞机的前端视向最受重视。飞机的雷达截面积测量值表格中有时给出其前向、尾向和侧面三个数值。如果截面积的数值是在运动的状态下测量得到的，那这个值则通常都在某段时间内起伏数值的平均值；否则，就是特殊定制的视向上的静态值。212目标多散射中心的距离成像模型雷达日标散射的特性可用一组散射中心类似于以散射中心为基本的分布式目标模型表示为（21）式中，TM是对T的M阶微分；
17、P是目标上散射中心的个数；21,MM是各散射中心的微分或积分的阶数；PMX是各散射中心的幅度；PT各散射中心的时延；因为分布式目标的时延直接与分布式目标径向的距离成比例，因此，此“冲击响应”也即目标散射中心距离像，对应TX的目标频率响应为（22）注意到式（21）与式（22）中第一部分所斯能够以的目标冲击冲击响应的细微差异；在这通过在时域中引入TM，对应的在频域表明了分布式目标上各个散射中心伴随频率的变化。由此可以推出，分布式目标距离像模型能够描述分布式目标各个散射中心在一维距离上的分布特征，同时夜比较顺利的反应了单个散射中心和分布式目标作为一个整体的频率响应的特性。22分布式目标回波技术分布式
18、目标回波声成像是反潜声自导鱼雷检测、估计和识别的基本依据。传统的分布式目标模型则是把分布式目标当成一个球体的半径随照射角变化的刚性球,也就是说目标是个单亮点,在某特定条件下对分布式目标回波模型的简化。但是扩展目标潜艇进行这样的简化,121PMPPMMMPTTXTXMPPMMMPMPMJWTPEJWXWX121兰州理工大学毕业设计5对分布式目标回波的描述是不完善的。为能够充分的反映分布式目标体积散布、不同部位反射能力和各部位之间的隐蔽效应,采用了多亮点目标回波模型。这个模型把分布式目标回波看做是由一些个亮点组成的反射体,该分布式目标回波为各亮点子回波相干加权之和。和单亮点目标模型相比,多亮点分布
19、式目标模型深层的反映了目标的体积散布和前后隐蔽特性,使分布式目标回波的延伸脉宽和幅度起伏的特性得到了反馈。经过理论分析和实验证明,在高频的情况下,任一分布式目标回波都是由若干个子回波加权相加而成的,每个子回波都被看做是从某一个散射点发出的,这个散射点就是亮点。它可能为真实的亮点,也可能为某个等效的亮点。这样以来,任一个分布式目标回波都为这些亮点回波加在一起的结果。依照形成机理可把亮点分解成几何固定亮点和镜反射亮点。前者在目标上有着相对固定的位置,它的目标强度是入射声波的照射角和俯仰角的函数,它主要是由分布式目标的几何形状决定。而后者在目标表面上可移动,且目标强度是目标轴和鱼雷之间夹角的函数,它
20、是根据波传播的声程确定的等效亮点。目标的回波模型开始于点目标或有限数目的多亮点目标模型。认为目标的回波是点目标或目标上有限数目的散射点对入射信号作用的结果。随研究的深入，单一点目标亦或是多亮点分布式目标模型在相当多的场合下很难准确的描述有关目标的特性，考察者们只能试采用更复杂的信号模型。如果目标尺度很大或观察距离很近，目标的径向尺寸大于阵列的分辨率，这时的目标可以被看作连续散射型目标。连续散射型分布式目标的回波则是通过某些基于面元法的思想来描述判断的，将分布式目标表面上的面积微元对回波的贡献在整个表面进行面积分，随之便可以轻松的获得回波响应。当然，如果表面形状复杂的话，对表面积分就很繁琐，甚至
21、得不到想要的结果，在这种情况下，可以指着将表面按特定的方式网格化，再运用数值计算法便可轻而易举的得出模拟的结果，人们在实际中发现，这些数值方法得出的结果具有足够的精确度，可以满足工程应用的要求。用面元法分析连续散射型目标回波的最大障碍在于难于得到复杂形状目标回波的解析解，如果采用数值计算，则计算量通常很大。分布源目标模型及分布式目标估计是一项新的技术，用分布源目标模型描述体目标特性具有合理性，有利于精确估计目标方位。在点目标或多亮点目标模型中，假定在能观测到的空间内只有几个空间离散点源，在空间上属于稀疏分布。而对分布源而言，入射分量在观测空间内基本上可谓是连续分布，不可从空间上分离一切的能量分
22、量。而且不能同分布源的能量扩展在观测空间中经常有着重叠，使得分布源信号在空间中不能够满足稀疏分布的特点。而对分布源信号来说，模型的建立要有很多的参数信息，因为不具备点目标模型中那些很方便的假设条件，所以可能分布源信号问题会变得相当的复杂。兰州理工大学毕业设计6第3章目标声成像回波模型原理31线列阵的概述线列阵应该是每个阵元按规定好的间距在一条直线上排列好才组成的。虽然它的结构简单，但是它的用途非常多，在舰艇声纳、近海警戒、海底探测并且在舰艇澡声测量里面都有非常多的应用。这些年以来，伴着造船技术有着长足的进步，安静型潜艇辐射澡声改变的已经非常低了，要想探测到低频澡声线谱，越来越长的线列阵被创造，
23、随之而来的阵元数也日益增多，线列阵声纳技术的发展也得到了长足的进步。世人皆知两个耳朵或者眼睛能够比做是一个阵列，然而线列阵就是模仿人体构造的最大声压法而设计出来的。非常充分的利用了澡声和信号特性不一的统计特性，把水听器在一个时刻收到的形态各异的信号，通过相关方位时间的补偿处理，然后通过能量积累，从而得到输出能亮最大的方向才是分布式目标到达的最终方向。311线列阵的历史及发展拖曳线列阵声纳的历史可以追溯到第一次世界大战，1917年，美国人HAYESHC博士展示了一种拉拽舷侧线列阵的声纳结构，在英国海军JOUT军舰舷侧的前面装了两个14元线列阵，与此同时在其尾部扯出了两个14元线列阵，这套系统能探
24、测到的大概距离差不多是19公里，能够应对那时候的潜艇攻击。历经了90多年的演变，拉扯线列阵声纳的孔径从4倍波长发生到超过100倍波长，动态范围从60DB增加到120DB。据说上个的世纪50年代在美国东海附近的海域发现了海底石油，开采海底石油成了海洋开发中的一个重要的内容，然而地质勘探是唯一的方法，彼时英国的东方石油公司以及生物探测公司、日本的PLSEY公司和加拿大的SERA公司，前前后后的使用拉扯线列阵用去接收地震回波。现实中，找到了这种拉扯式线列阵声纳那具备着远程监听检测潜艇低频辐射澡声的功能，于是日本军队最开始将它移植于声纳使用的频率范围与地址勘探不同，并且两者要求的信号的处理方式不同，而
25、拉拽阵自身的定深、噪声，使拉扯式声纳的研制碰到了瓶颈。1995年5月，发明出了一种致力于警戒的型号为AN/SQR14的远程拉扯式线列阵声呐，并且使这个声纳装在用于警戒的船上。就在这时，全新的战术拉扯式声纳由EXO公司开发出，型号为AN18。该产品附加在变深声纳拖鱼后，被32个隔振水听器组成的90M长的线列阵。接着，EXO公司使用了新的技术，制造出了没有拖鱼的长拖。兰州理工大学毕业设计7312线列阵的原理声纳信号处理技术的好坏决定了声纳工作状态。在拖曳DOUBLE线列阵声纳出现使，SINGLE线列阵遇到了很多关于信号处理方面的问题，自从拖线DOUBLE列阵声纳问世，类似的问题依旧还是声纳重中之重
26、的研究，而且伴着战术区域从深海到浅海的转移，安静型潜艇横空出世，拖曳式线列阵探测弱目标信号有了更加严格的性能要求。声纳的作用距离、测向测距精度、多目标分辨能力、目标识别以及跟踪能力等，都是实现战术指标所必须的，这些指标是达到战术指标的技术保证，如果单纯的以提高测距精度为目的，推出了延时估计精度的概念为了保证可靠性的绝对要求，湿端换能器的寿命指标被推出。上文提到的那么多的关键技术都直接关系到了声纳系统性能的探测，于是就显得非常的重要。如果按工作方式分，可以将它分成两类,被动式和主动联合式，通常情况下安装在潜艇上的只是被动式声纳，被动式拖曳线列阵声纳一般情况下可以由一下这些组成（1）线列阵接收目标
27、的辐射噪声和目标回波信号，由若干个水听器按照一定的间距布放，并采取隔振措施,配备用于水下姿态监视的传感器模块；（2）拖曳收放系统拖曳布放和回收线列阵,包括拖缆和绞车；（3）深度航向监视系统监视线列阵在水下的深度和航向信息，以便得到水听器阵列的姿态；（4）信号处理系统处理线列阵接收的各种信息，以实现目标的探测和相关参数估计；（5）显示控制系统显示信号处理系统的输出结果，与其它系统进行信息交换,并将控制命令下发给有关系统；（6）数据记录系统存储线列阵接收的数据。主被动联合式拖曳线列阵声纳除包括以上部分外,还有发射换能器基阵相应的拖曳收放系统以及发射机，接收声信号的声学模块包括水听器，前置放大器等元
28、器件“如用数字方式传输号,必须有A/D转换器、多路数据混合系统和编码电路等除了用于接收声信号的声学模块之外,还必须有若干配套部件隔振模块和仪表模块“隔振模块是一个机械滤波器,它把拖缆的可能抖动过滤掉,尽量降低抖动引起的噪声；仪表模块内通常包括深度传感器、温度传感器和航向传感器,这些传感器的信号传输到工作平台上之后,可以显示声学模块的姿态；拖缆有两方面的作用,一是作为机械的拖拉缆,二是用于信号和电流的传输。兰州理工大学毕业设计832亮点模型的简介目标是分布在相对体积内的很多反射点的一个大集体，通过分布式目标光照模型的研究得到，这些点的反射是不匀称的，反而存在很多反射能力很强的点，我们叫做亮点，亮
29、点的分布和分布式目标表面的曲率，面积和入射波的方向有着千丝万缕的联系。从入射的距离区分可以把它为远长场和近场，在远场的情况下分布式目标能够看成亮点；但是如果从能量的角度来看的话，亮点对回波有很大的奉献，而另外那些漫反射对回波的影响可以不予考虑。在高频的前提下，构成回波的各成分都能够等效成一个散射中心也就是亮点的回波，整个分布的目标相当于一组空间分布的量点。分布的目标回波在空间内的分布情况近场目标回波很大一方面是由强反射量点构成，在空间域，分布的目标回波则是由空间上来自四面八方的子反射回波相加而成，于是从方位的角度来处理目标回波，可得到目标亮点的空间方位信息。目标回波的时间方位应呈现出一定的分布
30、规律，而不是一个恒定的方位值。体目标回波的短时方位基本上随时间呈线性变化，即目标的空间方位发生了方位扩展。在已有的研究结果中，按照分布源各个分量之间的空间、时间相关性，可以将其划分为相干分布源CD、非相干分布源ID和部分相干分布源PD。对ID信号源，信号空间占据整个观测空间，噪声子空间变为零子空间。这也就是传统的子空间方法不适用于分布式目标定位的原因。SVALAEE等提出了信号子空间有效维数的概念，它实质上也就是相关矩阵的大特征值个数。噪声子空间是由小特征值对应的特征矢量构成的，而这些小的特征值是源于信号子空间，因而信号子空间与噪声子空间不完全正交。后来，YMENG等提出了DISPARE算法，
31、提高了非相关式分布源的参数估计性能。综上所述，目标的回波模型可以分为如图所示的几类分布式目标模型法的回波模型基于菲涅耳半波带近似法的回波模型基于赫姆霍兹积分近似波模型基于表面网格划分的回模型连续散射型目标的回波型三维空间多亮点回波模型垂直尺度多亮点回波模型水平尺度多亮点回波模多亮点回波模型点目标回波模型图31目标回波模型的分类兰州理工大学毕业设计933单亮点模型当目标的尺度较小或探测距离很大时，目标可以看为一个点目标或单亮点目标。从声学角度看，发出的声波可看成球面波，假设传播介质是无限均匀介质，目标的距离为R，则目标表面的入射声波为22KRWTJRERAP（31）这里是角频率，K是波数，A是发
32、射信号的幅度。入射声波随距离衰减。因为是一个点目标，入射波在表面激起二次辐射，接收到的二次辐射声压为22CRTSRATHTSTE（32）是随方位而变的反射强度。按照目标反射的声学机理，可以建立点目标反射的回波模型。设发射信号为TS，目标脉冲响应为冲激，即TTH，则回波为KRWTJIERAP（33）式中，是反射强度，R是鱼雷与目标的距离，C是声波的传播速度。回波随距离平方衰减。点目标模型是一种简单的模型，在远距离、远场时，采用点目标是合理的。但在有些应用场合并不成立，当目标的辐射面大于阵列的分辨率时，就不能认为是点目标，而是多亮点目标。超宽带窄脉冲的分辨率高，用它作为鱼雷发射信号时，从目标回波中
33、可容易分辨出多亮点。34多亮点回波模型单点模型不能确定分布的目标的许多的细小的部分，非常有可能出现欺骗和误判的情况。像前面那样说的一样，构成分布的目标的反射回波的是强反射子回波，每个子回波的时间位置和实际可观测到的量点数目的确定受攻击目标舷角影响，大多数情况下分布的目标存在5个亮点。对于建立模拟回波，确定多亮点回波模型的重中之重是要确定潜艇的亮点数和亮点分布。分布的目标的纵向尺寸比于横向尺寸要大得多，水平上的多亮点能够被分辨出来，水平的多量点模型是属于鱼雷反对抗模型中的一种。水平方向上的多量点模型结构如图32所示。建立水平尺度分布式目标多量点的回波模型的同时，应该把分布的目标与鱼雷的位置和姿态
34、加以确认，每个亮点的位置坐标）（IIZX,或极坐标）（IIR,或许会选择几何结构的尺度目标进一步解析得到，再由位置坐标或极坐标求出量点的时延和方位角参数。设目标相对鱼雷的兰州理工大学毕业设计10位置为,0RO，其中0是目标方位角，R是到观测点的距离，敌舷角为，目标长度为L。以鱼雷为中心，建立平面坐标系XOY，Z轴指向鱼雷运动方向，X轴指向鱼雷右方。当给定了参数、L、R、目标中心对应的方位角0和亮点分布后，就可由几何关系求得各亮点的时延和方位角，从而给出目标多亮点的回波信号模型图32多亮点模型的集合图假定有N个亮点等间隔分布在长度为L的目标上，根据尺度目标几何结构及所建立的坐标系，可解得亮点距离
35、参数IR，再得到时延I和方位角I。由几何关系可得）（）（）（COSL2/LR2L2/LRRI2I22I（34）式中，IL是亮点分布位置NIINLLI,2,1,11（35）其中时延为ITCR2II（36）延长线的长度0L可由关系式SINRSINLL0II0解得，方位角I为SINRSIN2/LL000（37）各亮点位置的直角坐标）（IIZX,可由极坐标）（IIR,求得，为兰州理工大学毕业设计11IIIIIICOSSINRZRX（38）设发射信号为EXP0TJTATS，其中TA为实包络，如矩形，三角形，钟形，梯形等。在理想水声信道中，每一亮点对信号的响应为一个冲激函数，可用一个卷积模型来描述目标回波
36、，则回波信号为N1IIIN1IIICR2TSTTSTHTSTE（39）式中，IR是鱼雷与潜艇第I个亮点的距离，I是亮点反射强度，它与敌舷角和目标结构及材质有关，是敌舷角和位置坐标的函数，C是声波的传播速度11。第4章分布式目标估计方法41MUSIC类算法在方位估计的领域之中，最早的构成波束类的算法因为瑞利限的限制，在分布的目标相互之间的方位有着较小的差别的时候，不可以对分布的目标方位进行很准确的估计，大势所趋慢慢的推演出了高精度的方位估算方法，其中MUSIC算法是佼佼者。在点目标估计中，因为只关心分布的目标的位置，即需要求出的参数只有一个，MUSIC算法利用接收数据协方差矩阵特征分解后的噪声子
37、空间和信号子空间的正交特性，将方位估计转化为对目标方位角的搜索，对应的谱峰位置即目标的方位角。本节利用这一原理，将此方法推广到分布源方位估计领域对目标方位角和扩展角进行二维谱峰搜索，从而得到目标参数。42基于子空间的DSPE参数估计方法子空间方法是阵列信号处理中的重要方法，不仅适用于点信号源的情形，而且还可以推广到分布式信号源参数估计。常见的分布式参数估计子空间方法有DSPE和DISPARE算法。这两种算法把信号子空间和噪声子空间的概念推广到分布式信源，实质上是MUSIC类算法，因而需要阵列流形的测量和校准。我们令2/,2/2L表示定义在区间2/,2/上的所有复值均方积分函数的希尔伯特空间。在
38、2/,2/2L中的内积和范数定义为DSSSSJIJI,2/2/（41）兰州理工大学毕业设计12CJICISSS,（42）这里C表示该函数为连续函数。则阵列输出向量可以写为NYSLXQIII1,（43）这里L表示从2/,2/2L映射到P维复向量观测空间PH的线性算子，满足PHLL2/,2/2（44）2/2/DSALS（45）空间PH的内积和范数定义为JHIJIXXXX,（46）DJIDIXXX,（47）这里D表示该函数为离散函数。根据定义，L的逆运算L满足2/,2/2LHLP（48）CDSXLSXL,（49）我们有CHHHSDSXASXDASXLXL2/2/,（410）则XAXLH（411）首先
39、，我们信号子空间和噪声子空间的概念推广到分布式目标。对于固定的参数矢量IY，角密度函数IIYS是关于DOA参数的一个随机过程。所谓信源子空间，我们是指随机过程QIYSII,1（412）的所有实现张成的线性空间。为了简单起见，我们假定信源子空间S是2/,2/2L上的闭合子空间。继而，我们推广该理论。信源子空间S在线性算子L下的值域定义为信号子空间，表示如下SSLRS（413）兰州理工大学毕业设计13R的正交补空间定义为噪声子空间，用R上表示，S的正交补空间用S上表示，则有关系RSLRSLSRLSRL（414）我们所定义的上述信号与噪声子空间的概念与点目标传统概念也是相符合的。点目标DOA为IY，
40、的角密度函数表示如下IIIIYYS（415）这里I为第I个信号的随机复包络，那么，点目标情形时的信源子空间为,1QYYSPANS（416）对,1QYYSPAN13进行线性运算L得到信号子空间恰好符合点目标时的信号子空间的一般定义。现在我们使用新定义的信号和噪声子空间解释点目标经典MUSIC算法。假定我们知道R，比如QMIREI,1,，那么QMISEAELIHI,1（417）上述矢量属于S的正交补空间，正交于S中的任一向量2/2/0DSEAIH任意SS，而且QMI,1（418）即2/2/0IIHJIHEYADYEA对于QJQMI,1,1（419）定义,1QMNEEE，我们有0NIHEYA,对于Q
41、J,1（420）MUSIC算法估计多个点源的DOA是通过最大化下面的“频率探测器”，211NHHNNHMUSICEYAYAEEYAYP（421）这里,Y为待估参数集合。沿着类似的思路，我们可以推导出适用于分布式信号源的MUSIC类算法。现在假定R维数是QM，基为QMEE,1并令,1QMNEEE，且有兰州理工大学毕业设计14QMISEAELIHI,1,（422）因此，对于所有SS，我们有02/2/DYSEAIINH（423）由于S空间为IIYS张成的，故上式又可写为2/2/0DYSEAIINH（424）对于所有,1QIYSII的实现都成立既然IIYS是随机函数，我们还有022/2/DYSEAEI
42、INH（425）2/2/2/2/,DDAEYPEAHNINH（426）那么，参数矢量可由下式估计得到1MAXARGNHNYEYHETRY（427）其中称2/2/,DDAYPAYHHI为矩阵阵列流形。上述算法为分布式信号参数估计方法DSPE。43相干分布源的方位估计方法分布式目标可以看作一些满足一定分布规律的散射点的集合,这些散射点之间有所联系,所以和空间中的多目标是有所区别的。对于这样的目标源,接受阵列的观测矢量为QIIITNDTSATX1,（428）在上面的式子中,A是阵列流型向量,TSII为第I个分布式目标信号源在T时刻的角信号密度函数,它包含了分布式目标的所有信息,分布式目标波达方向就是
43、指中心波达方向I。对于源内角度相关源来自同一源的不同角度的来波信号成分相关来说，角信号密度函数可以表示为IIIIIGTSTS,其中TSI为目标反射信号,IIG为第I个分布式目标的角信号分布函数,则T时刻阵列的观测数据TX可以表示为QIIIITNBTSTX1（429）其中IIB是第I个分布式目标的方向矢量,对比点目标的情况,也可以把它理解为广义流型兰州理工大学毕业设计15向量,可以表示为QIDGABIIII,2,1（430）Q为分布式目标个数,IIB由分布式目标信号源的角信号分布函数的数学形式及其未知的分布参数决定,是M阶方向矢量,M为阵元数。假设远场存在Q个分布式目标,信号矢量TQTSTSTS
44、TS,21和噪声矢量TN是相互独立的零均值随机矢量，其二阶原点矩分别为,2KTIKNTNEKTPKSTSEMNHSH（431）其中,KT为KRONNECKER函数,SP为信号协方差矩阵,MI为M阶单位矩阵,2N为噪声方差。由QIIIITNBTSTX1可以得出,观测数据矢量TX的协方差矩阵XRNIHHXTXTXNKXTXER11,2,1,22,21,2,12,11,1MMRMRMRMRRRMRRRXXXXXXXXX432其中N为观测数据长度,在目标源之间不相干的情况下,均匀线列阵信号的协方差矩阵是TOEPLITZ14矩阵,矩阵的秩等于目标源个数,通过对XR特征分解后的子空间的分析,利用方向矢量与
45、噪声子空间的正交特性来构造算法的空间谱,由谱峰的位置求出目标方位。但是当多个目标同时入射进入阵列并且目标信号相干时,矩阵将不再是TOEPLITZ矩阵,并且由于目标源之间相干,导致XR的秩亏损,从而可能导致目标漏报。为了解决这一问题,提出对该算法的改进方法,构造一个TOEPLITZ矩阵1,1,1,1,2,2,1RMRMRMRRRMRRRRT（433）其中MMXMIIMMRIMIR1,1,1,11，即元素IR是协方差矩阵XR下三角部分与主对角线平行的第I条线上元素的平均。其实TOEPLITZ用来解相关时有三种途径,其一是相位TOEPLITZ方法,其二是幅度TOEPLITZ方法,在这里用到的是幅度和
46、相位同时TOEPLITZ化。此矩阵可以理解为将采样数据的协方差矩阵进行时间平均和空间平均后的数据矩阵。这是因为兰州理工大学毕业设计16对于矩阵XR可以理解为其每一元素是不同的接收阵元采样数据的时间平均,而对其采取TOEPLITZ化则意味着在此基础上将不同的时间平均值再进行一次阵元之间的平均,即空间平均,这样做的目的是恢复矩阵的秩。此时的协方差矩阵RT称之为修正协方差矩阵，TOEPLITZ化解相干的原理不是通过降低阵元自由度而获得的,而是通过改变协方差矩阵的数据结构获得的,所以这与平滑类解相干算法不同,阵列的孔径得到了有效的利用。修正后的协方差矩阵的奇异值分解为HSSSHSSSTUUUUR,对应
47、信号子空间及噪声子空间有不同的奇异值理想情况下有NQQ121为了避免由于目标之间的相干性带来的矩阵秩的亏损从而导致对于目标源个数估计的不准确,这里使用对于目标源欠估计下的处理方法,结合针对相干源的TOEPLITZ方法,可以更有效的处理相干源的方位估计问题。一般情况下,由矩阵方向矢量及噪声子空间的正交性导出的MUSIC类方法用到的噪声子空间是由特征向量张成的子空间,并未考虑特征值在其中的影响,这样做的前提必须是对于目标源的个数估计十分准确的前提下,如果对于目标源的先验知识不足,或者说估计的不够准确,则此类方法往往不能做到准确估计每一个源的方位。为了解决这一问题,考虑到特征值的作用,在此对于噪声子
48、空间的重新定义为加权噪声子空间,如下1,1,12211MMNVVVU（434）其中M,21为修正协方差矩阵奇异值分解后的各个奇异值,此式表明噪声子空间的构成用到了全部奇异值及对应的奇异向量,相干信号导致的结果是信号子空间扩散到噪声子空间,用整个子空间的加权来构造噪声子空间可以充分利用所用的信息。对于信号子空间部分,由于奇异值较大,所以对于噪声子空间的贡献就小,这样一来,如果出现了对与源数目的漏报,就可以避免由于对源数目估计的不准确导致的估计结果错误,此时的噪声子空间称为加权的噪声子空间。利用方向矢量和加权后噪声子空间的正交性,即0IIHNBU16,可以得出DOA估计方法的多维谱峰搜索法1BUUBFHNNH（435）这里假定分布式目标信号源的角信号分布函数形式已知,具体表达式由未知参数I决定,I为含有
