CN110187310A - 一种雷达设备lfm信号瞬时频率曲线拟合方法及*** - Google Patents

一种雷达设备lfm信号瞬时频率曲线拟合方法及*** Download PDF

Info

Publication number
CN110187310A
CN110187310A CN201910396983.8A CN201910396983A CN110187310A CN 110187310 A CN110187310 A CN 110187310A CN 201910396983 A CN201910396983 A CN 201910396983A CN 110187310 A CN110187310 A CN 110187310A
Authority
CN
China
Prior art keywords
data
lfm signal
discrete
radar equipment
cached
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.)
Pending
Application number
CN201910396983.8A
Other languages
English (en)
Inventor
段云鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201910396983.8A priority Critical patent/CN110187310A/zh
Publication of CN110187310A publication Critical patent/CN110187310A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00

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

本发明公开了一种雷达设备LFM信号瞬时频率曲线拟合方法及***,通过数据缓存模块、相位信息提取模块和瞬时频率曲线拟合模块实现。数据缓存模块对LFM信号的离散数据进行缓存,形成缓存数据;相位信息提取模块对缓存数据进行离散希尔伯特变换,并提取出相位信息;瞬时频率曲线拟合模块对瞬时频率数据进行多项式曲线拟合。本发明解决了传统LFM信号瞬时频率曲线拟合方法拟合精度不高的问题。

Description

一种雷达设备LFM信号瞬时频率曲线拟合方法及***
技术领域
本发明涉及一种雷达设备LFM信号瞬时频率曲线拟合方法及***。
背景技术
雷达设备工作时,需要对输入的线性调频LFM信号的瞬时频率曲线进行拟合。传统的LFM信号瞬时频率曲线拟合方法使用数字示波器对LFM信号进行采集、量化,提取出频率信息,采用逐点方式对频率值进行曲线拟合,此方法实现简单,但由于点与点之间采用直线的方式进行连接,拟合的精度不高。
发明内容
本发明目的在于提供一种雷达设备LFM信号瞬时频率曲线拟合方法,解决传统LFM信号瞬时频率曲线拟合方法拟合精度不高的问题。
有鉴于此,本发明提供的技术方案是:一种雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,包括:
对LFM信号的离散数据进行缓存,形成缓存数据;
对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;
对所述瞬时频率数据进行多项式曲线拟合。
进一步地,所述对LFM信号的离散数据进行缓存的步骤包括:在一个脉冲宽度τ的时间范围内对LFM信号的离散数据进行缓存,形成缓存数据x(n),数据速率为fs,其中n为时域点索引值,n=1,2,...,N,N为缓存数据x(n)的长度,N=fs·τ;LFM信号的调频带宽为B,调频斜率μ=B/τ。
进一步地,所述离散希尔伯特变换依据算法:xh(n)=hilbert[x(n)],得到离散希尔伯特变换结果xh(n),xh(n)为复数,hilbert[·]表示离散希尔伯特变换。
进一步地,所述提取出相位信息依据算法:
式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算。
进一步地,对所述获得瞬时频率数据依据算法:
得到瞬时频率数据f(n),式中为缓存数据x(n+1)的相位信息。
进一步地,所述多项式的次数设置为10。
本发明的另一目的还在于提供一种雷达设备LFM信号瞬时频率曲线拟合***,其特征在于,包括:
数据缓存模块,用于对LFM信号的离散数据进行缓存,形成缓存数据;
相位信息提取模块,用于对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;
瞬时频率曲线拟合模块,用于对所述瞬时频率数据进行多项式曲线拟合。
进一步地,所述提取出相位信息依据算法:
式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算。
本发明实现了以下显著的有益效果:
实现简单,包括:对LFM信号的离散数据进行缓存,形成缓存数据;对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;对所述瞬时频率数据进行多项式曲线拟合。解决传统LFM信号瞬时频率曲线拟合方法拟合精度不高的问题,经过多次测试,此种方法有效、可行,能够对雷达设备输入的LFM信号的瞬时频率曲线进行拟合,目前本方法已成功应用在雷达设备中,并且在各项试验中,雷达设备所拟合的瞬时频率曲线精度满足***测试要求。
附图说明
图1为本发明的雷达设备LFM信号瞬时频率曲线拟合方法的流程图;
图2为本发明的雷达设备LFM信号瞬时频率曲线拟合***的结构示意图。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明,根据下面说明和权利要求书,本发明的优点和特征将更清楚。需要说明的是,附图均采用非常简化的形式且均适用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
需要说明的是,为了清楚地说明本发明的内容,本发明特举多个实施例以进一步阐释本发明的不同实现方式,其中,该多个实施例是列举式而非穷举式。此外,为了说明的简洁,前实施例中已提及的内容往往在后实施例中予以省略,因此,后实施例中未提及的内容可相应参考前实施例。
虽然该发明可以以多种形式的修改和替换来扩展,说明书中也列出了一些具体的实施图例并进行详细阐述。应当理解的是,发明者的出发点不是将该发明限于所阐述的特定实施例,正相反,发明者的出发点在于保护所有给予由本权利声明定义的精神或范围内进行的改进、等效替换和修改。同样的元器件号码可能被用于所有附图以代表相同的或类似的部分。
请参照图1,本发明的一种雷达设备LFM信号瞬时频率曲线拟合方法,包括:步骤S101,对LFM信号的离散数据进行缓存,形成缓存数据;步骤S102,对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;步骤S103,对所述瞬时频率数据进行多项式曲线拟合。
优选地,所述对LFM信号的离散数据进行缓存的步骤包括:在一个脉冲宽度τ的时间范围内对LFM信号的离散数据进行缓存,形成缓存数据x(n),数据速率为fs,其中n为时域点索引值,n=1,2,...,N,N为缓存数据x(n)的长度,N=fs·τ;LFM信号的调频带宽为B,调频斜率μ=B/τ。
优选地,所述离散希尔伯特变换依据算法:xh(n)=hilbert[x(n)],得到离散希尔伯特变换结果xh(n),xh(n)为复数,hilbert[·]表示离散希尔伯特变换。
优选地,所述提取出相位信息依据算法:
式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算。
优选地,对所述获得瞬时频率数据依据算法:
得到瞬时频率数据f(n),式中为缓存数据x(n+1)的相位信息。
优选地,所述多项式的次数设置为10。
请参照图2,本发明的另一目的还在于提供一种雷达设备LFM信号瞬时频率曲线拟合***,包括:数据缓存模块,用于对LFM信号的离散数据进行缓存,形成缓存数据;相位信息提取模块,用于对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;瞬时频率曲线拟合模块,用于对所述瞬时频率数据进行多项式曲线拟合。
优选地,所述提取出相位信息依据算法:
式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算。
在一个实施例中,本发明的一种雷达设备LFM信号瞬时频率曲线拟合方法,其具体步骤为:
第一步 搭建瞬时频率曲线拟合***
瞬时频率曲线拟合***,包括:数据缓存模块、相位信息提取模块和瞬时频率曲线拟合模块;所述数据缓存模块的功能为:对LFM信号的离散数据进行缓存,形成缓存数据;相位信息提取模块的功能为:对缓存数据进行离散希尔伯特变换,并提取出相位信息;瞬时频率曲线拟合模块的功能为:对瞬时频率数据进行多项式曲线拟合。
第二步 数据缓存模块对LFM信号的离散数据进行缓存
数据缓存模块在一个脉冲宽度τ的时间范围内对LFM信号的离散数据进行缓存,形成缓存数据x(n),数据速率为fs,其中n为时域点索引值,n=1,2,...,N,N为缓存数据x(n)的长度,N=fs·τ;LFM信号的调频带宽为B,调频斜率μ=B/τ。
第三步 相位信息提取模块对缓存数据进行离散希尔伯特变换,并提取出相位信息
相位信息提取模块对缓存数据x(n)进行离散希尔伯特变换:xh(n)=hilbert[x(n)],得到离散希尔伯特变换结果xh(n),xh(n)为复数,hilbert[·]表示离散希尔伯特变换;
使用公式:
提取出相位信息式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算。
第四步 瞬时频率曲线拟合模块对瞬时频率数据进行多项式曲线拟合
瞬时频率曲线拟合模块使用公式:
得到瞬时频率数据f(n),式中为缓存数据x(n+1)的相位信息;
对瞬时频率数据f(n)进行多项式曲线拟合,多项式的次数设置为10,得到多项式拟合曲线f′(n)。
至此,实现了雷达设备LFM信号瞬时频率曲线的拟合。
本发明实现了以下显著的有益效果:
实现简单,包括:对LFM信号的离散数据进行缓存,形成缓存数据;对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;对所述瞬时频率数据进行多项式曲线拟合。解决传统LFM信号瞬时频率曲线拟合方法拟合精度不高的问题,经过多次测试,此种方法有效、可行,能够对雷达设备输入的LFM信号的瞬时频率曲线进行拟合,目前本方法已成功应用在雷达设备中,并且在各项试验中,雷达设备所拟合的瞬时频率曲线精度满足***测试要求。
根据本发明技术方案和构思,还可以有其他任何合适的改动。对于本领域普通技术人员来说,所有这些替换、调整和改进都应属于本发明所附权利要求的保护范围。

Claims (8)

1.一种雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,包括:
对LFM信号的离散数据进行缓存,形成缓存数据;
对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;
对所述瞬时频率数据进行多项式曲线拟合。
2.根据权利要求1所述的雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,所述对LFM信号的离散数据进行缓存的步骤包括:在一个脉冲宽度τ的时间范围内对LFM信号的离散数据进行缓存,形成缓存数据x(n),数据速率为fs,其中n为时域点索引值,n=1,2,...,N,N为缓存数据x(n)的长度,N=fs·τ;LFM信号的调频带宽为B,调频斜率μ=B/τ。
3.根据权利要求2所述的雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,所述离散希尔伯特变换依据算法:xh(n)=hilbert[x(n)],得到离散希尔伯特变换结果xh(n),xh(n)为复数,hilbert[·]表示离散希尔伯特变换。
4.根据权利要求3所述的雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,所述提取出相位信息依据算法:
式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算,xh(n)=hilbert[x(n)],得到离散希尔伯特变换结果xh(n),xh(n)为复数,hilbert[·]表示离散希尔伯特变换。
5.根据权利要求4所述的雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,对所述获得瞬时频率数据依据算法:
得到瞬时频率数据f(n),式中为缓存数据x(n+1)的相位信息。
6.根据权利要求5所述的雷达设备LFM信号瞬时频率曲线拟合方法,其特征在于,所述多项式的次数设置为10。
7.一种雷达设备LFM信号瞬时频率曲线拟合***,其特征在于,包括:
数据缓存模块,用于对LFM信号的离散数据进行缓存,形成缓存数据;
相位信息提取模块,用于对所述缓存数据进行离散希尔伯特变换,并提取出相位信息,获得瞬时频率数据;
瞬时频率曲线拟合模块,用于对所述瞬时频率数据进行多项式曲线拟合。
8.根据权利要求7所述的雷达设备LFM信号瞬时频率曲线拟合***,其特征在于,所述提取出相位信息依据算法:
式中real{·}表示复数取实部运算,imag{·}表示复数取虚部运算,xh(n)=hilbert[x(n)],得到离散希尔伯特变换结果xh(n),xh(n)为复数,hilbert[·]表示离散希尔伯特变换。
CN201910396983.8A 2019-05-13 2019-05-13 一种雷达设备lfm信号瞬时频率曲线拟合方法及*** Pending CN110187310A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910396983.8A CN110187310A (zh) 2019-05-13 2019-05-13 一种雷达设备lfm信号瞬时频率曲线拟合方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910396983.8A CN110187310A (zh) 2019-05-13 2019-05-13 一种雷达设备lfm信号瞬时频率曲线拟合方法及***

Publications (1)

Publication Number Publication Date
CN110187310A true CN110187310A (zh) 2019-08-30

Family

ID=67716187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910396983.8A Pending CN110187310A (zh) 2019-05-13 2019-05-13 一种雷达设备lfm信号瞬时频率曲线拟合方法及***

Country Status (1)

Country Link
CN (1) CN110187310A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147590A (zh) * 2020-08-25 2020-12-29 西安电子科技大学 一种基于响应估计频域拟合的通道均衡方法
CN112526430A (zh) * 2020-12-09 2021-03-19 中国航空工业集团公司北京长城计量测试技术研究所 一种飞机电源供电特性频率瞬变参数校准方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5574639A (en) * 1994-10-12 1996-11-12 National Instruments Corporation System and method for constructing filters for detecting signals whose frequency content varies with time
CN102508206A (zh) * 2011-10-31 2012-06-20 哈尔滨工程大学 基于小波包去噪和功率谱熵的线性调频信号参数估计方法
CN107526064A (zh) * 2017-07-10 2017-12-29 上海电机学院 基于二维特征的自适应lfm信号参数估计方法
CN108519511A (zh) * 2018-03-28 2018-09-11 电子科技大学 一种线性调频信号频率特征参数的时域测量方法
CN109188370A (zh) * 2018-10-29 2019-01-11 北京遥感设备研究所 一种雷达设备lfm脉冲信号包络曲线拟合方法及***
CN109283498A (zh) * 2017-07-21 2019-01-29 北京遥感设备研究所 一种线性调频脉冲信号相位误差曲线生成方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5574639A (en) * 1994-10-12 1996-11-12 National Instruments Corporation System and method for constructing filters for detecting signals whose frequency content varies with time
CN102508206A (zh) * 2011-10-31 2012-06-20 哈尔滨工程大学 基于小波包去噪和功率谱熵的线性调频信号参数估计方法
CN107526064A (zh) * 2017-07-10 2017-12-29 上海电机学院 基于二维特征的自适应lfm信号参数估计方法
CN109283498A (zh) * 2017-07-21 2019-01-29 北京遥感设备研究所 一种线性调频脉冲信号相位误差曲线生成方法
CN108519511A (zh) * 2018-03-28 2018-09-11 电子科技大学 一种线性调频信号频率特征参数的时域测量方法
CN109188370A (zh) * 2018-10-29 2019-01-11 北京遥感设备研究所 一种雷达设备lfm脉冲信号包络曲线拟合方法及***

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张西托,饶伟,黄信荣,陈珊珊: "《LFM 信号参数估计算法研究》", 《雷达科学与技术》 *
李风民: "《希尔伯特变换及其应用研究》", 《中国无线电》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147590A (zh) * 2020-08-25 2020-12-29 西安电子科技大学 一种基于响应估计频域拟合的通道均衡方法
CN112147590B (zh) * 2020-08-25 2023-04-07 西安电子科技大学 一种基于响应估计频域拟合的通道均衡方法
CN112526430A (zh) * 2020-12-09 2021-03-19 中国航空工业集团公司北京长城计量测试技术研究所 一种飞机电源供电特性频率瞬变参数校准方法

Similar Documents

Publication Publication Date Title
CN110187310A (zh) 一种雷达设备lfm信号瞬时频率曲线拟合方法及***
Butz Fourier transformation for pedestrians
Zhang et al. Extraction of a source signal whose kurtosis value lies in a specific range
DE102016118734B4 (de) Skalierung von Festkomma-Schnelle-Fourier-Transformationen bei Radar- und Sonaranwendungen
Abbiss et al. Regularized iterative and non-iterative procedures for object restoration from experimental data
EP1540844B1 (en) Methods for transmitting a waveform having a controllable attenuation and propagation velocity
CN109283498A (zh) 一种线性调频脉冲信号相位误差曲线生成方法
CN105509771B (zh) 一种发动机滑油金属颗粒在线监测的信号降噪方法
CN110554362A (zh) 一种雷达设备大时宽带宽积lfm信号相位谱合成***及方法
CN106341101B (zh) 群时延滤波器系数的计算方法以及滤波器的模拟***
Richards The discrete-time Fourier transform and discrete Fourier transform of windowed stationary white noise
Saneva et al. Some Abelian and Tauberian results for the short-time Fourier transform
DE112014006818B4 (de) Frequenzgang-Kompensation in einem Digital-Analog-Wandler
Holá et al. Relations between minimal usco and minimal cusco maps
Gillberg et al. Frequency domain identification of continuous-time output error models, Part II: Non-uniformly sampled data and B-spline output approximation
Song et al. An approximation for the inverse first passage time problem
Shah et al. Integration of the Indian Stock Market: at the angle of Time-Frequency
JP2023500363A (ja) ユニバーサルノッチフィルタ
CN109188370B (zh) 一种雷达设备lfm脉冲信号包络曲线拟合方法及***
Tanahashi et al. Fine scale eddy cluster and energy cascade in homogeneous isotropic turbulence
CN109343421A (zh) 一种线性调频信号的幅相补偿***及方法
Zhu et al. Anti-jamming method based on equivalent sampling for UWB pulse train signal
Trigub The exact order of approximation to periodic functions by Bernstein-Stechkin polynomials
Illing Fourier analysis
Azimi et al. Bayesian estimation for the burr type III distribution under type II Doubly censored data

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190830

RJ01 Rejection of invention patent application after publication