雷达通信一体化波形设计常见优化指标及约束条件

参考的主要是刘凡大佬等人一系列的文章。

MUI干扰

多用户干扰指的是发射机给不同的用户发送不同的信息,发出去的信号是多个用户信息的叠加,我们希望每个用户只接受到自己的信息,所以我们会对发送信息S作一个类似预编码的处理,变成X,而我们对CSI已知,所以如下所示,
Y = S + ( H X − S ) ⏟ MUI  + W , \mathbf{Y}=\mathbf{S}+\underbrace{(\mathbf{H X}-\mathbf{S})}_{\text {MUI }}+\mathbf{W}, Y=S+MUI  (HXS)+W,
min ⁡ X ∥ H X − S ∥ F 2 \begin{aligned} &\min _{\mathbf{X}}\|\mathbf{H X}-\mathbf{S}\|_{F}^{2} \\ \end{aligned} XminHXSF2

雷达信噪比

雷达系统和通信系统的一个区别就是雷达杂波,一般雷达把不感兴趣的目标都认为是杂波,需要抑制,而通信对所有环境杂波都感兴趣,这也是我觉得雷达通信相互有些冲突的地方,如何把雷达探测和通信的信道估计结合,相互促进性能是一个可以解决的问题,就这点我感觉OTFS波形会比较好。
回到这里,我们假设雷达感兴趣的目标是位于 θ 0 \theta_0 θ0处的位置,那么其他方向的都是杂波,这里又应用了一个先验知识,就是雷达知道目标所在位置,这些约束都还挺强的。
q = α 0 A ( θ 0 ) x + ∑ k = 1 K α k A ( θ k ) x + u \mathbf{q}=\alpha_{0} \mathbf{A}\left(\theta_{0}\right) \mathbf{x}+\sum_{k=1}^{K} \alpha_{k} \mathbf{A}\left(\theta_{k}\right) \mathbf{x}+\mathbf{u} q=α0A(θ0)x+k=1KαkA(θk)x+u
其中 A ( θ ) \mathbf{A}(\theta) A(θ) 是发射天线和接收天线的导向矢量。
A ( θ ) = I N ⊗ [ a r ( θ ) a t T ( θ ) ] \mathbf{A}(\theta)=\mathbf{I}_{N} \otimes\left[\mathbf{a}_{r}(\theta) \mathbf{a}_{t}^{T}(\theta)\right] A(θ)=IN[ar(θ)atT(θ)]
同时可以将MIMO雷达处理技术结合进来,比如接收天线处再加个滤波器系数来实现噪声、干扰和干扰抑制,同时保持所需信号不失真 。
以下的文章这样做过,
Joint Transmit Waveform and Receive Filter Design for Dual-Function Radar-Communication Systems
c = w H q = α 0 w H A ( θ 0 ) x + w H ∑ k = 1 K α k A ( θ k ) x + w H u \begin{aligned} c &=\mathbf{w}^{H} \mathbf{q} \\ &=\alpha_{0} \mathbf{w}^{H} \mathbf{A}\left(\theta_{0}\right) \mathbf{x}+\mathbf{w}^{H} \sum_{k=1}^{K} \alpha_{k} \mathbf{A}\left(\theta_{k}\right) \mathbf{x}+\mathbf{w}^{H} \mathbf{u} \end{aligned} c=wHq=α0wHA(θ0)x+wHk=1KαkA(θk)x+wHu
输出信噪比如下
γ ( x , w ) = σ 0 2 ∣ w H A ( θ 0 ) x ∣ 2 w H [ ∑ k = 1 K σ k 2 A ( θ k ) x x H A H ( θ k ) ] w + σ u 2 w H w \gamma(\mathbf{x}, \mathbf{w})=\frac{\sigma_{0}^{2}\left|\mathbf{w}^{H} \mathbf{A}\left(\theta_{0}\right) \mathbf{x}\right|^{2}}{\mathbf{w}^{H}\left[\sum_{k=1}^{K} \sigma_{k}^{2} \mathbf{A}\left(\theta_{k}\right) \mathbf{x} \mathbf{x}^{H} \mathbf{A}^{H}\left(\theta_{k}\right)\right] \mathbf{w}+\sigma_{u}^{2} \mathbf{w}^{H} \mathbf{w}} γ(x,w)=wH[k=1Kσk2A(θk)xxHAH(θk)]w+σu2wHwσ02wHA(θ0)x2

刘凡还有一个不同的优化式子,直接去设杂波目标离感兴趣目标间隔P个距离门,
Y R = α 0 a ( θ 0 ) a H ( θ 0 ) X + ∑ p = − P , p ≠ 0 P ∑ m = 1 M α p , m a ( θ m ) a H ( θ m ) X J p + W R , \begin{aligned} \mathbf{Y}_{R} &=\alpha_{0} \mathbf{a}\left(\theta_{0}\right) \mathbf{a}^{H}\left(\theta_{0}\right) \mathbf{X} \\ &+\sum_{p=-P, p \neq 0}^{P} \sum_{m=1}^{M} \alpha_{p, m} \mathbf{a}\left(\theta_{m}\right) \mathbf{a}^{H}\left(\theta_{m}\right) \mathbf{X} \mathbf{J}_{p}+\mathbf{W}_{R}, \end{aligned} YR=α0a(θ0)aH(θ0)X+p=P,p=0Pm=1Mαp,ma(θm)aH(θm)XJp+WR,
J p = [ 0 … 0 ⏟ p zeros  1 … 0 0 … 0 0 ⋱ ⋮ ⋮ ⋮ ⋮ 1 0 … 0 0 … 0 ] ∈ R L × L \mathbf{J}_{p}=\left[\begin{array}{cccc} \underbrace{0 \ldots 0}_{\text {p zeros }} & 1 & \ldots & 0 \\ 0 \ldots 0 & 0 & \ddots & \vdots \\ \vdots & \vdots & \vdots & 1 \\ 0 \ldots 0 & 0 & \ldots & 0 \end{array}\right] \in \mathbb{R}^{L \times L} Jp=p zeros  000000100010RL×L
只能说J这个矩阵的设计确实很秒,喵喵喵。
之后就是雷达中的匹配滤波处理
D = 1 L Y R X H = 1 L α 0 a ( θ 0 ) a H ( θ 0 ) X X H + 1 L ∑ p = − P , p ≠ 0 P ∑ m = 1 M α p , m a ( θ m ) a H ( θ m ) X J p X H + 1 L W R X H , \begin{aligned} \mathbf{D} &=\frac{1}{L} \mathbf{Y}_{R} \mathbf{X}^{H}=\frac{1}{L} \alpha_{0} \mathbf{a}\left(\theta_{0}\right) \mathbf{a}^{H}\left(\theta_{0}\right) \mathbf{X} \mathbf{X}^{H} \\ &+\frac{1}{L} \sum_{p=-P, p \neq 0}^{P} \sum_{m=1}^{M} \alpha_{p, m} \mathbf{a}\left(\theta_{m}\right) \mathbf{a}^{H}\left(\theta_{m}\right) \mathbf{X} \mathbf{J}_{p} \mathbf{X}^{H}+\frac{1}{L} \mathbf{W}_{R} \mathbf{X}^{H}, \end{aligned} D=L1YRXH=L1α0a(θ0)aH(θ0)XXH+L1p=P,p=0Pm=1Mαp,ma(θm)aH(θm)XJpXH+L1WRXH,
将问题直接转成下面这个,感觉一瞬间就降低优化式子的难度了。
P I S L = ∑ p = − P , p ≠ 0 P ∥ X p J p X H ∥ F 2 P_{\mathrm{ISL}}=\sum_{p=-P, p \neq 0}^{P}\left\|\mathbf{X}_{p} \mathbf{J}_{p} \mathbf{X}^{H}\right\|_{F}^{2} PISL=p=P,p=0PXpJpXHF2

天线导向矢量约束

雷达在执行不同的任务时需要的波束是不同的,比如搜索目标的时候,肯定是需要全向天线,在跟踪目标的时候是定向波束,保证目标一直在强波束照射下。所以采用定向的时候,雷达大致是知道目标方向的,这里很难处理的是,通信用户和雷达探测目标应该不会在同一个方向吧,那么雷达采用定向天线的话就很容易造成通信质量降低,这又是一个Trade-off的问题,能不能集成来提升性能呢?

MIMO天线发射信号为 a ∗ ( θ ) X \mathbf{a}^{*}(\theta) \mathbf{X} a(θ)X, 在 θ \theta θ处的发射功率为 P ( θ ) = a ∗ ( θ ) R X a ( θ ) P(\theta)=\mathbf{a}^{*}(\theta) \mathbf{R}_{X} \mathbf{a}(\theta) P(θ)=a(θ)RXa(θ).
R X = 1 L X X H . \mathbf{R}_{X}=\frac{1}{L} \mathbf{X X}^{H}. RX=L1XXH.

其中 a ( θ ) = \mathbf{a}(\theta)= a(θ)= [ 1 , e j 2 π Δ sin ⁡ ( θ ) , … , e j 2 π ( N − 1 ) Δ sin ⁡ ( θ ) ] T ∈ C N × 1 . \left[1, e^{j 2 \pi \Delta \sin (\theta)}, \ldots, e^{j 2 \pi(N-1) \Delta \sin (\theta)}\right]^{T} \in \mathbb{C}^{N \times 1} \quad. [1,ej2πΔsin(θ),,ej2π(N1)Δsin(θ)]TCN×1.

全向天线约束

min ⁡ X ∥ H X − S ∥ F 2 \min _{\mathbf{X}}\|\mathbf{H X}-\mathbf{S}\|_{F}^{2} XminHXSF2
s.t. 1 L X X H = P T N I N \frac{1}{L} \mathbf{X} \mathbf{X}^{H}=\frac{P_{T}}{N} \mathbf{I}_{N} L1XXH=NPTIN,

这是可以直接解出来的,用奇异值分解加最小二乘。可惜我不是很明白,希望有知道的小伙伴能跟我讲一讲。
X = L P T N U I N × L V H \mathbf{X}=\sqrt{\frac{L P_{T}}{N}} \mathbf{U} \mathbf{I}_{N \times L} \mathbf{V}^{H} X=NLPT UIN×LVH
其中 U Σ V H = H H S \mathbf{U} \mathbf{\Sigma} \mathbf{V}^{H}=\mathbf{H}^{H} \mathbf{S} UΣVH=HHS H H S \mathbf{H}^{H} \mathbf{S} HHS S V D \mathrm{SVD} SVD

定向天线约束

定向类似于全向的方法,也是可以解出来的
min ⁡ X ∥ H X − S ∥ F 2  s.t.  1 L X X X H = R d , \begin{aligned} &\min _{\mathbf{X}}\|\mathbf{H X}-\mathbf{S}\|_{F}^{2} \\ &\text { s.t. } \frac{1}{L} \mathbf{X X} \mathbf{X}^{H}=\mathbf{R}_{d}, \end{aligned} XminHXSF2 s.t. L1XXXH=Rd,
where R d \mathbf{R}_{d} Rd is the desired Hermitian positive semidefinite covariance matrix. We consider its Cholesky decomposition, which is
R d = F F H \mathbf{R}_{d}=\mathbf{F} \mathbf{F}^{H} Rd=FFH
X = L F U ~ I N × L V ~ H \mathbf{X}=\sqrt{L} \mathbf{F} \tilde{\mathbf{U}} \mathbf{I}_{N \times L} \tilde{\mathbf{V}}^{H} X=L FU~IN×LV~H

相似性约束

由于天线波束方向约束的X可以直接解出来,而且采用X作为发射波形,其他的性能提高并不多,所以给X再添加一个约束。
∥ x − x 0 ∥ 2 2 ( ∥ x − x 0 ∥ ∞ ) ≤ ϵ \left\|\mathbf{x}-\mathbf{x}_{0}\right\|_{2}^{2}\left(\left\|\mathbf{x}-\mathbf{x}_{0}\right\|_{\infty}\right) \leq \epsilon xx022(xx0)ϵ
这个 X 0 X_0 X0就是波束约束解出来的,但这里同时也可以换成其他好的波形,那么具有抗多普勒的波形是否可以呢?

天线功率限制

天线发射的总功率肯定不是无限的,所以
 s.t.  1 L ∥ X ∥ F 2 = P T \text { s.t. } \frac{1}{L}\|\mathbf{X}\|_{F}^{2}=P_{T}  s.t. L1XF2=PT

恒模限制

不管是雷达发射机还是通信发射机,对PAPR的要求都很高,所以希望每一个时刻发射波形的功率恒定
 s.t.  ∣ x n ∣ = P max ⁡ T N , 1 ≤ n ≤ T N \text { s.t. }\left|x_{n}\right|=\sqrt{\frac{P_{\max }}{T N}}, 1 \leq n \leq T N  s.t. xn=TNPmax ,1nTN

结语

对于MIMO-OFDM雷达肯定还有一些其他的约束,而且以上优化都会存在先验知识很强,比如信道状态信息,雷达目标个数,虽然这些在目前通信处理和雷达处理中很容易解决,但能在发射波形上做一些优化吗?

你可能感兴趣的:(雷达通信一体化,其他)