基于双共轭梯度法的频控阵多波束方向调制方法

1.本发明属于天线技术领域，具体涉及一种基于双共轭梯度法的频控阵多波束方向调制方法。


背景技术：

2.方向调制(directional modulation,dm)技术作为一种具有巨大应用价值的物理层安全传输策略，已经受到了学术界的广泛关注。dm技术可以在自由空间沿着特定的方向为目标用户生成一个标准的符号星座，同时扭曲窃听者的接收信号。
3.目前，dm的实现方案主要包括在射频端通过优化的方式和在基带通过人工噪声辅助的方式。对于两种方式，均可使用相控阵(phased array,pa)或者频控阵(frequency diverse array,fda)来实现，而这两种阵列的主要区别在于，基于pa的dm方案只能实现角度维的物理层安全传输，而基于fda的dm方案可以同时实现角度维和距离维的物理层安全传输。
4.基于射频端优化方式的dm技术的基本思路是通过优化阵列天线的权值，从而在目标用户处产生标准的符号星座，同时扭曲窃听者位置处的接收信号。然而，由于该方法的计算复杂度一般比较高，因此在基带基于人工噪声(artificial noise,an)辅助的dm技术被提出，此方案的基本思路是对阵列发射的信号添加正交于用户信道的人工噪声，从而恶化窃听者的信道。但是，此方法需要计算复矩阵的伪逆，并且根据定义直接计算伪逆的计算复杂度仍较大，为用户数的三次方，因此本发明引入双共轭梯度法对矩阵的伪逆进行计算，其基本思路是将矩阵的伪逆计算问题转化成线性方程组的求解问题，然后通过双共轭梯度法对其迭代求解，从而降低计算复杂度。


技术实现要素：

5.针对现有基于人工噪声辅助的方向调制方法需要计算复矩阵的伪逆，存在根据定义直接计算伪逆的复杂度较高的缺陷，本发明提供了一种基于双共轭梯度法的频控阵多波束方向调制方法。本发明方法引入双共轭梯度法对矩阵的伪逆进行计算，将矩阵的伪逆计算问题转化成线性方程组的求解问题，并通过双共轭梯度法对其迭代求解，实现了对频控阵多波束的方向调制，并极大的降低了计算复杂度。
6.本发明采用的技术方案如下：
7.一种基于双共轭梯度法的频控阵多波束方向调制方法，包括以下步骤：
8.s1:根据fda发射机与接收机结构，基于人工噪声辅助的方向调制方法建立信号传输模型。
9.s2:将矩阵的伪逆计算问题转化成线性方程组的求解问题。
10.s3:通过双共轭梯度法对线性方程组迭代求解，实现对频控阵多波束的方向调制。
11.各步骤具体处理流程如下：
12.s1:根据fda发射机与接收机结构，基于人工噪声辅助的方向调制方法建立信号传
输模型，其阵列加权矢量形式为：
[0013][0014]
其中，ps表示总的发射功率，p1表示归一化矩阵，p2表示正交矩阵，β1表示信号的功率分配因子，β2表示人工噪声的功率分配因子，且满足表示人工噪声的功率分配因子，且满足表示人工噪声归一化因子，tr(
·
)表示矩阵的迹，z表示人工噪声矢量，且满足in表示n阶单位矩阵，s＝[s1,s2,...,sk]
t
表示发送给k个用户的保密信息，s1,s2,...,sk分别表示发送给第1,2,...,k个用户的保密信息。
[0015]
为了在目标用户处收到正确的信号，则矩阵p1与矩阵p2满足：
[0016]hh
(θd,rd)p1＝ikꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0017]hh
(θd,rd)p2＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0018]
其中(
·
)h表示矩阵的共轭转置，h(θd,rd)表示k个用户的导向矢量矩阵，ik表示k阶单位矩阵。
[0019]
由式(2)与式(3)可得：
[0020]
p1＝(hh(θd,rd))

ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0021]
p2＝i
n-(hh(θd,rd))
hh
(θd,rd)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0022]
其中(
·
)

表示矩阵的伪逆，即：
[0023]
(hh(θd,rd))

＝h(θd,rd)[hh(θd,rd)h(θd,rd)]-1
ꢀꢀ
(6)
[0024]
其中(
·
)-1
表示矩阵的逆。
[0025]
s2:将矩阵的伪逆计算问题转化成线性方程组的求解问题。
[0026]
由于(hh(θd,rd))

＝(h

(θd,rd))h，则将式(6)转化为：
[0027]h
(θd,rd)＝(hh(θd,rd)h(θd,rd))-1
h(θd,rd)hꢀꢀ
(7)
[0028]
令gt、qt分别表示h

(θd,rd)与hh(θd,rd)的第t列数据，1≤t≤n，则计算g
t
＝(hh(θd,rd)h(θd,rd))-1qt
，等价于求解方程组：
[0029]
h(θd,rd)hh(θd,rd)g
t
＝q
t
ꢀꢀ
(8)
[0030]
s3:通过双共轭梯度法对线性方程组迭代求解。
[0031]
令a＝hh(θd,rd)h(θd,rd)，x＝g
t
，b＝q
t
，则式(8)转化为ax＝b；
[0032]
双共轭梯度法的计算过程如下：
[0033]
step 1:取x的初始值x0，x0通常取零向量或者单位向量，容差δ，通常满足δ≤10-6
，以及最大迭代次数w，通常满足1≤w≤50；
[0034]
step 2:令r0＝b-ax0，取p0＝r0，其中r0、p0、均为中间变量；
[0035]
step 3:计算：
[0036]
x
i 1
＝xi αipi，r
i 1
＝r
i-αiapi，
[0037]
p
i 1
＝r
i 1
βipi，
[0038]
其中，a
*
表示a的共轭矩阵，表示ri的共轭向量，和分别表示αi与βi的共轭复数，(
·
，
·
)表示复内积，αi、βi、ri、pi、均为中间变量，下标i＝0,1,2,...表示迭代次数。
[0039]
step 4:如果||ri||/||b||＜δ或者i＝w，则输出xi作为h

(θd,rd)的第t列数据，结束；否则令i＝i 1，转step 3；其中||
·
||表示向量的二范数。
[0040]
通过n次迭代计算即求得h

(θd,rd)，最后对h

(θd,rd)共轭转置，计算得到(hh(θd,rd))

。
[0041]
由于基于人工噪声辅助的方向调制方法需要计算复矩阵的伪逆，其单列数据的计算复杂度为o(k3)，因此本发明将伪逆的计算问题转换成线性方程组的求解问题，然后通过引入双共轭梯度法对其进行迭代计算。在本发明方法中，进行一次迭代过程中的复数乘法量为2k2 8k 7，复数加法量为2k2 7k-4，则进行一次迭代的复杂度为o(k2)，那么进行u次迭代的复杂度为o(uk2)，从而大大降低了单列数据的计算复杂度。
附图说明
[0042]
图1为本发明的一种基于双共轭梯度法的频控阵多波束方向调制方法的流程示意图；
[0043]
图2为基于频控阵的发射机；
[0044]
图3为用户位置处的接收机；
[0045]
图4为误码率与信噪比的关系图；
[0046]
图5为当角度取用户所在位置时，误码率与空间距离的关系图；
[0047]
图6为当距离取用户所在位置时，误码率与空间角度的关系图。
具体实施方式
[0048]
下面根据附图和实施例详细阐述本发明的技术方案，此处所选的实施例仅用于解释该发明，并不能够限定此发明。
[0049]
s1:根据fda发射机与接收机结构，基于人工噪声辅助的方向调制方法建立信号传输模型。
[0050]
设频控阵发射机有n个阵元，阵元间距为d，其结构如图2所示，第n个阵元的发射信号为：
[0051][0052]
其中表示第n个阵元的加权系数，an与φn均是实数，分别表示幅度与相位加权，fn＝f0 δfn表示发射频率，f0表示发射的基本载波频率，δfn表示第n个阵元的频偏，n＝1,2,..n，是复基带信号，t表示时间，满足：
[0053][0054]
t表示波形持续时间，(
·
)
*
表示共轭算子。
[0055]
设接收机在空间远场区用户位置为(r,θ)，r表示距离，θ表示角度，其结构如图3所示，则接收到的信号为：
[0056][0057]
其中，rn＝r-(n-1)dsinθ表示第n个阵元到用户位置的距离，r表示第1个阵元到用户位置的距离，c表示光速。
[0058]
假设是窄带信号，则有(3)式改写为：
[0059][0060]
接收机对接收到的信号首先进行相位加权并分成m路信号分别进行下变频，其中m＝n,然后与冲击响应进行匹配滤波，最后将m路信号进行求和得到信号x。
[0061]
设满足如下正交关系：
[0062][0063]
其中，τ表示任意时延，δfm表示第m路信号的频偏，m＝1,2,...,m，δ(
·
)表示kronecker delta函数，满足：
[0064][0065]
则经过匹配滤波后的第m路信号xm为：
[0066][0067]
那么m路信号求和得到的信号x为：
[0068][0069]
则在用户位置(r,θ)处的归一化导向矢量为：
[0070][0071]
其中(
·
)
t
表示矩阵转置，那么k个用户的
导向矢量矩阵为：
[0072]
h(θd,rd)＝[h(θ1,r1),h(θ2,r2),...h(θk,rk)...,h(θk,rk)]
ꢀꢀ
(10)
[0073]
其中(θd,rd)＝[(θ1,r1),(θ2,r2),...(θk,rk),...,(θk,rk)]，(θk,rk)表示第k个用户的位置，k＝1,2,...,k。
[0074]
假设s＝[s1,s2,...,sk]
t
是发送给k个用户的保密信息，则fda的加权矢量为：
[0075][0076]
其中，ps表示总的发射功率，p1表示归一化矩阵，p2表示正交矩阵，β1表示信号的功率分配因子，β2表示人工噪声的功率分配因子，且满足表示人工噪声的功率分配因子，且满足表示人工噪声归一化因子，tr(
·
)表示矩阵的迹，z表示人工噪声矢量，且满足in表示n阶单位矩阵。
[0077]
为了在目标用户处收到正确的信号，则矩阵p1与矩阵p2满足：
[0078]hh
(θd,rd)p1＝ikꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0079]hh
(θd,rd)p2＝0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0080]
其中(
·
)h表示矩阵的共轭转置，ik表示k阶单位矩阵。
[0081]
由式(12)与式(13)得：
[0082]
p1＝(hh(θd,rd))

ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0083]
p2＝i
n-(h(θd,rd)h)
hh
(θd,rd)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(15)
[0084]
其中(
·
)

表示矩阵的伪逆，即：
[0085]
(hh(θd,rd))

＝h(θd,rd)[hh(θd,rd)h(θd,rd)]-1
ꢀꢀ
(16)
[0086]
(
·
)-1
表示矩阵的逆，则k个用户的接收信号矢量r为：
[0087][0088]
其中表示复噪声矢量，是目标用户的信道噪声功率。
[0089]
同理，对于窃听用户位置(θe,re)≠(θk,rk)收到的信号为：
[0090][0091]
其中表示复高斯噪声，是窃听用户的信道噪声功率。
[0092]
由式(17)与式(18)可知，用户位置收到的信号仅受到信道噪声的影响，而窃听者位置收到的信号不仅受到信道噪声的影响，同时受到人工噪声的影响。
[0093]
s2:将矩阵的伪逆计算问题转化成线性方程组的求解问题。
[0094]
由式(14)与式(15)可知，在计算p1与p2之前，需要计算式(16)，由矩阵的性质(hh(θd,rd))

＝(h

(θd,rd))h，式(16)可转化为计算：
[0095]h
(θd,rd)＝(hh(θd,rd)h(θd,rd))-1
h(θd,rd)hꢀꢀ
(19)
[0096]
令g
t
，q
t
分别是h

(θd,rd)与hh(θd,rd)的第t列数据，1≤t≤n，则计算g
t
＝(hh(θd,rd)h(θd,rd))-1qt
，等价于求解方程组
[0097]
h(θd,rd)hh(θd,rd)g
t
＝q
t
ꢀꢀ
(20)
[0098]
s3:通过双共轭梯度法对线性方程组迭代求解，实现对频控阵多波束的方向调制。
[0099]
令a＝hh(θd,rd)h(θd,rd)，x＝g
t
，b＝q
t
，则(20)式可转化为ax＝b；
[0100]
双共轭梯度法的计算过程如下：
[0101]
step 1:取x的初始值x0，x0通常取零向量或者单位向量，容差δ，通常满足δ≤10-6
，以及最大迭代次数w，通常满足1≤w≤50；
[0102]
step 2:令r0＝b-ax0，取p0＝r0，其中r0、p0、均为中间变量；
[0103]
step 3:对i＝0,1,2,...，计算：
[0104]
x
i 1
＝xi αipi，r
i 1
＝r
i-αiapi，
[0105]
p
i 1
＝r
i 1
βipi，
[0106]
其中，a
*
表示a的共轭矩阵，r
*
表示r的共轭向量，α
*
和β
*
分别表示α与β的共轭复数，(
·
，
·
)表示复内积，αi、βi、ri、pi、均为中间变量，下标i＝0,1,2,...表示迭代次数。
[0107]
step 4:如果||ri||/||b||＜δ或者i＝w，则输出xi作为h

(θd,rd)的第t列数据，结束；否则令i＝i 1，转step 3；其中||
·
||表示向量的二范数。
[0108]
通过n次迭代计算即可求得h

(θd,rd)，最后对h

(θd,rd)共轭转置，计算得到(hh(θd,rd))

。
[0109]
由于基于人工噪声辅助的方向调制方法需要计算复矩阵的伪逆，其单列数据的计算复杂度为o(k3)，则本发明将伪逆的计算问题转换成线性方程组的求解问题，然后通过引入双共轭梯度法对其进行迭代计算。在本发明的方法中，进行一次迭代过程中的复数乘法量为2k2 8k 7，复数加法量为2k2 7k-4，则进行一次迭代的复杂度为o(k2)，那么进行u次迭代的复杂度为o(uk2)，从而降低单列数据的计算复杂度。
[0110]
实施例
[0111]
对于基于双共轭梯度法的频控阵多波束方向调制方法的实验验证的参数设置为f0＝10ghz，d＝c/2f0，n＝512，δfn＝ηnδf，ηn是随机整数且满足0≤ηn≤2n-1，并且η1,η2,...,ηn互不相同，δf＝1khz，ps＝1，k个用户的保密信息s随机生成，调制方式为qpsk，用户个数k＝4，其位置分别为(-40
°
,60km)，(-10
°
,110km)，(20
°
,140km)，(70
°
,210km)。
[0112]
图4表示直接求解法与基于本发明不同迭代次数仿真的用户误码率随信噪比的变化情况，从中可以发现，当迭代次数为二次时，本发明提出的方法基本与直接求逆法的误码性能相同，从而证明了本发明的可行性。图5表示当空间角度分别取用户所在的位置时，误码率随空间距离变化情况，图6表示当空间距离取目标用户所在的位置时，误码率随空间角度的变化情况，结合图5与图6可以看到，在空间中，只有用户位置处的误码率才存在较低的
值，其他位置处的误码率均较高。呈现图5与图6现象的原因是，在用户处的人工噪声可以被完全消除，而窃听者却无法消除人工噪声，这就导致了窃听者位置处接收到的保密信号的质量非常低，因此会呈现比较高的误码率，并且结合图5与图6可知，本发明经过一次迭代后的误码性能可以有效逼近直接求逆法的误码性能。
[0113]
综上所述，本发明提出的基于双共轭梯度法的频控阵多波束方向调制方法能够在降低单列数据计算复杂度的情况下，实现信号调制的方向性。