CN110361697B - 一种基于协方差矩阵混合重构的稳健波束形成方法 - Google Patents
一种基于协方差矩阵混合重构的稳健波束形成方法 Download PDFInfo
- Publication number
- CN110361697B CN110361697B CN201910616737.9A CN201910616737A CN110361697B CN 110361697 B CN110361697 B CN 110361697B CN 201910616737 A CN201910616737 A CN 201910616737A CN 110361697 B CN110361697 B CN 110361697B
- Authority
- CN
- China
- Prior art keywords
- covariance matrix
- reconstruction
- desired signal
- vector
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000013598 vector Substances 0.000 claims abstract description 64
- 238000001228 spectrum Methods 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 9
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 abstract description 5
- 230000008859 change Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 18
- 239000000654 additive Substances 0.000 description 5
- 230000000996 additive effect Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radio Transmission System (AREA)
Abstract
本发明涉及一种基于协方差矩阵混合重构的稳健波束形成方法,主要针对阵列流型误差和期望信号误差,包括:获取天线阵列的接收数据协方差矩阵,对所述接收数据协防差矩阵进行去期望信号重构,得到第一重构协防差矩阵;利用Capon谱对所述接收数据协防差矩阵进行除期望信号范围外的区域进行重构,得到第二重构协防差矩阵;对两个重构协防差矩阵进行加权求和,得到了混合重构的协方差矩阵,再利用期望信号导向矢量校正方法对期望信号导向矢量进行校正,计算得到加权矢量。本发明的方法,通过协方差矩阵混合重构得到适应于SNR变化的协方差矩阵,并且对期望导向矢量误差有校正性能,具有较好的稳健性。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种基于协方差矩阵混合重构的稳健波束形成方法。
背景技术
阵列信号处理技术在雷达,通信等很多领域都得到了广泛的应用,它的最大优势在于能够对感兴趣的方向进行波束形成从而提高信号的能量,并且在干扰方向形成增益零陷。但是,很多阵列信号处理方法都是基于理想的阵列流型,而这种假设的理想模型往往无法满足实际中的情况。实际阵列雷达***中存在很大误差,包括阵列幅相误差、位置误差以及阵元间互耦等。由于阵列信号处理算法对阵列误差的敏感性,波束形成技术性能也会受到很大影响。常规波束形成器的优点是结构简单,处理方法容易实现,但是它不能自适应的调整权值,对干扰的抑制性能较差。自适应波束形成器能够根据接收信号数据自适应的调整权值,从而利用一定的优化准则得到加权向量,在空间干扰的来波方向形成零陷,在期望信号方向形成大的增益。
目前关于阵列流型误差下的稳健波束形成算法是研究热点。近三十年来,涌现了很多稳健的波束形成算法,但是高性能的稳健自适应波束形成技术仍有很多问题有待进一步的研究。关于期望信号导向矢量失配主要有三类校正法,一类是通过一定准则对期望信号导向矢量进行优化校正,另一种是对期望信号来波区域的幅度响应进行约束,还有一类是基于协方差矩阵重构的波束形成方法。基于协方差矩阵重构的波束形成方法,尤其在高SNR(SIGNAL-NOISE RATIO,信噪比)情况下,具有较好的性能。
虽然,关于抗阵列流型误差的稳健波束算法也有很多创新,但是均存在一些问题,例如,对于接受信号有要求,或者也是基于输入SNR情况对协方差矩阵的构造有所不同,还有利用不确定集上的Capon功率谱(传统波束形成功率谱)进行积分得到精确的IPN矩阵(协方差矩阵重构矩阵),再利用凸优化技术对期望信号进行校正,但是需要较多的先验信息且计算量较大。而且,对协方差矩阵进行可变加载值的稳健波束形成算法虽然有一定的稳健性,但是输出SINR(Signal to Interference plus Noise Ratio,信号与干扰加噪声比)与理想值相差仍然较大。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于协方差矩阵混合重构的稳健波束形成方法。本发明要解决的技术问题通过以下技术方案实现:
本发明提供了一种一种基于协方差矩阵混合重构的稳健波束形成方法,包括:
在本发明的一个实施例中,所述步骤S2包括:
λN≥λN-1≥...≥λ1,
在本发明的一个实施例中,所述步骤S4包括:
S41:对去除期望信号对应的特征值λN的特征值λi重新排序,将其中的最大值定义为最大干扰值λganrao,将最后M个较小的值相加,相当于噪声的估计值,M为阵元数的1/2,再对剩余的N-1个特征值λi与所述噪声的估计值进行比较,将大于所述噪声估计值的特征值λi去除,剩余的特征值λi进行求和,得到普通噪声值λzao;
S42:根据最大干扰值λganrao、期望信号对应的特征值λN和普通噪声值λzao,定义得到,
β=λ1+λ2+...+λN-1,
a=(10log10(λNα))α,
b=(10log10(β(1-α)))(1-α),
其中,若a<0,则取a为0.00001;
在本发明的一个实施例中,所述步骤S5包括:
S51:计算得到估计期望信号导向矢量误差e⊥,
其中,subjectto表示约束条件,-1表示求逆矩阵;
在本发明的一个实施例中,所述步骤S6包括:
与现有技术相比,本发明的有益效果在于:
1、本发明的基于协方差矩阵混合重构的稳健波束形成方法,通过协方差矩阵混合重构得到适应于SNR变化的协方差矩阵,并且对期望导向矢量误差有校正性能,具有较好的稳健性。
2、本发明的基于协方差矩阵混合重构的稳健波束形成方法,所选加权值能够适应干噪比不同或差异大的情况,在有阵列流型误差时,在整个变化SNR范围有较好的输出SINR。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂,以下特举较佳实施例,并配合附图,详细说明如下。
附图说明
图1是本发明实施例提供的一种基于协方差矩阵混合重构的稳健波束形成方法的流程图;
图2是本发明实施例提供的一种阵列幅相误差下的不同算法输出SINR随SNR变化的曲线图;
图3是本发明实施例提供的另一种阵列幅相误差下的不同算法输出SINR随SNR变化的曲线图;
图4是本发明实施例提供的一种任意阵列流型误差(加性)下的不同算法输出SINR随SNR变化的曲线图;
图5是本发明实施例提供的另一种任意阵列流型误差(加性)下的不同算法输出SINR随SNR变化的曲线图;
图6是本发明实施例提供的一种无阵列流型误差下的不同算法输出SINR随SNR变化的曲线图;
图7是本发明实施例提供的一种理想阵列流型下的不同算法输出SINR随快拍数变化的曲线图;
图8是本发明实施例提供的一种阵列幅相误差下的不同算法输出SINR随快拍数变化的曲线图;
图9是本发明实施例提供的另一种阵列幅相误差下的不同算法输出SINR随快拍数变化的曲线图;
图10是本发明实施例提供的一种任意阵列流型误差(加性)下的不同算法输出SINR随快拍数变化的曲线图;
图11是本发明实施例提供的另一种任意阵列流型误差(加性)下的不同算法输出SINR随快拍数变化的曲线图。
具体实施方式
为了进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及具体实施方式,对依据本发明提出的一种基于协方差矩阵混合重构的稳健波束形成方法进行详细说明。
有关本发明的前述及其他技术内容、特点及功效,在以下配合附图的具体实施方式详细说明中即可清楚地呈现。通过具体实施方式的说明,可对本发明为达成预定目的所采取的技术手段及功效进行更加深入且具体地了解,然而所附附图仅是提供参考与说明之用,并非用来对本发明的技术方案加以限制。
实施例一
请参见图1,图1是本发明实施例提供的一种基于协方差矩阵混合重构的稳健波束形成方法的流程图,如图所示,本实施例的一种基于协方差矩阵混合重构的稳健波束形成方法,包括:
具体地,所述步骤S2包括:
λN≥λN-1≥...≥λ1,
具体地,所述步骤S4包括:
S41:对去除期望信号对应的特征值λN的特征值λi重新排序,将其中的最大值定义为最大干扰值λganrao,将最后M个较小的值相加,相当于噪声的估计值,M为阵元数的1/2,再对剩余的N-1个特征值λi与所述噪声的估计值进行比较,将大于所述噪声估计值的特征值λi去除,剩余的特征值λi进行求和,得到普通噪声值λzao;
S42:根据最大干扰值λganrao、期望信号对应的特征值λN和普通噪声值λzao,定义得到,
β=λ1+λ2+...+λN-1,
a=(10log10(λNα))α,
b=(10log10(β(1-α)))(1-α),
其中,若a<0,则取a为0.00001;
具体地,所述步骤S5包括:
S51:计算得到估计期望信号导向矢量误差e⊥,
其中,subjectto表示约束条件,-1表示求逆矩阵;
具体地,所述权矢量w为,
根据所述权矢量w得到方向图,所述权矢量w能够适应干噪比不同或差异大的情况,在有阵列流型误差时,在整个变化SNR范围有较好的输出SINR。
本实施例的基于协方差矩阵混合重构的稳健波束形成方法,通过协方差矩阵混合重构得到适应于SNR变化的协方差矩阵,并且对期望导向矢量误差有校正性能,具有较好的稳健性。
实施例二
本实施例是对实施例一中的基于协方差矩阵混合重构的稳健波束形成方法通过仿真实验进一步地说明。在本实施例的仿真实验过程中,天线阵列结构采用均匀线性阵列,阵元数N为20,阵元间距d均为λ/2,干扰方向为-30和40°,干噪比分别为25dB和35dB,理想目标方向是0°,SNR取值变化范围为[-20:40],间隔取5dB。
在本实施例的仿真输出SINR随着输入SNR变化的实验中,快拍数L均为200,天线阵列的接收数据协方差矩阵采用估计值如下:
其中,N表示阵元数,L表示快拍数。
在本实施例的仿真实验过程中,对比的算法包括Capon、CMT、多参数约束、协方差矩阵重构以及投影和对角加载结合。
对于阵列幅相误差情况,请参见图2和图3,其中,图2中的期望信号在[-4:4]范围内随机变化,Proposed表示本发明的基于协方差矩阵混合重构的稳健波束形成方法。从图中可以看出,本发明的方法在整个SNR变化范围内的输出SINR都比较高,在低SNR情况下性能更好,高SNR时性能接近协方差矩阵重构,这与设计的初衷是符合的,验证了所述加权矢量的合理性,而且在期望信号随机抖动的情况下能够保证期望信号方向增益。
对于阵列流型误差情况,请参见图4-图6,其中,图4中的期望信号在[-4:4]范围内随机变化,从图中可以看出本发明方法的性能很稳定,依然保持了整个SNR范围的高SINR。
对于快拍数变化对输出SINR性能影响,请参见图7-图11,其中,图9和图11中的期望信号在[-4:4]范围内随机变化,从图中可以看出,在幅相误差情况下,快拍数大于20后的输出SINR性能更稳定,当期望信号抖动时候和无抖动时候,本发明方法的输出SINR高于其他算法且稳定性好。在加性误差下,无论在期望信号抖动和不抖动的情况,本发明方法的性能也都略高于其他算法且稳定性强,再次验证了本发明方法的稳健性。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (5)
1.一种基于协方差矩阵混合重构的稳健波束形成方法,其特征在于,包括:
具体地,包括:
S41:对去除期望信号对应的特征值λN的特征值λi重新排序,将其中的最大值定义为最大干扰值λganrao,将最后M个较小的值相加,相当于噪声的估计值,M为阵元数的1/2,再对剩余的N-1个特征值λi与所述噪声的估计值进行比较,将大于所述噪声估计值的特征值λi去除,剩余的特征值λi进行求和,得到普通噪声值λzao;
S42:根据最大干扰值λganrao、期望信号对应的特征值λN和普通噪声值λzao,定义得到,
β=λ1+λ2+...+λN-1,
a=(10log10(λNα))α,
b=(10log10(β(1-α)))(1-α),
其中,若a<0,则取a为0.00001;
2.根据权利要求1所述的方法,其特征在于,所述步骤S2包括:
λN≥λN-1≥...≥λ1,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910616737.9A CN110361697B (zh) | 2019-07-09 | 2019-07-09 | 一种基于协方差矩阵混合重构的稳健波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910616737.9A CN110361697B (zh) | 2019-07-09 | 2019-07-09 | 一种基于协方差矩阵混合重构的稳健波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110361697A CN110361697A (zh) | 2019-10-22 |
CN110361697B true CN110361697B (zh) | 2021-08-13 |
Family
ID=68218258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910616737.9A Active CN110361697B (zh) | 2019-07-09 | 2019-07-09 | 一种基于协方差矩阵混合重构的稳健波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110361697B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112904321B (zh) * | 2019-12-03 | 2023-01-17 | 复旦大学 | 一种校准均匀线性阵列互耦和阵元增益-相位误差的方法 |
CN111988077B (zh) * | 2020-08-20 | 2023-01-31 | 中国人民解放军空军工程大学 | 信息处理方法、装置、电子设备及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105204006A (zh) * | 2015-10-19 | 2015-12-30 | 电子科技大学 | 基于子空间的干扰噪声协方差矩阵重构的波束形成方法 |
CN105306123A (zh) * | 2015-10-15 | 2016-02-03 | 哈尔滨工程大学 | 一种抗阵列***误差的稳健波束形成方法 |
CN106569181A (zh) * | 2016-11-15 | 2017-04-19 | 大连大学 | 基于协方差矩阵重构稳健Capon波束形成的算法 |
CN106788653A (zh) * | 2016-11-25 | 2017-05-31 | 西安电子科技大学 | 一种基于协方差矩阵重构的自适应波束形成方法 |
CN107167778A (zh) * | 2017-07-03 | 2017-09-15 | 电子科技大学 | 一种基于协方差矩阵重构和导向矢量估计的稳健波束形成方法 |
CN108663668A (zh) * | 2018-05-18 | 2018-10-16 | 西安电子科技大学 | 基于iaa的干扰加噪声协方差矩阵重构稳健波束形成方法 |
CN108809398A (zh) * | 2018-05-29 | 2018-11-13 | 中国科学技术大学 | 一种基于信源个数约束的稳健自适应波束形成方法 |
CN109639333A (zh) * | 2018-12-05 | 2019-04-16 | 上海神添实业有限公司 | 一种基于有效重构协方差矩阵的波束形成方法 |
CN109959899A (zh) * | 2018-04-12 | 2019-07-02 | 哈尔滨工程大学 | 特征投影预处理及协方差矩阵稀疏重构主瓣干扰抑制算法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105302936B (zh) * | 2015-08-31 | 2018-10-02 | 中国科学院声学研究所 | 基于相关计算和协方差矩阵重构的自适应波束形成方法 |
CN108181507B (zh) * | 2017-12-25 | 2020-01-03 | 中国科学技术大学 | 一种稳健自适应波束形成方法 |
CN109600152A (zh) * | 2018-12-17 | 2019-04-09 | 西北工业大学 | 一种基于子空间基变换的自适应波束形成方法 |
CN109541552B (zh) * | 2018-12-27 | 2021-04-27 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | 雷达天线阵列的自适应波束形成方法和*** |
-
2019
- 2019-07-09 CN CN201910616737.9A patent/CN110361697B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105306123A (zh) * | 2015-10-15 | 2016-02-03 | 哈尔滨工程大学 | 一种抗阵列***误差的稳健波束形成方法 |
CN105204006A (zh) * | 2015-10-19 | 2015-12-30 | 电子科技大学 | 基于子空间的干扰噪声协方差矩阵重构的波束形成方法 |
CN106569181A (zh) * | 2016-11-15 | 2017-04-19 | 大连大学 | 基于协方差矩阵重构稳健Capon波束形成的算法 |
CN106788653A (zh) * | 2016-11-25 | 2017-05-31 | 西安电子科技大学 | 一种基于协方差矩阵重构的自适应波束形成方法 |
CN107167778A (zh) * | 2017-07-03 | 2017-09-15 | 电子科技大学 | 一种基于协方差矩阵重构和导向矢量估计的稳健波束形成方法 |
CN109959899A (zh) * | 2018-04-12 | 2019-07-02 | 哈尔滨工程大学 | 特征投影预处理及协方差矩阵稀疏重构主瓣干扰抑制算法 |
CN108663668A (zh) * | 2018-05-18 | 2018-10-16 | 西安电子科技大学 | 基于iaa的干扰加噪声协方差矩阵重构稳健波束形成方法 |
CN108809398A (zh) * | 2018-05-29 | 2018-11-13 | 中国科学技术大学 | 一种基于信源个数约束的稳健自适应波束形成方法 |
CN109639333A (zh) * | 2018-12-05 | 2019-04-16 | 上海神添实业有限公司 | 一种基于有效重构协方差矩阵的波束形成方法 |
Non-Patent Citations (2)
Title |
---|
Robust Adaptive Beamforming Based on Steering Vector Estimation and Covariance Matrix Reconstruction;Feng Shen et al.;《IEEE COMMUNICATIONS LETTERS》;20150930;第19卷(第9期);全文 * |
基于协方差矩阵重构的稳健自适应波束形成算法;王昊 等;《探测与控制学报》;20180228;第40卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110361697A (zh) | 2019-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107167778B (zh) | 一种基于协方差矩阵重构和导向矢量估计的稳健波束形成方法 | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN109946664B (zh) | 一种主瓣干扰下的阵列雷达导引头单脉冲测角方法 | |
CN108919199A (zh) | 多波束成像声纳稀疏阵列的旁瓣抑制方法、阵列稀疏方法 | |
CN107276658B (zh) | 色噪声下基于协方差矩阵重构的波束形成方法 | |
CN109407055A (zh) | 基于多径利用的波束形成方法 | |
CN106707250B (zh) | 基于互耦校正的雷达阵列自适应波束形成方法 | |
CN110196410B (zh) | 一种阵列天线主瓣干扰抑制方法及*** | |
CN105306123A (zh) | 一种抗阵列***误差的稳健波束形成方法 | |
CN109639333B (zh) | 一种基于有效重构协方差矩阵的波束形成方法 | |
CN105354171B (zh) | 一种改进特征矢量的投影子空间估计自适应波束合成方法 | |
CN110361697B (zh) | 一种基于协方差矩阵混合重构的稳健波束形成方法 | |
CN107728112B (zh) | 在目标导向矢量严重失配情况下的稳健波束形成方法 | |
CN107302391A (zh) | 基于互质阵列的自适应波束成形方法 | |
CN110824414A (zh) | 一种到达角估计的装置及方法 | |
CN111880198B (zh) | 基于交替极化敏感阵列的空时极化抗干扰方法 | |
CN110708103A (zh) | 一种无需预延迟的宽带波束形成方法 | |
CN107342836B (zh) | 脉冲噪声下的加权稀疏约束稳健波束形成方法及装置 | |
CN110208757B (zh) | 一种抑制主瓣干扰的稳健自适应波束形成方法及装置 | |
Shuai et al. | Fast and accurate covariance matrix reconstruction for adaptive beamforming using Gauss-Legendre quadrature | |
CN115372925A (zh) | 一种基于深度学习的阵列稳健自适应波束形成方法 | |
CN112668155B (zh) | 一种基于二次重构的稳健波束形成方法及*** | |
Zhou et al. | Robust superdirective beamforming for HF circular receive antenna arrays | |
CN109655093B (zh) | 传感器阵列的幅相误差有源校正方法 | |
CN111257863B (zh) | 一种高精度多点线性约束的自适应单脉冲测向方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |