CN110441746A - 一种时域门变换方法和装置 - Google Patents
一种时域门变换方法和装置 Download PDFInfo
- Publication number
- CN110441746A CN110441746A CN201910767035.0A CN201910767035A CN110441746A CN 110441746 A CN110441746 A CN 110441746A CN 201910767035 A CN201910767035 A CN 201910767035A CN 110441746 A CN110441746 A CN 110441746A
- Authority
- CN
- China
- Prior art keywords
- domain
- time domain
- data
- frequency
- time
- 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
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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种时域门变换方法和装置,所述方法的一个实施方式包括:采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据;获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据;依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。该实施方式能够增强数据处理算法在目标图像局部分析环节上的灵活性,并避免数据截断造成的数据失真问题。
Description
技术领域
本发明涉及综合试验测量技术领域,尤其涉及一种时域门变换方法和装置。
背景技术
在扫频雷达散射截面RCS(Radar Cross Section)测量中,为了去除目标区域之外的杂散影响,经常会用到软件时域加门的方法对目标时域响应进行选取,在上述过程中一般采用快速傅里叶变换FFT(Fast Fourier Transform)及其逆变换进行时域频域之间的变换,采用FFT算法可以快速计算出全部点的离散傅里叶变换DFT(Discrete FourierTransform)值,即z变换(针对离散时间信号的一种变换方法,z表示复频率)在z平面单位圆上的全部等间隔取样值。但时域门变换过程中往往不需要计算整个单位圆上z变换的取样,只需要对不模糊窗口中的一段进行分析,这时希望图像的采样集中在这一窗口内,以获得较高的分辨率,而窗口以外的部分可不考虑。
因此,针对以上不足,需要提供一种更为灵活的时域门变换方法。
发明内容
本发明要解决的技术问题是:如何设计一种更为灵活、适应性更强的时域门变换方法。
为了解决上述技术问题,本发明提供了一种时域门变换方法。
本发明实施例的时域门变换方法包括:采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据;获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据;依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。
优选地,所述方法进一步包括:使用扫频测量方法采集待测目标的频域RCS数据;以及,在扫频测量过程中的不模糊窗口长度和待测目标区域长度符合预设条件时:所述频域时域变换方法为CZT逆变换,所述时域频域变换方法为CZT变换;在所述不模糊窗口长度和待测目标区域长度不符合预设条件时:所述频域时域变换方法为快速傅里叶逆变换,所述时域频域变换方法为快速傅里叶变换。
优选地,所述预设条件为:所述不模糊窗口长度与待测目标区域长度的比值不小于预设阈值。
优选地,所述获取用于时域选通的时域门函数包括:在时域构建矩形门函数,将矩形门函数变换到频域;将频域的矩形门函数与预先确定的窗函数相乘,并将相乘的结果变换到时域得到所述时域门函数。
优选地,所述对所述频域数据进行重归一化处理包括:在采集的所述频域RCS数据的每一频点建立单位点目标信号,形成单位点目标频域数据;其中,单位点目标信号的幅度为单位数值;将单位点目标频域数据使用所述频域时域变换方法转换为单位点目标时域数据,将所述时域门函数施加在所述单位点目标时域数据得到时域选通的单位点目标数据;将时域选通的单位点目标数据执行所述时域频域变换方法,得到重归一化函数;将所述时域选通数据转换得到的频域数据除以所述重归一化函数,从而去除该频域数据中的数据走样。
本发明还提供了一种时域门变换装置。
本发明实施例的时域门变换装置可包括:数据采集单元,用于采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据;时域选通单元,用于获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据;重归一化单元,用于依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。
优选地,数据采集单元可进一步用于:使用扫频测量方法采集待测目标的频域RCS数据;时域选通单元可进一步用于:在时域构建矩形门函数,将矩形门函数变换到频域;将频域的矩形门函数与预先确定的窗函数相乘,并将相乘的结果变换到时域得到所述时域门函数;以及,在扫频测量过程中的不模糊窗口长度与待测目标区域长度的比值不小于预设阈值时:所述频域时域变换方法为CZT逆变换,所述时域频域变换方法为CZT变换;在所述不模糊窗口长度与待测目标区域长度的比值小于预设阈值时:所述频域时域变换方法为快速傅里叶逆变换,所述时域频域变换方法为快速傅里叶变换。
优选地,重归一化单元可进一步用于:在采集的所述频域RCS数据的每一频点建立单位点目标信号,形成单位点目标频域数据;其中,单位点目标信号的幅度为单位数值;将单位点目标频域数据使用所述频域时域变换方法转换为单位点目标时域数据,将所述时域门函数施加在所述单位点目标时域数据得到时域选通的单位点目标数据;将时域选通的单位点目标数据执行所述时域频域变换方法,得到重归一化函数;将所述时域选通数据转换得到的频域数据除以所述重归一化函数,从而去除该频域数据中的数据走样。
本发明的上述技术方案具有如下优点:
其一,在扫频测量过程中的不模糊窗口长度远大于待测目标区域长度(即不模糊窗口长度与待测目标区域长度的比值大于或等于预设阈值)时,本发明采用CZT变换(Chirp-Z变换,也称为线性调频Z变换)及其逆变换方法取代传统的FFT方法。FFT是一种全局变换方法,其需要增加FFT点数来提高频谱的分辨率,增加FFT点数的实质是在整个单位圆上均匀增加采样点数,从而造成运算量的成倍增加。CZT可以只做局部变换,其能够对回波频谱或者时间窗口中的某段进行局部细化,可在采样点数增加不多的情况下提高雷达的测量精度。在不模糊窗口长度远大于待测目标区域长度时,如果分辨率相同,CZT的运算量远小于FFT的运算量。通过以上设置,本发明能够增强数据处理算法在目标图像局部分析环节上的灵活性。
其二,在本发明实施例中,将预先确定的窗函数增加到矩形门函数来构建用于时域选通的时域门函数并对变换后的频域数据进行重归一化处理,避免数据截断造成的数据失真问题。
附图说明
图1是本发明实施例中时域门变换方法的主要步骤示意图;
图2是本发明实施例中采集的待测目标频域RCS曲线示意图;
图3是本发明实施例中由图2信号变换得到的时域RCS曲线示意图;
图4是本发明实施例的时域门函数示意图;
图5是本发明实施例中使用时域门函数进行时域选通的示意图;
图6是本发明实施例的重归一化函数示意图;
图7是本发明实施例中经过时域门变换完整流程的待测目标频域RCS曲线示意图;
图8是根据本发明实施例中时域门变换装置的组成部分示意图;
图9是根据本发明实施例可以应用于其中的示例性系统架构图;
图10是用来实现本发明实施例中时域门变换方法的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1是本发明实施例中时域门变换方法的主要步骤示意图,如图1所示,本发明实施例的时域门变换方法具体执行以下步骤:
步骤S101:采集待测目标的频域RCS数据,依据频域时域变换方法将频域RCS数据转换为时域RCS数据。
具体地,实际场景中,可使用扫频测量方法采集待测目标的频域RCS数据。可以理解,RCS数据可以是回波电压数据,也可以是通过回波电压计算得到的RCS数据。一般地,RCS数据为复数,包含幅值信息和相位信息。
在扫频测量过程中,频率采样间隔能够确定不模糊窗口的长度。在电磁散射领域,如果待测目标的长度(一般为通过待测目标中心的最长距离)不大于不模糊窗口长度,则成像不会发生重叠或者模糊;如果待测目标长度大于不模糊窗口长度,则一般情况下成像存在重叠或者模糊。在本发明实施例中,如果不模糊窗口长度和待测目标区域长度(即待测目标的适宜测量距离,一般地,待测目标区域长度略大于待测目标的长度,例如,若待测目标长度为5米,待测目标区域长度可以是6米或7米)符合预设条件,上述频域时域变换方法为CZT逆变换,下文将要介绍的时域频域变换方法为CZT变换;如果不模糊窗口长度和待测目标区域长度不符合预设条件时:频域时域变换方法为快速傅里叶逆变换,时域频域变换方法为快速傅里叶变换。
可以理解,上述预设条件可根据需求灵活设置,例如可设置为:不模糊窗口长度与待测目标区域长度的比值不小于预设阈值,预设阈值可根据应用环境设置,例如可设置为10或者20等。由于CZT变换以及FFT变换属于已知技术,此处不再详述。通过以上设置,能够在不模糊窗口长度远大于待测目标区域长度时采用CZT变换及其逆变换,从而提升计算速度和局部细节分析灵活性。
例如,在某次微波暗室扫频RCS测试中,频率为12.4~18GHz,采样1601点,约等于3.5MHz步进,对应的不模糊窗口长度约42米,而暗室静区只有2米。这样,如果分辨率相同,FFT需要的点数是CZT的约21倍,如果CZT使用1601点,FFT就需要补零到34307点,二者的运算量分别是94332和258440,可见CZT运算量只相当于FFT运算量的约二点七分之一。
图2是本发明实施例中采集的待测目标频域RCS曲线示意图,图3是本发明实施例中由图2信号变换得到的时域RCS曲线示意图。在图2中,横坐标为测量的纵向距离,单位是米(m),纵坐标为频域RCS信号的幅度,单位是dB。可以理解,上述纵向距离与测量时间具有对应关系。在图3中,横坐标为频率,单位是GHz,纵坐标为频域RCS信号的幅度,单位是dB。
步骤S102:获取用于时域选通的时域门函数,将时域门函数施加在时域RCS数据得到时域选通数据。
作为一个优选方案,获取时域门函数的具体步骤如下:首先在时域构建矩形门函数,并将矩形门函数变换到频域。之后将频域的矩形门函数与预先确定的频域内的窗函数相乘,并将相乘的结果变换到时域得到上述时域门函数。可以理解,上述时域门函数等于矩形门函数和时域内窗函数的卷积。在本发明实施例中,窗函数用于对矩形门函数进行平滑处理,从而消除数据截断过程中的吉布斯效应,避免频谱泄露。实际应用中,窗函数可根据工程经验确定,例如可选用现有的汉宁窗函数、海明窗函数或者凯塞窗函数等。具体应用中,由于时域门函数对应于连续空间,因此需要将其离散化。
确定时域窗函数之后,可将时域门函数施加在步骤S101中得到的时域RCS数据从而获得时域选通数据。一般地,上述施加动作具体为:将时域RCS数据与时域门函数相乘。
图4是本发明实施例的时域门函数示意图,图5是本发明实施例中使用时域门函数进行时域选通的示意图,二图中横纵坐标的含义均与图3相同,此处不再重复。
步骤S103:依据时域频域变换方法将时域选通数据转换为频域数据,对频域数据进行重归一化处理以去除数据走样。
在本步骤中,首先将步骤S102中得到的时域选通数据转换为频域数据,之后利用以下步骤进行重归一化处理以去除门变换过程中的数据失真。可以理解,针对目标数据的重归一化指的是:建立单位数据并执行与目标数据完全相同的处理过程从而获得数据失真情况、由此对目标数据进行修正以去除失真的过程。
首先,在采集的上述频域RCS数据的每一频点(每一频点对应于一个具体频率)建立单位点目标信号,形成单位点目标频域数据。其中,单位点目标信号的幅度为单位数值,例如1dB。之后,将单位点目标频域数据使用上述频域时域变换方法转换为单位点目标时域数据,并将上述时域门函数施加在上述单位点目标时域数据得到时域选通的单位点目标数据。接着,将时域选通的单位点目标数据执行上述时域频域变换方法,即可得到重归一化函数。最后,将步骤S103中依据时域频域变换方法由时域选通数据转换得到的待测目标频域数据除以重归一化函数从而得到最终的处理结果,即时域门变换后的数据,也即去除杂散后的目标频率响应,同时该数据已经消除了门变换过程引入的失真。
图6是本发明实施例的重归一化函数示意图,图7是本发明实施例中经过时域门变换完整流程的待测目标频域RCS曲线示意图。图6中坐标的含义与图2相同,此处不再重复。图7的横坐标为频率,单位是GHz,纵坐标为RCS幅度,单位是dBsm(平方米分贝数)。在图7中,灰度较大的曲线为理论曲线,灰度较小的为测试曲线,可以看到二者非常接近。
需要说明的是,步骤S102中的“获取用于时域选通的时域门函数”可执行于步骤S101之前、同时或者之后,本发明对此不作限制。另外,步骤S103中的重归一化步骤可以在采集到待测目标的频域RCS数据后即执行,也可在其它时间点执行,本发明也不对此进行限制。此外,实际应用中,如果采集的多组数据施加的时域门函数一致,可首先针对多组数据执行重归一化步骤,这样可缩短计算时间。
综上所述,在本发明实施例的技术方案中,优化了微波暗室扫频RCS测量中常使用的时域门变换方法,通过构建时域门函数、针对变换后的频域数据进行重归一化处理等,可使频率数据更为真实地反映目标的频域响应特性。本发明核心部分可采用CZT变换实现,提高了局部细节分析的灵活性,处理速度和局部分析灵活性得到双重保证。此外,通过构造低旁瓣、合理上升下降过程的时域门函数,有效降低原始数据采用门变换后的频谱泄漏,增强RCS测量后数据处理方法的适用性和数据的准确性。
需要说明的是,对于前述的各方法实施例,为了便于描述,将其表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明并不受所描述的动作顺序的限制,某些步骤事实上可以采用其它顺序进行或者同时进行。此外,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定是实现本发明所必须的。
为便于更好的实施本发明实施例的上述方案,下面还提供用于实施上述方案的相关装置。
请参阅图8所示,本发明实施例提供的时域门变换装置800包括:数据采集单元801、时域选通单元802和重归一化单元803。
其中,数据采集单元801可用于采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据。时域选通单元802可用于获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据。重归一化单元803可用于依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。
在本发明实施例中,数据采集单元801可进一步用于:使用扫频测量方法采集待测目标的频域RCS数据。时域选通单元802可进一步用于:在时域构建矩形门函数,将矩形门函数变换到频域;将频域的矩形门函数与预先确定的窗函数相乘,并将相乘的结果变换到时域得到所述时域门函数。同时,在扫频测量过程中的不模糊窗口长度与待测目标区域长度的比值不小于预设阈值时:所述频域时域变换方法为CZT逆变换,所述时域频域变换方法为CZT变换;在所述不模糊窗口长度与待测目标区域长度的比值小于预设阈值时:所述频域时域变换方法为快速傅里叶逆变换,所述时域频域变换方法为快速傅里叶变换。
较佳地,在本发明实施例中,重归一化单元803可进一步用于:在采集的所述频域RCS数据的每一频点建立单位点目标信号,形成单位点目标频域数据;其中,单位点目标信号的幅度为单位数值;将单位点目标频域数据使用所述频域时域变换方法转换为单位点目标时域数据,将所述时域门函数施加在所述单位点目标时域数据得到时域选通的单位点目标数据;将时域选通的单位点目标数据执行所述时域频域变换方法,得到重归一化函数;将所述时域选通数据转换得到的频域数据除以所述重归一化函数,从而去除该频域数据中的数据走样。
图9示出了可以应用本发明实施例中时域门变换方法或时域门变换装置的示例性系统架构900。
如图9所示,系统架构900可以包括第一终端设备901、第二终端设备902、第三终端设备903、网络904和服务器905(此架构仅仅是示例,具体架构中包含的组件可以根据申请具体情况调整)。网络904用以在任一终端设备和服务器905之间提供通信链路的介质。网络904可以包括各种连接类型,例如有线、无线通信链路或者光纤电缆等等。
用户可以使用任一终端设备通过网络904与服务器905交互,以接收或发送消息等。任一终端设备上可以安装有各种通讯客户端应用,例如时域门变换计算应用等(仅为示例)。
任一终端设备可以是具有显示屏并且支持网页浏览的各种电子设备,包括但不限于智能手机、平板电脑、膝上型便携计算机和台式计算机等等。
服务器905可以是提供各种服务的服务器,例如对用户利用任一终端设备所操作的时域门变换计算应用提供支持的后台服务器(仅为示例)。后台服务器可以对接收到的数据变换请求等进行处理,并将处理结果(变换后的数据--仅为示例)反馈给终端设备。
需要说明的是,本发明实施例所提供的时域门变换方法一般由服务器905执行,相应地,时域门变换装置一般设置于服务器905中。
应该理解,图9中的终端设备、网络和服务器的数目仅仅是示意性的。根据实现需要,可以具有任意数目的终端设备、网络和服务器。
本发明还提供了一种电子设备。本发明实施例的电子设备包括:一个或多个处理器;存储装置,用于存储一个或多个程序,当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现本发明所提供的时域门变换方法。
下面参考图10,其示出了适于用来实现本发明实施例的电子设备的计算机系统1000的结构示意图。图10示出的电子设备仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。
如图10所示,计算机系统1000包括中央处理单元(CPU)1001,其可以根据存储在只读存储器(ROM)1002中的程序或者从存储部分1008加载到随机访问存储器(RAM)1003中的程序而执行各种适当的动作和处理。在RAM 1003中,还存储有计算机系统1000操作所需的各种程序和数据。CPU 1001、ROM 1002以及RAM 1003通过总线1004彼此相连。输入/输出(I/O)接口1005也连接至总线1004。
以下部件连接至I/O接口1005:包括键盘、鼠标等的输入部分1006;包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分1007;包括硬盘等的存储部分1008;以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分1009。通信部分1009经由诸如因特网的网络执行通信处理。驱动器1010也根据需要连接至I/O接口1005。可拆卸介质1011,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器1010上,以便从其上读出的计算机程序根据需要被安装入存储部分1008。
特别地,根据本发明公开的实施例,上文的主要步骤图描述的过程可以被实现为计算机软件程序。例如,本发明实施例包括一种计算机程序产品,其包括承载在计算机可读介质上的计算机程序,该计算机程序包含用于执行主要步骤图所示的方法的程序代码。在上述实施例中,该计算机程序可以通过通信部分1009从网络上被下载和安装,和/或从可拆卸介质1011被安装。在该计算机程序被中央处理单元1001执行时,执行本发明的系统中限定的上述功能。
需要说明的是,本发明所示的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。在本发明中,计算机可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述任意合适的组合。计算机可读信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括但不限于:无线、电线、光缆、RF等等,或者上述的任意合适的组合。
附图中的流程图和框图,图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,上述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这根据所涉及的功能而定。也要注意的是,框图或流程图中的每个方框、以及框图或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
描述于本发明实施例中所涉及到的单元可以通过软件的方式实现,也可以通过硬件的方式来实现。所描述的单元也可以设置在处理器中,例如,可以描述为:一种处理器包括数据采集单元、时域选通单元和重归一化单元。其中,这些单元的名称在某种情况下并不构成对该单元本身的限定,例如,数据采集单元还可以被描述为“向时域选通单元提供待测目标的时域RCS数据的单元”。
作为另一方面,本发明还提供了一种计算机可读介质,该计算机可读介质可以是上述实施例中描述的设备中所包含的;也可以是单独存在,而未装配入该设备中的。上述计算机可读介质承载有一个或者多个程序,当上述一个或者多个程序被该设备执行时,使得该设备执行的步骤包括:采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据;获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据;依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。
在本发明实施例的技术方案中,优化了微波暗室扫频RCS测量中常使用的时域门变换方法,通过构建时域门函数、针对变换后的频域数据进行重归一化处理等,可使频率数据更为真实地反映目标的频域响应特性。本发明核心部分可采用CZT变换实现,提高了局部细节分析的灵活性,处理速度和局部分析灵活性得到双重保证。此外,通过构造低旁瓣、合理上升下降过程的时域门函数,有效降低原始数据采用门变换后的频谱泄漏,增强RCS测量后数据处理方法的适用性和数据的准确性。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种时域门变换方法,其特征在于,包括:
采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据;
获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据;以及
依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。
2.根据权利要求1所述的时域门变换方法,其特征在于,所述方法进一步包括:使用扫频测量方法采集待测目标的频域RCS数据;以及
在扫频测量过程中的不模糊窗口长度和待测目标区域长度符合预设条件时:所述频域时域变换方法为CZT逆变换,所述时域频域变换方法为CZT变换;
在所述不模糊窗口长度和待测目标区域长度不符合预设条件时:所述频域时域变换方法为快速傅里叶逆变换,所述时域频域变换方法为快速傅里叶变换。
3.根据权利要求2所述的时域门变换方法,其特征在于,所述预设条件为:所述不模糊窗口长度与待测目标区域长度的比值不小于预设阈值。
4.根据权利要求1所述的时域门变换方法,其特征在于,所述获取用于时域选通的时域门函数,包括:
在时域构建矩形门函数,将矩形门函数变换到频域;
将频域的矩形门函数与预先确定的窗函数相乘,并将相乘的结果变换到时域得到所述时域门函数。
5.根据权利要求1-4任一所述的时域门变换方法,其特征在于,所述对所述频域数据进行重归一化处理,包括:
在采集的所述频域RCS数据的每一频点建立单位点目标信号,形成单位点目标频域数据;其中,单位点目标信号的幅度为单位数值;
将单位点目标频域数据使用所述频域时域变换方法转换为单位点目标时域数据,将所述时域门函数施加在所述单位点目标时域数据得到时域选通的单位点目标数据;
将时域选通的单位点目标数据执行所述时域频域变换方法,得到重归一化函数;以及
将所述时域选通数据转换得到的频域数据除以所述重归一化函数,从而去除该频域数据中的数据走样。
6.一种时域门变换装置,其特征在于,包括:
数据采集单元,用于采集待测目标的频域RCS数据,依据频域时域变换方法将所述频域RCS数据转换为时域RCS数据;
时域选通单元,用于获取用于时域选通的时域门函数,将所述时域门函数施加在所述时域RCS数据得到时域选通数据;以及
重归一化单元,用于依据时域频域变换方法将所述时域选通数据转换为频域数据,对所述频域数据进行重归一化处理以去除数据走样。
7.根据权利要求6所述的装置,其特征在于,
数据采集单元进一步用于:使用扫频测量方法采集待测目标的频域RCS数据;
时域选通单元进一步用于:在时域构建矩形门函数,将矩形门函数变换到频域;将频域的矩形门函数与预先确定的窗函数相乘,并将相乘的结果变换到时域得到所述时域门函数;以及
在扫频测量过程中的不模糊窗口长度与待测目标区域长度的比值不小于预设阈值时:所述频域时域变换方法为CZT逆变换,所述时域频域变换方法为CZT变换;在所述不模糊窗口长度与待测目标区域长度的比值小于预设阈值时:所述频域时域变换方法为快速傅里叶逆变换,所述时域频域变换方法为快速傅里叶变换。
8.根据权利要求6或7所述的装置,其特征在于,重归一化单元进一步用于:
在采集的所述频域RCS数据的每一频点建立单位点目标信号,形成单位点目标频域数据;其中,单位点目标信号的幅度为单位数值;将单位点目标频域数据使用所述频域时域变换方法转换为单位点目标时域数据,将所述时域门函数施加在所述单位点目标时域数据得到时域选通的单位点目标数据;将时域选通的单位点目标数据执行所述时域频域变换方法,得到重归一化函数;将所述时域选通数据转换得到的频域数据除以所述重归一化函数,从而去除该频域数据中的数据走样。
9.一种电子设备,其特征在于,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序,
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-5中任一所述的时域门变换方法。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现如权利要求1-5中任一所述的时域门变换方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910767035.0A CN110441746B (zh) | 2019-08-20 | 2019-08-20 | 一种时域门变换方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910767035.0A CN110441746B (zh) | 2019-08-20 | 2019-08-20 | 一种时域门变换方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110441746A true CN110441746A (zh) | 2019-11-12 |
CN110441746B CN110441746B (zh) | 2021-07-09 |
Family
ID=68436640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910767035.0A Active CN110441746B (zh) | 2019-08-20 | 2019-08-20 | 一种时域门变换方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110441746B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113983994A (zh) * | 2021-10-25 | 2022-01-28 | 北京环境特性研究所 | 样品材料参数确定方法及装置 |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61212781A (ja) * | 1985-03-18 | 1986-09-20 | Nec Corp | パルスドツプラレ−ダ方式 |
US20050099330A1 (en) * | 1999-05-25 | 2005-05-12 | Safe Zone Systems, Inc. | Object detection method and apparatus |
CN102967855A (zh) * | 2012-12-12 | 2013-03-13 | 北京航空航天大学 | 异地连续定标rcs测量中改进背景相减技术的信号处理方法 |
CN103646011A (zh) * | 2013-12-09 | 2014-03-19 | 中国电子科技集团公司第四十一研究所 | 一种基于线性调频z变换的信号频谱细化方法 |
CN103777185A (zh) * | 2014-01-15 | 2014-05-07 | 北京环境特性研究所 | 基于暗室扫频获取目标体雷达散射截面的方法及装置 |
CN103823215A (zh) * | 2014-03-03 | 2014-05-28 | 中国科学院电子学研究所 | 线性调频连续波雷达测距方法 |
CN203881941U (zh) * | 2014-04-14 | 2014-10-15 | 西安阿尔特测控技术有限公司 | 一种高铁路基沉降监测雷达装置 |
CN106324579A (zh) * | 2016-10-31 | 2017-01-11 | 上海无线电设备研究所 | 一种目标强散射源rcs提取算法 |
CN106546817A (zh) * | 2016-11-08 | 2017-03-29 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
CN106569191A (zh) * | 2016-10-17 | 2017-04-19 | 西北工业大学 | 一种利用高分辨率成像获取目标rcs的方法 |
US20170149590A1 (en) * | 2015-11-25 | 2017-05-25 | Korea Aerospace Research Institute | Method and system for inverse chirp-z transformation |
CN106872951A (zh) * | 2017-01-03 | 2017-06-20 | 北京环境特性研究所 | 一种基于压缩感知的暗室宽带rcs测量方法 |
CN107037416A (zh) * | 2017-04-13 | 2017-08-11 | 北京环境特性研究所 | 基于扫频rcs测量的时域门变换方法 |
CN107861117A (zh) * | 2017-10-17 | 2018-03-30 | 电子科技大学 | 一种适用于连续波周界监视雷达的多目标参数测量方法 |
CN109932719A (zh) * | 2019-03-18 | 2019-06-25 | 西安电子科技大学 | 基于sar成像的rcs高精度测量方法 |
CN110061790A (zh) * | 2019-04-03 | 2019-07-26 | 北京环境特性研究所 | 一种扫频rcs测试中抑制射频干扰的方法 |
-
2019
- 2019-08-20 CN CN201910767035.0A patent/CN110441746B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61212781A (ja) * | 1985-03-18 | 1986-09-20 | Nec Corp | パルスドツプラレ−ダ方式 |
US20050099330A1 (en) * | 1999-05-25 | 2005-05-12 | Safe Zone Systems, Inc. | Object detection method and apparatus |
CN102967855A (zh) * | 2012-12-12 | 2013-03-13 | 北京航空航天大学 | 异地连续定标rcs测量中改进背景相减技术的信号处理方法 |
CN103646011A (zh) * | 2013-12-09 | 2014-03-19 | 中国电子科技集团公司第四十一研究所 | 一种基于线性调频z变换的信号频谱细化方法 |
CN103777185A (zh) * | 2014-01-15 | 2014-05-07 | 北京环境特性研究所 | 基于暗室扫频获取目标体雷达散射截面的方法及装置 |
CN103823215A (zh) * | 2014-03-03 | 2014-05-28 | 中国科学院电子学研究所 | 线性调频连续波雷达测距方法 |
CN203881941U (zh) * | 2014-04-14 | 2014-10-15 | 西安阿尔特测控技术有限公司 | 一种高铁路基沉降监测雷达装置 |
US20170149590A1 (en) * | 2015-11-25 | 2017-05-25 | Korea Aerospace Research Institute | Method and system for inverse chirp-z transformation |
CN106569191A (zh) * | 2016-10-17 | 2017-04-19 | 西北工业大学 | 一种利用高分辨率成像获取目标rcs的方法 |
CN106324579A (zh) * | 2016-10-31 | 2017-01-11 | 上海无线电设备研究所 | 一种目标强散射源rcs提取算法 |
CN106546817A (zh) * | 2016-11-08 | 2017-03-29 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
CN106872951A (zh) * | 2017-01-03 | 2017-06-20 | 北京环境特性研究所 | 一种基于压缩感知的暗室宽带rcs测量方法 |
CN107037416A (zh) * | 2017-04-13 | 2017-08-11 | 北京环境特性研究所 | 基于扫频rcs测量的时域门变换方法 |
CN107861117A (zh) * | 2017-10-17 | 2018-03-30 | 电子科技大学 | 一种适用于连续波周界监视雷达的多目标参数测量方法 |
CN109932719A (zh) * | 2019-03-18 | 2019-06-25 | 西安电子科技大学 | 基于sar成像的rcs高精度测量方法 |
CN110061790A (zh) * | 2019-04-03 | 2019-07-26 | 北京环境特性研究所 | 一种扫频rcs测试中抑制射频干扰的方法 |
Non-Patent Citations (2)
Title |
---|
张淑溢: ""雷达散射截面积计量技术研究"", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
王万富: ""雷达目标RCS宽带测量技术和应用"", 《制导与引信》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113983994A (zh) * | 2021-10-25 | 2022-01-28 | 北京环境特性研究所 | 样品材料参数确定方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110441746B (zh) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schuet et al. | A model-based probabilistic inversion framework for characterizing wire fault detection using TDR | |
CN109697522A (zh) | 一种数据预测的方法和装置 | |
CN108650048B (zh) | 一种高精度数字阵列多通道延时补偿方法 | |
EP4314784A1 (en) | Single data set calibration and imaging with uncooperative electromagnetic inversion | |
CN109799550A (zh) | 用于预测降雨强度的方法和装置 | |
CN109710507A (zh) | 一种自动化测试的方法和装置 | |
US20230280286A1 (en) | Stored grain inventory management neural network | |
CN110208760A (zh) | 一种基于时域上采样的雷达回波仿真方法 | |
CN109002389A (zh) | 页面自动化测试的方法和装置 | |
CN110441746A (zh) | 一种时域门变换方法和装置 | |
CN118348529A (zh) | 一种弱散射体的定量衍射断层扫描方法及系统 | |
Olencki et al. | A low-cost system for far-field non-anechoic measurements of antenna performance figures | |
Hao et al. | Novel method for non-stationary signals via high-concentration time–frequency analysis using SSTFrFT | |
CN117056709A (zh) | 时序预测模型的训练方法、装置和存储介质及电子设备 | |
JPH11295365A (ja) | 電磁界解析方法 | |
Weber et al. | Cable fault characterization by time-domain analysis from S-parameter measurement and sparse inverse Chirp-Z transform | |
Xu et al. | Harmonic parameter online estimation in power system based on Hann self-convolving window and equidistant two-point interpolated DFT | |
Xie et al. | Finite word length optimization for spaceborne SAR imaging systems | |
CN113820565A (zh) | 暂态电能质量事件分析方法、装置和计算机设备 | |
Liu et al. | Energy and spectrum efficient radio frequency fingerprint intelligent blind identification | |
CN115685100A (zh) | 光学相控阵天线的相位校准方法及装置 | |
WO2016181549A1 (ja) | 性能評価装置及び性能評価プログラム | |
Zyari et al. | Identifying multiple reflections in distributed-lumped high-frequency structures | |
Selva | Efficient type‐4 and type‐5 non‐uniform FFT methods in the one‐dimensional case | |
CN112751634B (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 | ||
GR01 | Patent grant |