CN107632292A - 一种对雷达信号进行调频傅立叶变换的方法 - Google Patents
一种对雷达信号进行调频傅立叶变换的方法 Download PDFInfo
- Publication number
- CN107632292A CN107632292A CN201710862287.2A CN201710862287A CN107632292A CN 107632292 A CN107632292 A CN 107632292A CN 201710862287 A CN201710862287 A CN 201710862287A CN 107632292 A CN107632292 A CN 107632292A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- mfrac
- signal
- frequency modulation
- 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
Abstract
本发明公开了一种对雷达信号进行调频傅立叶变换的方法,首先将连续调频傅立叶变换离散化,得到离散调频傅立叶变换的表达式,再将表达式中的指数因子展开成对一个复数进行离散傅立叶变换的形式。然后依据离散傅立叶变换的线性性质,对展开后的复数的实部和虚部分别求傅立叶变换,这里为了减少计算量,采用快速傅立叶变换来代替离散傅立叶变换。之后将原虚部经过快速傅立叶变换的结果乘上j,再和原实部经过快速傅立叶变换的结果相加。最后整理算术式,得到了一个复数形式的离散调频傅立叶变换的结果。本算法主要目的是估计信号的参数,分析信号成分。
Description
技术领域
本方法涉及雷达信号处理,主要是在离散调频傅立叶变换的基础上,提出一种利用快速傅立叶变换来实现离散调频傅立叶变换的算法。
背景技术
线性调频信号指持续期间频率连续线性变化的信号,是一种常用的雷达信号。由于其优点显著,应用十分广泛。线性调频信号具有峰值功率低,工作频率范围宽,距离分辨力高,作用距离远,信号处理起来相对简单等优点。
傅立叶变换是常用的处理信号的方法,其是将信号从时域搬到频域上,经过傅立叶变换可以得到信号的频谱图,进而分析出原始信号的频率以及其频率所对的幅值。但是对于调频信号来说,他是变频信号,如果继续使用傅立叶变换分析其频谱,存在局限性。
本发明提出来调频傅立叶变换的方法,该方法通过匹配傅立叶变换的一种极坐标的形式。主要是对待估计的参数取合适的范围,逐点计算进行搜索,然后找到信号原设定的参数值。
发明内容
鉴于线性调频信号的优点,现有技术中都提出了调频傅立叶变换的方法来处理调频信号。本发明的目的在于在现有的基础上提出一种对雷达信号进行调频傅立叶变换的方法首先是将连续调频信号离散化,然后将离散化后的信号进行展开推导,利用快速傅立叶变换的来实现该信号的离散调频傅立叶变换。
为了实现上述目的,本发明所采用的技术方案为一种对雷达信号进行调频傅立叶变换的方法,本方法对雷达信号进行离散调频傅立叶变换处理,对雷达信号进行离散调频傅立叶变换处理的目的是估计出雷达信号的参数。用fs表示雷达信号的采样周期,N表示雷达信号的采样序列的长度,k表示采样序列号,初始频率表示为离散调频信号表示为l′,k和l′是要估计的参数。
由于要处理的雷达信号不只包含一对含有参数k和l′的调频信号,所以要根据实际情况设定阈值,超过阈值的值所对应的参数都是所估计出得雷达信号的成分。
第一步将连续调频傅立叶变换的定义式离散化,得到离散调频傅立叶变换的表达式。
第二步将得到离散调频傅立叶变换的表达式中的指数因子展开,展成对一个复数进行离散傅立叶变换的形式。
第三步将第二步中得到的复数的实部和虚部分开。
第四步依据离散傅立叶变换的线性性质,对上述展开后的复数的实部和虚部分别求傅立叶变换,为减少计算量,采用快速傅立叶变换代替离散傅立叶变换。
第五步实部的快速傅立叶变换的结果与虚部快速傅立叶变换的结果乘上j之后相加。
第六步整理上一步的表达式,将表达式展开实部和虚部分别合并,得到一个复数形式。
第七步对上一步得到的复数求幅值。
第八步对参数k和l′取合适的范围,在所取得范围内分别取N个点,逐点进行前七步的运算,需要计算N×N次。
第九步根据具体情况设定阈值,超过阈值的幅值对应的参数k和l‘是所需要的估计值,估计值可能为多个。
与现有技术相比较,本发明具有如下有益效果。
(1)经过化简推导后的方法,又利用了快速傅立叶变换,所以大大地减少了计算量。
(2)由于是基于快速傅立叶变换来实现离散调频傅立叶变换,所以在很多开发软件上可以直接调用已有的快速傅立叶变换函数,提高了算法的可行性。
附图说明
图1离散调频傅立叶变换算法流程。
具体实施方式
下面结合推导过程对本发明作进一步说明。线性调频信号是常用的雷达信号,本发明主要以线性调频信号为例讲述具体的实施方式。
首先给出调频傅立叶变换的定义式:
其中l代表调频率,即频率线性变换的斜率;j是虚数;f是初始频率。F(l,f)表示调频傅立叶变换式;t是时间变量。
线性调频信号的表达式为:
f0是要估计参数值的线性调频信号的初始频率,l0是要估计参数值的线性调频信号的调频率。
对线性调频信号离散化,dt→Ts,Ts是采样周期;t→nTs,n是f(t)的函数离散f(n)的序列号;
fs是采样频率,N表示序列长度,取2的指数次幂。任意一点k代表的频率表示为 l′是离散调频率。调频信号离散化后的表达式为:
k0是线性调频信号的待估计的参数值f0离散化后的关系量n是f(t)的函数离散f(n)的序列号、l’0是离散调频率。
离散化后的调频傅里叶变换式为:
化简式(4)后:
对指数项展开:
令
得:
由式(7)可知,将离散调频傅立叶变换转化成对f‘(n)求傅立叶变换。
将(3)式代入(7)式得:
对上式进一步展开得:
令 然后对f’1(n),f’2(n)分别求快速傅立叶变换得到F‘1(k)、F‘2(k)。f’1(n)是f′(n)的实数部分,f’2(n)是f′(n)的虚数部分、F‘1(k)是f’1(n)快速傅立叶变换后的结果、F‘2(k)是f’2(n)快速傅立叶变换后的结果。最后F(l’,k)=F‘1(k)+jF‘2(k)。
调频傅立叶变换是匹配傅立叶变换的一种极坐标形式形式,主要是通过对l‘和k这两个调频参数进行扫描,当l’=l‘0,k=k0时,l’-k的幅度谱会出现尖峰。
对l‘和k这两个参数的扫描范围和步长的选取很关键,下面举例说明一下,在仿真中假设l‘0=25,k0=5,搜索范围l‘0∈[0,31],k0∈[0,31],需要计算32×32次。调频率为负的时候同样适用,假设l‘0=-23,k0=5,搜索范围l‘0∈[-31,0],k0∈[0,31],需要计算32×32次。搜索范围的选取,经过多次实验发现选取的搜索数应为大于l‘0的最小的2的指数次幂,这样的尖峰最明显,效果最佳。
如图1所示,其中×表示乘法器;+表示加法器;N是2的指数次幂,是变换点数;l‘表示离散调频率;f(n)是线性调频信号;n是序列的序号;经过推倒是对进行调频傅立叶变换
Claims (2)
1.一种对雷达信号进行调频傅立叶变换的方法,其特征在于:本方法对雷达信号进行离散调频傅立叶变换处理,对雷达信号进行离散调频傅立叶变换处理的目的是估计出雷达信号的参数;用fs表示雷达信号的采样周期,N表示雷达信号的采样序列的长度,k表示采样序列号,初始频率表示为离散调频信号表示为l′,k和l′是要估计的参数;
由于要处理的雷达信号不只包含一对含有参数k和l′的调频信号,所以要根据实际情况设定阈值,超过阈值的值所对应的参数都是所估计出得雷达信号的成分;
第一步将连续调频傅立叶变换的定义式离散化,得到离散调频傅立叶变换的表达式;
第二步将得到离散调频傅立叶变换的表达式中的指数因子展开,展成对一个复数进行离散傅立叶变换的形式;
第三步将第二步中得到的复数的实部和虚部分开;
第四步依据离散傅立叶变换的线性性质,对上述展开后的复数的实部和虚部分别求傅立叶变换,为减少计算量,采用快速傅立叶变换代替离散傅立叶变换;
第五步实部的快速傅立叶变换的结果与虚部快速傅立叶变换的结果乘上j之后相加;
第六步整理上一步的表达式,将表达式展开实部和虚部分别合并,得到一个复数形式;
第七步对上一步得到的复数求幅值;
第八步对参数k和l′取合适的范围,在所取得范围内分别取N个点,逐点进行前七步的运算,需要计算N×N次;
第九步根据具体情况设定阈值,超过阈值的幅值对应的参数k和l‘是所需要的估计值,估计值可能为多个。
2.根据权利要求1所述的一种对雷达信号进行调频傅立叶变换的方法,其特征在于:首先给出调频傅立叶变换的定义式:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>lt</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中l代表调频率,即频率线性变换的斜率;j是虚数;f是初始频率;F(l,f)表示调频傅立叶变换式;t是时间变量;
线性调频信号的表达式为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>l</mi>
<mn>0</mn>
</msub>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
f0是要估计参数值的线性调频信号的初始频率,l0是要估计参数值的线性调频信号的调频率;
对线性调频信号离散化,dt→Ts,Ts是采样周期;t→nTs,n是f(t)的函数离散f(n)的序列号;
任意一点k代表的频率表示为 l′是离散调频率;调频信号离散化后的表达式为:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mi>n</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>l</mi>
<mn>0</mn>
<mo>,</mo>
</msubsup>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
k0是线性调频信号的待估计的参数值f0离散化后的关系量n是f(t)的函数离散f(n)的序列号、l’0是离散调频率;
离散化后的调频傅里叶变换式为:
化简式(4)后:
对指数项展开:
令
得:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>l</mi>
<mo>,</mo>
</msup>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msup>
<mi>f</mi>
<mo>,</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo>&times;</mo>
<mi>k</mi>
<mo>&times;</mo>
<mi>n</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
由式(7)可知,将离散调频傅立叶变换转化成对;
将(3)式代入(7)式得:
对上式进一步展开得:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>f</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>T</mi>
<mi>s</mi>
<mo>&times;</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mi>n</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>l</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mo>&times;</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>&times;</mo>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
</mrow>
<mi>N</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mo>&times;</mo>
<mi>T</mi>
<mi>s</mi>
<mo>&times;</mo>
<mi>sin</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo>(</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mi>n</mi>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>l</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&pi;</mi>
<mo>&times;</mo>
<msup>
<mi>l</mi>
<mo>&prime;</mo>
</msup>
<mo>&times;</mo>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
</mrow>
<mi>N</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
令 然后对f’1(n),f’2(n)分别求快速傅立叶变换得到F‘1(k)、F‘2(k);f’1(n)是f′(n)的实数部分,f’2(n)是f′(n)的虚数部分、F‘1(k)是f’1(n)快速傅立叶变换后的结果、F‘2(k)是f’2(n)快速傅立叶变换后的结果;最后F(l’,k)=F‘1(k)+jF‘2(k);
调频傅立叶变换是匹配傅立叶变换的一种极坐标形式形式,主要是通过对l‘和k这两个调频参数进行扫描,当l’=l‘0,k=k0时,l′-k的幅度谱会出现尖峰。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710862287.2A CN107632292B (zh) | 2017-09-21 | 2017-09-21 | 一种对雷达信号进行调频傅立叶变换的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710862287.2A CN107632292B (zh) | 2017-09-21 | 2017-09-21 | 一种对雷达信号进行调频傅立叶变换的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107632292A true CN107632292A (zh) | 2018-01-26 |
CN107632292B CN107632292B (zh) | 2021-07-30 |
Family
ID=61101587
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710862287.2A Active CN107632292B (zh) | 2017-09-21 | 2017-09-21 | 一种对雷达信号进行调频傅立叶变换的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107632292B (zh) |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1268232A (zh) * | 1997-07-02 | 2000-09-27 | 艾利森电话股份有限公司 | 离散傅立叶变换(dft)及其反变换(idft)高效计算的方法与设备 |
EP1176516A2 (en) * | 2000-07-18 | 2002-01-30 | Victor Company Of Japan Limited | Recursive discrete fourier transformation apparatus |
CN101388001A (zh) * | 2008-06-25 | 2009-03-18 | 天津大学 | 基于全相位fft的高精度瞬间相位估计方法 |
CN101706770A (zh) * | 2009-11-03 | 2010-05-12 | 北京天童芯源科技有限公司 | 包含四条指令的支持快速傅里叶变换运算的方法 |
CN101833035A (zh) * | 2010-04-19 | 2010-09-15 | 天津大学 | 线性调频信号参数估计方法及其实施装置 |
CN102510363A (zh) * | 2011-09-30 | 2012-06-20 | 哈尔滨工程大学 | 一种强干扰源环境下的lfm信号检测方法 |
CN102707275A (zh) * | 2012-05-25 | 2012-10-03 | 北京航空航天大学 | 一种线性调频连续波雷达高度表的数字化处理方法 |
CN103262067A (zh) * | 2011-10-31 | 2013-08-21 | 深圳市英威腾电气股份有限公司 | 一种数据处理方法、数据处理装置和通信系统 |
CN104062655A (zh) * | 2014-06-27 | 2014-09-24 | 中国农业科学院农业资源与农业区划研究所 | 一种从雷达原始数据中提取运动参数做干涉测量的方法 |
US20150160374A1 (en) * | 2013-04-18 | 2015-06-11 | Weatherlink Co., Ltd | Weather information signal processing module |
US20160047891A1 (en) * | 2014-08-15 | 2016-02-18 | U.S.A. As Represented By The Administrator Of The National Aeronautics And Space Administration | HYPERFINE INTERPOLATED RANGE FINDING FOR CW LIDAR, RADAR, and SONAR USING REPEATING WAVEFORMS AND FOURIER TRANSFORM REORDERING |
US20160054438A1 (en) * | 2014-08-25 | 2016-02-25 | Texas Instruments Incorporated | Vibration parameters monitoring using fmcw radar |
DE102014116452A1 (de) * | 2014-11-11 | 2016-05-12 | Infineon Technologies Ag | Verfahren und Vorrichtung zur Verarbeitung von Radarsignalen |
CN105783974A (zh) * | 2014-12-25 | 2016-07-20 | 中国科学院声学研究所 | 一种线性调频信号的检测、参数估计方法及系统 |
CN105874351A (zh) * | 2013-12-06 | 2016-08-17 | 美国休斯研究所 | 编码孔径雷达(car)信号的处理方法及装置 |
CN106546949A (zh) * | 2016-10-28 | 2017-03-29 | 东南大学 | 一种基于频率预估计的双阵元正弦信号来波方向估计方法 |
US20170097405A1 (en) * | 2015-10-05 | 2017-04-06 | Analog Devices, Inc. | Scaling fixed-point fast fourier transforms in radar and sonar applications |
CN106840310A (zh) * | 2017-02-17 | 2017-06-13 | 水利部南京水利水文自动化研究所 | 一种连续调频波制雷达水位计测量方法 |
CN106842163A (zh) * | 2017-03-14 | 2017-06-13 | 西安电子科技大学 | 一种弹道目标回波信号时频特性估计方法 |
CN107153189A (zh) * | 2017-04-18 | 2017-09-12 | 上海交通大学 | 线性调频连续波雷达测距的信号处理方法 |
-
2017
- 2017-09-21 CN CN201710862287.2A patent/CN107632292B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1268232A (zh) * | 1997-07-02 | 2000-09-27 | 艾利森电话股份有限公司 | 离散傅立叶变换(dft)及其反变换(idft)高效计算的方法与设备 |
EP1176516A2 (en) * | 2000-07-18 | 2002-01-30 | Victor Company Of Japan Limited | Recursive discrete fourier transformation apparatus |
CN101388001A (zh) * | 2008-06-25 | 2009-03-18 | 天津大学 | 基于全相位fft的高精度瞬间相位估计方法 |
CN101706770A (zh) * | 2009-11-03 | 2010-05-12 | 北京天童芯源科技有限公司 | 包含四条指令的支持快速傅里叶变换运算的方法 |
CN101833035A (zh) * | 2010-04-19 | 2010-09-15 | 天津大学 | 线性调频信号参数估计方法及其实施装置 |
CN102510363A (zh) * | 2011-09-30 | 2012-06-20 | 哈尔滨工程大学 | 一种强干扰源环境下的lfm信号检测方法 |
CN103262067A (zh) * | 2011-10-31 | 2013-08-21 | 深圳市英威腾电气股份有限公司 | 一种数据处理方法、数据处理装置和通信系统 |
CN102707275A (zh) * | 2012-05-25 | 2012-10-03 | 北京航空航天大学 | 一种线性调频连续波雷达高度表的数字化处理方法 |
US20150160374A1 (en) * | 2013-04-18 | 2015-06-11 | Weatherlink Co., Ltd | Weather information signal processing module |
CN105874351A (zh) * | 2013-12-06 | 2016-08-17 | 美国休斯研究所 | 编码孔径雷达(car)信号的处理方法及装置 |
CN104062655A (zh) * | 2014-06-27 | 2014-09-24 | 中国农业科学院农业资源与农业区划研究所 | 一种从雷达原始数据中提取运动参数做干涉测量的方法 |
US20160047891A1 (en) * | 2014-08-15 | 2016-02-18 | U.S.A. As Represented By The Administrator Of The National Aeronautics And Space Administration | HYPERFINE INTERPOLATED RANGE FINDING FOR CW LIDAR, RADAR, and SONAR USING REPEATING WAVEFORMS AND FOURIER TRANSFORM REORDERING |
US20160054438A1 (en) * | 2014-08-25 | 2016-02-25 | Texas Instruments Incorporated | Vibration parameters monitoring using fmcw radar |
DE102014116452A1 (de) * | 2014-11-11 | 2016-05-12 | Infineon Technologies Ag | Verfahren und Vorrichtung zur Verarbeitung von Radarsignalen |
CN105783974A (zh) * | 2014-12-25 | 2016-07-20 | 中国科学院声学研究所 | 一种线性调频信号的检测、参数估计方法及系统 |
US20170097405A1 (en) * | 2015-10-05 | 2017-04-06 | Analog Devices, Inc. | Scaling fixed-point fast fourier transforms in radar and sonar applications |
CN106546949A (zh) * | 2016-10-28 | 2017-03-29 | 东南大学 | 一种基于频率预估计的双阵元正弦信号来波方向估计方法 |
CN106840310A (zh) * | 2017-02-17 | 2017-06-13 | 水利部南京水利水文自动化研究所 | 一种连续调频波制雷达水位计测量方法 |
CN106842163A (zh) * | 2017-03-14 | 2017-06-13 | 西安电子科技大学 | 一种弹道目标回波信号时频特性估计方法 |
CN107153189A (zh) * | 2017-04-18 | 2017-09-12 | 上海交通大学 | 线性调频连续波雷达测距的信号处理方法 |
Non-Patent Citations (3)
Title |
---|
MANKUN XU AND XIJIAN PING: "An improved fast Fourier transform in polar coordinate system", 《PROCEEDINGS OF THE 2003 IEEE INTERNATIONAL CONFERENCE ON INTELLIGENT TRANSPORTATION SYSTEMS, SHANGHAI, CHINA》 * |
徐楚林: "线性调频信号检测的快速调频变换算法", 《声学技术》 * |
陈晶: "线性调频信号与正弦调频信号参数估计方法", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN107632292B (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107103908B (zh) | 复调音乐多音高估计方法及伪双谱在多音高估计中的应用 | |
Shen et al. | Low-complexity sinusoidal-assisted EMD (SAEMD) algorithms for solving mode-mixing problems in HHT | |
CN104297740B (zh) | 基于相位分析的雷达目标多普勒谱估计方法 | |
CN104502698B (zh) | 电力信号的频率测量方法和系统 | |
CN103245937A (zh) | 基于微多普勒效应的微动目标特征提取方法 | |
CN106059973B (zh) | 频偏估计方法和系统 | |
CN105991492A (zh) | 跳频信号识别方法 | |
CN107886111A (zh) | 一种改进型局部均值分解方法 | |
Golmankhaneh et al. | Comparison of iterative methods by solving nonlinear Sturm-Liouville, Burgers and Navier-Stokes equations | |
CN108009347A (zh) | 基于同步压缩联合改进广义s变换的时频分析方法 | |
CN110632596A (zh) | 一种太赫兹sar多频振动误差补偿方法 | |
Batiha | Numerical solution of Bratu-type equations by the variational iteration method | |
Xiong et al. | Singularity power spectrum distribution | |
KR20150094332A (ko) | 스펙트럼 센싱 방법 및 스펙트럼 센싱 장치 | |
CN106646403B (zh) | K分布雷达杂波实时模拟方法及系统 | |
CN103744066A (zh) | 数字正交鉴相与匹配滤波的优化方法 | |
CN107632292A (zh) | 一种对雷达信号进行调频傅立叶变换的方法 | |
CN106770192A (zh) | 基于插值法的激光诱导击穿光谱连续背景校正方法 | |
US11126693B2 (en) | Method and device for reducing noise in a modulated signal | |
CN107315714A (zh) | 一种去卷积功率谱估计方法 | |
KR101167066B1 (ko) | 해수면 온도 극값을 이용한 기후 민감도 추정 방법 | |
Rudolf et al. | Error bounds of MCMC for functions with unbounded stationary variance | |
Harikrishnan et al. | Computing the multifractal spectrum from time series: an algorithmic approach | |
Geng et al. | Coupled Camassa–Holm equations, N-peakons and infinitely many conservation laws | |
JP5958453B2 (ja) | 無線局識別装置、無線局識別方法および無線局識別プログラム |
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 |