技术领域
[0001] 本
发明属于阵列
信号处理中的波束合成领域,主要涉及宽带波束形成的
硬件实现资源消耗分析及其降低方法。
背景技术
[0002] 在传统的阵列
信号处理技术中,以均匀线阵为例可知在M个阵元的阵列每个阵元接收的信号为:
[0003]
[0004] 对于窄带信号s(n)≈s(n-τm),每个阵元接收信号可变为 定义空间
相位φ=ωsτ=2πf0dsinθ/c=2πdsinθ/λ。将阵列信号写成矩阵形式:
[0005]
[0006] 定义列向量x(n)=[x0(n) x1(n) … xM-1(n)]Τ,a(θ)=[1 e-jφ … e-j(M-1)φ]Τ。阵列接收公式可变换为
[0007] 窄带波束形成的过程可以理解为,把信号到每一个阵元的载波引入的相位补偿好,就可以得到高增益的波束合成。引入权向量w=a(θ)Η,最后得到波束合成的结果[0008] 然而,对于宽带波束形成的情况,s(n)≠s(n-τm)。仅仅补偿载波的相位是不够的,需要对信号本身的延时进行补偿。通过延迟
滤波器h(n)对s(n-τm)补偿延时,使s(n)=h(n)*s(n-τm)。将τm分为整数部分和小数部分τm=τi_m+τf_m。延时滤波器也分为整数延迟滤波器和分数延迟滤波器h(n)=hi(n)*hf(n)。分数延迟滤波器通过FARROW结构实现。FARROW分数延迟滤波器在补偿
相位延迟的过程中确定了波束的指向。这样权向量w=[1 1 … 1]T,即仅需将每个阵列补偿后的信号相加,可得到很高的幅度增益。
[0009] 在实际工程应用中,宽带波束形成方法如下:将收到的实信号xrx_m(n)进行数字下变频变化为I/Q两路的基带信号s(n-τ)=si(n)+i·sq(n),对应为复数信号;将得到的复数信号进行相位补偿和时间延迟,时间延迟分为整数延迟和分数延迟,分数延迟利用分数延迟滤波器实现;通过分数延迟后恢复了基带信号s(n),将M个阵列的信号通过权向量加成,增加了波束增益。
[0010] 在阵列信号处理中,数字接收机做波束合成时,硬件实现上首先要求通过数字下变频的方法,把接收到的信号转化为I/Q两路的基带信号。I/Q两路代表
算法上复数的
实部和
虚部。实现上要通过许多复数运算,复数运算的加法比实数运算的加法需要多一个加法器,复数运算的乘法比实数运算的乘法多用至少两个乘法器。以一个9×9阶的FARROW延迟滤波器为例子,FARROW滤波器的系数为实数,输入为复数。一个复数乘以实数运算需要两个乘法器,一个需要9×9×2=162个乘法器,9×8×2=144个加法器。
[0011] 宽带波束形成中使用的分数延迟FARROW结构本身会耗费大量硬件资源,接收的阵列信号,变为I/Q两路后代表复数信号,再通过FARROW,将耗费更多的资源。对处理方式进行改变,以降低资源消耗十分有必要。
发明内容
[0012] 本发明的目的在于提供一种降低硬件资源消耗的宽带波束形成设计方法,完全针对于硬件工程实现,在实现宽带波束形成时,可以降低系统所消耗的资源。具体表现在,将舍去传统宽带波束形成中使用的复数运算处理,对实数进行运算。将接受到的实信号不首先进行数字下变频变成I/Q两路,而是直接对接收的实数信号进行处理。这样可以减少宽带波束形成时的FARROW结构所消耗的资源。最后再将实信号转化为复信号。这样也不会引起实信号的相位模糊问题,不会对后续的信号处理造成影响。
[0013] 本发明的思路是:对于有M个阵元接收的阵列信号接收系统,我们接收到的射频实信号记为 接收机第一步是通过模拟变频,将接
射频信号搬移到较低频,经AD变换为
数字信号记为 不通过数字DDC将
信号恢复成 的复数信号,而是直接通过延时补偿滤波器。延时补偿
滤波器将时间延迟τm补偿掉的同时,会在较低频引入新的相位变化量。此时信号表达式记为 在这一步中,由于进行实数运算,FARROW结构实现的硬
件资源将减少一半。在完成延时补偿后,我们再进行数字下变频,将信号变为I/Q两路恢复成复数基带信号 最后,我们引入权向量
其中φm=(ωs-ω0)τm,权向量w可将FARROW引入的相位一起补偿掉。得到宽带波束合成的H
结果y(n)=wx(n)=Ms(n)。
[0014] 一种降低硬件资源消耗的宽带波束形成设计方法,具体步骤如下:
[0015] S1、对于有M个阵元任意共形阵阵列,建立信号接收模型。第m个阵元的坐标表示pm=[xm ym zm],选择坐标原点为基准点。某信号s(t)以入射方向 入射方向向量记为第m个阵元接收到的射频实信号为其中τm是 第m个阵元相对于参考点的延时 ,
[0016] S2、通过模拟前端将射频端信号变为较低频信号为将中频信号通过高速AD
采样,较低频
模拟信号变为数字中
频实信号,信号表达为:
[0017] S3、延时τm分为两部分,整数部分和小数部分,即τm=τi_m+τf_m,其中τi_m为整数部分,τi_m为小数部分。对其进行延时补偿,分为整数延时补偿和小数延时补偿(分数延迟)。如下S31与S32所示,两个步骤可以互换。
[0018] S31、对信号进行整数延时补偿,通过对信号延迟τi_m个采样点即可补偿。整数延时补偿后信号变为
[0019] S32、对整数延时补偿后的信号,进行小数部分延时补偿。通过FARROW结构分数延时滤波器,信号变为:
[0020]
[0021] S4、对延时补偿后的信号进行数字下变频,将实信号转化为I/Q两路表示的复数基带信号。I/Q两路的信号综合表示为:
[0022]
[0023] 其中, 为I路, 为Q路。
[0024] S5、设置波束合成的方向为 引入权向量 其中φm=(ωs-ω0)τ1_m, 当波束合成方向
与信号入射方向相同时(即θ=θ1, 时),权向量与数字下变频后的信号相乘得到具有高增益的接收信号y(n)=wHx(n)=Ms(n),完成宽带波束合成。
[0025] 本发明的有益效果是:
[0026] 对于接收到的实数信号直接进行处理,可以在延时补偿时,将使用的FARROW结构滤波器减少一半,硬件上将节约一半的资源。同时,在通过延时补偿后,再进行数字DDC变成I/Q两路的复数信号,保留了信号的全部信息,不会造成相位模糊。
附图说明
[0028] 图2是一种均匀线阵示意图
[0029] 图3是一种线阵宽带波束形成的
能量累积方向图。
[0030] 图4是一种均匀面阵示意图
[0031] 图5是一种面阵宽带波束形成的能量累积方向图。
具体实施方式
[0032] 下面结合
实施例和附图,详细说明本发明的技术方案。
[0033] 如图1所示:
[0034] S1、对于有M个阵元任意共形阵阵列,建立信号接收模型。第m个阵元的坐标表示pm=[xm ym zm],选择坐标原点为基准点。某信号s(t)以入射方向 入射方向向量记为第m个阵元接收到的射频实信号为其中τm是第m个阵元 相对于参考点的 延时 ,
[0035] S2、通过模拟前端将射频端信号变为较低频信号为将中频信号通过高速AD采样,较低频模拟信号变为数字中
频实信号,信号表达为:
[0036] S3、延时τm分为两部分,整数部分和小数部分,即τm=τi_m+τf_m,其中τi_m为整数部分,τi_m为小数部分。对其进行延时补偿,分为整数延时补偿和小数延时补偿(分数延迟)。如下S31与S32所示,两个步骤可以互换。
[0037] S31、对信号进行整数延时补偿,通过对信号延迟τi_m个采样点即可补偿。整数延时补偿后信号变为
[0038] S32、对整数延时补偿后的信号,进行小数部分延时补偿。通过FARROW结构分数延时滤波器,信号变为:
[0039]
[0040] S4、对延时补偿后的信号进行数字下变频,将实信号转化为I/Q两路表示的复数基带信号。I/Q两路的信号综合表示为:
[0041]
[0042] 其中, 为I路, 为Q路。
[0043] S5、设置波束合成的方向为 引入权向量 其中φm=(ωs-ω0)τ1_m, 当波束合成方向
与信号入射方向相同时(即θ=θ1, 时),权向量与数字下变频后的信号相乘得到具有高增益的接收信号y(n)=wHx(n)=Ms(n),完成宽带波束合成。
[0044] 实施例1
[0045] 本发明是通过新的方法达到降低硬件资源消耗,同时保证宽带波束形成与传统方法结果一致。如附图2所示为均匀线阵示意图,由M=36个阵元构成的均匀线阵接收1个远场信源发射的宽带信号。远场信源为LFM(线性调频)信号,记为其中, ωs=2πfs。记波束合成的指向为
θ1=30°,
[0046] 具体如下:
[0047] ①、对于有36个阵元的均匀线阵,建立信号接收模型。第m个阵元的坐标表示pm=[xm 0 0]=[(m-1)λ/2 0 0],λ为信号
波长,且λ=cf0。阵元间距设置为信号半波长,所以xm=(m-1)λ/2。选择坐标原点为基准点,同时也是第一个阵元
位置。某信号s(t)以入射方向(θ,0),θ=30°,入射方向向量记为r=-[sinθ 0 cosθ]。第m个阵元接收到的射频实信号为其中τm是第m个阵元相对于参考点的延时,τm=-1/c(xm sinθ)=-(m-1)f0/4;
[0048] ②、通过模拟前端将射频端信号变为较低频信号为将较低频信号通过高速AD采样,较低频模拟信号变为数字
较低频实信号,信号表达为:
[0049] ③、首先对信号进行整数延时补偿,通过对信号延迟τi_m个采样点即可补偿。整数延时补偿后信号变为 然后对整数延时补偿后的信号,进行小数部分延时补偿。通过FARROW结构分数延时滤波器,此处取11阶FARROW滤波器效果较好,信号变为
[0050] ④、对延时补偿后的信号进行数字下变频,将实信号转化为I/Q两路表示的复数基带信号。I/Q两路的信号综合表示为:
[0051]
[0052] ⑤、设置波束合成的方向为 引入权向量 其中φm=(ωs-ω0)τ1_m, 当波束合成方向与信号入
射方向相同时(即θ=θ1, 时),权向量与数字下变频后的信号相乘得到具有高增益的接收信号y(n)=wHx(n)=36s(n),完成宽带波束合成。
[0053] 如附图3所示,通过扫描设置不同的波束方向
角,得到了本实施例宽带波束形成的输出能量累积方向图。通过此图,可以看出宽带波束合成的方向指向。图中所显示的方向主瓣与我们预设值的θ=30°一致,证明了本方法的可行性。
[0054] 同时,在此实施例中,实施步骤S3将实数信号直接通过11阶的FARROW滤波器。此处消耗的乘法器为11×11=121,加法器个数为11×10=110。而用传统方法,将接收信号变为基带复数信号,复数信号通过11阶FARROW滤波器,将消耗乘法器为11×11×2=242,加法器为11×10×2=220。本方法将比传统宽带波束形成技术节约一半的资源。
[0055] 实施例2
[0056] 本发明的又一个实施例,是通过均匀面阵来验证本发明能在得到同样波束形成效果的情况下减少一半资源。如附图4示意图表示均匀面阵,该面阵由M×N=144个阵元构成,其中M=12,N=12。面阵接收1个远场信源发射的宽带信号。远场信源为LFM(线性调频)信号,记为 其中, ωs=2πfs。记波束合成的指向为θ1=20°,
[0057] 具体如下:
[0058] ①、对于有M×N个阵元的均匀面阵建立信号接收模型,其中M=12,N=12。第mn个阵元的坐标表示pm,n=[xm yn 0]=[(m-1)λ/2 (n-1)λ/2 0],λ为信号波长,且λ=cf0。行和列之间,阵元间距均设置为信号半波长,xm=(m-1)λ/2,yn=(n-1)λ/2。选择坐标原点为基准点,同时也是第一个阵元位置。某信号s(t)以入射方向 θ=20°, 入射方向向量记为 第mn个阵元接收到的射频实信号为
其中τm是第mn个阵元的阵元相对于参考点的延时,
[0059] ②、通过模拟前端将射频端信号变为较低频信号为将较低频信号通过高速AD采样,较低频模拟信号变为数字
较低频实信号,信号表达为:
[0060] ③、首先对信号进行整数延时补偿,通过对信号延迟τi_m个采样点即可补偿,整数延时补偿后信号变为 然后对整数延时补偿后的信号,进行小数部分延时补偿。通过FARROW结构分数延时滤波器,信号变为
此处选用15阶FARROW滤波器即可得到良好的效果。
[0061] ④、对延时补偿后的信号进行数字下变频,将实信号转化为I/Q两路表示的复数基带信号。I/Q两路的信号综合表示为
[0062] ⑤、设置波束合成的方向为 引入权向量 其中φm=(ωs-ω0)τ1_m, 当波束合成方向与信号入
射方向相同时(即θ=θ1, 时),权向量与数字下变频后的信号相乘得到具有高增益的接收信号y(n)=wHx(n)=144s(n),完成宽带波束合成。
[0063] 如附图5所示,通过扫描设置不同的波束方向角,得到了本实施例宽带波束形成的输出能量累积方向图。通过此图,可以看出宽带波束合成的方向指向。图中所显示
俯仰角和方位角的方向主瓣都与我们预设值的θ=20°, 一致,证明了本方法的可行性。
[0064] 同时,在此实施例中,实施步骤S3将实数信号直接通过15阶的FARROW滤波器。此处消耗的乘法器为15×15=225,加法器个数为15×14=210。而用传统方法,将接收信号变为基带复数信号,复数信号通过15阶FARROW滤波器,将消耗乘法器为15×15×2=450,加法器为15×14×2=420。本方法将比传统宽带波束形成技术节约一半的资源。