CN108196237A - 一种抑制fmcw雷达回波信号中寄生调幅的方法 - Google Patents
一种抑制fmcw雷达回波信号中寄生调幅的方法 Download PDFInfo
- Publication number
- CN108196237A CN108196237A CN201711428715.7A CN201711428715A CN108196237A CN 108196237 A CN108196237 A CN 108196237A CN 201711428715 A CN201711428715 A CN 201711428715A CN 108196237 A CN108196237 A CN 108196237A
- Authority
- CN
- China
- Prior art keywords
- sampling
- sequence
- points
- point
- signal
- 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 41
- 230000002401 inhibitory effect Effects 0.000 title abstract description 7
- 238000005070 sampling Methods 0.000 claims abstract description 121
- 230000003071 parasitic effect Effects 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 abstract description 23
- 238000001914 filtration Methods 0.000 abstract description 10
- 238000013461 design Methods 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 abstract description 4
- 230000009466 transformation Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000007704 transition 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
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
本发明公开了一种抑制FMCW雷达回波信号中寄生调幅的方法。首先,提取雷达回波信号x(n)的所有极值点,得到一次极大值序列和一次极小值序列,然后提取一次极大值序列中的局部极大值点,得到二次极大值序列;提取一次极小值序列中的局部极小值点,得到二次极小值序列。在整个时间采样轴上,采用三次样条插值函数对二次极大值序列和二次极小值序列进行拟合,得到上包络和下包络。计算上下包络的均值,作为寄生信号的骨架。再计算雷达回波信号与寄生信号骨架之差,作为滤除寄生调幅后的信号。本方法无需复杂的滤波器设计和傅里叶变换计算,效率更高。
Description
技术领域
本发明属于调频连续波(简称FMCW)雷达探测与应用技术领域,具体涉及一种抑制FMCW雷达回波信号中寄生调幅的方法。
背景技术
FMCW雷达是由调制信号发生器、压控振荡器(VCO)、混频器、收发天线等模块组成,主要应用于汽车防撞、无人机等方面,可以获得被测目标与雷达之间的距离和相对速度。FMCW雷达具有结构简单、发射功率低、分辨率和灵敏度高、天线部件尺寸小等优点。调制信号发生器产生的调制信号经过压控振荡器输出调频连续波,通常调制信号为锯齿波或三角波。FMCW从天线发射出去,遇到目标后反射回来与本振信号混频得到差频信号。对差频信号进行频率提取,经典方法是采用傅里叶变换对差频信号直接估计。通过提取的频率计算出雷达和目标间的相对速度和距离。FMCW雷达在工作时,调制电路产生三角波或锯齿波形式的调制信号,调制信号输入至压控振荡器,产生调频信号。通过功分器分为两路信号,一路经天线发射出去,一路传输给接收机用于混频。发射天线将调频信号发射出去,遇到目标后产生散射信号被接收天线接收。接收机将回波信号与发射信号混频,得到差频信号,再通过数据处理方法提取目标的距离和速度参数。理想情况下,压控振荡器输出等幅的调频信号。实际上,由于压控振荡器的非线性,使得调频信号的幅度不是恒定不变的,而是随频率变化的,即压控振荡器输出的信号是一个调频调幅信号。混频接收的时候,差频信号相当于叠加了一个频率与调制信号相同的寄生信号。由于该信号是FMCW雷达内部电路系统产生的,不随目标的距离速度变化而改变,称该信号为寄生调幅。
针对FMCW雷达探测系统中产生的寄生调幅,现有的方法是采用高通滤波器滤除系统中的低频部分。但是,利用双线性变化法及四阶贝塞尔滤波器滤除寄生信号时,由于模拟角频率和数字角频率的非线性关系,滤波效果有限【参考文献:陈丽.面向大型结构的调频连续波雷达位移测量系统VCO及回波的影响.重庆大学,2015】;利用八阶巴特沃斯高通滤波器滤除低频部分,过渡带较长【参考文献:齐国清.FMCW雷达近距离反射性能研究.大连海事大学,2016】;利用CIC-HB-FIR滤波器结构顺序滤波,3种滤波器的设计和运算过程较为复杂【参考文献:李朋.LFMCW汽车防撞雷达信号处理及硬件实现.电子科技大学,2016】。总而言之,现有的滤波处理方法,用硬件实现时,滤波器性能通常受外界温度及环境影响,精度难以提高。用软件实现时,计算过程复杂,通常的数字滤波器在四阶以上,运算复杂度较大。
鉴于滤波器去除FMCW回波信号中寄生调幅的缺点,有必要设计一种基于FMCW雷达回波信号的时域特征,有效滤除寄生调幅信号的方法。
发明内容
本发明针对现有技术的缺陷,提供一种抑制FMCW雷达回波信号中寄生调幅的方法,可以实现对调频连续波雷达回波信号中寄生调幅的抑制,与传统的基于频域滤波的寄生调幅抑制方法相比,本方法无需复杂的滤波器设计和傅里叶变换计算,计算速度快,效率更高;无需额外硬件条件。
一种抑制FMCW雷达回波信号中寄生调幅的方法,包括以下步骤:
步骤1:提取原始FMCW雷达回波信号x(n),n=0,1,2,…,N-1的所有局部极值点,由其中的局部极大值点构成一次极大值序列,其中的局部极小值点构成一次极小值序列;其中,N为原始FMCW雷达回波信号中的采样点个数(原始FMCW雷达回波信号的长度);
步骤2:提取一次极大值序列中的所有局部极大值点,得到二次极大值序列;提取一次极小值序列中的所有局部极小值点,得到二次极小值序列;
步骤3:求出二次极大值序列的三次样条插值函数,以拟合出上包络u(n);求出二次极小值序列的三次样条插值函数,以拟合出下包络l(n);
步骤4:计算上下包络的均值即为寄生信号的骨架;
步骤5:计算原始FMCW雷达回波信号x(n)与寄生信号的骨架之差y(n)=x(n)-h(n),得到抑制掉寄生调幅后的信号,实现对寄生调幅的抑制。
以下对各个步骤具体说明:
进一步地,所述步骤1中,通过差值比较法提取原始FMCW雷达回波信号x(n),n=0,1,2,…,N-1的所有局部极值点,具体处理过程为:遍历原始FMCW雷达回波信号中各个采样点,将每个采样点的采样值分别与该信号中与其相邻的两个采样点的采样值进行比较,以确定该采样点是否为局部极值点;若x(n)-x(n-1)>0且x(n)-x(n+1)≥0,则原始FMCW雷达回波信号中第n个采样点为局部极大值点;若x(n)-x(n-1)<0且x(n)-x(n+1)≤0,则第n个采样点为局部极小值点,其中n=1…,N-2。因为步骤3是根据离散的局部极大值/极小值进行三次样条插值,两头的采样值x(0)和x(N-1)不影响插值后的曲线,所以不需要对x(0)、x(N-1)进行判断。
进一步地,所述步骤2中,通过差值比较法提取一次极大值序列的所有局部极大值点,具体处理过程为:遍历一次极大值序列中的各个采样点,将每个采样点的采样值分别与该序列中与其相邻的两个采样点的采样值进行比较,以确定该采样点是否为局部极大值点;若某一采样点的采样值大于该序列中与其相邻的两个采样点的采样值,则该采样点为局部极大值点。
通过差值比较法提取一次极小值序列的所有局部极小值点,具体处理过程为:遍历一次极小值序列中的各个采样点,将每个采样点的采样值分别与该序列中与其相邻的两个采样点的采样值进行比较,以确定该采样点是否为局部极小值点;若某一采样点的采样值小于该序列中与其相邻的两个采样点的采样值,则该采样点为局部极小值点。
进一步地,所述步骤3中拟合出上包络u(n)具体包括以下步骤:
1)求解二次极大值序列的三次样条插值函数;
设二次极大值序列中采样点的个数为L,L个采样点对应的采样时间分别为x1,x2,…,xL,其中x1<x2…<xL,对应的采样值分别为:y1,y2,…,yL;
设三次样条插值函数的表达式为:
其中,Si(x)表示xi到xi+1的曲线,hi=xi+1-xi;Mi=Si″(x),Si″(x)为Si(x)的二阶导数,Mi的求解方法为:
利用Si(x)的二阶导数可以导出方程:
μiMi-1+2Mi+λiMi+1=di
其中,μi、λi和di均为中间参数;以上共得到L-2个方程;
根据Si(x)的边界条件:S1″′(x)=S2″′(x),SL-1″′(x)=SL″′(x),其中Si″′(x)为Si(x)的三阶导数,可得:
将上述条件带入方程中可得方程组:
其中,a=μ1/λ1,b=λL-1/μL-1;由此求解得出M1,M2,…,ML。
将所得解带入三次样条插值函数表达式,即可求得二次极大值序列的三次样条插值函数;
2)根据二次极大值序列的三次样条插值函数计算原始FMCW雷达回波信号中除二次极大值序列中的L个采样点以外的N-L个采样点对应的函数值;
对于每一个采样点,首先判断其对应的采样时间x所在区间,然后将x代入相应区间对应的三次样条插值函数计算该采样点对应的函数值,具体地:
若x≤x2,则将其代入S1(x)进行计算;
若x∈(xi,xi+1],则将其代入Si(x)进行计算;i∈[2,L-2];
若x>xL-1,则将其代入SL-1(x)进行计算;
3)将计算得到的N-L个采样点对应的函数值与二次极大值序列中L个采样点对应的采样值,按照采样时间进行排序,得到上包络u(n),n=0,1,2,…,N-1。
同样地,所述步骤3中拟合出下包络l(n)具体包括以下步骤:
ⅰ)求解二次极小值序列的三次样条插值函数;
设二次极小值序列中采样点的个数为K,K个采样点对应的采样时间分别为x1,x2,…,xK,其中x1<x2…<xK,对应的采样值分别为:y1,y2,…,yK;
设三次样条插值函数的表达式为:
其中,hi=xi+1-xi;Mi=Si″(x),Si″(x)为Si(x)的二阶导数,Mi的求解方法为:
利用Si(x)的二阶导数可以导出方程:
μiMi-1+2Mi+λiMi+1=di
其中,μi、λi和di均为中间参数;以上共得到K-2个方程;
根据Si(x)的边界条件:S1″′(x)=S2″′(x),SK-1″′(x)=SK″′(x),其中Si″′(x)为Si(x)的三阶导数,得到:
将上述条件带入方程中可得方程组:
其中,a=μ1/λ1,b=λK-1/μK-1;
由此求解得出M1,M2,…,MK,将所得解带入三次样条插值函数表达式,即可求得二次极小值序列的三次样条插值函数;
ⅱ)根据二次极小值序列的三次样条插值函数计算原始FMCW雷达回波信号中除二次极小值序列中的K个采样点以外的N-K个采样点对应的函数值;
对于每一个采样点,首先判断其对应的采样时间x所在区间,然后将x代入相应区间对应的三次样条插值函数计算该采样点对应的函数值,具体地:
若x≤x2,则将其代入S1(x)进行计算;
若x∈(xi,xi+1],则将其代入Si(x)进行计算;i∈[2,K-2];
若x>xK-1,则将其代入SK-1(x)进行计算;
ⅲ)将计算得到的N-K个采样点对应的函数值与二次极小值序列中K个采样点对应的采样值,按照采样时间进行排序,得到下包络l(n),n=0,1,2,…,N-1。
有益效果:
本发明提出一种抑制FMCW雷达回波中寄生调幅的方法,该方法基于FMCW雷达中频信号的时域表征形态,对寄生信号进行骨架提取,能有效抑制寄生调幅。相对于传统模拟滤波器受环境影响、精度不高、不稳定、高阶数字滤波器设计较为复杂、运算复杂度较大等缺点,本发明无需额外的硬件需求,计算速度快,抑制效果好,可应用于车载雷达,无人机等测速测距雷达系统的寄生调幅抑制中。
附图说明
为说明本发明的技术方案,下面对本发明的一个实施例加以描述。
图1示出了本发明方法的流程图。
图2示出了FMCW雷达系统回波信号采样后的序列。
图3示出了回波信号x(n)经傅里叶变换后的频谱图。
图4示出了一次极大值序列f1(i)分布图。
图5示出了一次极小值序列f2(j)分布图。
图6示出了二次极大值序列g1(p)分布图。
图7示出了二次极小值序列g2(q)分布图。
图8示出了提取的寄生信号骨架h(n)的波形。
图9示出了滤除寄生信号后的回波信号y(n)。
图10示出了经过本方法后的回波信号频谱图。
具体实施方式
以下将结合附图与具体实例对本发明的具体细节进行描述。
本发明通过对寄生信号的骨架提取从而抑制寄生调幅,如图1所示,FMCW雷达探测前方目标,将采样后的回波信号进行两次极值搜索,再求出二次极值数据,在整个时间采样轴上,采用三次样条插值函数对二次极大值序列和二次极小值序列进行拟合,得到上包络和下包络;然后求取上下包络的均值,将原始信号减去该均值,得到滤除寄生调幅后的信号。
实施例1:
本实例中,FMCW雷达探测前方3m目标。雷达发射信号的频率从24GHz到24.4GHz,带宽B=400MHz。发射信号采用三角波调制,调制周期T=4ms。回波信号与发射信号混频得到差频信号,通过模数转换器对差频信号采样,采样频率Fs=781.25KHz,采样点数N为3125点。从而,频率分辨率寄生信号的频率为如图2所示,采样后的信号有明显的寄生调幅现象。对采样后的信号进行傅里叶变换加以分析,傅里叶变换后的频谱图见图3。从图3可见,寄生信号的幅值非常大,以至于有用的中频信号被淹没。
采用本发明的方法,提取回波信号的极大值序列和极小值序列。由于杂波和噪声影响,一次极大值序列f1(i)如图4中的红点所示,一次极小值序列f2(j)如图5中的黑点所示。再提取一次极大值序列f1(i)的极大值序列,得到二次极大值序列g1(p),如图6中的黑点所示。提取一次极小值序列f2(j)的极小值序列,得到二次极小值序列g2(q),如图7中的黑点所示。对二次极大值序列g1(p)和二次极小值序列g2(q),求出三次样条插值函数并拟合出上下包络。对上下包络求均值,得到寄生信号骨架如图8中的黑线所示。将FMCW回波信号x(n)减去寄生信号骨架h(n),得到滤除寄生信号后的信号y(n),如图9所示。从图9中可以看出,由调制信号引起的寄生调幅现象已经被抑制。对y(n)进行傅里叶变换,其频谱图见图10。从图10中可以看出,寄生调幅信号的频谱分量已经被抑制。在CPU为Pentium(R)Dual-Core2.00GHz,内存为2.00GB的PC上运行,本发明所提的算法耗时200ms。采用传统的双线性法数字滤波进行寄生调幅抑制,算法耗时1.5s。
Claims (5)
1.一种抑制FMCW雷达回波信号中寄生调幅的方法,其特征在于,包括以下步骤:
步骤1:提取原始FMCW雷达回波信号x(n),n=0,1,2,…,N-1的所有局部极值点,由其中的局部极大值点构成一次极大值序列,其中的局部极小值点构成一次极小值序列;其中,N为原始FMCW雷达回波信号中的采样点个数(原始FMCW雷达回波信号的长度);
步骤2:提取一次极大值序列中的所有局部极大值点,得到二次极大值序列;提取一次极小值序列中的所有局部极小值点,得到二次极小值序列;
步骤3:求出二次极大值序列的三次样条插值函数,以拟合出上包络u(n);求出二次极小值序列的三次样条插值函数,以拟合出下包络l(n);
步骤4:计算上下包络的均值即为寄生信号的骨架;
步骤5:计算原始FMCW雷达回波信号x(n)与寄生信号的骨架之差y(n)=x(n)-h(n),得到抑制掉寄生调幅后的信号,实现对寄生调幅的抑制。
2.根据权利要求1所述的抑制FMCW雷达回波信号中寄生调幅的方法,其特征在于,所述步骤1中,通过差值比较法提取原始FMCW雷达回波信号x(n),n=0,1,2,…,N-1的所有局部极值点,具体处理过程为:遍历原始FMCW雷达回波信号中各个采样点,将每个采样点的采样值分别与该信号中与其相邻的两个采样点的采样值进行比较,以确定该采样点是否为局部极值点;若x(n)-x(n-1)>0且x(n)-x(n+1)≥0,则原始FMCW雷达回波信号中第n个采样点为局部极大值点;若x(n)-x(n-1)<0且x(n)-x(n+1)≤0,则第n个采样点为局部极小值点,其中n=1…,N-2。
3.根据权利要求1所述的抑制FMCW雷达回波信号中寄生调幅的方法,其特征在于,所述步骤2中,通过差值比较法提取一次极大值序列的所有局部极大值点,具体处理过程为:遍历一次极大值序列中的各个采样点,将每个采样点的采样值分别与该序列中与其相邻的两个采样点的采样值进行比较,以确定该采样点是否为局部极大值点;若某一采样点的采样值大于该序列中与其相邻的两个采样点的采样值,则该采样点为局部极大值点;
通过差值比较法提取一次极小值序列的所有局部极小值点,具体处理过程为:遍历一次极小值序列中的各个采样点,将每个采样点的采样值分别与该序列中与其相邻的两个采样点的采样值进行比较,以确定该采样点是否为局部极小值点;若某一采样点的采样值小于该序列中与其相邻的两个采样点的采样值,则该采样点为局部极小值点。
4.根据权利要求1所述的抑制FMCW雷达回波信号中寄生调幅的方法,其特征在于,所述步骤3中拟合出上包络u(n)具体包括以下步骤:
1)求解二次极大值序列的三次样条插值函数;
设二次极大值序列中采样点的个数为L,L个采样点对应的采样时间分别为x1,x2,…,xL,其中x1<x2…<xL,对应的采样值分别为:y1,y2,…,yL;
设三次样条插值函数的表达式为:
其中,hi=xi+1-xi;Mi=Si″(x),Si″(x)为Si(x)的二阶导数,Mi的求解方法为:
首先,利用Si(x)的二阶导数导出方程:
μiMi-1+2Mi+λiMi+1=di
其中,
然后,根据Si(x)的边界条件:S1″′(x)=S2″′(x),SL-1″′(x)=SL″′(x),其中Si″′(x)为Si(x)的三阶导数,得到:
联立上述方程,求解得出M1,M2,…,ML,将所得解带入三次样条插值函数表达式,即可求得二次极大值序列的三次样条插值函数;
2)根据二次极大值序列的三次样条插值函数计算原始FMCW雷达回波信号中除二次极大值序列中的L个采样点以外的N-L个采样点对应的函数值;
对于每一个采样点,首先判断其对应的采样时间x所在区间,然后将x代入相应区间对应的三次样条插值函数计算该采样点对应的函数值,具体地:
若x≤x2,则将其代入S1(x)进行计算;
若x∈(xi,xi+1],则将其代入Si(x)进行计算;i∈[2,L-2];
若x>xL-1,则将其代入SL-1(x)进行计算;
3)将计算得到的N-L个采样点对应的函数值与二次极大值序列中L个采样点对应的采样值,按照采样时间进行排序,得到上包络u(n),n=0,1,2,…,N-1。
5.根据权利要求1所述的抑制FMCW雷达回波信号中寄生调幅的方法,其特征在于,所述步骤3中拟合出下包络l(n)具体包括以下步骤:
ⅰ)求解二次极小值序列的三次样条插值函数;
设二次极小值序列中采样点的个数为K,K个采样点对应的采样时间分别为x1,x2,…,xK,其中x1<x2…<xK,对应的采样值分别为:y1,y2,…,yK;
设三次样条插值函数的表达式为:
其中,hi=xi+1-xi;Mi=Si″(x),Si″(x)为Si(x)的二阶导数,Mi的求解方法为:
首先,利用Si(x)的二阶导数导出方程:
μiMi-1+2Mi+λiMi+1=di
其中,
然后,根据Si(x)的边界条件:S1″′(x)=S2″′(x),SK-1″′(x)=SK″′(x),其中Si″′(x)为Si(x)的三阶导数,得到:
联立上述方程,求解得出M1,M2,…,MK,将所得解带入三次样条插值函数表达式,即可求得二次极小值序列的三次样条插值函数;
ⅱ)根据二次极小值序列的三次样条插值函数计算原始FMCW雷达回波信号中除二次极小值序列中的K个采样点以外的N-K个采样点对应的函数值;
对于每一个采样点,首先判断其对应的采样时间x所在区间,然后将x代入相应区间对应的三次样条插值函数计算该采样点对应的函数值,具体地:
若x≤x2,则将其代入S1(x)进行计算;
若x∈(xi,xi+1],则将其代入Si(x)进行计算;i∈[2,K-2];
若x>xK-1,则将其代入SK-1(x)进行计算;
ⅲ)将计算得到的N-K个采样点对应的函数值与二次极小值序列中K个采样点对应的采样值,按照采样时间进行排序,得到下包络l(n),n=0,1,2,…,N-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711428715.7A CN108196237B (zh) | 2017-12-26 | 2017-12-26 | 一种抑制fmcw雷达回波信号中寄生调幅的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711428715.7A CN108196237B (zh) | 2017-12-26 | 2017-12-26 | 一种抑制fmcw雷达回波信号中寄生调幅的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108196237A true CN108196237A (zh) | 2018-06-22 |
CN108196237B CN108196237B (zh) | 2021-06-25 |
Family
ID=62584047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711428715.7A Active CN108196237B (zh) | 2017-12-26 | 2017-12-26 | 一种抑制fmcw雷达回波信号中寄生调幅的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108196237B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110275142A (zh) * | 2019-05-16 | 2019-09-24 | 中国人民解放军63921部队 | 一种宽带雷达直采数据预处理方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101957446A (zh) * | 2010-09-26 | 2011-01-26 | 深圳市汉华安道科技有限责任公司 | 一种fmcw雷达测距的方法和装置 |
CN103226194A (zh) * | 2013-03-26 | 2013-07-31 | 中国科学院电子学研究所 | 一种基于经验模式分解的InSAR干涉相位滤波方法 |
CN104392427A (zh) * | 2014-12-09 | 2015-03-04 | 哈尔滨工业大学 | 经验模态分解和稀疏表示相结合的sar图像去噪方法 |
US20150066486A1 (en) * | 2013-08-28 | 2015-03-05 | Accusonus S.A. | Methods and systems for improved signal decomposition |
CN104618286A (zh) * | 2014-12-30 | 2015-05-13 | 苏州东奇信息科技股份有限公司 | 基于冲击滤波器成型的严格带限高效调制系统 |
CN104700151A (zh) * | 2014-05-26 | 2015-06-10 | 国网辽宁省电力有限公司 | 一种基于三次样条插值曲线拟合的风电功率考核方法 |
CN106646012A (zh) * | 2016-09-09 | 2017-05-10 | 中国舰船研究设计中心 | 一种用频装备发射谱包络提取方法 |
-
2017
- 2017-12-26 CN CN201711428715.7A patent/CN108196237B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101957446A (zh) * | 2010-09-26 | 2011-01-26 | 深圳市汉华安道科技有限责任公司 | 一种fmcw雷达测距的方法和装置 |
CN103226194A (zh) * | 2013-03-26 | 2013-07-31 | 中国科学院电子学研究所 | 一种基于经验模式分解的InSAR干涉相位滤波方法 |
US20150066486A1 (en) * | 2013-08-28 | 2015-03-05 | Accusonus S.A. | Methods and systems for improved signal decomposition |
CN104700151A (zh) * | 2014-05-26 | 2015-06-10 | 国网辽宁省电力有限公司 | 一种基于三次样条插值曲线拟合的风电功率考核方法 |
CN104392427A (zh) * | 2014-12-09 | 2015-03-04 | 哈尔滨工业大学 | 经验模态分解和稀疏表示相结合的sar图像去噪方法 |
CN104618286A (zh) * | 2014-12-30 | 2015-05-13 | 苏州东奇信息科技股份有限公司 | 基于冲击滤波器成型的严格带限高效调制系统 |
CN106646012A (zh) * | 2016-09-09 | 2017-05-10 | 中国舰船研究设计中心 | 一种用频装备发射谱包络提取方法 |
Non-Patent Citations (4)
Title |
---|
HUI LI等: ""Ensemble Empirical Mode Decomposition and Hilbert-Huang Transform Applied to Bearing Fault Diagnosis"", 《2010 3RD INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING (CISP2010)》 * |
S.A. ELGAMEL等: ""Empirical mode decomposition-based monopulse processor for enhanced radar tracking in the presence of high-power interference"", 《IET RADAR, SONAR AND NAVIGATION》 * |
盖强: "《局部波分析理论及工程应用》", 31 October 2015, 国防工业出版社 * |
盛怀茂: ""汽车防撞毫米波雷达FMCW雷达前端集成关键技术研究"", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110275142A (zh) * | 2019-05-16 | 2019-09-24 | 中国人民解放军63921部队 | 一种宽带雷达直采数据预处理方法 |
CN110275142B (zh) * | 2019-05-16 | 2021-04-13 | 中国人民解放军63921部队 | 一种宽带雷达直采数据预处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108196237B (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111352102B (zh) | 一种基于调频连续波雷达的多目标个数检测方法及装置 | |
JP6564472B2 (ja) | 車両レーダーシステム | |
CN107861117B (zh) | 一种适用于连续波周界监视雷达的多目标参数测量方法 | |
CN109407071B (zh) | 雷达测距方法、雷达测距装置、无人机和存储介质 | |
EP3324205B1 (en) | Decentralised radar system | |
JP2020067455A (ja) | 妨害信号抑圧を行うfmcwレーダー | |
CN103760558B (zh) | 一种太赫兹雷达isar成像方法 | |
CN109425854B (zh) | 多分辨率多普勒处理 | |
EP3771921B1 (en) | Method for processing non-interrupted phase synchronization signal of bistatic sar based on coded signal | |
CN103412309B (zh) | 移不变双基地前视合成孔径雷达nlcs成像方法 | |
CN102169176A (zh) | 观测信号处理装置 | |
CN110770598A (zh) | 微波雷达及其数据处理方法、无人机 | |
CN103983974A (zh) | 双站调频连续波合成孔径雷达成像方法 | |
CN110857975B (zh) | 雷达范围精确度改进方法 | |
CN106054154B (zh) | 一种基于icpf的机动目标的步进频率信号高分辨率成像方法 | |
Hyun et al. | Development of short-range ground surveillance radar for moving target detection | |
CN111060909B (zh) | 一种机载雷达斜前视超分辨成像方法 | |
Gardill et al. | In-situ time-frequency analysis of the 77 GHz bands using a commercial chirp-sequence automotive FMCW radar sensor | |
CN108196237B (zh) | 一种抑制fmcw雷达回波信号中寄生调幅的方法 | |
US20220268880A1 (en) | Radar device | |
RU2296346C2 (ru) | Способ измерения дальности в импульсно-доплеровских радиолокационных станциях | |
Tran et al. | An experimental study of radar tomographic imaging in a multi-bistatic scenario | |
CN109061596B (zh) | 一种基于子空间分析的调频间断连续波雷达杂散频谱消除方法 | |
Wang et al. | Velocity estimation of moving targets using SAR | |
CN112965067A (zh) | 一种适用于fmcw汽车雷达目标速度扩展方法 |
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 |