CN103780294A - 一种最大信噪比输出相控阵天线加权方法 - Google Patents

一种最大信噪比输出相控阵天线加权方法 Download PDF

Info

Publication number
CN103780294A
CN103780294A CN201410016170.9A CN201410016170A CN103780294A CN 103780294 A CN103780294 A CN 103780294A CN 201410016170 A CN201410016170 A CN 201410016170A CN 103780294 A CN103780294 A CN 103780294A
Authority
CN
China
Prior art keywords
signal
power
array antenna
weighting
phased array
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.)
Granted
Application number
CN201410016170.9A
Other languages
English (en)
Other versions
CN103780294B (zh
Inventor
曾浩
杨清昆
幸倩
刘玲
张云
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN201410016170.9A priority Critical patent/CN103780294B/zh
Publication of CN103780294A publication Critical patent/CN103780294A/zh
Application granted granted Critical
Publication of CN103780294B publication Critical patent/CN103780294B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

一种最大信噪比输出相控阵天线加权方法,其实现步骤为:1)把阵列接收射频信号转换为数字中频信号y(k);2)根据传统相控阵天线理论计算第一组权矢量w;3)加权输出为h(k);4)估计各个通道接收信号的期望信号功率和噪声功率计算第二组权矢量a,其中第i路对应权值为6)相控阵的最终输出为相对于传统相控阵天线,本发明的相控阵天线加权方法,不会增加硬件成本,无需通道校正,能够确保阵列输出具有最大信噪比,满足工程应用需求。

Description

一种最大信噪比输出相控阵天线加权方法
技术领域
本发明涉及到相控阵天线领域,具体涉及到一种相控阵天线的加权方法。
背景技术
相控阵天线通过对接收信号进行相位加权,使得天线主瓣指向期望信号方向,从而提高了天线G/T值,它可以用于雷达、声纳、通信系统中。传统相控阵天线假设各个阵元接收信号具有相同的信噪比(Signal to Noise Ratio,SNR),权矢量仅仅对信号的相位进行校正。但由于实际每个阵元对应的通道传输特性不一致,导致仅仅是相位加权的传统相控阵天线输出信号无法达到最大SNR。虽然通过各种误差校正方法可以克服通道一致性误差,但是,一方面,校正方法本身设计和实现非常复杂,另一方面,随着相控阵天线使用时间的增加,通道特性会发生变化,导致早期的校正权值失效。总之,相控阵天线通道的不一致性和校正方法的失效,会导致各个阵元接收信号具有不同的SNR输入,此时按照传统的相控阵天线进行相位加权,无法得到最大SNR输出。
发明内容
本发明所要解决的技术问题是:由于通道不一致性误差,各个阵元间接收信号的SNR并不相同,仅仅对接收信号进行传统相控阵天线的相位校正加权,导致阵列输出无法达到最大SNR。
解决上述问题的方法是一种最大信噪比输出相控阵天线加权方法,其实现步骤是:
第一步,相控阵天线接收到的射频信号为x(t)=[x1(t)…xN(t)]T,N为天线阵元数,符号“T”表示求转置,xi(t)表示第i个阵元接收的射频信号,i=1,...,N;通过正交下变频把x(t)转换到中频,并进行模数变换,得到数字中频复信号y(k)=[y1(k)…yN(k)]T,y(k)表示第k个快拍,K是快拍数,k=1,...,K;
第二步,在已知期望信号角度情况下,根据相控阵天线基本原理,计算第一次加权采用的权矢量w=[w1...wN]T,wi表示第i个通道对应权值,i=1,...,N;
第三步,利用权矢量w分别对每个通道的中频复信号yi(k)进行加权,得到第一次加权后的输出信号h(k)=[h1(k)…hN(k)]T,其中
Figure BDA0000456916320000021
表示第i个通道第一次加权后的输出信号,符号“*”表示求共轭,i=1,...,N;
第四步,利用谱估计法,分别估计信号hi(k)中期望信号功率
Figure BDA0000456916320000022
和噪声功率
Figure BDA0000456916320000023
对第i个通道,i=1,...,N,其估计步骤为:
1)计算长度为K的信号序列hi(k)的离散傅里叶变换Hi(k),即信号频谱Hi(k),并计算其模的平方值|Hi(k)|2,从而求得接收信号总功率
Figure BDA0000456916320000024
2)将长度为K的频谱Hi(k)按序号平均分为M段,M为正整数,每段的序列长度
Figure BDA0000456916320000025
按照式计算第j段的功率,j=1,...,M;
3)对于M段的各个功率Eij,取功率最小值
Figure BDA0000456916320000027
符号min{}表示取集合元素最小值,并认为其只包含了噪声分量,从而计算得到总的噪声功率 P i n = ME i min ;
4)由接收信号总功率Pi与噪声功率可以求得期望信号功率
Figure BDA00004569163200000210
第五步,根据估计出的第i个通道信号功率
Figure BDA00004569163200000211
和噪声功率
Figure BDA00004569163200000212
计算相控阵天线第二次加权采用的权矢量a=[a1...aN]T,其中
Figure BDA0000456916320000031
i=1,...,N;
第六步,计算相控阵天线的第二次加权输出,即相控阵天线最终的输出为 z out ( k ) = Σ i = 1 N a i h i ( k ) .
本发明的有益效果是:相对传统相控阵天线,本发明的相控阵天线加权方法,无需通道校正,也不会增加硬件成本,能够确保阵列输出具有最大信噪比,满足工程应用需求。
本发明可以应用到相控阵天线领域。
附图说明
图1是均匀线阵布局和信号入射关系示意图;
图2是最大信噪比输出相控阵天线系统结构框图;
图3是最大信噪比输出相控阵天线加权方法流程图。
具体实施方法
如图1所示均匀线阵布局和信号入射示意图,对于一个均匀线阵,其阵元数为N,阵元间距为信号载波波长的一半。期望信号入射到阵列的角度θ为波达方向(Direction of Arrival,DOA),其范围
Figure BDA0000456916320000033
设第1个阵元位于原点,其余阵元在横轴正向,则阵列接收到的模拟射频信号可以用一个列矢量表示,即x(t)=[x1(t)…xN(t)]T,符号“T”表示求转置。
图2所示为最大信噪比输出相控阵天线系统结构框图。阵列接收到的模拟射频信号x(t),经过射频前端实现正交下变频变换,输出模拟中频信号,模拟中频信号经过ADC采样,输出中频数字复信号y(k),y(k)表示第k个快拍,K是快拍数,k=1,...,K。数字中频复信号先后进行两次加权,得到最终输出zout(k),该输出具有最大信噪比。第一次加权的权矢量为w=[w1...wN]T,输出信号为h(k)=[h1(k)...hN(k)]T,第二次加权的权矢量为a=[a1...aN]T,输出信号为zout(k)。
图3为最大信噪比输出相控阵天线加权方法流程图,本发明方法分为以下几个步骤:
第一步,相控阵天线接收到的射频信号为x(t)=[x1(t)…xN(t)]T,通过射频前端和ADC,得到数字中频复信号y(k)=[y1(k)…yN(k)]T。射频前端通过正交下变频,产生两路中频信号,分别代表信号实部和虚部。而每个通道分别对实部和虚部模拟信号进行ADC,得到数字复中频信号y(k)。
第二步,在已知期望信号波达方向θ情况下,根据传统相控阵天线理论,计算均匀线阵的权矢量w=[1e-jπsinθ...e-jπ(N-1)sinθ]T,用于第一次加权。
第三步,利用权矢量w分别对每个通道的中频信号yi(k)进行第一次加权,得到加权后的输出信号h(k)=[h1(k)…hN(k)]T,其中
Figure BDA0000456916320000041
表示第i个通道第一次加权后的输出信号,符号“*”表示求共轭,i=1,...,N。该权矢量仅仅是对接收信号的相位进行校正。在各个通道SNR相同情况下,相位加权输出也是具有最大SNR,但如果通道SNR不相等,则加权输出不具有最大SNR。
第四步,利用谱估计法,分别估计信号hi(k)中,期望信号功率
Figure BDA0000456916320000042
和噪声功率
Figure BDA0000456916320000043
对第i个通道,其估计步骤为:
1)计算长度为K的信号序列hi(k)的离散傅里叶变换Hi(k),即信号频谱Hi(k),Hi(k)长度也是K,k=1,...,K,然后,计算其模的平方值|Hi(k)|2,从而根据帕塞瓦尔定理求得接收信号总功率
2)将长度为K的频谱Hi(k)按序号平均分为M段,M为正整数,每段的序列长度M必须足够大,从而保证期望信号功率不能存在于所有段内,即保证有些段内仅仅存在噪声信号频谱而不存在期望信号频谱。计算每段频谱对应的功率,其中,第j段的功率计算方法为
Figure BDA0000456916320000053
j=1,...,M;
3)对于所有M段的功率,取功率最小值认为其只包含了噪声分量,并且认为每段内的噪声功率是相等的,从而计算得到总的噪声功率 P i n = ME i min ;
4)由接收信号总功率Pi与噪声功率
Figure BDA0000456916320000056
可以求得期望信号功率
Figure BDA0000456916320000057
第五步,根据估计出的第i个通道信号功率
Figure BDA0000456916320000058
和噪声功率
Figure BDA0000456916320000059
计算相控阵天线第二次加权采用的权矢量a=[a1...aN]T,其中i=1,...,N。该权矢量是通过建立信号模型,求解许瓦兹不等式,得出的理论结果,只要按照该参数加权,阵列输出就可以保证是最大SNR。
第六步,计算相控阵天线的第二次加权输出,即相控阵天线最终的输出为
Figure BDA00004569163200000511
而且输出zout(k)具有最大SNR,该最大SNR理论上等于所有阵元接收信号SNR之和。
本发明不仅可以应用于均匀线阵,也可以用于其他平面阵列。
相对于传统相控阵天线,本发明的相控阵天线加权方法,在数字域实现,不会增加硬件成本,无需通道校正,能够确保阵列输出具有最大信噪比,满足工程应用需求。

Claims (1)

1.一种最大信噪比输出相控阵天线加权方法,其实现步骤是:
第一步,相控阵天线接收到的射频信号为x(t)=[x1(t)…xN(t)]T,N为天线阵元数,符号“T”表示求转置,xi(t)表示第i个阵元接收的射频信号,i=1,...,N;通过正交下变频把x(t)转换到中频,并进行模数变换,得到数字中频复信号y(k)=[y1(k)…yN(k)]T,y(k)表示第k个快拍,K是快拍数,k=1,...,K;
第二步,在已知期望信号角度情况下,根据相控阵天线基本原理,计算第一次加权采用的权矢量w=[w1...wN]T,wi表示第i个通道对应权值,i=1,...,N;
第三步,利用权矢量w分别对每个通道的中频复信号yi(k)进行加权,得到第一次加权后的输出信号h(k)=[h1(k)…hN(k)]T,其中
Figure FDA0000456916310000011
表示第i个通道第一次加权后的输出信号,符号“*”表示求共轭,i=1,...,N;
第四步,利用谱估计法,分别估计信号hi(k)中期望信号功率和噪声功率
Figure FDA0000456916310000013
对第i个通道,i=1,...,N,其估计步骤为:
1)计算长度为K的信号序列hi(k)的离散傅里叶变换Hi(k),即信号频谱Hi(k),并计算其模的平方值|Hi(k)|2,从而求得接收信号总功率
Figure FDA0000456916310000014
2)将长度为K的频谱Hi(k)按序号平均分为M段,M为正整数,每段的序列长度
Figure FDA0000456916310000015
按照式
Figure FDA0000456916310000016
计算第j段的功率,j=1,...,M;
3)对于M段的各个功率Eij,取功率最小值
Figure FDA0000456916310000017
符号min{}表示取集合元素最小值,并认为其只包含了噪声分量,从而计算得到总的噪声功率 P i n = ME i min ;
4)由接收信号总功率Pi与噪声功率
Figure FDA0000456916310000019
可以求得期望信号功率
Figure FDA00004569163100000110
第五步,根据估计出的第i个通道信号功率和噪声功率
Figure FDA0000456916310000022
计算相控阵天线第二次加权采用的权矢量a=[a1...aN]T,其中
Figure FDA0000456916310000023
i=1,...,N;
第六步,计算相控阵天线的第二次加权输出,即相控阵天线最终的输出为 z out ( k ) = Σ i = 1 N a i h i ( k ) .
CN201410016170.9A 2014-01-14 2014-01-14 一种最大信噪比输出相控阵天线加权方法 Expired - Fee Related CN103780294B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410016170.9A CN103780294B (zh) 2014-01-14 2014-01-14 一种最大信噪比输出相控阵天线加权方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410016170.9A CN103780294B (zh) 2014-01-14 2014-01-14 一种最大信噪比输出相控阵天线加权方法

Publications (2)

Publication Number Publication Date
CN103780294A true CN103780294A (zh) 2014-05-07
CN103780294B CN103780294B (zh) 2016-11-16

Family

ID=50572188

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410016170.9A Expired - Fee Related CN103780294B (zh) 2014-01-14 2014-01-14 一种最大信噪比输出相控阵天线加权方法

Country Status (1)

Country Link
CN (1) CN103780294B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104702327A (zh) * 2015-03-24 2015-06-10 重庆大学 一种分布式相控阵接收天线的信号合成方法
CN107135022A (zh) * 2017-05-04 2017-09-05 成都航天通信设备有限责任公司 一种快速多天线加权合并方法
CN107579788A (zh) * 2017-07-29 2018-01-12 西安电子科技大学 underlay频谱共享时频重叠信号信噪比估计方法
CN109495419A (zh) * 2019-01-14 2019-03-19 哈尔滨工业大学 加权分数傅里叶变换与分组交织编码协同传输系统及方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102544751A (zh) * 2011-09-28 2012-07-04 中国电子科技集团公司第十研究所 多目标中频数字相控阵天线
CN103245934B (zh) * 2013-06-03 2014-10-29 重庆大学 一种大动态期望信号下的波束合成方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104702327A (zh) * 2015-03-24 2015-06-10 重庆大学 一种分布式相控阵接收天线的信号合成方法
CN107135022A (zh) * 2017-05-04 2017-09-05 成都航天通信设备有限责任公司 一种快速多天线加权合并方法
CN107135022B (zh) * 2017-05-04 2020-06-09 成都航天通信设备有限责任公司 一种快速多天线加权合并方法
CN107579788A (zh) * 2017-07-29 2018-01-12 西安电子科技大学 underlay频谱共享时频重叠信号信噪比估计方法
CN107579788B (zh) * 2017-07-29 2019-11-05 西安电子科技大学 underlay频谱共享时频重叠信号信噪比估计方法
CN109495419A (zh) * 2019-01-14 2019-03-19 哈尔滨工业大学 加权分数傅里叶变换与分组交织编码协同传输系统及方法
CN109495419B (zh) * 2019-01-14 2021-03-30 哈尔滨工业大学 加权分数傅里叶变换与分组交织编码协同传输系统及方法

Also Published As

Publication number Publication date
CN103780294B (zh) 2016-11-16

Similar Documents

Publication Publication Date Title
CN106227701B (zh) 一种阵列信号的接收通道的幅相误差的自校正方法
CN104506259B (zh) 宽带多天线通信系统接收通道间的时延差估计及校正方法
CN102608565B (zh) 一种基于均匀圆阵列的波达方向估计方法
JP6148622B2 (ja) レーダ装置
CN110099016B (zh) 一种基于深度学习网络的毫米波稀疏阵面信道估计方法
CN109765519B (zh) 一种模数混合天线阵列的角度估计方法
CN105319545B (zh) 提高stap检测性能的mimo-ofdm雷达波形设计方法
CN101926104B (zh) 用于使用模拟波束控制来传送信号的方法、发射站、接收站和前同步码结构
CN103780294A (zh) 一种最大信噪比输出相控阵天线加权方法
CN109787671B (zh) 一种混合波束成形装置及方法
CN102082591A (zh) 一种圆阵列天线数字波束形成方法
CN103983952A (zh) 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法
CN109581276A (zh) 一种基于求和求差嵌套阵的doa估计方法
CN105891771A (zh) 一种提高估计精度的基于连续分布的角度估计方法与设备
CN104849711A (zh) 基于频域的i-ofdm mimo雷达信号的多普勒补偿方法
CN113162670A (zh) 数字多波束校正与合成方法
CN103323827A (zh) 基于快速傅里叶变换的mimo雷达系统角度估计方法
CN105158735A (zh) 基于压缩采样阵列的空频二维谱估计方法
CN104535987A (zh) 适用于均匀圆阵列声纳系统的幅相误差自校正方法
CN107302391A (zh) 基于互质阵列的自适应波束成形方法
CN108306667A (zh) 大规模mimo中基于线性搜索的混合到达角估计
CN102544751A (zh) 多目标中频数字相控阵天线
CN106707250A (zh) 基于互耦校正的雷达阵列自适应波束形成方法
CN104168047A (zh) 利用压缩感知的单端时域波束搜索方法
CN111880198A (zh) 基于交替极化敏感阵列的空时极化抗干扰方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20161116

Termination date: 20180114

CF01 Termination of patent right due to non-payment of annual fee