Multi-carrier Index Modulation Based on Prolate Spheroidal Wave Functions with Dual-mode
-
摘要: 该文在基于信号分组优化的椭圆球面波函数(PSWFs)多载波调制的基础上,引入双模索引调制思想,提出了双模PSWFs多载波索引调制解调方法(DM-MCM-PSWFs)。该方法利用未被激活的剩余部分子载波加载第2星座图产生的调制符号,以传输额外的信息比特,实现了基于信号分组优化的PSWFs多载波调制中频谱资源的进一步利用,有效提高了系统频带利用率以及误码性能。理论和仿真分析表明,相较于基于信号分组优化的PSWFs多载波调制,所提方法以适当牺牲系统复杂度为代价,具有更高的系统频带利用率和更优的系统误码性能,当误比特率为1×10–5, n=7, k=3时,所提方法系统频带利用率、误码性能可分别提升9.15%, 2.4 dB。Abstract: With the foundation of multi-carrier modulation based on Prolate Spheroidal Wave Functions (PSWFs) with signal grouping optimization, the idea of dual-mode index modulation is introduced, and proposes a method of Multi-Carrier index Modulation based on PSWFs with Dual-Mode (DM-MCM-PSWFs) is introduced. The remaining sub-carriers that are not activated to load the modulation symbols generated by the second constellation to transmit additional information bits are used in this method, resulting in the improvements of efficiency of spectrum resources in multi-carrier modulation based on PSWFs with signal grouping optimization, and it improves effectively the system spectral efficiency and error performance. Theoretical and simulation analysis indicate that compared with multi-carrier modulation based on PSWFs with signal grouping optimization, the proposed method can attain higher system spectral efficiency and better system error performance at the expense of acceptable sacrifice of system complexity. When the bit rate is 1×10-5, n=7, and k=3, the system spectral efficiency and error performance of the proposed method can increase by 9.15% and 2.4 dB respectively.
-
1. 引言
相比较于传统的相控阵雷达,多输入多输出(MIMO)雷达由于其独特的优异特性正在引起越来越多的关注[1–3]。MIMO雷达发射正交波形,根据其阵列结构可以分为统计MIMO雷达[2]和集中式MIMO雷达[3]。统计MIMO雷达天线阵列各阵元之间有较大间隔,利用空间分集技术使雷达反射起伏对MIMO雷达造成的影响大大降低。集中式MIMO雷达发射阵元和接收阵元都紧密排列,利用波形分集技术和多个通道相干处理,可以改善弱目标和低速目标检测性能,并且能对更多的目标进行角度估计。集中式MIMO雷达根据收发基地是否分置分为单基地MIMO雷达和双基地MIMO雷达。本文主要研究的是双基地MIMO雷达的测角问题。
文献[4]提出一种基于2维搜索的Capon算法对双基地MIMO雷达的波离角(Direction Of Departure, DOD)和波达角(Direction Of Arrival, DOA)进行估计。由于需要进行2维搜索,所以Capon算法的计算量很大。为了解决这个问题,文献[5]利用均匀线阵的旋转不变特性,提出ESPRIT算法对目标角度进行估计。文献[6]提出的传播因子算法(Propagator Method, PM)不需要对接收信号的自相关矩阵进行特征分解,计算量小,但是估计性能有所损失。文献[7]提出了基于十字型阵列的MIMO雷达降维ESPRIT算法,虽然该算法是基于单基地MIMO雷达,但同样是对2维角度进行估计。文献[8]提出的多目标跟踪算法解决了联合对角化方法低信噪比信号子空间扩展的问题。文献[9]提出了一种基于实值处理的酉ESPRIT算法。该算法采用前后平滑技术[10]将接收信号从复数域转换到实数域,能够有效降低运算量。由于采用前后平滑技术,酉ESPRIT算法最多可以处理两个相干信源。文献[11]提出了色噪声条件下的双基地MIMO雷达角度估计方法,能够有效抑制色噪声影响,提高测角精度。文献[12]提出的半实值MUSIC算法对自相关矩阵的实部进行特征分解,可对任意阵列的双基地MIMO雷达角度进行估计。文献[13]提出了一种基于波束域的ESPRIT算法。波束域ESPRIT算法将接收信号从阵元域转换到波束域,降低信号维度,同时提高角度估计性能。
近年来,基于张量的双基地MIMO雷达测角算法得到了广泛的关注。张量算法将接收信号构造成一个多维张量模型,通过张量分解提取角度信息。文献[14]提出了一种高阶奇异值分解(Higher-Order Singular Value Decomposition, HOSVD)张量算法。该算法首次将前后平滑技术应用到张量模型中用于构建实值张量模型,并对该实值张量进行HOSVD得到信号子空间。通过HOSVD方法得到的信号子空间相比于传统的矩阵方法[4–13]有更高的精确度,因此测角性能会得到提升。文献[15]提出了基于平行因子分解(PARallel FACtor decomposition, PARAFAC)的张量算法。PARAFAC算法不对信号子空间进行估计,目标的角度信息直接从分解得到的加载矩阵中提取。文献[16]提出的基于实值处理的酉PARAFAC算法证明了文献[14]中的实值张量模型仍然满足PARAFAC分解条件。文献[17]将实值PARAFAC算法用于双基地MIMO雷达互耦条件下的角度估计。张量算法需要对信号进行多维处理,计算复杂度极高。
本文提出一种基于实值处理的联合波束域双基地MIMO雷达测角方法。本文算法的发射和接收波束矩阵通过凸优化方法进行设计,在保证旋转不变结构的同时,滤波器的带宽可以灵活控制且旁瓣电平可以被有效抑制,达到提高主副瓣比的目的。更重要的是,本文通过凸优化方法对波束矩阵的结构进行设计,这种特殊结构是将滤波后信号构造成实值模型的必要条件。最后,本文通过建立映射关系的方法对插值误差进行补偿,提高测角精度。仿真实验验证了本文算法的有效性。
这里给出文中常用符号的定义:
(⋅)T ,(⋅)H 和(⋅)∗ 分别表示矩阵或矢量的转置,共轭转置以及共轭。diag(⋅) 表示对角化操作。vec(⋅) 表示矢量化操作。⊗ 和⊙ 分别表示Kronecker积和Khatri-Rao积。Re{⋅} 和Im{⋅} 分别表示取矩阵或者矢量的实部和虚部。IM 表示M×M 维的单位矩阵。2. 信号模型
假设双基地MIMO雷达发射和接收阵列为均匀线阵,其中发射阵元数为
M ,接收阵元数为N 。假设发射和接收阵列的阵元间距均为d ,为了避免相位模糊,需要d≤λ/2 ,λ 表示信号波长,本文采用d=λ/2 。M 个发射阵元发射正交信号S=[s1(t), s2(t),···,sM(t)]T∈CM×L ,L 表示单个脉冲内的采样点数。由于发射信号彼此正交,因此满足SSH/L=IM 。假设远场目标的个数为K ,第k 个目标的DOD和DOA分别用θk 和φk 表示。一个相干处理间隔内的脉冲数为Q ,则接收阵列接收到的第q 个脉冲信号可以表示为Xq = BΛ qATS+Nq (1) 其中,
A=[a(θ1),a(θ2),···,a(θK)] 和B=[b(φ1), b(φ2),···,b(φK)] 分别表示发射和接收导向矩阵。第k 个目标的发射导向矢量为a(θk)=[e−jπ(M−1)sinθk/2,···,1,···,ejπ(M−1)sinθk/2]T (2) 第
k 个目标的接收导向矢量为b(φk)=[e−jπ(N−1)sinφk/2,···,1,···,ejπ(N−1)sinφk/2]T (3) 矩阵
Λ q 的形式为Λ q=diag(cq) ,其中矢量cq=[α1,q,α2,q,···,αK,q]T 。{αk,q}Qq=1 表示第k 个目标的反射系数。Nq 表示均值为零的高斯白噪声矩阵。对式(1)右乘SH/L 进行匹配滤波得到Yq=BΛ qAT+¯Nq (4) 其中,
¯Nq=NqSH/L 表示匹配滤波后的噪声项,Yq 为N×M 维矩阵。将Yq 按列重排为列向量yq=(A⊙B)cq+¯nq (5) 式中,
yq=vec(Yq) ,¯nq=vec(¯Nq) 。将Q 个脉冲信号集中写成MN×Q 的矩阵形式Y=(A⊙B)C+˜N (6) 其中,
C=[c1,c2,···,cQ]∈CK×Q ,˜N=[¯n1,¯n2,···, ¯nQ]∈CMN×Q 。3. 实值联合波束域测角算法
3.1 空域滤波器设计
定义一个
¯M×1 (¯M<M) 维的虚拟发射导向矢量¯a(θ)=[e−jπ(¯M−1)sinθ/2,···,1,···,ejπ(¯M−1)sinθ/2]T (7) ¯M 称为发射波束空间的维度。为了保证波束域转换后的旋转不变结构,本小节目的之一是设计一个MׯM 的矩阵Wt 使WHta(θ) 在感兴趣的目标区域尽可能等于¯a(θ) 。同时,要求空域滤波器的旁瓣电平要尽可能地低。更重要的是,本文对Wt 的矩阵结构有特殊的要求,下一小节将证明这种结构是后续实值处理提取角度信息的必要条件。当M 是偶数时,Wt 可以通过求解式(8)最小最大化问题得到minWt,Γ tmaxj||WHta(θj)||,θj∈¯Θt,j=1,2,···,Jt,s.t.||WHta(θi)−¯a(θi)||≤α,θi∈Θ t,i=1,2,···,It,Wt=[Γ Tt,Π ¯MΓ HtΠ M/2]T (8) 式(8)中,
Γ t∈C(M/2)ׯM 。Π n 表示一个n×n 维的交换矩阵,仅在反对角线上具有元素1,其它元素等于0。Θ t 表示感兴趣的目标区域(相对于发射端),¯Θt 表示Θ t 的补集。It 和Jt 表示在Θ t 和¯Θt 的采样点数。正参数α 限制了在Θ t 区域最大可接受的WHta(θ) 和¯a(θ) 之间的偏差。如果M 是奇数,Wt 可以通过求解式(9)得到minWt,Γ t,gtmaxj||WHta(θj)||,θj∈¯Θt,j=1,2,···,Jt,s.t.||WHta(θi)−¯a(θi)||≤α,θi∈Θ t,i=1,2,···,It,Wt=[Γ Tt,gt,Π ¯MΓ HtΠ (M−1)/2]T,g∗t=Π ¯Mgt (9) 其中,
gt 是一个¯M×1 维矢量,Γ t∈C(M−1)/2ׯM 。式(8)和式(9)里的等式约束是实值处理提取角度信息的必要条件,在下一小节中进行证明。凸优化问题式(8)和式(9)可以用CVX工具箱[18]进行求解。同样,当
N 是偶数时,NׯN (¯N<N) 维的矩阵Wr 可以通过求解式(10)最小最大化问题得到:minWr,Γ rmaxj||WHrb(φj)||,φj∈¯Θr,j=1,2,···,Jr,s.t.||WHrb(φi)−¯b(φi)||≤β,φi∈Θ r,i=1,2,···,Ir,Wr=[Γ Tr,Π ¯NΓ HrΠ N/2]T (10) 式(10)中,
Γ r∈C(N/2)ׯN ,Θ r 表示感兴趣的目标区域(相对于接收端),¯b(φ)=[e−jπ(¯N−1)sinφ/2,···, 1,···,ejπ(¯N−1)sinφ/2]T 表示¯N×1 维的虚拟接收导向矢量。¯N 表示接收波束空间的维度。正参数β 限制了在Θ r 区域最大可接受的WHrb(φ) 和¯b(φ) 之间的偏差。同理,当N 是奇数时,Wr 可以用式(9)方法求得。在上述优化问题中,
Θ t 和Θ r 分别表示发射和接收空域滤波器的带宽,其大小可以按照需求自由设置,只要通过对优化问题求解就可以得到以Θ t 和Θ r 为带宽的发射和接收空域滤波器。因而相对于传统DFT空域滤波器,所提方法设计的空域滤波器带宽控制更加灵活。在得到
Wt 和Wr 后,发射和接收波束矩阵分别定义为¯Wt=WtU¯M (11) ¯Wr=WrU¯N (12) 式中,
U 为酉矩阵,其定义在文献[13]中给出。发射和接收空域滤波器响应可以分别表示为
Pt(θ)=||¯WHta(θ)||2/||a(θ)||2 (13) Pr(φ)=||¯WHrb(φ)||2/||b(φ)||2 (14) 发射波束空间的维度
¯M 要根据发射空域滤波器的带宽Θ t 来设定。为此,定义一个正定矩阵Ξ =∫Θ ta(θ)aH(θ)dθ (15) 本文要求矩阵
Ξ 的¯M 个最大特征值之和应该超过所有特征值之和一定的比例,本文中这个标准设定为99%。在文献[19]中这种方法被称为基于离散长球序列的方法,其基本原理是最大化滤波器主瓣能量相对于总能量的比例。同理,可以用相同的方法确定接收波束空间的维度。3.2 DOD和DOA估计
为了降低运算量和提高测角精度,本文对接收信号进行实值处理。首先,采用前后平滑技术构造一个中心厄米特矩阵[10]
¯Y=[Y,Π MNY∗Π Q] (16) 可以发现,矩阵
¯Y 满足条件¯Y∗=Π MN¯YΠ 2Q 。接下来,利用上一小节得到的发射和接收波束矩阵将¯Y 从阵元域转换到波束域。然而,由于Wt 和Wr 与DFT矩阵的结构不同,对¯Y 滤波后的信号结构也因此与文献[13]不同。这里需要强调的是,转换后的矩阵必须是实值的才能进行后续的实值处理。接下来,将会利用Wt 和Wr 的矩阵结构信息证明波束域转换后的矩阵是实值的。定理 1 如果矩阵
Wt 和Wr 满足式(8)—式(10)的等式约束条件,则¯M¯N×2Q 维的矩阵Z=(¯WHt⊗¯WHr)¯YU2Q (17) 是实值的。其中,
¯Wt 和¯Wr 的形式分别由式(11)和式(12)定义。证明 要证明矩阵
Z 是一个实值矩阵,只需证明Z∗=Z 。将等式¯Y∗=Π MN¯YΠ 2Q 代入到Z∗ 的表达式中得到Z∗=(¯WTt⊗¯WTr)¯Y∗U∗2Q=(¯WTt⊗¯WTr)Π MN¯YΠ 2QU∗2Q=(UT¯MWTt⊗UT¯NWTr)(Π M⊗Π N)¯YU2Q=(UT¯MWTtΠ M⊗UT¯NWTrΠ N)¯YU2Q (18) 在式(18)的推导过程中,利用
Π MN=Π M⊗ Π N 和Π 2QU∗2Q=U2Q 的等式条件。如果M 是偶数,由式(8)可知Wt=[Γ Tt,Π ¯MΓ HtΠ M/2]T 。如果M 为奇数,根据式(9)Wt=[Γ Tt,gt,Π ¯MΓ HtΠ (M−1)/2]T 。此时,g∗t=Π ¯Mgt 。不失一般性,以M 为奇数为例,将Wt 表达式代入式(14)UT¯MWTtΠ M=UT¯M[Π ¯MΓ Ht,gt,Γ TtΠ (M−1)/2]=UT¯MΠ ¯M[Γ Ht,Π ¯Mgt,Π ¯MΓ TtΠ (M−1)/2]=UT¯MΠ ¯M[Γ Ht,g∗t,Π ¯MΓ TtΠ (M−1)/2]=UH¯MWHt (19) 同样,用相同的方法可以证明
UT¯NWTrΠ N=UH¯NWHr (20) 基于式(19)和式(20),式(18)可以重写为
Z∗=(UT¯MWTtΠ M⊗UT¯NWTrΠ N)¯YU2Q=(UH¯MWHt⊗UH¯NWHr)¯YU2Q=(¯WHt⊗¯WHr)¯YU2Q=Z (21) 证毕
从上面的证明过程可以看出,如果不对
Wt 和Wr 的矩阵结构进行事先设计就无法构造实值信号模型。为了从实值矩阵Z 中提取角度信息,需要对矩阵Z 的结构有一个更直观的认识。首先,将式(6)代入式(16),矩阵¯Y 重写为¯Y=[Y,Π MNY∗Π Q]=[(A⊙B)C,(Π MA∗⊙Π NB∗)C∗Π Q]+⌢N=(A⊙B)[C,C∗Π Q]+⌢N (22) 式中,
⌢N=[˜N,Π MN˜N∗Π Q] 。由于导向矩阵A 和B 中心共轭对称,存在A=Π MA∗ 和B=Π NB∗ 。将式(22)代入到式(17)中去,得到Z=(¯WHt⊗¯WHr)¯YU2Q=(UH¯MWHtA⊙UH¯NWHrB)[C,C∗Π Q]U2Q+ˆN (23) 式中
ˆN=(¯WHt⊗¯WHr)⌢NU2Q 表示噪声项。如果不考虑求解式(8)—式(10)时插值误差的影响,则存在¯A=WHtA 和¯B=WHrB 。在上面的等式中,¯A=[¯a(θ1),¯a(θ2),···,¯a(θK)] 表示虚拟发射导向矩阵,¯B=[¯b(φ1),¯b(φ2),···,¯b(φK)] 表示虚拟接收导向矩阵。此时,式(23)可以重写为Z=(UH¯M¯A⊙UH¯N¯B)[C,C∗Π Q]U2Q+ˆN (24) 可以发现,此时
UH¯M¯A 和UH¯N¯B 分别包含K 个目标的DOD和DOA信息。构造Z 的自相关矩阵R=(1/2Q)ZZH ,并对其进行特征分解得到信号子空间Es 。Es 的维度为¯M¯N×K ,由矩阵R 最大的K 个特征值对应的特征矢量构成。另外,Es 和UH¯M¯A⊙UH¯N¯B 张成的子空间相同,可以表示为Es=(UH¯M¯A⊙UH¯N¯B)T ,其中,T 是一个K×K 维的非奇异矩阵。已知Es=(UH¯M¯A⊙UH¯N¯B)T ,利用阵列的旋转不变特性,式(25)和式(26)关系成立(Kt,1⊗I¯N)EsΨ t=(Kt,2⊗I¯N)Es (25) (I¯M⊗Kr,1)EsΨ r=(I¯M⊗Kr,2)Es (26) 式中,
Kt,1=Re{UH¯M−1JtU¯M} ;Kt,2=Im{UH¯M−1 JtU¯M} ;Kr,1=Re{UH¯N−1JrU¯N} ;Kr,2=Im{UH¯N−1 JrU¯N} 。在上面的等式关系中,Jt=[0,I¯M−1] 和Jr=[0,I¯N−1] 分别是(¯M−1)ׯM 和(¯N−1)ׯN 维的选择矩阵。式(25)和式(26)中,Ψ t=T−1Δ tT 以及Ψ r=T−1Δ rT 。Δ t 和Δ r 具有式(27)和式(28)的形式Δ t=diag{tan(πsinθ1/2),tan(πsinθ2/2),···,tan(πsinθK/2)} (27) Δ r=diag{tan(πsinφ1/2),tan(πsinφ2/2),···,tan(πsinφK/2)} (28) 式(25)和式(26)中的
Ψ t 和Ψ r 可以用最小二乘方法得到,之后对Ψ t 和Ψ r 进行特征分解,其第k 个特征值分别用γk 和μk 表示。本文采用同时Schur分解方法[20]进行角度配对,由于不是本文研究重点,这里不再赘述。此时,第k 个目标的DOD和DOA估计值分别为ˆθk=arcsin[2arctan(γk)/π] (29) ˆφk=arcsin[2arctan(μk)/π] (30) 3.3 插值误差补偿
在上一小节中,目标的角度值已经在式(24)成立的基础上估计出来。然而,式(24)是在不考虑插值误差的情况下才成立的。实际情况中,
WHtA 和WHrB 不可能完全等于¯A 和¯B ,因此若直接将式(29)和式(30)作为最终的角度估计值,测角精度必然会受到插值误差的影响,精度下降。因此,必须采取措施对插值误差进行补偿。首先,建立式(23)的无噪声信号模型˜Z=(¯WHtA⊙¯WHrB)[C,C∗Π Q]U2Q (31) 可以发现,基于此模型的角度估计精度仅受插值误差的影响。观察式(31)的结构,根据式(32)和式(33)分别确定DOD和DOA的映射关系
ηiKt,1[¯WHta(θi)]=Kt,2[¯WHta(θi)] (32) ζiKr,1[¯WHrb(φi)]=Kr,2[¯WHrb(φi)] (33) 式中,
θi 和φi 分别表示在Θ t 和Θ r 区域采样的真实角度值。ηi 和ζi 的表达式分别为ηi=tan(πsin¯θi/2) 和ζi=tan(πsin¯φi/2) 。¯θi 和¯φi 分别表示θi 和φi 的映射值。值得注意是,采样间隔要尽可能的小以满足雷达对测角精度的要求。θi 和φi 对应的映射DOD和DOA表示为¯θi=arcsin[2arctan(ηi)/π] 和¯φi= arcsin[2arctan(ζi)/π] 。据此,θi 和¯θi 的映射关系可以一一确定,同样φi 和¯φi 的映射关系也可以得到确定。当目标的角度估计值在有噪声模信号模型下根据式(29)和式(30)得到后,更精确的角度估计值可以根据映射操作得到。4. 仿真实验
假设双基地MIMO雷达的发射和接收阵列均为半波长间隔等距线阵。发射阵元个数
M=10 ,接收阵元数N=8 。定义均方根误差(Root Mean Square Error, RMSE)RMSE=√12PKK∑k=1P∑p=1[(ˆθk,p−θk)2+(ˆφk,p−φk)2] (34)
其中,
P=400 表示总的蒙特卡洛实验次数。ˆθk,p 和ˆφk,p 分别表示第p 次蒙特卡洛实验第k 个目标的DOD和DOA估计值(映射后)。在所有的仿真实验中,感兴趣的空间区域Θ t=[10∘,40∘] ,Θ r=[0∘,30∘] ,参数α=0.04 ,β=0.05 ,发射和接收波束矩阵通过求解式(8)—式(10)得到。图1(a)和图1(b)分别为发射和接收空域滤波器的归一化幅频响应,由式(13)和式(14)给出。文献[13]中的DFT空域滤波器被用于比较。发射和接收波束空间的维度分别为¯M=8 和¯N=6 。从图1(a)和图1(b)中可以看出,本文方法设计的发射和接收空域滤波器仅在感兴趣空间区域有高增益,其主副瓣比达到20 dB,而传统的DFT空域滤波器主副瓣比只有8 dB。相比于传统的DFT空域滤波器,本文方法设计的空域滤波器能够按需求灵活设计带宽,且旁瓣电平能够被有效抑制,会极大提高低信噪比下的测角性能。实验 1 插值误差补偿的必要性验证
假设3个不相关信源,分别处于
(θ1,φ1)= (15∘,5∘) ,(θ2,φ2)=(25∘,10∘) ,(θ3,φ3)=(35∘,25∘) ,信噪比均为SNR=5 dB,脉冲数Q=50 。图2(a)给出本文算法映射前的角度估计结果。从图2(a)的仿真结果中可以看出,测角结果在不进行误差补偿的情况下,测角精度严重降低。这种损失不仅受噪声影响,而且也受到插值误差的影响。进行误差补偿后,图2(b)给出本文算法映射后的角度估计结果。经过统计,插值误差补偿前的RMSE为0.24∘ ,补偿后为0.09∘ 。可以看出,插值误差补偿后测角精度得到明显提高。实验 2 非相关信源的测角性能估计
假设两个非相关信源分别位于
(θ1,φ1)=(15∘,5∘) 和(θ2,φ2)=(35∘,25∘) 。本文算法与PM算法,ESPRIT算法,酉ESPRIT算法,波束域ESPRIT算法,酉PARAFAC张量算法以及文献[21]中的克拉美罗界(Cramer-Rao Bound, CRB)进行比较。图3给出了RMSE随着信噪比的变化情况,此时的脉冲数固定为50。图4给出了RMSE随着脉冲数的变化情况,此时的信噪比固定为0 dB。可以看出,本文算法的测角性能明显优于其他算法。这种优异性在低信噪比情况下尤其明显。实验 3 相关信源的测角性能估计
两个相关信源分别位于
(θ1,φ1)=(15∘,5∘) 和(θ2,φ2)=(35∘,25∘) ,相关系数为0.99。图5给出了RMSE随着信噪比的变化情况,此时的脉冲数固定为50。图6给出检测概率随信噪比的变化情况,此时的脉冲数依然固定为50。当估计的目标角度值与真实目标角度值之差的绝对值小于等于0.1∘ 时则认为目标检测成功。可以看出,PM算法和ESPRIT算法在处理相关信源时存在较大的性能损失。由于文献[13]提出的传统波束域ESPRIT算法也采用前后平滑技术进行实值处理,因此对两个相关信源也可以高精度估计。本文算法的测角精度和检测概率要明显的高于其他算法。5. 结论
本文提出了一种改进的基于实值波束域的双基地MIMO雷达测角方法。本文算法适用于发射波形为全正交波形的双基地MIMO雷达,且要求目标所在的大体空间位置已知。本文算法的空域滤波器采用凸优化方法设计,能够灵活控制带宽和抑制旁瓣电平,达到提高主副瓣比的目的。发射和接收波束空间的维度要根据空域滤波器带宽来确定。同时,为了对滤波之后的信号进行实值处理,必须对发射和接收波束矩阵的结构进行设定。本文证明了这种结构是进行实值处理的必要条件。另外,本文算法在处理相关信源时也有良好的性能。最后,仿真实验验证了本文算法的有效性。
-
表 1 n=4,k=2时DM-MCM-PSWFs的一种映射方案
比特信号 信号索引 子载波映射 [0,0] {0,0,1,1} {sAI(1),sAI(2),sBI(1),sBI(2)} [0,1] {0,1,0,1} {sAI(1),sBI(1),sAI(2),sBI(2)} [1,0] {1,0,0,1} {sBI(1),sAI(1),sAI(2),sBI(2)} [1,1] {0,1,1,0} {sAI(1),sBI(1),sBI(2),sAI(2)} 表 2 不同多载波调制方法系统频带利用率
调制方法 g n k SE(bit/s/Hz) Eb/N0(dB) ρ(%) DM-MCM-PSWFs 15 6 3 3.09 11.98 / MCM-PSWFs-SGO-2PAM 9 10 7 2.41 11.05 28.2 MCM-PSWFs-SGO-4PAM 23 4 1 1.90 13.46 62.6 23 4 2 2.85 14.95 8.4 MCOM-PSWFs-8QAM 1 92 92 2.85 12.45 8.4 表 3 信号索引检测乘法运算量
调制方式 运算量 n k 乘法次数(B=1.44 MHz) DM-MCM-PSWFs-ML O(ng2⌊Ckn⌋) 4 1/2/3 368 MCM-PSWFs-SGO-ML O(2kg2⌊Ckn⌋) 4 1 184 4 2 368 4 3 552 MCM-PSWFs-SGO-OS O(gnlog2n) 4 2 16 -
[1] CAI Yunlong, QIN Zhijin, CUI Fangyu, et al. Modulation and multiple access for 5G networks[J]. IEEE Communications Surveys & Tutorials, 2018, 20(1): 629–646. [2] ZHANG Lin, LIANG Yingchang, and XIAO Ming. Spectrum sharing for internet of things: A survey[J]. IEEE Wireless Communications, 2019, 26(3): 132–139. doi: 10.1109/MWC.2018.1800259 [3] HOU Zhanwei, SHE Changyang, LI Yonghui, et al. Prediction and communication co-design for ultra-reliable and low-latency communications[J]. IEEE Transactions on Wireless Communications, 2020, 19(2): 1196–1209. doi: 10.1109/TWC.2019.2951660 [4] BUZZI S, D’ANDREA C, LI Dejian, et al. MIMO-UFMC transceiver schemes for millimeter-wave wireless communications[J]. IEEE Transactions on Communications, 2019, 67(5): 3323–3336. doi: 10.1109/TCOMM.2019.2896122 [5] MICHAILOW N, MATTHÉ M, GASPAR I S, et al. Generalized frequency division multiplexing for 5th generation cellular network[J]. IEEE Transactions on Communications, 2014, 62(9): 3045–3061. doi: 10.1109/TCOMM.2014.2345566 [6] WANG Hongxing, LU Faping, LIU Chuanhui, et al. Frequency domain multi-carrier modulation based on prolate spheroidal wave functions[J]. IEEE Access, 2020, 8: 99665–99680. doi: 10.1109/ACCESS.2020.2997679 [7] 陆发平, 王红星, 刘传辉, 等. 椭圆球面函数频域调制解调方法[J]. 电子与信息学报, 2020, 42(8): 1888–1895. doi: 10.11999/JEIT190642LU Faping, WANG Hongxing, LIU Chuanhui, et al. PSWFs Frequency Domain Modulation and Demodulation Method[J]. Journal of Electronics &Information Technology, 2020, 42(8): 1888–1895. doi: 10.11999/JEIT190642 [8] KHARE K and GEORGE N. Sampling theory approach to prolate spheroidal wavefunctions[J]. Journal of Physics A: Mathematical and General, 2003, 36(39): 10011–10021. doi: 10.1088/0305-4470/36/39/303 [9] 王红星, 陆发平, 刘传辉, 等. 严格奇偶对称的椭圆球面波函数信号构建与低复杂度检测方法[J]. 中国科学: 信息科学, 2020, 50(5): 766–776. doi: 10.1360/SSI-2019-0121WANG Hongxing, LU Faping, LIU Chuanhui, et al. Strict parity symmetric prolate spheroidal wave functions signal construction and low complexity detection method[J]. Scientia Sinica Informationis, 2020, 50(5): 766–776. doi: 10.1360/SSI-2019-0121 [10] CAI Biao, LIU Aijun, and LIANG Xiaohu. Low-complexity selective mapping methods for multicarrier faster-than-nyquist signaling[J]. IEEE Access, 2020, 8: 31420–31431. doi: 10.1109/ACCESS.2020.2973394 [11] HU Yue, WANG Yafeng, and WANG Haocheng. A decoding method based on RNN for OvTDM[J]. China Communications, 2020, 17(4): 1–10. [12] LI Qiang, WEN Miaowen, CLERCKX B, et al. Subcarrier index modulation for future wireless networks: Principles, applications, and challenges[J]. IEEE Wireless Communications, 2020, 27(3): 118–125. doi: 10.1109/MWC.001.1900335 [13] DANG Shuping, GUO Shuaishuai, COON J P, et al. Enhanced Huffman coded OFDM with index modulation[J]. IEEE Transactions on Wireless Communications, 2020, 19(4): 2489–2503. doi: 10.1109/TWC.2020.2965524 [14] LI Qiang, WEN Miaowen, BASAR E, et al. Index modulated OFDM spread spectrum[J]. IEEE Transactions on Wireless Communications, 2018, 17(4): 2360–2374. doi: 10.1109/TWC.2018.2793238 [15] 王红星, 陆发平, 刘传辉, 等. 基于信号分组优化的椭圆球面波多载波调制解调方法[J]. 中国科学: 信息科学, 2021, 51(7): 1168–1182.WANG Hongxing, LU Faping, LIU Chuanhui, et al. Multi-carrier modulation scheme based on prolate spheroidal wave functions with signal grouping optimization[J]. Scientia Sinence Informationis, 2021, 51(7): 1168–1182. [16] MAO Tianqi, JIANG Rui, and BAI Ruowen. Optical dual-mode index modulation aided OFDM for visible light communications[J]. Optics Communications, 2017, 391: 37–41. doi: 10.1016/j.optcom.2017.01.005 [17] 陆发平, 刘传辉, 康家方. 基于顺序统计量的索引调制信号检测方法[J]. 无线电通信技术, 2019, 45(1): 57–61. doi: 10.3969/j.issn.1003-3114.2019.01.11LU Faping, LIU Chuanhui, and KANG Jiafang. Index modulation signal detection method based on order statistic[J]. Radio Communications Technology, 2019, 45(1): 57–61. doi: 10.3969/j.issn.1003-3114.2019.01.11 期刊类型引用(7)
1. 黄中瑞,史英春,唐波,秦立龙. 阵元等功率约束下的MIMO雷达发射加权矩阵优化算法. 电子与信息学报. 2022(05): 1856-1864 . 本站查看
2. 李洪兵,欧阳婉,杨明磊,叶舟,陈业伟,陈伯孝. L型分布式面阵二维高精度角度估计方法研究. 火控雷达技术. 2022(04): 24-31 . 百度学术
3. 谢前朋,潘小义,陈吉源,肖顺平. 基于新型阵列的双基地电磁矢量传感器MIMO雷达高分辨角度参数估计. 电子与信息学报. 2021(02): 270-276 . 本站查看
4. 师俊朋,文方青,艾林,张弓,龚政辉. 空域色噪声背景下双基地MIMO雷达角度估计. 系统工程与电子技术. 2021(06): 1477-1485 . 百度学术
5. 郑志东. 互耦条件下基于4阶累积量的双基地MIMO雷达角度估计. 电子与信息学报. 2021(09): 2720-2727 . 本站查看
6. 刘东贺,赵永波,庞晓娇,曹成虎,陈胜. 基于降维波束空间的实值ESPRIT单基地MIMO雷达测角算法. 电子与信息学报. 2021(12): 3639-3646 . 本站查看
7. 徐保庆,剡熠琛,同非,任哲毅. 基于波束域实值HOSVD的双基地MIMO雷达测角算法. 火控雷达技术. 2021(04): 6-12 . 百度学术
其他类型引用(3)
-