留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于小波时间延迟估计的宽带信号波束形成算法研究

刘影 谢驰

刘影, 谢驰. 基于小波时间延迟估计的宽带信号波束形成算法研究[J]. 电子科技大学学报, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
引用本文: 刘影, 谢驰. 基于小波时间延迟估计的宽带信号波束形成算法研究[J]. 电子科技大学学报, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
LIU Ying, XIE Chi. Broadband Signal Beamforming Algorithm Based on the Wavelet and Time Delay Estimation[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
Citation: LIU Ying, XIE Chi. Broadband Signal Beamforming Algorithm Based on the Wavelet and Time Delay Estimation[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004

基于小波时间延迟估计的宽带信号波束形成算法研究

doi: 10.3969/j.issn.1001-0548.2018.02.004
基金项目: 

国家自然科学基金 61401075

详细信息
    作者简介:

    刘影(1984-), 女, 博士生, 主要从事信号与信息处理方面的研究

  • 中图分类号: TN911.6

Broadband Signal Beamforming Algorithm Based on the Wavelet and Time Delay Estimation

图(4)
计量
  • 文章访问数:  7054
  • HTML全文浏览量:  2265
  • PDF下载量:  172
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-05-04
  • 修回日期:  2017-11-01
  • 刊出日期:  2018-03-30

基于小波时间延迟估计的宽带信号波束形成算法研究

doi: 10.3969/j.issn.1001-0548.2018.02.004
    基金项目:

    国家自然科学基金 61401075

    作者简介:

    刘影(1984-), 女, 博士生, 主要从事信号与信息处理方面的研究

  • 中图分类号: TN911.6

摘要: 在阵列信号接收中,由于信号多径效应,阵列天线接收到不同时间延迟的期望信号,这使得高精度的波束形成技术成为信号处理中的一个难点。该文提出了一种考虑时延的波束形成算法。该算法首先通过小波算法对宽带接收信号进行接收信号时间延迟估计,将时延估计值与预存时延估计值进行误差比较,并对接近期望信号方向的空间多波束进行迭代优化,最后实现在空间中形成整体的自适应波束。仿真结果验证了改进算法的有效性。

English Abstract

刘影, 谢驰. 基于小波时间延迟估计的宽带信号波束形成算法研究[J]. 电子科技大学学报, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
引用本文: 刘影, 谢驰. 基于小波时间延迟估计的宽带信号波束形成算法研究[J]. 电子科技大学学报, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
LIU Ying, XIE Chi. Broadband Signal Beamforming Algorithm Based on the Wavelet and Time Delay Estimation[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
Citation: LIU Ying, XIE Chi. Broadband Signal Beamforming Algorithm Based on the Wavelet and Time Delay Estimation[J]. Journal of University of Electronic Science and Technology of China, 2018, 47(2): 183-188. doi: 10.3969/j.issn.1001-0548.2018.02.004
  • 目前自适应阵列信号处理技术在雷达、声纳和移动通信系统中有着越来越多的重要应用[1]。宽带信号波束形成算法可以采用时域和频域分析方法。文献[2]采用基于均匀一致性旁瓣级的设计方法,将宽带信号用一组有限脉冲响应(finite impulse response, FIR)滤波器来拟合,从而实现恒定束宽的波束形成。但该方法存在迭代收敛步长难以确定的问题,因而算法精度不高,收敛速度慢。FIR抽头系数较多,虽然可以采用无限脉冲响应(infinite impulse response, IIR)减轻矩阵运算的高复杂性,但IIR算法的不稳定性也使得IIR算法的收敛速度变慢[3]。时域的自适应滤波处理虽然对窄带干扰具有较好的抑制能力,但对滤波算法的稳定性要求较高,且时域收敛速度较慢,不适合用于干扰快速变化的情况[4-5]

    频域方法对阵列接收数据按照频率进行分块,并通过离散傅里叶变换(discrete Fourier transform, DFT)将时域数据转换成频域数据,由于在子块中信号频宽较窄,可以通过独立的窄带波束形成器实现波束成形,最后通过合成实现宽带信号波束成形[6]。文献[7-8]分别采用频率矩和滑动窗函数的方法,在有效地进行宽带波束形成的同时也抑制了宽带干扰。而小波算法可以根据频率的高低自动调节窗口大小,是一种自适应的时频分析方法,可以实现数据实时性处理能力的提高。

    对于阵列信号处理技术中的波束形成方法,目前最广泛采用的是最小方差无失真响应(minimum variance distortionless response, MVDR)波束合成、对角加载方法、旁瓣对消器(generalized sidelobe canceller, GSC)[9-12]。目前对波束形成技术的使用往往是在已知接收信号频率带宽的情况下进行的,然而在实际通讯中随着城市中建筑物的增多,在无线通讯中接收端不可避免地接收到多径信号。由于多径效应而产生的不同时间延迟扩展信号,使波束形成系统的通信质量下降。因此精确的时间延迟估计(time delay estimation, TDE)是提高波束形成精确度的重要内容。文献[13-14]提出了基于强噪声干扰的广义相关时间延迟估计测量算法,文献[15-16]提出了高阶与分数阶累积量方法,但是对于宽带信号时域算法的抗噪声性能降低,而频域算法无法提供足够多的有效频率采样点,现有的时间延迟算法效果不理想[17-18]。因此,本文提出了一种基于时间延迟估计的宽带信号波束形成方法,该方法通过引入加权函数因子,并对每次波束图主瓣迭代所获得的时间延迟估计值和实际期望信号的时延估计值之间误差与系统预存阈值相比较,从而不断迭代加权函数得到期望波束图的主瓣方向,形成整体的自适应波束。该算法在降低计算低复杂度和提高收敛速度方面具有一定优势。

    • 时间延迟估计是利用对不同的接收器接收到信号的时间差进行估计,从而确定信源的距离、方位、速度和移动方向等相关参量。设发射信号为x(t),采样信号时间间隔为Ts,经过阵元个数为M的阵列系统采样后的接收信号为:

      $$ y(t + \tau ) = \sum\limits_{m = 1}^M {\sum\limits_{i = 1}^I {{a_m}x(t + {\tau _i})} } + n(t) $$ (1)

      式中,I为多径的数目;τiτ分别为第i个多径路线所产生的时间延迟和系统产生的总时间延迟量;x(t)为发射信号的采样值;n(t)为零均值噪声;ai为第i个路线的幅度。将接收到的多径信号进行泰勒公式展开,得:

      $$ \sum\limits_{i = 1}^I {{a_i}(t + {\tau _i})} = x(t) + {a_1}{x^{(1)}}(t){\tau _1} + \frac{1}{2}{a_2}{x^{(2)}}(t)\tau _2^2 + \cdots $$ (2)

      y(t+τ)进行泰勒公式展开,有:

      $$ y(t + \tau ) = y(t) + {y^{(1)}}(t)\tau + \frac{1}{2}{y^{(2)}}(t){\tau ^2} + \cdots $$ (3)

      阵列系统接收到的时间延迟误差τ近似为:

      $$ \tau = \frac{{\sum\limits_{i = 1}^I {{\tau _i}} }}{I} $$ (4)
    • 设小波母函数ψ(t)为一平方可积函数,假定所要分析的函数是任意函数f(t),设pjf(t)代表f(t)在子空间Vj上的投影,该投影可以用小波变换的尺度空间$\{ {\phi _{j, k}}(t) = {2^{-j/2}}\phi ({2^{-j}}t-k), k \in Z\} $的线性组合来表示,有:

      $$ {p_j}f(t) = \sum\limits_k {{c_{j, k}}{\phi _{j, k}}(t)} $$ (5)

      式中,cj, k是各线性组合的权值。设pj+1f(t)为在子空间Vj+1上的投影,则:

      $$ {p_{j + 1}}f(t) = \sum\limits_k {{c_{j + 1, k}}} {\phi _{j + 1, k}}(t) $$ (6)

      式中,pj+1f(t)是在分辨率为j+1下所观察到的f(t)概貌,而cj+1, k是其离散逼近。同理,在空间Wj中,如果有带通函数ψ(t),其整数位移的集合{ψj+1, k(t), kZ}必构成Wj+1中的一组正交归一基,即:

      $$ < {\psi _{j + 1, k}}(t), {\psi _{j + 1, m}}(t) > = \delta (k-m) $$ (7)
      $$ {D_{j + 1}}f(t) = \sum\limits_k {{d_{j + 1, k}}} {\psi _{j + 1, k}}(t) $$ (8)

      式中,Dj+1f(t)表示在分辨率jj+1下对f(t)逼近的细节差异,称Dj+1f(t)是在分辨率j+1下的细节函数,dj+1, k是其离散细节,其中:

      $$ {d_{j, k}} = \sum\limits_m {h(m-2k){c_{j-1, k}}} $$ (9)
      $$ {c_{j, k}} = \sum\limits_m {g(m-2k)} {c_{j-1, k}} $$ (10)

      式(9)和式(10)表明,经j-1尺度空间的尺度系数和抽头长度为2kh(k)滤波器和g(k)滤波器的加权求和可以得到j尺度空间的尺度系数和小波系数。第j级小波分解对应的矩阵${2^{J-j + 1}} \times {2^{J-j + 1}}$(J为分解层数)为:

      $$ {W_i} = \left[{\begin{array}{*{20}{c}} {{h_k}}& \cdots &{{h_{2k-1}}}&0& \cdots &0&{{h_0}}& \cdots &{{h_{k-1}}}\\ {{h_{k-2}}}& \cdots &{{h_{2k - 3}}}&{{h_{2k - 2}}}& \cdots &0&0& \cdots &{{h_{k - 3}}}\\ \vdots&\cdots&\vdots&\vdots&\cdots&\vdots&\vdots&\cdots&\vdots \\ {{h_{k + 2}}}& \cdots &0&0& \cdots &{{h_1}}&{{h_2}}& \cdots &{{h_{k + 1}}}\\ {{g_k}}& \cdots &{{g_{2k - 1}}}&0& \cdots &0&{{g_0}}& \cdots &{{g_{k - 1}}}\\ {{g_{k - 2}}}& \cdots &{{g_{2k - 3}}}&{{g_{2j - 2}}}& \cdots &0&0& \cdots &{{g_{k - 3}}}\\ \vdots&\cdots&\vdots&\vdots&\cdots&\vdots&\vdots&\cdots&\vdots \\ {{g_{k + 2}}}& \cdots &0&0& \cdots &{{g_1}}&{{g_2}}& \cdots &{{g_{k + 1}}} \end{array}} \right] $$ (11)

      小波分解的矩阵表达式为:

      $$ {{\boldsymbol{W}}_T} = \left[{\begin{array}{*{20}{c}} {{{\boldsymbol{W}}_j}}&0& \cdots &0\\ 0&{\boldsymbol{I}}& \cdots &0\\ \vdots&\vdots&\ddots&\vdots \\ 0&0& \cdots &{\boldsymbol{I}} \end{array}} \right] \times \cdots \times \left[{\begin{array}{*{20}{c}} {{{\boldsymbol{W}}_2}}&0\\ 0&{\boldsymbol{\rm{I}}} \end{array}} \right] \times {{\boldsymbol{W}}_1} $$ (12)

      因此,小波滤波器h(k)和尺度滤波器g(k)将信号分解为低频分量和高频分量,从而减小了信号的相干性,提高自适应信号处理的速度。图 1为基于多通道小波算法的时间延迟估计框图。

      图  1  基于多通道小波滤波器组的时间延迟估计结构图

      对基于多通道小波算法的多通道数据采集中,通道抽取数目设为P个,每一通道中的信号分解级数为J,则算法步骤为:

      1) 对发射信号x(t)先进行P次多速率抽取,得到P个子带信号,再对每个子带信号进行j级小波变换,获得逼近信号Cp, j, x和细节信息Dp, j, x,同时对接收信号y(t)进行P次多速率抽取,得到P个子带信号,然后对每个子带信号进行j级小波变换,得到逼近信号Cp, j, x和细节信息Dp, j, x,其中p∈[1, P],j∈[1, J]。

      2)在第p个子带的第j级上,对逼近信号Cp, j, x和细节信息Dp, j, x分别进行延时并与Cp, j, xDp, j, x比较,找到使平方误差ep(τp)最小的时延τp,其中:${e_p}({\tau _p}) = {\left| {{C_{p, j, x(t-\tau)}}-{C_{p, j, y}}} \right|^2} + {\left| {{D_{p, j, (x-\tau)}} -{D_{p, j, y}}} \right|^2}$。

      3)设τp为小波变换第j级得到的时间延迟估计值,计算实际信号的时间延迟估计值τ*=2P2jτp

    • 可以根据接收到期望用户的信号特征对系统中已经预存的多个波束进行补偿,达到在空间中形成整体自适应波束的目的,从而通过多波束形成技术有效增加波束的覆盖能力。波束切换方法在多波束中对期望信号的方向进行跟踪,因而新算法只需优化子波束,避免计算复杂的高维系统矩阵,降低系统成本。本文设计的基于小波时间延迟估计的波束形成算法可以分为二级优化处理过程。第一级是对均匀线阵进行多波束运算,在空间中覆盖整个输出波束的K个子波束;第二级是对子波束中最接近期望信号方向的子波束,通过实际的期望信号方向和估计期望信号方向到达阵元的迭代的延时估计值误差与设定的阈值进行比较,最终获得主瓣指向期望信号的多波束方向图。在第一系统级中,信号源X(f, θ)被具有M个阵元的均匀线阵接收,并进行K(KM)点的离散傅里叶变换形成多波束,其中第k个波束的表达式为:

      $$ {{\boldsymbol{Y}}}({f_k}, {\theta _k}) = {{\boldsymbol{w}}}_k^{\rm{H}}X(f, \theta ) = \sum\limits_{k = 1}^K {\sum\limits_{m = 1}^M {{x_m}({f_k}, \theta )} {{\rm{e}}^{{\rm{j}}(m-1){\varphi _k}}}}\\ \;\;\;\;\;k = 1, 2 \cdots, K,\;\;m = 1, 2, \cdots, M $$ (13)

      式中,${x_m}({f_k}, \theta)$为第m个阵元接收的第k个子频率带的信号;θ为期望信号方向;${\varphi _k} = k \cdot 2\pi /M$,为第k个波束的相位量,则该子波束的方向为:

      $$ {\theta _k} = \arcsin \frac{{k\lambda }}{{Md}} $$ (14)

      假设一个16阵元的均匀线阵,阵元间距为λ/2,离散傅里叶波束形成器输出16个波束,分别为k=-8, -7, -6, …, 6, 7,则通过式(14)可以得到16个波束中每个波束相应的方位角:90°, ±61.04°, ±48.59°, ±38.68°, ±30°, ±22.02°, ±14.48°, ±7.18°, 0°。因此,对于实际期望方向与多波束形成的方向有误差时,可以在阵元之间插入相移动时波束图的主瓣方向对准期望信号方向,但相移是频率函数,因而多用于对窄带信号源的波束形成处理中。对于宽带信号源,可以在数据处理中插入时间延迟估计值来补偿已知的权重,形成波束的主瓣方向与实际期望信号的方向一致。设第k个子波束的主瓣方向最接近期望信号的方向,该子波的波束图为Y(fk, θk),参考波束与期望波束的主瓣响应之间的相位差为ξ。对于已知的波束主瓣区域为[θNL, θNR],其中θNLθNR分别为主瓣区域的左角度和右角度。对于子波第i+1次迭代后的波束图为Y*(fk, θi),其中${\theta _i} \in [({{\theta} _{\rm{NR}}} + {{\theta }_{\rm{NL}}})/2, \theta]$,i=1, 2, …I为在该区间离散化值,波束图的迭代调整为:

      $$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\mathit{\boldsymbol{Y}}_{\mathit{i + }{\rm{1}}}^*\left( {{f_k},{\theta _\mathit{i}}} \right) = \\ \left\{ \begin{array}{l} \mathit{\boldsymbol{Y}}_\mathit{i}^*\left( {{f_k},{\theta _\mathit{i}}} \right)\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\left| {{\tau ^*} - \tau } \right| \le \varepsilon \;\;\;\;\;(主瓣区域)\\ {K_m}\left[ {\mathit{\boldsymbol{Y}}_i^*\left( {{f_k},{\theta _\mathit{i}}} \right) - \mathit{\boldsymbol{Y}}\left( {{f_k},{\theta _\mathit{i}}} \right)} \right]\;\;\;\;\;\;\;\varepsilon < \left| {{\tau ^*} - \tau } \right|\;\;\;\;\;\left( {旁瓣区域} \right) \end{array} \right. \end{array} $$ (15)

      式中,τ*为在子波速方向的信号到达阵列的延时估计值;τ为期望信号方向到达阵元的延时估计值;ε为阈值,由τ*τ之间的允许误差值确定;Km是波束主瓣区域的迭代增益。根据最小均方误差准则使输出估计误差ε的均方值最小,通过拉格朗日乘子法可得:

      $$ {{\boldsymbol{R}}_{Y_{i + 1}^*}}{{\boldsymbol{w}}_{i + 1}} = {{\boldsymbol{r}}_{{Y_{i + 1}}Y_{i + 1}^*}} $$ (16)

      其中形成的新的子波束图的协方差矩阵${{{\boldsymbol{R}}}_{Y_{i + 1}^*}}$和互相关向量${{{\boldsymbol{r}}}_{{Y_{i + 1}}Y_{i + 1}^*}}$分别为:

      $$ {{{\boldsymbol{R}}}_{Y_{i + 1}^*}} = \sum\limits_{i = 1}^I {\sum\limits_{k = 1}^K {{\boldsymbol{Y}}_{i + 1}^*({f_k}, {\theta _i}){{\boldsymbol{a}}}({f_k}, {\theta _i}){{{\boldsymbol{a}}}^{\rm{H}}}({f_k}, {\theta _i})} } $$ (17)
      $$ {{{\boldsymbol{r}}}_{{{\boldsymbol{Y}}_{i + 1}}Y_{i + 1}^*}} = \sum\limits_{i = 1}^I {\sum\limits_{k = 1}^K {{{\boldsymbol{Y}}}_{i + 1}^{*'}({f_k}, {\theta _i}){\boldsymbol{Y}}({f_k}, {\theta _i}){{\boldsymbol{a}}}({f_k}, {\theta _i})} } $$ (18)

      式中,a为方向导向矢量。

      根据式(16)解得第i+1次迭代的波束形成加权向量为:

      $$ {{\boldsymbol{w}}_{i + 1}} = {\boldsymbol{R}}_{{Y_{i + 1}}}^{-1}{{\boldsymbol{r}}_{{Y_{i + 1}}Y_{i + 1}^*}} $$ (19)

      本文提出的波束形成算法步骤为:

      1)通过多波束形成算法式(14)在空间形成多个子波束,寻找最接近期望信号方向的子波束,对相位误差区域离散化,初始化Km, ε

      2)根据式(15)调整子波束图Y*(fk, θi),通过式(17)和式(18)得到新的协方差矩阵RY*i+1和互相关向量${{{\boldsymbol{r}}}_{{Y_{i + 1}}Y_{i + 1}^*}}$;

      3) 通过式(19)计算新的加权向量和波束图;

      4) 计算实际的期望信号方向和估计期望信号方向到达阵元的时间延迟估计值之间的误差,并设定允许的误差值为ε,如果误差值小于ε则停止迭代,否则转到步骤2)。

    • 首先对宽带线性调频接收信号进行时间延迟估计,假设阵列接收到的信号带宽为B=200kHz,载波频率为f0=1.2MHz,系统的采样频率为fs=2B,采样点数为1 024点,信号函数表达式为:

      $$ x\left( t \right) = \frac{1}{{\sqrt T }}\sin c\left( {t/T} \right){{\rm{e}}^{{\rm{j2 \mathsf{ π} }}\left( {{f_0}\mathit{t + K}{\mathit{t}^\mathit{2}}/2} \right)}} $$ (20)

      式中,K=B/T,设系统接收到的信号是由3个多径信号相加构成:τ1=104Tsτ2=186Ts, τ3=240Ts每个信号的幅度分别为a1=1, a2=1.1, a3=0.5。当信噪比为15 dB时,小波滤波器分解次数为4次。图 2为基于小波函数的时间延迟估计误差曲线。

      图  2  时间延迟估计误差曲线

      图 2看以出,小波滤波器组在第4级时的估计时间延迟τp分别为10、18、24 ns,信号的时间延迟为τ*=2×P×2jτp(j为小波函数级数为4,P为滤波器的阶数8),可见基于多通道小波滤波组的时延估计为(0.131 2×2),(0.233 6×2), (0.300 8×2)μs由此可见基于多通道小波滤波器组的时延估计方法能正确估计出信号的实际时延。设多波束形成系统为一个8阵元的均匀线阵,阵元间距为λ/2,K=8,形成8个子波束,角度区域为[-90°, 90°],通过式(14)得到每个波束的方位角为90°, ±48.59°, ±30°, ±14.48°, 0°。图 3为该多波束算法的波束图。信噪比为SNR=30dB的期望信号入射方向q=40°。对子波束的主瓣方向与实际期望信号方向的误差角度区域以间隔1°进行离散化,得到N=8个离散角度,主瓣迭代增益Km=0.1,主瓣区域设计能容忍的误差最大值ε=0.05μs,根据给出的方法进行加权迭代。图 4a为基于小波时间延迟的波束形成算法的波束图,图 4b为波束形成误差迭代图。

      图  3  多波束算法的波束图

      图  4  基于小波时间延迟的波束形成算法

      图 3可知,由于经过离散傅里叶运算子波束呈sinx/x的形状,具有良好的正交性,但图中的子波束的指向性相对实际的期望信号具有一定的相位误差,需要进行修正。

      图 4a可知,经过引入加权函数因子,通过不断迭代加权函数进一步使子波束的主瓣方向优化拟合期望信号方向,通过迭代计算后选择的子波束能够对期望信号有很好的指向性,但其他子波束之间仍具有良好的正交性,从而可以获得整个信号空间的波束图。从图 4b看出经过80次迭代后时延误差小于给定阈值。

    • 本文基于小波多分辨方法的时间延迟估计算法,通过小波的分解与重构算法在不同的分辨率级别上进行时延估计,以提高时延估计的准确性和算法的速度。本文提出通过引入加权函数因子,对每次迭代产生的主波束的时间延迟估计值和预存期望方向主波束的时间延迟估计值进行误差比较的多波束形成算法,通过对加权函数的不断迭代得到期望的主瓣方向,从而在空间中形成整体的自适应波束图。

参考文献 (18)

目录

    /

    返回文章
    返回