一种自适应优化互质阵的排布方法及目标DOA估算方法技术

本发明专利技术公开了一种自适应优化互质阵的排布方法及目标DOA估算方法，包括以下步骤，步骤1：将互质阵中的任意一个子阵位置反转，即坐标位置从正轴反转到负轴；步骤2：将步骤1中反转后的子阵进行平移，平移因子通过优化方程得到为FA；本发明专利技术自适应优化互质阵的排布算法简便、易行，具有更强的普适性；不需要针对不同单元数的互质阵重新计算最优的排布方式；利用本发明专利技术得到的差分互质阵与传统差分互质阵相比具有更多连续的单元数，孔径数更大，用于通信、成像、定位等领域，可在不增加单元个数的前提下，获得更好的信息处理效果。

【技术实现步骤摘要】
一种自适应优化互质阵的排布方法及目标DOA估算方法
本专利技术涉及天线领域，具体涉及一种自适应优化互质阵的排布方法及基于该互质阵所生成的差分互质阵的目标DOA估算方法。
技术介绍
天线阵作为通信、成像、定位等诸多领域中信号的收发设备而被广泛研究。如何在有限阵列单元的情况下实现尽可能大的孔径是天线领域中的技术难题。传统的N元天线阵由于孔径有限最多能够检测N-1个目标；为了解决这一问题，互质阵被提出来。虽然由互质阵产生的差分互质阵具有比传统互质阵或线性阵更大的孔径，即自由度较传统阵列要高，从而进一步明显的改善了信号的接收能力和检测能力。但是连续单元数还是没有被充分挖掘，即自由度没有达到可实现的最大值，洞的问题没有从根本上减少。此外，现有减小洞即提高自由度的方法操作复杂、计算量大、针对不同互质数组成的互质阵需要重新计算优化的排布方式，可操作性严重不足。
技术实现思路
本专利技术提供一种操作简单、自由度高，不需要针对不同单元数的互质阵重新计算最优排布的自适应优化互质阵的排布方法，及将该互质阵所得到的差分互质阵用于目标DOA估算的方法。本专利技术采用的技术方案是：一种自适应优化互质阵的排布方法，包括以下步骤：步骤1：将互质阵中的任意一个子阵位置反转，即坐标位置从正轴反转到负轴；步骤2：将步骤1中反转后的子阵进行平移得到优化后的互质阵，平移因子为FA。进一步的，所述的平移因子FA如下：式中：Adis为子阵列A的位置矩阵，Bdis为子阵列B的位置矩阵，函数g(a，b)为矩阵a和b的差中连续的最大值。一种基于差分互质阵目标DOA估算方法，包括以下步骤：S1：根据步骤2得到的优化后的互质阵得到优化的差分互质阵，将该差分互质阵作为收发天线阵；S2：采用步骤S1中的收发天线阵发射探测脉冲，经空间目标反射后接收目标的回波信号；S3：将步骤S2得到的回波信号进行矩阵化处理；S4：采用Capon方法计算来自不同路径目标的DOA，完成目标DOA的估算。本专利技术方法的有益效果：(1)本专利技术自适应优化互质阵的排布算法简便、易行，具有更强的普适性；不需要针对不同单元数的互质阵重新计算最优的排布方式；(2)本专利技术所生成的差分互质阵与传统互质阵得到的差分互质阵相比具有更多连续的单元数，孔径数更大，用于通信、成像、定位等领域，可在不增加单元个数的前提下，获得更好的信息处理效果；(3)本专利技术中的目标DOA估算方法，可明显提高DOA的定位精度和分辨率以及降低旁瓣。附图说明图1为传统互质阵及子阵的示意图。图2为当Mc＝3，Nc＝8时，由传统互质阵产生的传统差分互质阵的排布示意图。图3为本专利技术优化互质阵排布方法的流程图。图4为本专利技术中优化后的互质阵及其子阵的示意图。图5为当Mc＝3，Nc＝8时，由本专利技术的优化互质阵产生的优化差分互质阵的排布示意图。图6为基于传统差分互质阵和优化差分互质阵的包含DOA信息的功率谱。图7为图6的放大示意图。具体实施方式下面结合附图和具体实施例对本专利技术做进一步说明。一种自适应优化互质阵的排布方法，包括以下步骤：步骤1：将互质阵中的任意一个子阵位置反转，即坐标位置从正轴反转到负轴；步骤2：将步骤1中反转后的子阵进行平移得到优化后的互质阵，平移因子为FA。假设互质阵中的任意一个子阵A位置反转，得到优化后的互质阵子阵A；将优化后的互质子阵A与原互质阵中的子阵B组成新的互质阵，并利用新的互质阵得到优化的差分互质阵。平移因子FA如下：式中：Adis为子阵列A的位置矩阵，Bdis为子阵列B的位置矩阵，函数g(a，b)为矩阵a和b的差中连续的最大值。选取任意互质的两个自然数Mc和Nc，根据传统的互质阵排布可知阵列单元数为N＝2Mc+Nc-1；若不对该阵进行任何操作其有效孔径数为N，则所能检测到的目标数最大为N-1。为了得到更大的孔径，则需要将传统互质阵进行变形以实现虚拟均匀线阵，该线阵具有更多的连续单元。因此，可以得到更大的孔径。该变形操作所形成的阵列称为差分互质阵。相应的阵列单元位置为SDCA＝{±(ncMcd-mcNcd)}，其中(1≤nc≤Nc-1)，(1≤mc≤2Mc-1)。该种方式产生的差分互质阵产生的连续单元个数有限，即自由度有限或出现较多的洞。正是为了解决现有技术中存在的上述问题，将传统互质阵的任意一个子阵位置反转，即坐标位置从正轴反转到负轴；反转后的子阵进行平移，平移因子为FA。FA通过求解两个子阵的距离差距阵的最大连续值而获得，该距离差矩阵的最大连续值的物理意义为差分互质阵的最大连续单元个数。具体过程如下：取互质阵的两个子阵单元数分别为Mc和Nc，假设Mc＜Nc。此时，由传统方法得到的互质阵单元数为N＝2Mc+Nc-1。若不对该阵进行任何操作其有效孔径数为N，则所能检测到的目标数最大为N-1。相应的互质阵排布如图1所示，子阵A的位置为[0,Mc,2Mc,…,(Nc-1)Mc]d，子阵B的位置为[0,Nc,2Nc,…,(2Mc-1)Nc]d，d为天线工作中心频率的半波长。利用这种方式产生的差分互质阵的自由度最小为2McNc+1个，例如当Mc＝3，Nc＝8时，相应的差分互质阵单元的位置排布如图2所示。可见连续单元的个数为26×2+1＝53个，即自由度为53大于2McNc+1＝49个。该差分互质阵共有4McNc-2Nc+1＝81个孔径包含了14个洞即如图2所示深灰色的圆圈所示。因此，利用传统方式产生的差分互质阵所能实现的有效孔径数为81-14＝67个。为了在不增加阵列单元数目的前提下，得到尽可能多的连续单元个数即自由度和有效孔径数，提出本专利技术方法。将互质阵中的一个子阵例如A进行反转后得到其新的位置为[-(Nc-1)Mc,…,-2Mc,-Mc,0]d；通过变换后单元最大位置从传统的max((Nc-1)Mcd,(2Mc-1)Ncd)即从(Nc-1)Mcd，(2Mc-1)Ncd中取最大值变为[(Nc-1)Mc+(2Mc-1)Nc]d。可以看出[(Nc-1)Mc+(2Mc-1)Nc]d始终大于max((Nc-1)Mcd,(2Mc-1)Ncd)。因此，反转后的互质阵所产生的差分互质阵单元最大位置恒大于传统方法所产生的差分互质阵的单元最大位置。但是单纯的反转并不能解决连续单元个数最优的问题，因此提出优化平移因子FA的普适方法。具体方法如下：式中：Adis为子阵列A的位置矩阵，Bdis为子阵列B的位置矩阵，函数g(a，b)为矩阵a和b的差中连续的最大值；本专利技术中a＝-Adis+FAd，将子阵列A进行反转平移变换后替换原互质子阵A，而子阵B不做任何操作。因此，此方法只需要对一个子阵进行变换，较为简便。另外，优化方程得出的结果为得到连续单元数最多的差分互质阵所对应的互质阵的最优平移因子FA。FA的优化过程如图3所示，针对已设定好单元数目的互质阵，设置平移因子FA的取值范围，该范围可以很大，步长可以很小；FA取值范围和步长根据实际应用场景来确定；把FA的初始值用来求得差分互质阵，所得到的连续单元个数为Q1，FA的第二个取值用同样的方法得到差分互质阵，所得到的连续单元个数为Q2；比较Q1和Q2的大小，若Q1＜Q2成立，说明初始值不是最优解即不是最优的平移因子数FA，则排除FA取值范围的初始值，将第二个值设置为初始值，并得到相应的差分互质阵，计算并得到本文档来自技高网...

【技术保护点】
1.一种自适应优化互质阵的排布方法，其特征在于，包括以下步骤：步骤1：将互质阵中的任意一个子阵位置反转，即坐标位置从正轴反转到负轴；步骤2：将步骤1中反转后的子阵进行平移得到优化后的互质阵，平移因子为FA。

【技术特征摘要】
1.一种自适应优化互质阵的排布方法，其特征在于，包括以下步骤：步骤1：将互质阵中的任意一个子阵位置反转，即坐标位置从正轴反转到负轴；步骤2：将步骤1中反转后的子阵进行平移得到优化后的互质阵，平移因子为FA。2.根据权利要求1所述的一种自适应优化互质阵的排布方法，其特征在于，所述平移因子FA的优化方程如下：式中：Adis为子阵列A的位置矩阵，Bdis为子阵列B的位置矩阵，函数g(a，b)为矩阵a和b...

【专利技术属性】
技术研发人员：李冰，刘仕奇，赵德双，胡斌杰，孙柏，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：四川,51