数字波束形成递推最小二乘算法研究.docx

数字波束形成递推最小二乘算法研究摘要：数字波束形成(DBF)是一门新兴的阵列天线技术，自适应波束形成(ADBF) 针对复杂信号环境对波束形成的一个很好的补充自适应技术的核心思想是优化 理论，常用的最佳准则有最小均方误差(MMSE),最大信噪比(MaxSNR),线性约 束最小方差(LCMv),最大似然(ML),最小二乘法(LS)等；本问主要对递推最小二 乘算法(RLS)进行了包括了流程设计、性能分析和仿真评价关键字：DBF、ADBF、RLS、一数字波束形成简介数字波束形成(DBF)⑴是随着数字信号处理方法的发展而建立起来的一门新 技术，这种技术不仅能充分保留阵列天线上收集的信息，而且能利用复杂的数字 信号信号处理方法对信息进行处理，因此它具有--般雷达不具有的超低分辨率和 低副瓣的性能，波束扫描，自适应波朿控制等由于采用了先进的数字信号处理 方法和自适应技术，DBF雷达不仅性能优越，血且非常灵活，被认为是下〜代雷 达系统中广泛使用的一项新技术相对于一般模拟波朿形成技术，数字波束形成技术具有非常多的优点，比如 超低分辨率，超低副瓣，抗电子干扰，扫描迅速，多H标处理，高性能并行数字 处理等自适应波束形成(ADBF)针对复杂信号环境对波束形成的-•个很好的补充。

它的核心思想是采用各种优化算法和自适应算法，根据阵列对不同环境的响应， 产牛不同的加权系数，己取得最好的信号特征自适应波朿形成可以实现自适应 的空一时处理，进行各种非线性处等，波束灵活，抗干扰强，具有较好的自校正 和自适应能力阵列信号处理是现代信号处理〔2］的一个分支，其本质是利用空间分散排列的传 感器阵列和多通道接收机来获取信号的吋域和空域等多维信息，以达到检测信号 和提取其参数的H的阵列信号处理的主要内容可分为波束形成技术，空间谱估 计等方面技术，他们都是基于对信号进行空间采样的的数据进行处理，因此这些 技术是相互渗透和相互关联的二数字波束形成原理和典型波束形成准则2.1数字波束形成原理假设接收天线为N阵元均匀直线阵⑶，FI标的来波方向为乩第k个波束指向 为％,相邻阵元间距为d,信号波长为儿且各阵元都是各向同性的,对K阵元的 加权分别为W“，W2K,…,WNK，信号是窄带信号简单地讲，数字多波束形成器就 是一个乘加器，如图2-1所示右表示第i阵的接收信号经过A/D变换和数字正交 后的复数字信号加权系数“认=a^exp (-/△①J其中△①k = Yrising表示提供 的阵内相位补偿值，©表示降低天线副瓣所需的幅度加权。

进行相位和幅度补偿后, 并对各阵列的输出信号相加，即可得：FkW = ^wikXi (2-1)对式(2・1)进行求模，就可以得到第K个波束的天线方向图函数天线的方向图不 是唯一的，根据要求，相同的数据可以用不同方法加权(改变权因子)，以便形成不 同形式的波束和任意多的波束通过数字波束合成后输出的信号(如信号确实投射 在该波束内)可以大幅度的提高信号的SNRo最后得到的复输出信号再直接送入后 续处理单元图2・1自适应数字波束形成原理图表示成矩阵形式如下：环⑼=叽X认其中为权矢量，Xk为数据矢量，T表示转置，它们 都为列矢量2.2自适应数字波束形成的原理DBF技术是当代雷达的热点技术之一，对于阵列信号处理来说，它往往和自 适应技术相互联系因为而对非常复杂的环境时候，恶化的工作环境将会使不具 备自我适应，自我调整能力的DBF雷达的性能急剧下降因此，为了提高DBF 系统的抗干扰能力，必须求助于自适应数字波朿形成(ADBF)技术，自适应阵列是 实现ADBF的基础口适应阵列的结构框图如图2-1示从图中可以看出自适应阵列是由按一定空间排列的多个阵单元构成，是一种 在实际环境下自行控制其方向图的天线系统。

它能实时地对外界未知的干扰环境 作出反应，在干扰的到达方向形成零点或降低此方向的副瓣电平，这样便可以保 证接收所需信号与干扰噪声的信噪比有一个最佳值口适应阵列正是利用这种空 间特性，改善了阵列输出SNR,抑制了强干扰自适应技术的核心思想是优化理论〔4】，我们熟悉的优化理论有拉格朗H求极值 函数的变分法，最陡下降法，最小二乘法等H前常用的最佳准则有最小均方误差(MMSE),最大信噪比(MaxSNR),线性 约束最小方差(LCMv),最大似然(ML),最小二乘法(LS)等，下面简单介绍几种较 常使用的准则2.3最小均方误差准则(MMSE)最小均方误差准则就是使阵列输出压⑹=wQGJK代表K时刻)与参考信号 d(k)的均方误差最小，均方误差为：E(f2(/c)) = E[d伙)F - + W?RWk (2-2)其• |• r=E(d(k)Xk),R=E(XXH) ~般地将R称为互相关矩阵[役 将式(2-4)对于权向 量求梯度，得到梯度算子：00(£(严伙)))=-2厂+ 2/?%令梯度算子为零，得到 最小均方误差准则下的最佳权向量：woptR~rr2.4最大信噪比原则(MaxSNR)最大信噪比原则是基于期望信号的功率与噪声功率之比的最大的准则，假设 期望信号为S,且Rt=E (ssH), Rw=E(uuh),其中u表示噪声，则有：0-/ = E(\Whs\2) = WhRsW 此=E(\Whu\2) = WhRuW这时，输出的信噪比SNR为：於 WhRsW经计算，使得输出信噪比最大的最佳权向量是对应于矩阵尺补心的最大特征值的特 征值的特征向量，得到的最佳权系数满足：R^Rs^opt = ^max^opt2.5线性约束最小方差准则(LCMV)在已知期待信号的来波方向和参考信号的条件下，最小方差准则是通过最小 化阵列输出的噪声方差，来取得对信号S的较好的增益。

经权重后的波束形成器 的输出为:y(t) = WTX{t) = WTs + WSz为保证波束形成对信号S增益，必须对波 束形成器的权向量加以限制，使其在信号S的方向产牛一定的增益，即：WHa0 = g其中，如为期望信号的方向矢量，则最佳权重可以表示为：W°pt =顽甌三递推最小二乘算法(RLS)3.1自适应算法简介口适应算法所采用的最优准则有最小均方误差(LMS)准则⑹，最小二乘(LS)准 则、最大信嗓比准则和统计检测准则等，其中最小均方误差(LMS)准则和递推最小 二乘法(LS)准则是H前最为流行的自适应算法准则LMS算法和RLS算法由于采 用的最优准则不同，I大I此这两种算法在性能，复杂度等方面均有许多差别对于某个口适应算法，衡量其算法性能的指标主要为以下儿个：(1) 算法的收敛速度用算法达到稳定状态(即与最优值的接近程度达到 一定值)需要迭代的次数表示2) 算法的跟踪性能它指信道发牛变化时，算法自适应跟踪信道的能力3) 算法的稳健性只当输入病态的情况下算法能否正常工作4) 算法计算的复朵程度自适应波束形成算法可有多种分类方式，一般的分类方法根据算法是否需要 显式的训练序列，可分为非盲算法和盲算法两类，另外，根据权值作用于阵元或 波束，还可分为阵元空间和波束空间两类。

盲算法一般利用调制信号本身固有的、与具体承载的信息比特无关的一些特 征，如恒模、子空问、有限符号集和高阶统计量等与非盲算法相比，它无需发 送已知符号，故可节省系统频谱资源其缺点是-般收敛速度较慢，且存在相位 模糊问题非盲算法相对盲算法而言，通常误差较小，收敛速度也较快，但需要 占用一定的系统资源用于传送参考信号、训练序列等即使非盲算法收敛速度快， 但仍然跟不上快衰落变化的速率要求对于通信而言，算法的收敛速度是很重要 的一个性质所以本文主要介绍非盲算法中的RLS算法非盲算法主要包括最小均方误差算法LMS、采样矩阵求逆法SMI和迭代最小 二乘算法RLSLMS算法由于其收敛速度很慢，因此在信号环境变化很快时其性 能不是很好，且算法性能对阵列信号协方差矩阵的特征值散布度很敏感，当散布 度很大吋，很难收敛SMI算法收敛速度比LMS快，但是算法复杂了不少，而 且在硬件实现中的有限字长效应会给求逆运算带来数值上的不稳定RLS算法是 基于使每一快拍的阵列输出平方和最小的准则，即最小二乘（LS）准则它利用 了从算法初始化后得到的所有阵列数据信息，用递推方法来完成矩阵的求逆运算, 因而收敛速度快，对特征值的散布度不敏感，且能实现收敛速度与计算复杂性之 间的折衷。

一般在大信噪比的情况下，RLS比LMS的收敛速度快一个数量级3.2递推最小二乘算法（RLS）RLS算法是SMI算法的迭代过程，该算法的性能准则是使参考信号d（〃）和阵列输出信号只“）之间的加权均方误差最小与LMS不同的是，RLS算法的性能量 度弁•则是使信号的加权平方误差和最小化RLS算法在权值迭代调整过程中用一 个与信号相关的附加矩阵代替LMS算法中的步长，从而提高了收敛速度RLS算 法的权向量迭代式为：/ 、 (、 戶(刃)X (刀 +1)W5 + 1) = W5) + ——4 ； x2 + X(刃 + \)卩5) X (/7 + 1) (3・17)[d(/7 + 1) - ^(/7)7(/7 + 1)]P(n + 1)二=扛应—PS*； 5 + 1）/57"）（3-18） 久 Q +才（刀+ 1）"（刀）X （门+ 1）其中“⑷表示参考信号，p（"）为迭代过程中的过渡矩阵，其中2称为遗忘因 子（2 < 1）0引入遗忘因子的作用是让离n吋刻较近的误差有较大的权重，距离较 远的拥有降低的权重，确保以前观测到的数据被渐渐“遗忘”，从而使滤波器工作 在一个平稳状态下下面是算法推导：R5)=工2 = 0=久* £ /Tf (i)u H (i) + u (71)* W(H)i = 0=A* R(n -1) + uH (n) * u(n)设.4 = 2? (n ) 3 1 = —1} C = w (h ) D = 1则 r=r+ci)w.rl = E — BC{D 十 CBBCyCg8才 IT： 5 一 55加％膚5 一 1) ~1 一才％ 日 5)/T：(n-LM5)-令：P(7i) = fi- (n) /?(n) = AR(n-l) + x(n)x,y(n)则：尸何=2[£5-1)-弘対山尸少一叨只口 -1血何W0 =-——v 其中 A¥U Oi)P(n・l)M(/i)为增益向量，又：P(n)//(n)=〃(□)2 R(n-1) + ”(n)/z〃 (n)R '(n-Q/zCn)/l + //(n)Rj(n—l)〃H(n)=心)nr(n) = S a"—'u(i)d*(i) = 2r(n— 1) + d^(n)u(n) i=0P(n)=犷[P(n -1)-k(n)uH(n)P(n -1)]所以w(n) = R_1(n) r(n) = P(n) r(n)=P(n — 1)厂(71 —1) + (n)[P(n — l)u(n)-k(n)uH (n)P(n — l)w(n)] - k(n)uH (n)P(n — l)r(n 一 1)=w(n 一 1) + /(“)£(〃) 一 k(n)u H (n)vv(/7 -1)化简得：w(/i) = w(n-1) + k(n)e* (n) (*)式中:e(n) = d(n) 一 wH (n 一 l)w(n)先验误差3.3总结RLS算法的步骤1、 初始化：w(0)=0, R(0)= o I,2、 更新：对于“1、2 ••-计算:滤波:y(n) = wH (n 一 l)w(n)估计误差:e(n) = d(n)~ y(n)更新R(n) =P(n-l)u(n)a + u" (n)P(n - l)u。