![阵列信号处理及MATLAB实现(第2版)](https://wfqqreader-1252317822.image.myqcloud.com/cover/854/41202854/b_41202854.jpg)
3.5 基于投影和斜投影的波束形成算法
3.5.1 基于投影的波束形成算法
1. EBS波束形成算法
假设有l个期望信号,J个干扰,对有限次快拍下的协方差矩阵进行特征分解:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-86-13.jpg?sign=1738897358-f1r6Nji3v4ZhFGv87f0xF0rEvcqhUPle-0-bf84e916cb76bbf0f3707eacf1fb3f9c)
式中,是相应的M个特征值,其对应的特征向量为
,记
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-2.jpg?sign=1738897358-LCfZ1lpdTAyaotfgtJI8IRwr97KcEZB9-0-e6692f77d6b22193b56b8744fc7f7b61)
Us的列向量张成信号子空间,而Un的列向量张成噪声子空间。在SMI算法中,自适应权为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-3.jpg?sign=1738897358-BpnYCuFzaXKQx9e3NKQXrcwN2fV67vu4-0-2ae700ceeed12c72f4c02a8c77a40f99)
上式表明权向量由信号子空间分量和噪声子空间分量构成。在理想情况下,期望信号位于信号子空间,有,因此,权向量
仅为信号子空间的分量,噪声子空间的分量为零。ESB算法就基于这种原理,摒弃权向量在噪声子空间中的分量而仅保留在信号子空间中的分量,成为基于特征结构的自适应波束形成算法或投影算法,即
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-6.jpg?sign=1738897358-hKjzA00W4lkCMDg3hNrRXbeGAPKA3W2t-0-c620faab4b24966e26bf2039e69d9afd)
当数据协方差矩阵中含有较强的期望信号时,该算法较为有效。当期望信号功率较小时,直接摒弃权向量在噪声子空间中的分量将会有较大的误差。一种极端情况是,
中不含期望信号,即在理想情况下,
不成立。权向量在噪声子空间的分量不为零,此时
不是最优权了,它将导致输出SINR性能下降。另外,由于噪声子空间的扰动,使自适应方向图发生畸变。所以ESB算法不适用于小期望信号。
2. EBS改进算法(IESB)
(1)进行特征分解后,特征值从大到小排列,计算第J+1和J+2两个特征值之比,当
大于某个门限值时,则构成
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-13.jpg?sign=1738897358-yHhu546yBGcZRrQMHqN7lRwcCt8LX5PJ-0-36935f1d1e6b538bb52ca3e09e937827)
否则
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-14.jpg?sign=1738897358-g210uOx7vxYyAr15Jmz5OE2NV3y7fRgH-0-c99cdcb9131875baf1b00cadc9052f43)
(2)对进行奇异值分解
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-16.jpg?sign=1738897358-7ZmMDYdol9uoNU3ZqmBIkFlGAANN5IiZ-0-b1678817653d8c780027efaba2b37775)
(3)将SMI方法求得的权向量向
的大特征值对应的左奇异向量列空间
投影,即
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-87-20.jpg?sign=1738897358-PeXAOQtWmWAOIBHzayfzG4iMxCZNfvty-0-3255b461d5b8c8d1ba72b3b148063832)
由于引入了期望信号导向向量,并且在期望信号功率与噪声功率相当或更弱时,去除了干扰较大的特征向量,该方法能在输入信号较大时保持基于特征结构的自适应波束形成算法性能,又能在期望信号较小(甚至为零)时具有较好的波束保形能力。但是,该方法计算量较大,需要进行一次特征分解和一次奇异值分解。
3.5.2 基于斜投影的波束形成算法
采样矩阵求逆(SMI)算法是最常用的自适应波束形成算法,该算法具有较快收敛速度。但是SMI在少快拍数、高信噪比和相干信源情况下会导致副瓣电平升高,主瓣偏移,波束畸变,输出信干噪比(SINR)下降。文献[20,21]提出了对角线加载的波束形成算法来抑制方向图畸变,分析了加载量对自适应阵列信干噪比的影响。对角线加载技术能减弱小特征值对应的噪声波束的影响,改善了方向图畸变。但是加载量的确定一直以来是一个比较困难的问题,至今仍没有很好地解决,从而限制了对角加载技术的实际应用。文献[25,26]中利用投影算子改善了波束形成的稳健性,但投影算法在相干信源情况下性能下降,而且投影算子需要知道期望信号和干扰信号的方向向量,这在实际系统中很难满足。斜投影算子是投影算子扩展,文献[27,28]中将斜投影成功应用于单输入多输出系统的盲辨识,文献[29]中将斜投影应用到DOA估计,文献[30]中将斜投影应用到卷积混合信号的盲分离。文中将研究基于斜投影的波束形成算法,对接收信号进行斜投影可有效消除干扰,进而提高波束形成的稳健性,而且该算法在少快拍数和相干信源情况下仍具有较好的波束形成性能。
1. 斜投影基础原理
矩阵(F≤N),矩阵A的投影矩阵为
;
为矩阵A的正交投影矩阵。斜投影是正交投影扩展,斜投影算子
为沿着与子空间Rang(B)平行的方向,到子空间Rang(A)上的投影算子:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-88-5.jpg?sign=1738897358-KiTKmbUwGfRCPvSBv00C8P0rl41jCfGY-0-2311aa41dd1ad877a91a9db7e65072c6)
其中,EAB是方阵,且为非对称矩阵,它有以下特性:EABA=A;EABB=O。
定义A=[H,S],则
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-88-6.jpg?sign=1738897358-8EoCOJ3ORk1yzyqFy42HewZJn1WH0rxS-0-4e820585aaa42c1ab2f5d3961ae2905f)
矩阵A的投影矩阵PA
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-88-7.jpg?sign=1738897358-MtJo1jdmOCAsA1LxrRyvy1koBvyvaSAW-0-5e4e768a64b035313e01ba9431c6d620)
可见,矩阵A的投影矩阵PA可表示为两个斜投影矩阵之和。
2. 阵列信号的模型
假设有K个信源,信源i的信号为期望信号,考虑M元均匀线阵接收信号可表示为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-89-1.jpg?sign=1738897358-y7hEkLNHQKs7KKiB0yBWhjCgSJNv8M7M-0-83d7c83a1ec05981e710c652d65dd44d)
其中,是信源k的发射信号;
是信源k的波达方向;
是信源i的归一化方向向量;n(t)为均值为0,方差是
的白噪声;K为信源数。方向矩阵
;信源矩阵
;矩阵
包含K-1个信源方向向量,为干扰信源的方向矩阵;
包含K-1个信源信号,为干扰信源矩阵。
接收信号的协方差矩阵表示为
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-89-10.jpg?sign=1738897358-oH563BnOB98EnsHaO8OUluAPIHR7UdRu-0-77354369d4adc9382d485c7a98cca057)
其中,是期望信号的功率,
是干扰源信源协方差矩阵。
对接收信号的协方差矩阵进行特征分解,有
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-89-13.jpg?sign=1738897358-PKz3BoCQEBOukhSJ1cAWKNtLWRUaZWdX-0-c030d7dd559828d2a47ce6494914f352)
其中,US是信号子空间,UN是噪声子空间。
定义矩阵RA:
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-89-14.jpg?sign=1738897358-Ln8LAMLZsIuAdkXUJ56jDZSnbVLMnMN0-0-3babde6521bc43b59ddc975ccca6402f)
3. 基于斜投影的波束形成算法
对阵列接收信号进行斜投影,则
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-89-16.jpg?sign=1738897358-FzUygsXu9pUOJzXpEDTzB5G3pf4zvlJI-0-abb68e36b731d80d09e350715b9c7d65)
其中,矩阵是沿着与子空间
平行的方向,到子空间
上的投影算子。
斜投影后对信号进行空域匹配滤波,得到
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-90-1.jpg?sign=1738897358-cpZbYUF0JnNaFdNEWfuuyYosdeyxA541-0-b8b97ba7f823806b15180e8eb5c34238)
斜投影矩阵需要知道期望信号和所有干扰源方向向量,在实际的系统中很难满足。
引理3.5.1 斜投影矩阵,也可写成
![](https://epubservercos.yuewen.com/4102F9/21440186201518106/epubprivate/OEBPS/Images/37342-00-90-4.jpg?sign=1738897358-FjFfuWi5Uklbo8Nbc7ZXrayN8pXmET5E-0-fc3bc228f947be85981d9db2161b1d88)
其中,是伪逆矩阵,
斜投影的波束形成算法的具体步骤如下。
步骤1:计算接收信号的协方差矩阵,并进行特征值分解,根据式(3-53)计算出,进而根据公式(3-56)计算出斜投影矩阵
步骤2:对阵列接收到的信号进行斜投影,如式(3-54)所示。
步骤3:对斜投影后的信号进行空域匹配滤波,如式(3-55)所示。这样就实现了斜投影的波束形成。
总之,基于斜投影的波束形成算法在不同SNR情况下皆具有较好的波束形成性能,且在少快拍数情况下仍具有较好的波束形成性能。而且,该算法只需要期望信号的方向和接收信号,是一种稳健且性能优越的波束形成算法。