关于Capon最优加权系数和MVDR加权系数等价的证明

首先给出Capon的功率谱,主要思想是最大化输出信噪比,即保证阵列对期望信号方向相应一定时,最小化输出功率,可以形成下面的优化问题
min ⁡ w w H R w  s.t.  w H a 0 = 1 \min _{\mathbf{w}} \mathbf{w}^{\mathbf{H}} \mathbf{R} \mathbf{w} \text { s.t. } \mathbf{w}^{\mathrm{H}} \mathbf{a}_{0}=1 wminwHRw s.t. wHa0=1
为了求解上述优化问题,可以用拉格朗日乘数法进行计算,得到如下函数
J ( w , λ ) = w H R w + λ ( w H a 0 − 1 ) J(\mathbf{w}, \lambda)=\mathbf{w}^{\mathrm{H}} \mathbf{R} \mathbf{w}+\lambda\left(\mathbf{w}^{\mathrm{H}} \mathbf{a}_{0}-1\right) J(w,λ)=wHRw+λ(wHa01)
其中 λ ≥ 0 \lambda\ge0 λ0为拉格朗日算子,对 w ∗ \mathbf{w}^{*} w求导,并令其等于零,有
R w + λ a 0 = 0 \mathbf{R} \mathbf{w}+\lambda \mathbf{a}_{0}=0 Rw+λa0=0
可以得到解 w = − λ R − 1 a 0 \mathbf{w}=-\lambda \mathbf{R}^{-1} \mathbf{a}_{0} w=λR1a0,代替约束条件 w H a 0 = 1 \mathbf{w}^{\mathrm{H}}\mathbf{a}_{0}=1 wHa0=1中,可以求得
λ = − 1 a 0 H R − 1 a 0 \lambda=-\frac{1}{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} λ=a0HR1a01
因此,最优解为
w = R − 1 a 0 a 0 H R − 1 a 0 \mathbf{w}=\frac{\mathbf{R}^{-1} \mathbf{a}_{0}}{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} w=a0HR1a0R1a0
上述过程就是标准的Capon波束形成(SCB),此时将最优解带入目标函数中,可以求得阵列的输出功率为
P = w H R w = 1 a 0 H R − 1 a 0 P=\mathbf{w}^{\mathrm{H}} \mathbf{R} \mathbf{w}=\frac{1}{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} P=wHRw=a0HR1a01
该功率又被称为Capon功率,它是在无失真接收期望信号时,阵列输出的最小功率,当导向矢量进行变化时,上述功率会形成在空间上的功率分布情况,又被称为Capon功率谱。
同时我们考虑到 R = σ s 2 a 0 a 0 H + R i + n \mathbf{R}=\sigma_{\mathrm{s}}^{2} \mathbf{a}_{0} \mathbf{a}_{0}^{H}+\mathbf{R}_{\mathrm{i}+n} R=σs2a0a0H+Ri+n,那么则有 w H R w = σ s 2 w H a 0 a 0 H w + w H R i + n w = σ s 2 + w H R i + n w \mathbf{w}^{\mathrm{H}} \mathbf{R} \mathbf{w}=\sigma_{\mathrm{s}}^{2} \mathbf{w}^{\mathrm{H}} \mathbf{a}_{0} \mathbf{a}_{0}^{\mathrm{H}} \mathbf{w}+\mathbf{w}^{\mathrm{H}} \mathbf{R}_{\mathrm{i}+\mathrm{n}} \mathbf{w}=\sigma_{\mathrm{s}}^{2}+\mathbf{w}^{\mathrm{H}} \mathbf{R}_{\mathrm{i}+\mathrm{n}} \mathbf{w} wHRw=σs2wHa0a0Hw+wHRi+nw=σs2+wHRi+nw,因此上述优化问题,可以等价为
min ⁡ w w H R i + n w  s.t.  w H a 0 = 1 \min _{\mathbf{w}} \mathbf{w}^{\mathrm{H}} \mathbf{R}_{\mathrm{i}+\mathrm{n}} \mathbf{w} \text { s.t. } \mathbf{w}^{\mathrm{H}} \mathbf{a}_{0}=1 wminwHRi+nw s.t. wHa0=1
该优化问题的物理意义是基于期望信号被无失真接收的条件下最小化干扰加噪声功率,也被称为最小方差无失真响应波束形成(minimum variance distortionless respond,MVDR),同样地,我们可以用拉格朗日乘数法对其进行求解,可以得到该优化问题的解为
w = R i + n − 1 a 0 a 0 H R i + n − 1 a 0 \mathbf{w}=\frac{\mathbf{R}_{\mathrm{i}+n}^{-1} \mathbf{a}_{0}}{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}_{\mathrm{i}+\mathbf{n}}^{-1} \mathbf{a}_{0}} w=a0HRi+n1a0Ri+n1a0
此时,阵列的输出为 w H R i + n w = E { ∣ w H x i + n ( t ) ∣ 2 } \mathbf{w}^{\mathrm{H}} \mathbf{R}_{\mathrm{i}+\mathrm{n}} \mathbf{w}=E\left\{\left|\mathbf{w}^{\mathrm{H}} \mathbf{x}_{\mathrm{i}+\mathrm{n}}(t)\right|^{2}\right\} wHRi+nw=E{wHxi+n(t)2}表示输出数据中干扰加噪声信号的方差,下面我们来证明,这两种最优加权系数是等价的。
根据矩阵求逆引理,可以得到
R i + n − 1 = ( R − σ s 2 a 0 a 0 H ) − 1 = R − 1 + σ s 2 R − 1 σ s 2 a 0 a 0 H R − 1 1 − σ s 2 a 0 H R − 1 a 0 \mathbf{R}_{\mathrm{i}+\mathrm{n}}^{-1}=\left(\mathbf{R}-\sigma_{s}^{2} \mathbf{a}_{0} \mathbf{a}_{0}^{\mathrm{H}}\right)^{-1}=\mathbf{R}^{-1}+\frac{\sigma_{s}^{2} \mathbf{R}^{-1} \sigma_{s}^{2} \mathbf{a}_{0} \mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1}}{1-\sigma_{s}^{2} \mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} Ri+n1=(Rσs2a0a0H)1=R1+1σs2a0HR1a0σs2R1σs2a0a0HR1
将其带入MVDR最优加权系数中可以得到,分子为
R i + n − 1 a 0 = R − 1 a 0 1 − σ s 2 a 0 H R − 1 a 0 \mathbf{R}_{\mathrm{i+n}}^{-1} \mathbf{a}_{0}=\frac{\mathbf{R}^{-1} \mathbf{a}_{0}}{1-\sigma_{s}^{2} \mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} Ri+n1a0=1σs2a0HR1a0R1a0
和分母为
a 0 H R i + n − 1 a 0 = a 0 H R − 1 a 0 1 − σ s 2 a 0 H R − 1 a 0 \mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}_{\mathrm{i+n}}^{-1} \mathbf{a}_{0}=\frac{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}}{1-\sigma_{s}^{2} \mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} a0HRi+n1a0=1σs2a0HR1a0a0HR1a0
因此可以得到
R i + n − 1 a 0 a 0 H R i + n − 1 a 0 = R − 1 a 0 a 0 H R − 1 a 0 \frac{\mathbf{R}_{\mathrm{i}+\mathbf{n}}^{-1} \mathbf{a}_{0}}{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}_{\mathrm{i}+\mathbf{n}}^{-1} \mathbf{a}_{0}}=\frac{\mathbf{R}^{-1} \mathbf{a}_{0}}{\mathbf{a}_{0}^{\mathrm{H}} \mathbf{R}^{-1} \mathbf{a}_{0}} a0HRi+n1a0Ri+n1a0=a0HR1a0R1a0
从上面可以看出,两种形式的最优加权系数是等价的。
当然,还可以通过仿真来验证两者的等价性,仿真参数设置如下:

参数名参数值
阵元数10
阵元间距λ/2
快拍数100
期望信号角度-5°
干扰信号角度-40°、30°
噪声功率1
信号和干扰功率0dBw

通过对阵列的方向图仿真可以得到如下的结果
在这里插入图片描述
可以看出,两种加权系数的方向图完全重合,同时在期望信号方向上达到最大增益,干扰角度方向上形成零陷,也就说明了两种加权系数的有效性和等价性
在这里插入图片描述

将两者画成上下两幅图也可以看出以上结论。
那么,画出两者的空间谱如下
在这里插入图片描述

从空间谱也可以看出,两种的等价性
在这里插入图片描述
两种方法的空间谱如上,从两幅图也可以看出以上结论。

  • 6
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值