CN101296014A - 一种宽带数字波束形成方法 - Google Patents
一种宽带数字波束形成方法 Download PDFInfo
- Publication number
- CN101296014A CN101296014A CNA2008100624796A CN200810062479A CN101296014A CN 101296014 A CN101296014 A CN 101296014A CN A2008100624796 A CNA2008100624796 A CN A2008100624796A CN 200810062479 A CN200810062479 A CN 200810062479A CN 101296014 A CN101296014 A CN 101296014A
- Authority
- CN
- China
- Prior art keywords
- signal
- output
- filtering
- sub
- filtration module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000001914 filtration Methods 0.000 claims abstract description 253
- 238000005070 sampling Methods 0.000 claims abstract description 81
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 33
- 230000006798 recombination Effects 0.000 claims description 16
- 238000005215 recombination Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000013461 design Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000036541 health Effects 0.000 claims description 6
- 230000008676 import Effects 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 230000000474 nursing effect Effects 0.000 claims description 6
- 230000008521 reorganization Effects 0.000 claims description 6
- 230000011664 signaling Effects 0.000 claims description 6
- 238000012856 packing Methods 0.000 claims description 5
- 230000003111 delayed effect Effects 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 37
- 239000000047 product Substances 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000013068 control sample Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种宽带数字波束形成方法。该方法包括时分复用中频正交采样、正交解调、低通滤波、内插延时与相位旋转以及波束求和等步骤。本发明的波束形成方法着重于用单个采样电路对多路带通信号进行时分复用中频正交采样,用FIR滤波器一次性完成与信号匹配的低通滤波、内插延时和相位旋转,实现真正的延时波束形成。在滤波的执行方式上,公开了一种级联结构的多通道滤波器,同一通道信号的滤波在不同的子滤波模块中级联完成,不同通道信号共享同一滤波模块,滤波输出的同时完成样本抽取,适合于在ASIC或FPGA中实现,较传统滤波器节省大量逻辑资源。
Description
技术领域
本发明涉及雷达、声纳及无线通信领域,特别涉及宽带数字波束形成方法。
背景技术
波束形成是广泛应用于雷达、声纳、通信、航天等领域的一项技术。波束形成器通常与声纳发射或接收阵列协同工作,同样它也可以与电磁波发射或接收阵列协同工作。
为了让传感器阵列指向一个或多个特定的方向,加强从某个或多个方向的信号,基本的波束形成方法是对每一路阵元接收的信号按该波束方向进行延时,然后再把补偿至同相位的各路信号相加,得到指向该方向的波束输出(参见图1)。近年来,现代声定位、高速水声通讯以及声成像等技术领域为了得到更高的性能,普遍采用宽带信号机制。用简单的相移法对宽带信号进行波束形成是不适用的,在时域可以采用延时的方法来处理,在频域主要应用宽带信号的FFT方法来实现多波束输出。本发明涉及的主要是用时域的延时方法来进行接收波束形成。
一个多波束形成***往往需要承担巨大的运算量,对计算速度和计算并行度提出了很高的要求。对于多波束形成技术的实现,一个关键的技术问题就是希望能够将波束精确地定向在期望的方向上,同时又希望输入信号的采样率不致过高,以减轻软、硬件的压力。在数字***中,常规的“延时求和”的波束形成方法通过样本序列的移位相加实现,而延时的精度与数据采样率成正比。为了得到较高的波束定向精度和良好的波束形成性能需要对阵元序列进行精确的延时控制,因此需要很高的采样频率,同时也大大增加了波束形成的计算量。随着信号频率的提高和阵列规模的扩大,波束形成***的复杂度(包括数据存储、传输及运算速率等)将急剧上升,这是制约数字波束形成***实现的主要障碍。
针对以上矛盾,发展了内插波束形成方法,这种方法虽然可以降低对采样速率的要求,但代价是增加了补零内插的操作,增加了额外的运算量。这种波束形成器的结构如图2所示,先对各路阵元信号的相邻两采样点之间等间隔***(D-1)个零值点,然后各自通过进行低通滤波,得到一个采样频率更高的信号,以上两步执行了内插增速采样。然后根据波束方向对内插后的各路信号进行相应的“延时求和”,此时的延时会比利用原始信号更精确,因为现在的最小延时单位是采样周期的1/D。
目前在波束形成***中一种比较常用的方法是正交采样波束形成,包括四个基本的步骤:
1)正交采样过程:将带通波形变换到复基带,得到复包络信号;
2)时延复包络信号,补偿各路复包络信号的相对延时;
3)相位旋转:对不同延时,补偿频谱搬移(正交采样)所产生的相位偏差;
4)迭加各通道信号,得到复包络波束输出。
这一波束形成方法可以将带通信号变换到基带进行处理,由于复基带信号的奈奎斯特频率等于信号的带宽,对于信号中心频率与带宽之比较大的***,可以显著降低波束形成***的采样频率和运算量。对于步骤2中的延时可以采用前述补零内插的方法来处理。
在当今大多数声纳***当中,载波频率一般在几千赫兹到几百千赫兹,而ADC器件的工作频率能达到几十兆赫兹,因此在波束形成器中可以采用时分复用采样的结构。
本发明中的宽带数字波束形成方法是针对声纳***的应用环境描述的,但其原理在雷达***中亦是相通的。
发明内容
本发明的目的是克服现有技术的不足,提供一种宽带数字波束形成方法。
一种宽带数字波束形成方法包括如下步骤:
1)对N路接收阵元的带通信号xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct)进行前置放大和抗混叠滤波等信号调理,式中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交信号;
2)用单个采样电路和ADC器件对经过调理后的N路信号进行时分复用中频正交采样、模/数转换,然后进行正交解调,合成一路数字信号,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列;
3)对步骤2)输出的序列进行周期性符号反转,使之补偿成为{Ii,Qi,Ii,Qi,Ii…}序列;
4)把从步骤3)输出的已经合成为一路数据流的各路信号同时输入级联结构的多通道FIR滤波器I和滤波器Q,分别进行同相和正交分量的滤波,一次性完成与信号匹配的低通滤波、内插延时和相位旋转,滤波输出的同时完成样本抽取,每路信号的延时不仅要补偿波束到达各个接收阵元产生的相对延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外相对延时,滤波器I和滤波器Q的每个通道滤波系数由内插系数 和相位旋转系数合成,滤波器I、Q的系数vi T、wi T以及滤波器I、Q的输出IDi(nTO)、QDi(nTO)分别为:
其中xi是第i路信号序列,τi是该路信号的延时,TO是复包络信号的输出周期,IDi(t)和QDi(t)分别是延时至波前对齐的基带复包络信号的同相和正交分量,ai T和bi T分别是同相分量和正交分量的延时内插系数组,该系数是补零内插方法的低通滤波系数中对应某一内插点的系数组;
5)对从步骤4)输出的经过延时补偿至波前对齐的N路基带复包络信号进行波束求和,序列输出基带上波束形成的结果。
所述的时分复用中频正交采样:在符合带通信号采样准则的前提下设置采样频率为接收带通信号中心频率的4N/(4k+1)倍或者4N/(4k+3)倍,其中k=0,1,2…,等效为每路信号的采样频率为接收带通信号中心频率的4/(4k+1)倍或者4/(4k+3)倍,其中k=0,1,2…,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列。
所述的级联结构的多通道FIR滤波器I和滤波器Q:由信号重组模块和若干个子滤波模块组成,每个子滤波模块阶数为2M,分别对每路信号完成2M阶的滤波计算,滤波器输出信号的采样频率等于原始信号采样频率的1/2M,按照原始信号的采样频率和降采样后所需的采样频率来确定子滤波器模块的阶数,子滤波模块个数由所需的滤波阶数和子滤波模块的阶数共同确定,滤波过程中采样信号从第一级子滤波模块依次流向最后一级子滤波模块,子滤波模块对每一个通道的滤波计算则从最后一级子滤波模块开始,并把每一个通道的滤波结果依次往前一级子滤波模块传递,在后一级子滤波模块输出某一通道信号的滤波结果到前一级子滤波模块之后,前一级子滤波模块开始该通道信号的滤波,最终在第一级子滤波模块处完成整个滤波计算,输出滤波结果。
所述的信号重组模块:把按采样所得的信号排列顺序转换成2M个同一通道的信号连续排列的顺序并依次输出,重组后的采样信号输入到第一级子滤波模块数据RAM的输入端;
所述的子滤波模块:包括一个存储采样信号的双端口数据RAM、存储系数的双端口系数RAM、乘法器、加法器以及二路选择器,乘法器的两个输入端连接数据RAM和系数RAM的输出端,完成采样信号和对应滤波系数的乘积,并把乘积结果输出至加法器,乘法器的输出端和二路选择器的输出端连接到加法器的两个输入端,加法器的输出端同时连接到本级与前一级二路选择器的输入端,本级二路选择器的两个输入端分别连接本级和后一级加法器的输出端,最后一级二路选择器的两个输入分别是本级加法器的输出和常数零,加法器和二路选择器协同工作,在每路通道信号滤波的第一个周期二路选择器选通后一级加法器的输出同本级乘法器的输出相加,如此完成最后一级子滤波模块至本级子滤波模块滤波结果的累加,在剩余周期内二路选择器全部选通本级加法器的输出同乘法器输出进行相加。
所述的数据RAM:被划分为N个连续的数据块,每个数据块连续存储2M个同一通道的采样信号,所有子滤波模块数据RAM的初始读取地址都为第1个数据块的起始处地址,初始写入地址都为第N个数据块的起始处地址,读取地址和写入地址在数据RAM空间内循环递增,本级数据RAM的输出端连接到本级乘法器的输入端和后一级数据RAM的输入端,第一级数据RAM的输入端连接信号重组模块的输出端。
所述的内插系数的求取方法:该系数是从补零内插方法的低通滤波系数中获取的,内插等价于对补零信号进行线性时不变滤波,通过内插将一个序列转换到一个较高的等效采样频率的内插序列,进行D倍的补零内插首先要在原始序列x(n)相邻两采样点之间等间隔***(D-1)个零值点,然后对补零后的序列xp(n)进行低通滤波,得到内插序列xD(n):
对应于某一个内插点xD(n)的插值系数只要选择计算该内插点时原始信号所对应的低通滤波系数h(n-rD)(-∞<r<+∞),即为对应该延时的一个内插系数组,若设计该低通滤波系数的幅频特性和复包络信号匹配,则在内插延时的同时完成了与信号匹配的低通滤波。
另一种宽带数字波束形成方法包括如下步骤:
1)对N路接收阵元的带通信号xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct)进行前置放大和抗混叠滤波等信号调理,式中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交信号;
2)用单个采样电路和ADC器件对经过调理后的N路信号进行时分复用中频正交采样、模/数转换,然后进行正交解调,合成一路数字信号,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列;
3)对步骤2)输出的序列进行周期性符号反转,使之补偿成为{Ii,Qi,Ii,Qi,Ii…}序列;
4)把从步骤3)输出的已经合成为一路数据流的各路信号输入级联结构的多通道FIR滤波器I,进行同相分量的滤波,一次性完成与信号匹配的低通滤波、内插延时、相位旋转和带通调制,滤波输出的同时完成样本抽取,每路信号的延时不仅要补偿波束到达各个接收阵元产生的相对延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外相对延时,滤波器I的每个通道滤波系数由内插系数 和相位旋转系数合成,滤波器I的系数vi T以及滤波器I的输出IDi(nTO)为:
上式中xi是第i路信号序列,τi是该路信号的延时,TO是通带信号的输出周期,IDi(t)是延时至波前对齐的基带复包络信号的同相分量,ai T和bi T分别是同相分量和正交分量的延时内插系数组,该系数是补零内插方法的低通滤波系数中对应某一内插点的系数组;
5)对从步骤4)输出的经过延时补偿至波前对齐的N路带通信号进行波束求和,序列输出通带上波束形成的结果。
所述的时分复用中频正交采样:在符合带通信号采样准则的前提下设置采样频率为接收带通信号中心频率的4N/(4k+1)倍或者4N/(4k+3)倍,其中k=0,1,2…,等效为每路信号的采样频率为接收带通信号中心频率的4/(4k+1)倍或者4/(4k+3)倍,其中k=0,1,2…,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列。
所述的级联结构的多通道FIR滤波器I:由信号重组模块和若干个子滤波模块组成,每个子滤波模块阶数为4M,分别对每路信号完成4M阶的滤波计算,滤波器输出信号的采样频率等于原始信号采样频率的1/4M,按照原始信号的采样频率和降采样后所需的采样频率来确定子滤波器模块的阶数,子滤波模块个数由所需的滤波阶数和子滤波模块的阶数共同确定,滤波过程中采样信号从第一级子滤波模块依次流向最后一级子滤波模块,子滤波模块对每一个通道的滤波计算则从最后一级子滤波模块开始,并把每一个通道的滤波结果依次往前一级子滤波模块传递,在后一级子滤波模块输出某一通道信号的滤波结果到前一级子滤波模块之后,前一级子滤波模块开始该通道信号的滤波,最终在第一级子滤波模块处完成整个滤波计算,输出滤波结果。
所述的信号重组模块:把按采样所得的信号排列顺序转换成4M个同一通道的信号连续排列的顺序并依次输出,重组后的采样信号输入到第一级子滤波模块数据RAM的输入端。
所述的子滤波模块:包括一个存储采样信号的双端口数据RAM、存储系数的双端口系数RAM、乘法器、加法器以及二路选择器。乘法器的两个输入端连接数据RAM和系数RAM的输出端,完成采样信号和对应滤波系数的乘积,并把乘积结果输出至加法器,乘法器的输出端和二路选择器的输出端连接到加法器的两个输入端,加法器的输出端同时连接到本级与前一级二路选择器的输入端,本级二路选择器的两个输入端分别连接本级和后一级加法器的输出端,最后一级二路选择器的两个输入分别是本级加法器的输出和常数零,加法器和二路选择器协同工作,在每路通道信号滤波的第一个周期二路选择器选通后一级加法器的输出同本级乘法器的输出相加,如此完成最后一级子滤波模块至本级子滤波模块滤波结果的累加,在剩余周期内二路选择器全部选通本级加法器的输出同乘法器输出进行相加;
所述的数据RAM:被划分为N个连续的数据块,每个数据块连续存储2M个同一通道的采样信号,所有子滤波模块数据RAM的初始读取地址都为第1个数据块的起始处地址,初始写入地址都为第N个数据块的起始处地址,读取地址和写入地址在数据RAM空间内循环递增,本级数据RAM的输出端连接到本级乘法器的输入端和后一级数据RAM的输入端,第一级数据RAM的输入端连接信号重组模块的输出端。
所述的内插系数的求取方法:该系数是从补零内插方法的低通滤波系数中获取的。内插等价于对补零信号进行线性时不变滤波,通过内插将一个序列转换到一个较高的等效采样频率的内插序列,进行D倍的补零内插首先要在原始序列x(n)相邻两采样点之间等间隔***(D-1)个零值点,然后对补零后的序列xp(n)进行低通滤波,得到内插序列xD(n):
对应于某一个特定内插点xD(n)的插值系数只要选择计算该内插点时原始信号所对应的低通滤波系数h(n-rD)(-∞<r<+∞),即为对应该延时的一个内插系数组,若设计该低通滤波系数的幅频特性和复包络信号匹配,则在内插延时的同时完成了与信号匹配的低通滤波。
本发明的波束形成方法着重于用单个采样电路对多路带通信号进行时分复用中频正交采样,用FIR滤波器一次性完成与信号匹配的低通滤波、内插延时和相位旋转,实现真正的延时波束形成。在滤波的执行方式上,公开了一种级联结构的多通道滤波器,同一通道信号的滤波在不同的子滤波模块中级联完成,不同通道信号共享同一滤波模块,滤波输出的同时完成样本抽取,适合于在ASIC或FPGA中实现,较传统滤波器节省大量逻辑资源。
附图说明
图1是“延时求和”波束形成原理示意图;
图2是内插延时的波束形成器结构示意图;
图3(a)是补零内插方法中原始信号的频谱示意图;
图3(b)是补零内插方法中补零信号的频谱示意图;
图3(c)是补零内插方法中内插信号的频谱示意图;
图4是本发明的原理示意图:由信号调理、时分复用中频正交采样、低通滤波、内插延时与相位旋转以及波束求和等步骤组成;
图5是级联的多通道滤波器结构框图;
图6是信号重组前后的信号排列方式示意图;
图7是子滤波模块结构和数据流示意图;
图8(a)是本发明一个优选实施例的波束形成流程图;
图8(b)是本发明另一个优选实施例的波束形成流程图。
具体实施方式
本发明在前述正交采样波束形成的基础上,结合时分复用中频正交采样、复解调以及内插延时等方法对带通信号进行波束形成,提出了具有以下特点的波束形成方法:1、用单个采样电路对多路带通信号进行时分复用中频正交采样,设置合适的采样频率使得在采样后能方便地进行复解调。2、采用FIR滤波器一次性完成与信号匹配的低通滤波、内插延时和相位旋转,实现真正的延时波束形成。3、内插延时的内插系数从补零内插方法的低通滤波系数中获取。4、提出了一种级联结构的多通道滤波器,同一通道信号的滤波在不同的子滤波模块中级联完成,不同通道信号共享同一滤波模块,滤波输出的同时完成样本抽取。
图8(a)和图8(b)示出了本发明的两个声纳波束形成优选实施例,该例中具有50个接收阵元,接收的水声信号为150KHz-200KHz的带通信号(中心频率为175KHz,解调后基带带宽25KHz),输出基带波束形成结果的宽带波束形成方法具体步骤如下所述:
一种宽带数字波束形成方法包括如下步骤:
1)对N路接收阵元的带通信号xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct)进行前置放大和抗混叠滤波等信号调理,式中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交信号;本例中对50路接收阵元的带通信号进行信号调理;
2)用单个采样电路和ADC器件对经过调理后的N路信号进行时分复用中频正交采样、模/数转换,然后进行正交解调,合成一路数字信号,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列;本例中用一个高速ADC器件取代50个低速ADC器件,采样电路的通道轮换选通开关和后续的ADC器件由采样脉冲控制协同工作,按照采样脉冲的频率对50路阵元信号依次轮流选通并进行模/数转换,节省了模/数转换的成本。
3)对步骤2)输出的序列进行周期性符号反转,使之补偿成为{Ii,Qi,Ii,Qi,Ii…}序列;
4)把从步骤3)输出的已经合成为一路数据流的各路信号同时输入级联结构的多通道FIR滤波器I和滤波器Q,分别进行同相和正交分量的滤波,一次性完成与信号匹配的低通滤波、内插延时和相位旋转,滤波输出的同时完成样本抽取,每路信号的延时不仅要补偿波束到达各个接收阵元产生的相对延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外相对延时,滤波器I和滤波器Q的每个通道滤波系数由内插系数 和相位旋转系数合成,滤波器I、Q的系数vi T、wi T以及滤波器I、Q的输出IDi(nTO)、QDi(nTO)分别为:
其中xi是第i路信号序列,τi是该路信号的延时,TO是复包络信号的输出周期,IDi(t)和QDi(t)分别是延时至波前对齐的基带复包络信号的同相和正交分量,ai T和bi T分别是同相分量和正交分量的延时内插系数组,该系数是补零内插方法的低通滤波系数中对应某一内插点的系数组;
5)对从步骤4)输出的经过延时补偿至波前对齐的N路基带复包络信号进行波束求和,序列输出基带上波束形成的结果。
所述的时分复用中频正交采样:在符合带通信号采样准则的前提下设置采样频率为接收带通信号中心频率的4N/(4k+1)倍或者4N/(4k+3)倍,其中k=0,1,2…,等效为每路信号的采样频率为接收带通信号中心频率的4/(4k+1)倍或者4/(4k+3)倍,其中k=0,1,2…,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列。在现有的正交采样技术中,可以通过模拟方式或数字方式实现,考虑到模拟方式在两个正交通道之间的相位通常有2°~3°的偏差,影响解调质量,故采用数字解调方式。带通信号可以表示为xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct),其中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交分量,ωc是信号的中心频率,xi(t)的复包络可表示为若控制采样频率fs=4fc/(4k+3),其中k=0,1,2…,则采样序列为{Ii,Qi,-Ii,-Qi,Ii…};若控制采样频率fs=4fc/(4k+1),其中k=0,1,2…,则采样序列为{Ii,-Qi,-Ii,Qi,Ii…}。在本例中控制采样频率fs=4N□fc=35MHz,则每路信号的等效采样频率是fs=4fc=700KHz,这样便可以进行时分复用中频正交采样,方便后续的正交解调,正交解调的精度与ADC器件的精度和采样脉冲时钟源精度成正比。
在步骤4)的滤波器I和滤波器Q系数的设计中,每路信号的延时不仅要补偿波束到达各个接收阵元产生的延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外的延时。如果把某路带通信号延时τi时间:
xi(t-τi)=Ii(t-τi)cos(ωc(t-τi))-Qi(t-τi)sin(ωc(t-τi))=IDi(t)cos(ωct)-QDi(t)sin(ωct)
其中:
延时至波前对齐的基带复包络信号记为IDi(t)和QDi(t)分别是延时至波前对齐的基带信号的同相和正交分量。每一个具体的延时τi都分为整数周期延时和分数周期延时,整数周期延时可以由序列的简单移位来得到,分数周期延时则由内插得到:
其中Ii和Qi分别代表由一定数量的Ii(t)和Qi(t)样本点组成的序列,ai T和bi T为内插系数组,从补零内插方法的低通滤波系数中获取。
所述的内插系数的求取方法:内插等价于对补零信号进行线性时不变滤波,内插的实质可以看作脉冲串采样的恢复***,通过内插将一个序列转换到一个较高的等效采样频率的内插序列。进行D倍的补零内插首先要在原始序列x(n)相邻两采样点之间等间隔***(D-1)个零值点,然后对补零后的序列xp(n)进行低通滤波,得到内插序列xD(n):
对应于某一个内插点xD(n)的插值系数只要选择计算该内插点时原始信号所对应的低通滤波系数h(n-rD)(-∞<r<+∞),即为对应该延时的一个内插系数组,如ai T和bi T即为两组内插系数组。若设计该低通滤波系数的截止频率和带外衰减等幅频特性和复包络信号匹配,则在内插延时的同时完成了与信号匹配的低通滤波。
进行宽带波束形成不仅其复包络需加以适当延时,还必须对其相位进行旋转,即Ii和Qi分别与cos(ωcτi)和sin(ωcτi)交叉相乘。经过延时至波前对齐的基带信号的同相和正交分量表示为:
内插延时和相位旋转这两步可以合起来简单地通过一个FIR滤波器来实现,避免了把带通信号分解为同相和正交分量,然后又要对经过延时的同相和正交分量交叉相乘。因为每路信号的数据是由Ii和Qi相间构成的,所以不需要把通道数据额外地分成Ii和Qi,直接让通道数据通过滤波器I和滤波器Q即可,每路信号需要FIR滤波器I和FIR滤波器Q这两个滤波器来分别处理同相和正交分量。滤波器I和滤波器Q的系数vi T和wi T由ai T和bi T和相位旋转系数合并而成:
滤波器I和滤波器Q的输出IDi(nTO)和QDi(nTO)如下:
步骤5)对从步骤4)输出的经过延时补偿至波前对齐的N路基带复包络信号进行波束求和,序列输出基带上波束形成的结果。
本例中根据接收信号的通带截止频率、过渡带宽及带外衰减等指标要求对每路信号的同相和正交分量分别进行60阶的低通滤波,同时根据波束形成的性能要求需要对采样信号进行10倍的补零内插。在补零内插方法中,进行10倍的补零内插在频谱上相当于对原始信号的频率进行了10倍压缩,如图3所示,因此需要设计600阶的低通滤波系数以供内插系数选取之用。由于本发明中采样后信号的同相和正交分量是合在一路数据流中的,因此相当于对一路信号进行120阶的滤波。考虑选用12个10阶子滤波模块进行120阶的滤波,滤波后每路信号的采样率为60KHz,符合基带信号的采样定理。显然,该滤波器输出也是各路信号交织在一起的一路数据流。
所述的级联结构的多通道FIR滤波器I和滤波器Q如图5所示:由信号重组模块和若干个子滤波模块组成,每个子滤波模块阶数为2M,分别对每路信号完成2M阶的滤波计算,由于中频正交采样后同相和正交分量相间排列,每个子滤波模块等效于对同相和正交分量分别进行M阶的滤波计算。该滤波器每隔2M个信号样本进行一次滤波计算,滤波器输出信号的采样频率等于原始信号采样频率的1/2M,由此可以按照原始信号的采样频率和降采样后所需的采样频率来确定子滤波器模块的阶数,子滤波模块个数由所需的滤波阶数和子滤波模块的阶数共同确定。本例中每个子滤波模块阶数2M等于10,即分别对每路信号完成10阶的滤波计算。
所述的信号重组模块:把按采样所得的信号排列顺序转换成2M个同一通道的信号连续排列的顺序并依次输出,参见图6,重组后的采样信号输入到第一级子滤波模块数据RAM的输入端;
所述的子滤波模块:如图7所示,包括一个存储采样信号的双端口数据RAM、存储系数的双端口系数RAM、乘法器、加法器以及二路选择器。
所述的数据RAM:如图7所示,被划分为N个连续的数据块,每个数据块连续存储2M个同一通道的采样信号,参见图6。所有子滤波模块数据RAM的读/写地址都是相同的,初始读取地址都为第1个数据块的起始处地址,初始写入地址都为第N个数据块的起始处地址,读取地址和写入地址在数据RAM空间内循环递增,本级数据RAM的输出端连接到本级乘法器的输入端和后一级数据RAM的输入端,第一级数据RAM的输入端连接信号重组模块的输出端。
根据数据RAM读/写地址的设置,每级数据RAM的读取地址比写入地址提前2M个地址,即一个连续的数据块。在读取本级数据RAM中的某个数据做乘法运算的同时把该数据按照写入地址写入下一级数据RAM,因此每个通道数据在数据RAM中的存储是每往后一级,要落后2M个地址。例如在第一级数据RAM中1通道的2M个数据存储在第N个数据块,2通道的2M个数据存储在第1个数据块,那么第二级数据RAM中1通道的2M个数据存储在第N-1个数据块,2通道的2M个数据存储在第N个数据块,依次类推。各级数据RAM同一通道信号为不同时刻的采样信号,第一级数据RAM为最近的采样信号,最后一级数据RAM则为最早的采样信号。按照数据RAM的存储特性,各级系数RAM的各通道系数亦要做对应配置,使得每个采样信号都能和正确对应的滤波系数相乘,从最后一级系数RAM到第一级系数RAM的同一通道的所有系数组成该通道的滤波系数。
乘法器的两个输入端连接数据RAM和系数RAM的输出端,完成采样信号和对应滤波系数的乘积,并把乘积结果输出至加法器,乘法器的输出端和二路选择器的输出端连接到加法器的两个输入端,加法器的输出端同时连接到本级与前一级二路选择器的输入端,本级二路选择器的两个输入端分别连接本级和后一级加法器的输出端,最后一级二路选择器的两个输入分别是本级加法器的输出和常数零,加法器和二路选择器协同工作,在每路通道信号滤波的第一个周期,上一级子滤波模块刚好完成该路信号的滤波运算并将结果输出至本级二路选择器,本级二路选择器此时选通后一级加法器的输出同本级乘法器的输出相加,如此完成最后一级子滤波模块至本级子滤波模块滤波结果的累加。在剩余周期内二路选择器全部选通本级加法器的输出同本级乘法器输出进行相加。
滤波过程中采样信号从第一级子滤波模块依次流向最后一级子滤波模块,子滤波模块对每一个通道的滤波计算则从最后一级子滤波模块开始,并把每一个通道的滤波结果依次往前一级子滤波模块传递,在后一级子滤波模块输出某一通道信号的滤波结果到前一级子滤波模块之后,前一级子滤波模块开始该通道信号的滤波,最终在第一级子滤波模块处完成整个滤波计算,输出滤波结果。
如前例中所述的一种宽带数字波束形成方法,如果其中FIR滤波器的子滤波模块阶数为4M,则滤波输出完成了4M倍的降采样,因此每路信号输出的时间间隔是采样信号时间间隔的4M倍,即TO=4MTs=4M/fs。对于fs=4fc/(2k+1),TO=(2k+1)M/fc,cos(ωcnTO)=1,sin(ωcnTO)=0,则在通带上的波束形成输出为:
1)对N路接收阵元的带通信号xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct)进行前置放大和抗混叠滤波等信号调理,式中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交信号;本例中对50路接收阵元的带通信号进行信号调理;
2)用单个采样电路和ADC器件对经过调理后的N路信号进行时分复用中频正交采样、模/数转换,然后进行正交解调,合成一路数字信号,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列;本例中用一个高速ADC器件取代50个低速ADC器件,采样电路的通道轮换选通开关和后续的ADC器件由采样脉冲控制协同工作,按照采样脉冲的频率对50路阵元信号依次轮流选通并进行模/数转换,节省了模/数转换的成本。
3)对步骤2)输出的序列进行周期性符号反转,使之补偿成为{Ii,Qi,Ii,Qi,Ii…}序列;
4)把从步骤3)输出的已经合成为一路数据流的各路信号输入级联结构的多通道FIR滤波器I,进行同相分量的滤波,一次性完成与信号匹配的低通滤波、内插延时、相位旋转和带通调制,滤波输出的同时完成样本抽取,每路信号的延时不仅要补偿波束到达各个接收阵元产生的相对延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外相对延时,滤波器I的每个通道滤波系数由内插系数 和相位旋转系数合成,滤波器I的系数vi T以及滤波器I的输出IDi(nTO)为:
上式中xi是第i路信号序列,τi是该路信号的延时,TO是通带信号的输出周期,IDi(t)是延时至波前对齐的基带复包络信号的同相分量,ai T和bi T分别是同相分量和正交分量的延时内插系数组,该系数是补零内插方法的低通滤波系数中对应某一内插点的系数组;
5)对从步骤4)输出的经过延时补偿至波前对齐的N路带通信号进行波束求和,序列输出通带上波束形成的结果。
所述的时分复用中频正交采样:在符合带通信号采样准则的前提下设置采样频率为接收带通信号中心频率的4N/(4k+1)倍或者4N/(4k+3)倍,其中k=0,1,2…,等效为每路信号的采样频率为接收带通信号中心频率的4/(4k+1)倍或者4/(4k+3)倍,其中k=0,1,2…,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列。在本例中控制采样频率fs=4N□fc=35MHz,则每路信号的等效采样频率是fs=4fc=700KHz,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}序列。
本例中根据接收信号的通带截止频率、过渡带宽及带外衰减等指标要求对每路信号的同相和正交分量分别进行60阶的低通滤波,同时根据波束形成的性能要求需要对采样信号进行10倍的补零内插。在补零内插方法中,进行10倍的补零内插在频谱上相当于对原始信号的频率进行了10倍压缩,如图3所示,因此需要设计600阶的低通滤波系数以供内插系数选取之用。由于本发明中采样后信号的同相和正交分量是合在一路数据流中的,因此相当于对一路信号进行120阶的滤波。考虑选用15个8阶子滤波模块进行120阶的滤波,滤波后每路信号的采样率为87.5KHz,符合带通信号的采样定理。显然,该滤波器输出也是各路信号交织在一起的一路数据流。
所述的级联结构的多通道FIR滤波器I:由信号重组模块和若干个子滤波模块组成,每个子滤波模块阶数为4M,分别对每路信号完成4M阶的滤波计算,其余结构同前例中所述的FIR滤波器I结构相同。本例中每个子滤波模块阶数4M等于8,即分别对每路信号完成8阶的滤波计算。
前文为了说明本发明对上述优选实施例进行披露。对于本领域中的技术人员来说变化和修改是显而易见的,例如采用若干个ADC器件对多路信号进行模/数转换等,所有这些变化和修改均意欲为后续的权利要求书所包含。
Claims (10)
1.一种宽带数字波束形成方法,其特征在于包括如下步骤:
1)对N路接收阵元的带通信号xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct)进行前置放大和抗混叠滤波等信号调理,式中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交信号;
2)用单个采样电路和ADC器件对经过调理后的N路信号进行时分复用中频正交采样、模/数转换,然后进行正交解调,合成一路数字信号,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列;
3)对步骤2)输出的序列进行周期性符号反转,使之补偿成为{Ii,Qi,Ii,Qi,Ii…}序列;
4)把从步骤3)输出的已经合成为一路数据流的各路信号同时输入级联结构的多通道FIR滤波器I和滤波器Q,分别进行同相和正交分量的滤波,一次性完成与信号匹配的低通滤波、内插延时和相位旋转,滤波输出的同时完成样本抽取,每路信号的延时不仅要补偿波束到达各个接收阵元产生的相对延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外相对延时,滤波器I和滤波器Q的每个通道滤波系数由内插系数 和相位旋转系数合成,滤波器I、Q的系数vi T、wi T以及滤波器I、Q的输出IDi(nTO)、QDi(nTO)分别为:
其中xi是第i路信号序列,τi是该路信号的延时,TO是复包络信号的输出周期,IDi(t)和QDi(t)分别是延时至波前对齐的基带复包络信号的同相和正交分量,ai T和bi T分别是同相分量和正交分量的延时内插系数组,该系数是补零内插方法的低通滤波系数中对应某一内插点的系数组;
5)对从步骤4)输出的经过延时补偿至波前对齐的N路基带复包络信号进行波束求和,序列输出基带上波束形成的结果。
2.根据权利要求1所述的一种宽带数字波束形成方法,其特征在于所述的时分复用中频正交采样:在符合带通信号采样准则的前提下设置采样频率为接收带通信号中心频率的4N/(4k+1)倍或者4N/(4k+3)倍,其中k=0,1,2…,等效为每路信号的采样频率为接收带通信号中心频率的4/(4k+1)倍或者4/(4k+3)倍,其中k=0,1,2…,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列。
3.根据权利要求1所述的一种宽带数字波束形成方法,其特征在于所述的级联结构的多通道FIR滤波器I和滤波器Q:由信号重组模块和若干个子滤波模块组成,每个子滤波模块阶数为2M,分别对每路信号完成2M阶的滤波计算,滤波器输出信号的采样频率等于原始信号采样频率的1/2M,按照原始信号的采样频率和降采样后所需的采样频率来确定子滤波器模块的阶数,子滤波模块个数由所需的滤波阶数和子滤波模块的阶数共同确定,滤波过程中采样信号从第一级子滤波模块依次流向最后一级子滤波模块,子滤波模块对每一个通道的滤波计算则从最后一级子滤波模块开始,并把每一个通道的滤波结果依次往前一级子滤波模块传递,在后一级子滤波模块输出某一通道信号的滤波结果到前一级子滤波模块之后,前一级子滤波模块开始该通道信号的滤波,最终在第一级子滤波模块处完成整个滤波计算,输出滤波结果。
4.根据权利要求3所述的一种宽带数字波束形成方法,其特征在于所述的信号重组模块:把按采样所得的信号排列顺序转换成2M个同一通道的信号连续排列的顺序并依次输出,重组后的采样信号输入到第一级子滤波模块数据RAM的输入端;
所述的子滤波模块:包括一个存储采样信号的双端口数据RAM、存储系数的双端口系数RAM、乘法器、加法器以及二路选择器,乘法器的两个输入端连接数据RAM和系数RAM的输出端,完成采样信号和对应滤波系数的乘积,并把乘积结果输出至加法器,乘法器的输出端和二路选择器的输出端连接到加法器的两个输入端,加法器的输出端同时连接到本级与前一级二路选择器的输入端,本级二路选择器的两个输入端分别连接本级和后一级加法器的输出端,最后一级二路选择器的两个输入分别是本级加法器的输出和常数零,加法器和二路选择器协同工作,在每路通道信号滤波的第一个周期二路选择器选通后一级加法器的输出同本级乘法器的输出相加,如此完成最后一级子滤波模块至本级子滤波模块滤波结果的累加,在剩余周期内二路选择器全部选通本级加法器的输出同乘法器输出进行相加。
所述的数据RAM:被划分为N个连续的数据块,每个数据块连续存储2M个同一通道的采样信号,所有子滤波模块数据RAM的初始读取地址都为第1个数据块的起始处地址,初始写入地址都为第N个数据块的起始处地址,读取地址和写入地址在数据RAM空间内循环递增,本级数据RAM的输出端连接到本级乘法器的输入端和后一级数据RAM的输入端,第一级数据RAM的输入端连接信号重组模块的输出端。
5.根据权利要求1所述的一种宽带数字波束形成方法,其特征在于所述的内插系数的求取方法:该系数是从补零内插方法的低通滤波系数中获取的,内插等价于对补零信号进行线性时不变滤波,通过内插将一个序列转换到一个较高的等效采样频率的内插序列,进行D倍的补零内插首先要在原始序列x(n)相邻两采样点之间等间隔***(D-1)个零值点,然后对补零后的序列xp(n)进行低通滤波,得到内插序列xD(n):
对应于某一个内插点xD(n)的插值系数只要选择计算该内插点时原始信号所对应的低通滤波系数h(n-rD)(-∞<r<+∞),即为对应该延时的一个内插系数组,若设计该低通滤波系数的幅频特性和复包络信号匹配,则在内插延时的同时完成了与信号匹配的低通滤波。
6.一种宽带数字波束形成方法,其特征在于包括如下步骤:
1)对N路接收阵元的带通信号xi(t)=Ii(t)cos(ωct)-Qi(t)sin(ωct)进行前置放大和抗混叠滤波等信号调理,式中Ii(t)和Qi(t)分别为带通信号xi(t)的基带同相和正交信号;
2)用单个采样电路和ADC器件对经过调理后的N路信号进行时分复用中频正交采样、模/数转换,然后进行正交解调,合成一路数字信号,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列;
3)对步骤2)输出的序列进行周期性符号反转,使之补偿成为{Ii,Qi,Ii,Qi,Ii…}序列;
4)把从步骤3)输出的已经合成为一路数据流的各路信号输入级联结构的多通道FIR滤波器I,进行同相分量的滤波,一次性完成与信号匹配的低通滤波、内插延时、相位旋转和带通调制,滤波输出的同时完成样本抽取,每路信号的延时不仅要补偿波束到达各个接收阵元产生的相对延时,还要补偿因对每个接收阵元在不同时刻采样所引入的额外相对延时,滤波器I的每个通道滤波系数由内插系数 和相位旋转系数合成,滤波器I的系数vi T以及滤波器I的输出IDi(nTO)为:
上式中xi是第i路信号序列,τi是该路信号的延时,TO是通带信号的输出周期,IDi(t)是延时至波前对齐的基带复包络信号的同相分量,ai T和bi T分别是同相分量和正交分量的延时内插系数组,该系数是补零内插方法的低通滤波系数中对应某一内插点的系数组;
5)对从步骤4)输出的经过延时补偿至波前对齐的N路带通信号进行波束求和,序列输出通带上波束形成的结果。
7.根据权利要求6所述的一种宽带数字波束形成方法,其特征在于所述的时分复用中频正交采样:在符合带通信号采样准则的前提下设置采样频率为接收带通信号中心频率的4N/(4k+1)倍或者4N/(4k+3)倍,其中k=0,1,2…,等效为每路信号的采样频率为接收带通信号中心频率的4/(4k+1)倍或者4/(4k+3)倍,其中k=0,1,2…,每路信号的采样输出为{Ii,-Qi,-Ii,Qi,Ii…}或{Ii,Qi,-Ii,-Qi,Ii…}序列。
8.根据权利要求6所述的一种宽带数字波束形成方法,其特征在于所述的级联结构的多通道FIR滤波器I:由信号重组模块和若干个子滤波模块组成,每个子滤波模块阶数为4M,分别对每路信号完成4M阶的滤波计算,滤波器输出信号的采样频率等于原始信号采样频率的1/4M,按照原始信号的采样频率和降采样后所需的采样频率来确定子滤波器模块的阶数,子滤波模块个数由所需的滤波阶数和子滤波模块的阶数共同确定,滤波过程中采样信号从第一级子滤波模块依次流向最后一级子滤波模块,子滤波模块对每一个通道的滤波计算则从最后一级子滤波模块开始,并把每一个通道的滤波结果依次往前一级子滤波模块传递,在后一级子滤波模块输出某一通道信号的滤波结果到前一级子滤波模块之后,前一级子滤波模块开始该通道信号的滤波,最终在第一级子滤波模块处完成整个滤波计算,输出滤波结果。
9.根据权利要求8所述的一种宽带数字波束形成方法,其特征在于所述的信号重组模块:把按采样所得的信号排列顺序转换成4M个同一通道的信号连续排列的顺序并依次输出,重组后的采样信号输入到第一级子滤波模块数据RAM的输入端。
所述的子滤波模块:包括一个存储采样信号的双端口数据RAM、存储系数的双端口系数RAM、乘法器、加法器以及二路选择器。乘法器的两个输入端连接数据RAM和系数RAM的输出端,完成采样信号和对应滤波系数的乘积,并把乘积结果输出至加法器,乘法器的输出端和二路选择器的输出端连接到加法器的两个输入端,加法器的输出端同时连接到本级与前一级二路选择器的输入端,本级二路选择器的两个输入端分别连接本级和后一级加法器的输出端,最后一级二路选择器的两个输入分别是本级加法器的输出和常数零,加法器和二路选择器协同工作,在每路通道信号滤波的第一个周期二路选择器选通后一级加法器的输出同本级乘法器的输出相加,如此完成最后一级子滤波模块至本级子滤波模块滤波结果的累加,在剩余周期内二路选择器全部选通本级加法器的输出同乘法器输出进行相加;
所述的数据RAM:被划分为N个连续的数据块,每个数据块连续存储2M个同一通道的采样信号,所有子滤波模块数据RAM的初始读取地址都为第1个数据块的起始处地址,初始写入地址都为第N个数据块的起始处地址,读取地址和写入地址在数据RAM空间内循环递增,本级数据RAM的输出端连接到本级乘法器的输入端和后一级数据RAM的输入端,第一级数据RAM的输入端连接信号重组模块的输出端。
10.根据权利要求6所述的一种宽带数字波束形成方法,其特征在于所述的内插系数的求取方法:该系数是从补零内插方法的低通滤波系数中获取的。内插等价于对补零信号进行线性时不变滤波,通过内插将一个序列转换到一个较高的等效采样频率的内插序列,进行D倍的补零内插首先要在原始序列x(n)相邻两采样点之间等间隔***(D-1)个零值点,然后对补零后的序列xp(n)进行低通滤波,得到内插序列xD(n):
对应于某一个特定内插点xD(n)的插值系数只要选择计算该内插点时原始信号所对应的低通滤波系数h(n-rD)(-∞<r<+∞),即为对应该延时的一个内插系数组,若设计该低通滤波系数的幅频特性和复包络信号匹配,则在内插延时的同时完成了与信号匹配的低通滤波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100624796A CN101296014B (zh) | 2008-06-12 | 2008-06-12 | 一种宽带数字波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100624796A CN101296014B (zh) | 2008-06-12 | 2008-06-12 | 一种宽带数字波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101296014A true CN101296014A (zh) | 2008-10-29 |
CN101296014B CN101296014B (zh) | 2011-12-07 |
Family
ID=40066066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100624796A Expired - Fee Related CN101296014B (zh) | 2008-06-12 | 2008-06-12 | 一种宽带数字波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101296014B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129069A (zh) * | 2010-12-17 | 2011-07-20 | 浙江大学 | 基于多fpga的波束形成器输出数据带宽控制装置 |
CN102156281A (zh) * | 2010-12-17 | 2011-08-17 | 浙江大学 | 基于单fpga的波束形成器输出数据带宽控制方法及装置 |
CN102313574A (zh) * | 2011-07-05 | 2012-01-11 | 刘继国 | 基于延迟线的信号测量通道合并方法及*** |
CN101644778B (zh) * | 2008-12-30 | 2012-03-14 | 中国科学院声学研究所 | 手持式成像声纳及其成像方法 |
CN102680946A (zh) * | 2012-06-06 | 2012-09-19 | 西安电子科技大学 | 基于微分算子的超宽带波束形成器及实现方法 |
CN102768358A (zh) * | 2011-05-05 | 2012-11-07 | 中国科学院声学研究所 | 一种基于fpga的水下实时成像方法及*** |
CN104199037A (zh) * | 2014-08-07 | 2014-12-10 | 哈尔滨工程大学 | 一种参量阵声纳瞬态宽带激励信号的递推滤波调制方法 |
CN105468882A (zh) * | 2014-07-28 | 2016-04-06 | 航天恒星科技有限公司 | 卫星自主定轨方法及*** |
CN106019289A (zh) * | 2016-05-19 | 2016-10-12 | 杭州电子科技大学 | 宽带相控声学多普勒测速***信号处理方法 |
CN106411379A (zh) * | 2016-09-29 | 2017-02-15 | 电子科技大学 | 一种降低硬件资源消耗的宽带波束形成设计方法 |
CN107797096A (zh) * | 2017-10-20 | 2018-03-13 | 电子科技大学 | 一种基于麦克风面阵的鸣笛检测定位方法 |
CN109361473A (zh) * | 2018-12-06 | 2019-02-19 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 高速大容量光子传输网络 |
CN110333488A (zh) * | 2019-07-23 | 2019-10-15 | 上海雪狸传感技术有限公司 | 一种毫米波宽角波束扫描雷达传感器 |
CN111917428A (zh) * | 2020-07-21 | 2020-11-10 | 湖南拉盖尔电子科技有限公司 | 一种数字中频信号快速基带化和降采样方法 |
CN116505994A (zh) * | 2023-06-26 | 2023-07-28 | 成都金支点科技有限公司 | 一种多波束形成方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101098179B (zh) * | 2006-06-30 | 2010-06-30 | 中国科学院声学研究所 | 一种宽带频域数字波束形成方法 |
-
2008
- 2008-06-12 CN CN2008100624796A patent/CN101296014B/zh not_active Expired - Fee Related
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644778B (zh) * | 2008-12-30 | 2012-03-14 | 中国科学院声学研究所 | 手持式成像声纳及其成像方法 |
CN102156281A (zh) * | 2010-12-17 | 2011-08-17 | 浙江大学 | 基于单fpga的波束形成器输出数据带宽控制方法及装置 |
CN102129069B (zh) * | 2010-12-17 | 2012-05-09 | 浙江大学 | 基于多fpga的波束形成器输出数据带宽控制装置 |
CN102129069A (zh) * | 2010-12-17 | 2011-07-20 | 浙江大学 | 基于多fpga的波束形成器输出数据带宽控制装置 |
CN102768358A (zh) * | 2011-05-05 | 2012-11-07 | 中国科学院声学研究所 | 一种基于fpga的水下实时成像方法及*** |
CN102768358B (zh) * | 2011-05-05 | 2014-03-26 | 中国科学院声学研究所 | 一种基于fpga的水下实时成像方法及*** |
CN102313574A (zh) * | 2011-07-05 | 2012-01-11 | 刘继国 | 基于延迟线的信号测量通道合并方法及*** |
CN102313574B (zh) * | 2011-07-05 | 2013-08-14 | 刘继国 | 基于延迟线的信号测量通道合并方法及*** |
CN102680946A (zh) * | 2012-06-06 | 2012-09-19 | 西安电子科技大学 | 基于微分算子的超宽带波束形成器及实现方法 |
CN102680946B (zh) * | 2012-06-06 | 2014-06-18 | 西安电子科技大学 | 基于微分算子的超宽带波束形成器及实现方法 |
CN105468882A (zh) * | 2014-07-28 | 2016-04-06 | 航天恒星科技有限公司 | 卫星自主定轨方法及*** |
CN104199037A (zh) * | 2014-08-07 | 2014-12-10 | 哈尔滨工程大学 | 一种参量阵声纳瞬态宽带激励信号的递推滤波调制方法 |
CN106019289A (zh) * | 2016-05-19 | 2016-10-12 | 杭州电子科技大学 | 宽带相控声学多普勒测速***信号处理方法 |
CN106019289B (zh) * | 2016-05-19 | 2018-06-26 | 杭州电子科技大学 | 宽带相控声学多普勒测速***信号处理方法 |
CN106411379A (zh) * | 2016-09-29 | 2017-02-15 | 电子科技大学 | 一种降低硬件资源消耗的宽带波束形成设计方法 |
CN106411379B (zh) * | 2016-09-29 | 2019-09-27 | 电子科技大学 | 一种降低硬件资源消耗的宽带波束形成设计方法 |
CN107797096A (zh) * | 2017-10-20 | 2018-03-13 | 电子科技大学 | 一种基于麦克风面阵的鸣笛检测定位方法 |
CN109361473A (zh) * | 2018-12-06 | 2019-02-19 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 高速大容量光子传输网络 |
CN110333488A (zh) * | 2019-07-23 | 2019-10-15 | 上海雪狸传感技术有限公司 | 一种毫米波宽角波束扫描雷达传感器 |
CN111917428A (zh) * | 2020-07-21 | 2020-11-10 | 湖南拉盖尔电子科技有限公司 | 一种数字中频信号快速基带化和降采样方法 |
CN116505994A (zh) * | 2023-06-26 | 2023-07-28 | 成都金支点科技有限公司 | 一种多波束形成方法和装置 |
CN116505994B (zh) * | 2023-06-26 | 2023-09-01 | 成都金支点科技有限公司 | 一种多波束形成方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101296014B (zh) | 2011-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101296014B (zh) | 一种宽带数字波束形成方法 | |
CN101304404B (zh) | 一种宽带时域数字波束形成方法 | |
CN102798840B (zh) | 外辐射源雷达宽带信道化接收***及fpga实现方法 | |
US4777612A (en) | Digital signal processing apparatus having a digital filter | |
US4866647A (en) | Continuously variable digital delay circuit | |
CN101098179B (zh) | 一种宽带频域数字波束形成方法 | |
US6014682A (en) | Methods and apparatus for variable-rate down-sampling filters for discrete-time sampled systems using a fixed sampling rate | |
US5592517A (en) | Cascaded comb integrator interpolating filters | |
CN1702986B (zh) | 基带多径衰落信道模拟器 | |
EP2904419B1 (en) | Improvements in and relating to radar receivers | |
GB2419482A (en) | Analog to digital signal conversion circuit and method | |
CN108051783A (zh) | 宽带相控阵雷达回波模拟器发射多通道延时补偿方法及*** | |
CN102680946B (zh) | 基于微分算子的超宽带波束形成器及实现方法 | |
US5831879A (en) | Digital transmit filter | |
CN101881821B (zh) | 一种分数阶傅里叶域信道化接收方法 | |
US6272441B1 (en) | Method for determining the pulse response of a broad band linear system and a measuring circuit for carrying out the method | |
CN107144821B (zh) | 宽带数字阵雷达中基于时延波束形成的高效接收通道 | |
US4779054A (en) | Digital inphase/quadrature product detector | |
CN103188186A (zh) | 重采样处理装置和方法、以及数字调制信号产生装置 | |
CN112526458A (zh) | 一种基于提取参数分数时延的宽带nlfm发射波束形成方法 | |
US4225832A (en) | Self-adapting equalizer | |
CN104242871B (zh) | 一种基于嵌套阵的多倍频程恒定束宽波束形成方法及*** | |
US5490511A (en) | Digital phase shifting apparatus | |
CN100463457C (zh) | 一种多个通道复接的级联积分梳状滤波器 | |
CN117579108A (zh) | 一种低复杂度宽带接收数字波束形成器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111207 Termination date: 20120612 |