CN117991243A - 一种基于fpga的声纳信号匹配滤波处理方法及其系统 - Google Patents
一种基于fpga的声纳信号匹配滤波处理方法及其系统 Download PDFInfo
- Publication number
- CN117991243A CN117991243A CN202410397142.XA CN202410397142A CN117991243A CN 117991243 A CN117991243 A CN 117991243A CN 202410397142 A CN202410397142 A CN 202410397142A CN 117991243 A CN117991243 A CN 117991243A
- Authority
- CN
- China
- Prior art keywords
- data
- matched filter
- output
- signal
- control 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
- 238000001914 filtration Methods 0.000 title claims abstract description 45
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 239000000872 buffer Substances 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000005070 sampling Methods 0.000 claims description 40
- 238000009825 accumulation Methods 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000000630 rising effect Effects 0.000 claims description 3
- 238000003491 array Methods 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F15/00—Digital computers in general; Data processing equipment in general
- G06F15/76—Architectures of general purpose stored program computers
- G06F15/78—Architectures of general purpose stored program computers comprising a single central processing unit
- G06F15/7867—Architectures of general purpose stored program computers comprising a single central processing unit with reconfigurable architecture
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0248—Filters characterised by a particular frequency response or filtering method
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于FPGA的声纳信号匹配滤波处理方法及其系统,涉及声纳信号处理技术领域,所述系统包括时序控制模块、数据缓存移位寄存器、滤波器系数存储RAM、乘累加器与求模模块,系统在每个采样输入至下一采样输入前,实时完成匹配滤波输出;所述复数乘法累加器包括一个DSP乘法器与若干加法器、寄存器;FPGA上实现一种时域乘累加的实时匹配滤波处理方法;本发明方法能够大幅降低对乘法器硬件资源的需求,且能够保持匹配滤波输出幅度的稳定。
Description
技术领域
本发明涉及声纳信号处理技术领域,尤其涉及一种基于FPGA的声纳信号匹配滤波处理方法及其系统。
背景技术
在雷达、声纳等系统中,往往利用发射宽带信号,并借助专门的脉冲压缩处理算法来实现所需的距离维高分辨力。匹配滤波常见的脉冲压缩处理算法,在各种宽带系统中得到普遍的使用。它本质上属于FIR数字滤波器的一种,但特殊之处在于其滤波器系数刚好为发射信号的共轭,这样对于一个刚好为发射信号拷贝的回波信号,经过匹配滤波器之后的输出恰好为一个窄脉冲,其脉冲宽度由发射信号的带宽决定,如此即实现了距离上的高分辨力(详见图1)。
由于匹配滤波处理是一种FIR数字滤波器,因此在一个异构嵌入式系统中,非常适合于在FPGA内实现,以减少对DSP或ARM处理器的算力占用。
目前FPGA内实现匹配滤波算法分为频域乘法实现和时域乘累加实现。频域乘法实现首先对输入信号进行FFT操作,然后将其与系统冲击响应进行相乘,最后对乘积进行IFFT操作得到输出。由于调用了FPGA的FFT核,计算速度快。但是由于需要缓存一段数据进行FFT操作,因此输出结果不是实时的,而且分段缓存之间可能存在拼接不连续问题。
相对而言,时域乘累加实现基于FIR架构(详见图2),即
其中为当前时刻缓存的前N个采样点的输入信号,为匹配滤波器系数,/>为匹配滤波输出。
可以直接调用FPGA的FIR滤波器的IP核来实现,更多的做法是将输入数据和滤波器系数分别存储在一片RAM内,在每个采样时刻,用N个乘法器分别实现每个输入采样信号与系数的相乘,最后通过一个加法器得到输出。这种实现方法结构简单,且能够每个采样点实时输出滤波结果。
然后,常规时域乘累加实现匹配滤波处理的架构存在的主要问题有以下几点:
1、对FPGA乘法器资源要求较大,假设匹配滤波器阶数为N,则每个采样点都需要做N点乘累加运算。当N为1024甚至更大时,显然常规时域乘累加做实现要大量的乘法器来实现,这可能突破了FPGA资源限制,或者说导致其它模块无法实现。在多通道数据处理系统如多波束声纳系统中,往往存在256甚至更多的并行通道需要同时完成匹配滤波,则需要的资源更多。
2、采用RAM做数据缓存,当一个新的采样输入时,要将RAM内所有数据全部更新一遍,存在大量数据搬移工作。
3、在发射脉宽变化范围比较大的系统中,为满足最大发射脉宽的需求,匹配滤波器的阶数会比较大,这样在发射短脉宽时,匹配滤波器系数绝大部分都为0,这会导致匹配滤波器输出值大幅降低,影响输出数据的动态范围。
由此可以看出,急需一种匹配滤波处理的FPGA实现方法来解决上述问题。
发明内容
本发明提供一种基于FPGA的声纳信号匹配滤波处理方法及其系统。
根据本发明的一个较佳实施例,本发明提供一种基于FPGA的声纳信号匹配滤波处理方法,包括下列步骤:
(a)移入采样信号数据Xi至数据缓存移位寄存器,其中i=1,2,……,N;
(b)数据缓存移位寄存器根据时序控制模块输出的数据选择信号D_SELECT,在不同的CLK_CTRL时钟沿依次输出采样信号数据Xi到乘累加器;
(c)滤波器系数存储RAM根据时序控制模块输出的读地址ADD_r,在不同的CLK_CTRL时钟沿依次输出匹配滤波器系数Ci到乘累加器,其中i=1,2,……,N;
(d)乘累加器根据时序控制模块输出的节拍控制信号,完成当前缓存器内采样信号数据Xi和匹配滤波器系数Ci的复数乘累加,得到累加输出结果;
(e)求模模块利用乘累加器的累加输出结果和匹配滤波器系数矩阵中非0值的个数K,计算得到匹配滤波器最终的匹配滤波输出。
作为优选的技术方案,在所述基于FPGA的声纳信号匹配滤波处理方法中,进一步包括步骤(f)当系统调整发射波形参数时,通过SPI接口将新的匹配滤波器系数Ci下发至系数更新模块,系数更新模块将写使能信号EN_w设置为高电平,并在N个控制时钟周期内,通过地址总线ADD_w和数据总线DATA_w,将匹配滤波器系数Ci依次写入系数更新模块,同时将匹配滤波器系数矩阵中非0值的个数K传递给求模模块。
作为优选的技术方案,在步骤(a)中,所述数据缓存移位寄存器中输入数据相匹配的采样率由采样时钟决定,每个采样时钟输入一个数据,每个数据位宽为2D,每个数据中高D位为实部,低D位为虚部。
作为优选的技术方案,在步骤(b)中,所述数据缓存移位寄存器输出数据相匹配的采样率与输入数据相同,每个采样时钟输出一个结果,每个结果位宽为D。
作为优选的技术方案,在上述方法中,每个CLK_DATA上升沿时刻,数据向右移动一个存储空间,并移入一个新的数据,数据移入移出为时钟驱动自动进行。
作为优选的技术方案,在步骤(c)中,每个匹配滤波器系数Ci为2D比特,其中高D位为实部,低D位为虚部,匹配滤波器系数Ci经一次加载后,无需随每个采样点进行移位。
作为优选的技术方案,在步骤(d)中,用一个DSP乘法器完成匹配滤波所需的所有乘法计算,时序控制模块根据内部计数值,在不同时刻输出不同选择信号,改变乘法器的两个输入,得到相应乘法输出,并对输出分组进行累加,具体做法步骤如下:
(1)利用D_SELECT信号和ADD_r,依次选择Xi和Ci进行相乘;
(2)每个数据可表示为,每个滤波器系数可表示为,两者的相乘可表示为:
,
利用时序控制模块输出的X_SELECT信号和C_SELECT信号,分别选择和的实部与虚
部,并经过四次乘法,分别得到;
(3)利用时序控制模块输出的Z_DISPAC信号,将Z_DISPAC信号分别输出至相应的
累加模块,由四个累加模块最终得到四个累加值,其中,
;
(4)在下一个采样数据输入之前,完成当前缓存器内采样信号数据Xi和匹配滤波器系数Ci的复数乘累加,得到累加输出结果:
;
。
作为优选的技术方案,在步骤(e)中,求模模块利用乘累加模块的四个累加值,以及匹配滤波器系数矩阵中非0值的个数K,得到匹配滤波器最终
的匹配滤波输出,步骤如下:
(1)分别得到实部与虚部累加结果:
(2)根据K值,对与分别左移位,得到与;K表示N个滤波器
系数中非0值个数;
(3)分别对与进行平方操作,得到;
(4)对求和并开方,得到匹配滤波器最终匹配滤波输出Y;
。
根据本发明的另一个较佳实施例,本发明提供一种基于FPGA的声纳信号匹配滤波处理系统,其包括:
时序控制模块;
数据缓存移位寄存器,其中所述数据缓存移位寄存器和所述时序控制模块可通信地连接,所述数据缓存移位寄存器可接收采样信号数据Xi,i=1,2,……,N;
滤波器系数存储RAM,其中所述滤波器系数存储RAM和所述时序控制模块可通信地连接;
乘累加器,其中所述乘累加器可通信地连接于数据缓存移位寄存器、所述滤波器系数存储RAM以及所述时序控制模块,所述数据缓存移位寄存器根据时序控制模块输出的数据选择信号D_SELECT,在不同的CLK_CTRL时钟沿依次输出采样信号数据Xi到所述乘累加器,滤波器系数存储RAM根据时序控制模块输出的读地址ADD_r,在不同的CLK_CTRL时钟沿依次输出匹配滤波器系数Ci到乘累加器,其中i=1,2,……,N;乘累加器根据时序控制模块输出的节拍控制信号,完成当前缓存器内采样信号数据Xi和匹配滤波器系数Ci的复数乘累加,得到累加输出结果;
求模模块,其中所述求模模块可通信地连接于所述乘累加器,所述求模模块利用所述乘累加器的累加输出结果和匹配滤波器系数矩阵中非0值的个数K,得到匹配滤波器最终的匹配滤波输出。
作为优选的技术方案,所述基于FPGA的声纳信号匹配滤波处理系统进一步包括系数更新模块和可通信地连接于所述系数更新模块的SPI接口,所述SPI接口包含时钟信号线、数据信号线与片选信号线,所述SPI接口用于修改匹配滤波器系数矩阵,所述系数更新模块可通信地连接于所述求模模块。
作为优选的技术方案,所述乘累加器包括一个DSP乘法器和连接于所述DSP乘法器的多个累加模块,其中所述DSP乘法器和所述累加模块可通信地连接于所述时序控制模块。
作为优选的技术方案,数据缓存移位寄存器为N阶移位寄存器数组,每个寄存器为2D比特,其中高D位为输入信号实部,低D位为输入信号虚部。
作为优选的技术方案,时序控制模块的控制时钟为内部逻辑运行基准时钟,时序控制模块由一个主计数器驱动,并在计数器的不同计数值上,控制其它各模块进行相应的操作,计数器的驱动时钟为主控时钟CLK_CTRL。
本发明的优点:
在FPGA上实现了一种声纳信号匹配滤波处理方法。与现有方案需要消耗大量DSP乘法器资源,且无法保证发射脉宽变化时匹配滤波输出幅度稳定性相比,本发明通过合理设计时序,分时复用DSP乘法器资源,仅使用一个DSP乘法器即可完成整个匹配滤波处理,且能实现发射脉宽变化时匹配滤波输出幅度稳定。本发明还通过移位寄存器数组来自动实现采样数据的移入,避免采用BRAM的方式需要进行采样数据手动移入。
本发明利用时序控制模块实现单个DSP乘法器完成匹配滤波处理。
本发明通过对乘累加实部虚部自适应移位实现匹配滤波输出幅度稳定。
本发明利用寄存器数组实现采样数据的自动移入。
附图说明
图1为现有技术中匹配滤波实现距离高分辨图;其中,上图为两个目标回波;下图为匹配滤波器输出;
图2为现有技术中常规时域乘累加实现匹配滤波处理架构图;
图3为本发明实施例1中系统 FPGA内匹配滤波模块流程图;
图4为本发明实施例1中数据缓存移位寄存器流程图;
图5为本发明实施例1中滤波器系数存储模块流程图;
图6为本发明实施例1中系数更新模块流程图;
图7为本发明实施例1中乘累加器流程图;
图8为本发明实施例1中求模模块流程图;
图9为本发明实施例1工作流程图。
具体实施方式
本发明提供了一种基于FPGA的声纳信号匹配滤波处理方法及其系统。
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施例,进一步阐述本发明。
实施例1:
开启系统后处于待机状态,系统检测到一个新的采样信号数据X被移入,从而触发匹配滤波的处理流程,输出数据;
所述系统包括时序控制模块、数据缓存移位寄存器、滤波器系数存储RAM、乘累加器与求模模块,系统在每个采样输入至下一采样输入前,实时完成匹配滤波输出;所述乘累加器包括一个DSP乘法器与若干加法器、寄存器;
所述数据缓存移位寄存器中输入数据相匹配的采样率由采样时钟决定,每个采样时钟输入一个数据;每个数据位宽为2D,每个数据中高D位为实部,低D位为虚部。2*D是位宽,D为非零正数,具体由数据类型决定。当D=16时,此时的数据就是32位二进制数,其中,高D位指的是第32-17位对应的部分,低D位指的是第16-1位对应的部分。
输出数据相匹配的采样率与输入数据相同,每个采样时钟输出一个结果;每个结果位宽为D;
时序控制模块的控制时钟为内部逻辑运行基准时钟;
还包括系数更新模块、SPI接口,所述SPI接口包含时钟信号线、数据信号线与片选信号线,所述SPI接口用于修改匹配滤波器系数矩阵;所述匹配滤波器系数矩阵包含N个匹配滤波器系数,其中前K个系数非0;每个系数位宽为2D,每个系数位宽中高D位为实部,低D位为虚部。
数据缓存移位寄存器为N阶移位寄存器数组,每个寄存器为2D比特,其中高D位为输入信号实部,低D位为输入信号虚部;
每个CLK_DATA上升沿时刻,数据向右移动一个存储空间,并移入一个新的数据;采用移位寄存器替代了RAM,数据移入移出为时钟驱动自动进行,无需手动搬运;
而在每个新的采样输入完成后,至下一个采样输入之前,根据时序控制模块输出的数据选择信号D_SELECT,在不同的CLK_CTRL时钟沿依次输出数据Xi(i=1,2,……,N)到乘累加器。
滤波器系数存储RAM存储匹配滤波器的N个匹配滤波器系数Ci(i=1,2,……,N);每个系数为2D比特,其中高D位为实部,低D位为虚部;由于匹配滤波器系数经一次加载后,无需随每个采样点进行移位,因此可以用FPGA内部的DP-RAM实现;
在每个新的采样输入完成后,至下一个采样输入之前,根据时序控制模块输出的读地址ADD_r,在不同的CLK_CTRL时钟沿依次输出数据Ci(i=1,2,……,N)到乘累加器。
所述系数更新模块通过SPI接口接受更高层的主控控制,当系统调整发射波形参数时,通过SPI接口将新的匹配滤波器系数下发至系数更新模块,系数更新模块将写使能信号EN_w设置为高电平,并在N个控制时钟周期内,通过地址总线ADD_w和数据总线DATA_w,将系数Ci(i=1,2,……,N)依次写入系数更新模块;同时将匹配滤波器系数矩阵中非0值的个数K传递给求模模块。
乘累加器根据时序控制模块输出的节拍控制信号,在下一个采样数据输入之前,
完成当前缓存器内数据Xi和匹配滤波器系数Ci的复数乘累加,得到4个累加值,累加输出结果如下:
。
由于需要用一个DSP乘法器完成匹配滤波所需的所有乘法计算,因此需要时序控制模块根据内部计数值,在不同时刻输出不同选择信号,改变乘法器的两个输入,得到相应乘法输出,并对输出分组进行累加,具体做法步骤如下:
1、利用D_SELECT信号和ADD_r,依次选择Xi和Ci(i=1,2,……,N)进行相乘;
2、每个数据可表示为,每个滤波器系数可表示为。
两者的相乘可表示为:
利用时序控制模块输出的X_SELECT信号和C_SELECT信号,分别选择和的实部与虚
部,并经过四次乘法,分别得到。
3、利用时序控制模块输出的Z_DISPAC信号,将Z_DISPAC信号分别输出至相应的累
加模块。乘累加器的每个累加模块由一个加法器和一个寄存器组成。四个累加模块最终得
到四个累加值。
,
求模模块利用乘累加器输出的四个累加值,得到匹配滤波
器最终匹配滤波输出Y,步骤如下:
1、分别得到实部与虚部累加结果:
2、根据K值,对与分别左移位,得到与;K表示N个滤波器
系数中非0值个数;
3、分别对与进行平方操作,得到;
4、对求和并开方,得到匹配滤波器最终匹配滤波输出Y:
。
平方操作可以继续复用乘累加器的DSP乘法器,此时乘累加操作已经结束,DSP乘法器处于空闲状态。开方运算调用CORDIC IP核进行。
时序控制模块为整个匹配滤波处理模块的主控部分,由一个主计数器驱动,并在计数器的不同计数值上,控制其它各模块进行相应的操作。计数器的驱动时钟为主控时钟CLK_CTRL。
本系统的工作流程如下
整个系统初始处于待机状态,当系统检测到一个新的采样信号X被移入时,马上触发匹配滤波的处理流程,整个工作流程的时序图如图5所示,其中方框中的数字代表计数器计数,即主控时钟CLK_CTRL周期数。
在计数值为3时完成初始化;
在计数值为4时开始第一个乘法计算;
在计数值为21时完成第一次乘累加计算;
在计数值为16i+5时完成第i次乘累加计算;
在计数值为16N+5时完成所有N次乘累加计算;
在计数值为16N+9时完成实部平方计算;
在计数值为16N+13时完成虚部平方计算;
在计数值为16N+14时完成实部虚部平方求和计算;
在计数值为16N+18时完成CORDIC开方计算;
在计数值为16N+19时准备数据输出(见图9)。
以上显示和描述了本发明的基本原理、主要特征及本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (13)
1.一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,包括下列步骤:
(a)移入采样信号数据Xi至数据缓存移位寄存器,其中i=1,2,……,N;
(b)数据缓存移位寄存器根据时序控制模块输出的数据选择信号D_SELECT,在不同的CLK_CTRL时钟沿依次输出采样信号数据Xi到乘累加器;
(c)滤波器系数存储RAM根据时序控制模块输出的读地址ADD_r,在不同的CLK_CTRL时钟沿依次输出匹配滤波器系数Ci到乘累加器,其中i=1,2,……,N;
(d)乘累加器根据时序控制模块输出的节拍控制信号,完成当前缓存器内采样信号数据Xi和匹配滤波器系数Ci的复数乘累加,得到累加输出结果;
(e)求模模块利用乘累加器的累加输出结果和匹配滤波器系数矩阵中非0值的个数K,计算得到匹配滤波器最终的匹配滤波输出。
2.根据权利要求1所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在所述基于FPGA的声纳信号匹配滤波处理方法中,进一步包括步骤(f)当系统调整发射波形参数时,通过SPI接口将新的匹配滤波器系数Ci下发至系数更新模块,系数更新模块将写使能信号EN_w设置为高电平,并在N个控制时钟周期内,通过地址总线ADD_w和数据总线DATA_w,将匹配滤波器系数Ci依次写入系数更新模块,同时将匹配滤波器系数矩阵中非0值的个数K传递给求模模块。
3.根据权利要求1所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在步骤(a)中,所述数据缓存移位寄存器中输入数据相匹配的采样率由采样时钟决定,每个采样时钟输入一个数据,每个数据位宽为2D,每个数据中高D位为实部,低D位为虚部。
4.根据权利要求3所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在步骤(b)中,所述数据缓存移位寄存器输出数据相匹配的采样率与输入数据相同,每个采样时钟输出一个结果,每个结果位宽为D。
5.根据权利要求4所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在上述方法中,每个CLK_DATA上升沿时刻,数据向右移动一个存储空间,并移入一个新的数据,数据移入移出为时钟驱动自动进行,无需手动搬运。
6.根据权利要求5所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在步骤(c)中,每个匹配滤波器系数Ci为2D比特,其中高D位为实部,低D位为虚部,匹配滤波器系数Ci经一次加载后,无需随每个采样点进行移位。
7.根据权利要求6所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在步骤(d)中,用一个DSP乘法器完成匹配滤波所需的所有乘法计算,时序控制模块根据内部计数值,在不同时刻输出不同选择信号,改变乘法器的两个输入,得到相应乘法输出,并对输出分组进行累加,具体做法步骤如下:
(1)利用D_SELECT信号和ADD_r,依次选择Xi和Ci进行相乘;
(2)每个数据可表示为,每个滤波器系数可表示为/>,两者的相乘可表示为:
,
,
利用时序控制模块输出的X_SELECT信号和C_SELECT信号,分别选择和的实部与虚部,并经过四次乘法,分别得到;
(3)利用时序控制模块输出的Z_DISPAC信号,将Z_DISPAC信号分别输出至相应的累加模块,由四个累加模块最终得到四个累加值,其中,
;
;
(4)在下一个采样数据输入之前,完成当前缓存器内采样信号数据Xi和匹配滤波器系数Ci的复数乘累加,得到累加输出结果:
;
。
8.根据权利要求7所述的一种基于FPGA的声纳信号匹配滤波处理方法,其特征在于,在步骤(e)中,求模模块利用乘累加模块的四个累加值,以及匹配滤波器系数矩阵中非0值的个数K,得到匹配滤波器最终的匹配滤波输出,步骤如下:
(1)分别得到实部与虚部累加结果:
;
(2)根据K值,对与/>分别左移/>位,得到/>与/>;K表示N个滤波器系数中非0值个数;
(3)分别对与/>进行平方操作,得到/>;
(4)对求和并开方,得到匹配滤波器最终匹配滤波输出Y;
。
9.一种基于FPGA的声纳信号匹配滤波处理系统,用于执行如权利要求1至9任一所述的基于FPGA的声纳信号匹配滤波处理方法,其特征在于,包括:
时序控制模块;
数据缓存移位寄存器,其中所述数据缓存移位寄存器和所述时序控制模块可通信地连接,所述数据缓存移位寄存器可接收采样信号数据Xi,i=1,2,……,N;
滤波器系数存储RAM,其中所述滤波器系数存储RAM和所述时序控制模块可通信地连接;
乘累加器,其中所述乘累加器可通信地连接于数据缓存移位寄存器、所述滤波器系数存储RAM以及所述时序控制模块,所述数据缓存移位寄存器根据时序控制模块输出的数据选择信号D_SELECT,在不同的CLK_CTRL时钟沿依次输出采样信号数据Xi到所述乘累加器,滤波器系数存储RAM根据时序控制模块输出的读地址ADD_r,在不同的CLK_CTRL时钟沿依次输出匹配滤波器系数Ci到乘累加器,其中i=1,2,……,N;乘累加器根据时序控制模块输出的节拍控制信号,完成当前缓存器内采样信号数据Xi和匹配滤波器系数Ci的复数乘累加,得到累加输出结果;
求模模块,其中所述求模模块可通信地连接于所述乘累加器,所述求模模块利用所述乘累加器的累加输出结果和匹配滤波器系数矩阵中非0值的个数K,得到匹配滤波器最终的匹配滤波输出。
10.根据权利要求9所述的一种基于FPGA的声纳信号匹配滤波处理系统,其特征在于,所述基于FPGA的声纳信号匹配滤波处理系统进一步包括系数更新模块和可通信地连接于所述系数更新模块的SPI接口,所述SPI接口包含时钟信号线、数据信号线与片选信号线,所述SPI接口用于修改匹配滤波器系数矩阵,所述系数更新模块可通信地连接于所述求模模块。
11.根据权利要求9所述的一种基于FPGA的声纳信号匹配滤波处理系统,其特征在于,所述乘累加器包括一个DSP乘法器和连接于所述DSP乘法器的多个累加模块,其中所述DSP乘法器和所述累加模块可通信地连接于所述时序控制模块。
12.根据权利要求9所述的一种基于FPGA的声纳信号匹配滤波处理系统,其特征在于,数据缓存移位寄存器为N阶移位寄存器数组,每个寄存器为2D比特,其中高D位为输入信号实部,低D位为输入信号虚部。
13.根据权利要求9所述的一种基于FPGA的声纳信号匹配滤波处理系统,其特征在于,时序控制模块的控制时钟为内部逻辑运行基准时钟,时序控制模块由一个主计数器驱动,并在计数器的不同计数值上,控制其它各模块进行相应的操作,计数器的驱动时钟为主控时钟CLK_CTRL。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410397142.XA CN117991243B (zh) | 2024-04-03 | 2024-04-03 | 一种基于fpga的声纳信号匹配滤波处理方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410397142.XA CN117991243B (zh) | 2024-04-03 | 2024-04-03 | 一种基于fpga的声纳信号匹配滤波处理方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117991243A true CN117991243A (zh) | 2024-05-07 |
CN117991243B CN117991243B (zh) | 2024-07-02 |
Family
ID=90895478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410397142.XA Active CN117991243B (zh) | 2024-04-03 | 2024-04-03 | 一种基于fpga的声纳信号匹配滤波处理方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117991243B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1323105A (zh) * | 2001-03-19 | 2001-11-21 | 深圳市中兴集成电路设计有限责任公司 | 一种相关器装置 |
CN1859038A (zh) * | 2006-01-11 | 2006-11-08 | 华为技术有限公司 | 一种时隙同步实现装置和方法 |
CN101262240A (zh) * | 2008-04-25 | 2008-09-10 | 浙江大学 | 一种易于硬件实现的全数字频率变换的方法及其装置 |
CN101278491A (zh) * | 2005-09-28 | 2008-10-01 | 日本电气株式会社 | 调制器、滤波器、滤波器增益控制方法和代码调制方法 |
CN101282322A (zh) * | 2008-03-05 | 2008-10-08 | 中科院嘉兴中心微系统所分中心 | 一种应用于无线中程传感网物理层的内插数字滤波器装置 |
JP2009153096A (ja) * | 2007-11-28 | 2009-07-09 | Sony Corp | 受信装置、受信方法、およびプログラム |
CN102739195A (zh) * | 2012-06-06 | 2012-10-17 | 华为技术有限公司 | 一种fir滤波器的处理方法、装置和系统 |
CN105741313A (zh) * | 2014-12-09 | 2016-07-06 | 南京理工大学 | 一种基于fpga的快速实时动目标检测系统 |
CN108900177A (zh) * | 2018-06-14 | 2018-11-27 | 福建星海通信科技有限公司 | 一种fir滤波器及其对数据进行滤波的方法 |
CN110058201A (zh) * | 2019-04-22 | 2019-07-26 | 西安电子工程研究所 | 一种基于fpga资源复用实现多波形多速率时域脉压的方法 |
CN113317801A (zh) * | 2021-04-23 | 2021-08-31 | 清华大学 | 单通道的脑电图信号处理装置及方法 |
CN114142829A (zh) * | 2021-11-30 | 2022-03-04 | 珠海一微半导体股份有限公司 | Fir滤波器复用系统 |
-
2024
- 2024-04-03 CN CN202410397142.XA patent/CN117991243B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1323105A (zh) * | 2001-03-19 | 2001-11-21 | 深圳市中兴集成电路设计有限责任公司 | 一种相关器装置 |
CN101278491A (zh) * | 2005-09-28 | 2008-10-01 | 日本电气株式会社 | 调制器、滤波器、滤波器增益控制方法和代码调制方法 |
CN1859038A (zh) * | 2006-01-11 | 2006-11-08 | 华为技术有限公司 | 一种时隙同步实现装置和方法 |
JP2009153096A (ja) * | 2007-11-28 | 2009-07-09 | Sony Corp | 受信装置、受信方法、およびプログラム |
CN101282322A (zh) * | 2008-03-05 | 2008-10-08 | 中科院嘉兴中心微系统所分中心 | 一种应用于无线中程传感网物理层的内插数字滤波器装置 |
CN101262240A (zh) * | 2008-04-25 | 2008-09-10 | 浙江大学 | 一种易于硬件实现的全数字频率变换的方法及其装置 |
CN102739195A (zh) * | 2012-06-06 | 2012-10-17 | 华为技术有限公司 | 一种fir滤波器的处理方法、装置和系统 |
CN105741313A (zh) * | 2014-12-09 | 2016-07-06 | 南京理工大学 | 一种基于fpga的快速实时动目标检测系统 |
CN108900177A (zh) * | 2018-06-14 | 2018-11-27 | 福建星海通信科技有限公司 | 一种fir滤波器及其对数据进行滤波的方法 |
CN110058201A (zh) * | 2019-04-22 | 2019-07-26 | 西安电子工程研究所 | 一种基于fpga资源复用实现多波形多速率时域脉压的方法 |
CN113317801A (zh) * | 2021-04-23 | 2021-08-31 | 清华大学 | 单通道的脑电图信号处理装置及方法 |
CN114142829A (zh) * | 2021-11-30 | 2022-03-04 | 珠海一微半导体股份有限公司 | Fir滤波器复用系统 |
Non-Patent Citations (7)
Title |
---|
吴太亮;刘峥;: "基于FPGA的时域脉冲压缩器研究", 制导与引信, no. 04, 15 December 2007 (2007-12-15) * |
王本明;赵前晟;丁海锋;罗丰;: "基于CPCI总线的通用FPGA信号处理板的设计", 电子科技, no. 10, 15 October 2009 (2009-10-15) * |
胡莉, 曾高荣, 陈红艳: "基于时分复用乘法器的FIR数字滤波器的设计", 西南科技大学学报, no. 02, 30 June 2004 (2004-06-30), pages 1 - 4 * |
陆聪;杨维明;王旭光;曾张帆;: "NLFM脉冲压缩及其FPGA时域实现", 湖北大学学报(自然科学版), no. 04, 5 July 2018 (2018-07-05) * |
陈进: "介绍一种实现可编程非递归数字滤波器(FIR)的研制方法", 计算机应用研究, no. 03, 30 May 1986 (1986-05-30), pages 1 - 4 * |
陈风波;董奇才;李念军;鲍振;: "乘法器复用技术在滑窗FIR滤波处理中的运用", 微计算机信息, no. 08, 15 March 2010 (2010-03-15) * |
靳鹏: "某雷达信号采集处理板的设计与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 4, 15 April 2013 (2013-04-15), pages 2 - 5 * |
Also Published As
Publication number | Publication date |
---|---|
CN117991243B (zh) | 2024-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108681984B (zh) | 一种3*3卷积算法的加速电路 | |
US5249578A (en) | Ultrasound imaging system using finite impulse response digital clutter filter with forward and reverse coefficients | |
CN110058201B (zh) | 一种基于fpga资源复用实现多波形多速率时域脉压的方法 | |
CN111337890B (zh) | 一种lfmcw雷达目标回波信号模拟方法 | |
EP1176718A2 (en) | Hardware accelerator for normal least-mean-square algorithm-based coefficient adaptation | |
CN102055435B (zh) | 一种窄带数字滤波器 | |
CN102349235A (zh) | 利用前向-后向滤波获得期望的非零相移 | |
CN101190136B (zh) | 实时产生滤波器系数的方法和装置 | |
JPH0418808A (ja) | 自動等化器及び半導体集積回路 | |
CN117991243B (zh) | 一种基于fpga的声纳信号匹配滤波处理方法及其系统 | |
CN115033293A (zh) | 零知识证明硬件加速器及生成方法、电子设备和存储介质 | |
CN109633613A (zh) | 一种高超声速平台联合脉冲压缩和弹速补偿的fpga实现方法 | |
CN101912277B (zh) | 一种基于流水线设计的实时数字正交解调的方法与装置 | |
CN102353940B (zh) | 基于fpga的脉冲压缩优化方法 | |
CN116111982A (zh) | 一种基于fpga的iir滤波方法、存储介质及装置 | |
CN114185014B (zh) | 一种应用于雷达信号处理的并行卷积方法及装置 | |
US9449597B2 (en) | Ultrasound receiving module, method and system | |
CN114327364A (zh) | 水声信号的高精度连续延时型fpga实时波束形成器 | |
CN210466088U (zh) | 微控制器和硬件模块 | |
US5271276A (en) | Phase regulating apparatus of ultrasonic measuring devices | |
US5864515A (en) | Sonar data processing | |
CN109921800A (zh) | 基于自定制浮点格式的频域脉冲压缩方法 | |
CN110808935B (zh) | 线性调频信号自相关运算的精确高效实现方法及装置 | |
CN100367915C (zh) | 基于反向初始化iir的彩色血流成像壁滤波方法和装置 | |
CN107911099B (zh) | 一种数字成形滤波方法及滤波装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |