CN109614579A - 相控阵天气雷达最大值提取信号方法 - Google Patents

相控阵天气雷达最大值提取信号方法 Download PDF

Info

Publication number
CN109614579A
CN109614579A CN201811581355.9A CN201811581355A CN109614579A CN 109614579 A CN109614579 A CN 109614579A CN 201811581355 A CN201811581355 A CN 201811581355A CN 109614579 A CN109614579 A CN 109614579A
Authority
CN
China
Prior art keywords
signal
pmax
phased array
maximum value
calculate
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
Application number
CN201811581355.9A
Other languages
English (en)
Other versions
CN109614579B (zh
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.)
Lei Xiang Technology (beijing) Co Ltd
Original Assignee
Lei Xiang Technology (beijing) Co Ltd
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 Lei Xiang Technology (beijing) Co Ltd filed Critical Lei Xiang Technology (beijing) Co Ltd
Priority to CN201811581355.9A priority Critical patent/CN109614579B/zh
Publication of CN109614579A publication Critical patent/CN109614579A/zh
Application granted granted Critical
Publication of CN109614579B publication Critical patent/CN109614579B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S2013/0236—Special technical features
    • G01S2013/0245—Radar with phased array antenna
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Electromagnetism (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Discrete Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种相控阵天气雷达最大值提取信号方法，其步骤如下，1、快速傅氏变换、2、提取Pmax（x,y,z）、3、计算V（x,y,z）、4、计算DV（x,y,z）、4、判断是否是信号。本发明采用提取谱中最大值，同时利用空间径向速度分布均匀性甄别是否是信号，在信号的谱形不确定的情况下，在提取小信号时，也不会容易出现错误。

Description

相控阵天气雷达最大值提取信号方法
技术领域
本发明涉及一种气象雷达的信号处理技术。
背景技术
在对相控雷达探测信号完成快速傅氏变换（FFT）后，就得到了回波信号和接收通道噪声的速度功率密度谱，接着的工作就是从这个谱中提取信号功率。通常采用的方式是根据信号的谱形特征在这个谱上找到信号的位置，测量信号谱宽，根据谱宽和信号功率谱密度计算信号功率。信号的谱形常常是不确定的，因此这些方法在提取小信号时，容易出现错误。
发明内容
本发明的目的是提供一种相控阵天气雷达最大值提取信号方法，采用提取谱中最大值，同时利用空间径向速度分布均匀性甄别是否是信号，以解决在信号的谱形常常不确定的情况下，提取小信号时，不容易出现错误等技术问题。
为了实现上述发明目的，本发明所采用的技术方案如下：
相控阵天气雷达最大值提取信号方法，
1）快速傅氏变换
在对相控雷达探测信号完成快速傅氏变换FFT后，得到
P（I），I=1,2，……128；取其中最大值Pmax，同时利用Pmax在谱中的位置得到多普勒速度V；其中P（I）是……
2）提取Pmax（x,y,z）
对雷达探测空间的信号完成上述处理后，就可得到Pmax（x,y,z）和V（x,y,z）；
3）计算V（x,y,z）、DV（x,y,z）
对于一个V（x,y,z），可以找到一组V（x+i,y+j,z+k）,i=1,2,……N, j=1,2,……N, k=1,2,……N,N不大于10；计算V（x+i,y+j,z+k）中相邻点的差值，并将这些差值的绝对值求平均，得DV（x,y,z）；其中DV是V（x+i,y+j,z+k）中相邻点的差值绝对值的平均值；
4）判断是否是信号
当DV（x,y,z）<B，Pmax（x,y,z）为信号，否则为噪声；其中B是判断是否为信号的阈值。
本发明的优点是采用提取谱中最大值，同时利用空间径向速度分布均匀性甄别是否是信号，在信号的谱形不确定的情况下，在提取小信号时，也不会容易出现错误。
具体实施方式
本发明采用提取谱中最大值，同时利用空间径向速度分布均匀性甄别是否是信号。
本发明的流程：
一、流程
1、快速傅氏变换
2、提取Pmax（x,y,z）
3、计算V（x,y,z）
4、计算DV（x,y,z）
5、判断是否是信号。
本发明上述流程包含的具体方法：
在对相控雷达探测信号完成快速傅氏变换（FFT）后，得到
P（I），I=1,2，……128
取其中最大值Pmax，同时利用Pmax在谱中的位置得到多普勒速度V。
对雷达探测空间的信号完成上述处理后，就可得到Pmax（x,y,z）和V（x,y,z）。对于一个Pmax（x,y,z），可以找到一组,i=1,2,……N, j=1,2,……N, k=1,2,……N,N一般不大于10。计算V（x+i,y+j,z+k）中相邻点的差值，并将这些差值的绝对值求平均，得DV（x,y,z），当DV（x,y,z）<B，如B=5，Pmax（x,y,z）为信号，否则为噪声。

Claims (1)

1.相控阵天气雷达最大值提取信号方法，其步骤如下：
1）快速傅氏变换
在对相控雷达探测信号完成快速傅氏变换FFT后，得到
P（I），I=1,2，……128；取其中最大值Pmax，同时利用Pmax在谱中的位置得到多普勒速度V；其中P（I）是功率谱密度；
2）提取Pmax（x,y,z）
对雷达探测空间的信号完成上述处理后，就可得到Pmax（x,y,z）和V（x,y,z）；
3）计算V（x,y,z）、DV（x,y,z）
对于一个V（x,y,z），可以找到一组V（x+i,y+j,z+k）,i=1,2,……N, j=1,2,……N, k=1,2,……N,N不大于10；计算V（x+i,y+j,z+k）中相邻点的差值，并将这些差值的绝对值求平均，得DV（x,y,z）；其中DV是V（x+i,y+j,z+k）中相邻点的差值绝对值的平均值；
4）判断是否是信号
当DV（x,y,z）<B，Pmax（x,y,z）为信号，否则为噪声；其中B是判断是否为信号的阈值。
CN201811581355.9A 2018-12-24 2018-12-24 相控阵天气雷达最大值提取信号方法 Active CN109614579B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811581355.9A CN109614579B (zh) 2018-12-24 2018-12-24 相控阵天气雷达最大值提取信号方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811581355.9A CN109614579B (zh) 2018-12-24 2018-12-24 相控阵天气雷达最大值提取信号方法

Publications (2)

Publication Number Publication Date
CN109614579A true CN109614579A (zh) 2019-04-12
CN109614579B CN109614579B (zh) 2022-05-17

Family

ID=66011404

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811581355.9A Active CN109614579B (zh) 2018-12-24 2018-12-24 相控阵天气雷达最大值提取信号方法

Country Status (1)

Country Link
CN (1) CN109614579B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080001808A1 (en) * 2004-12-30 2008-01-03 Passarelli Richard E Jr System and method for processing data in weather radar
CN102540156A (zh) * 2010-12-14 2012-07-04 中国航天科工集团第二研究院二十三所 一种基于双通道的风廓线雷达信号和数据处理方法
CN102866391A (zh) * 2012-09-05 2013-01-09 中北大学 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法
CN103983975A (zh) * 2014-06-05 2014-08-13 中国气象局气象探测中心 基于两种雷达的大气运动垂直速度检测方法及系统
CN104122538A (zh) * 2013-04-24 2014-10-29 何平 确定风廓线雷达噪声功率的方法
CN104316929A (zh) * 2014-09-12 2015-01-28 何平 一种功率谱识别方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080001808A1 (en) * 2004-12-30 2008-01-03 Passarelli Richard E Jr System and method for processing data in weather radar
CN102540156A (zh) * 2010-12-14 2012-07-04 中国航天科工集团第二研究院二十三所 一种基于双通道的风廓线雷达信号和数据处理方法
CN102866391A (zh) * 2012-09-05 2013-01-09 中北大学 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法
CN104122538A (zh) * 2013-04-24 2014-10-29 何平 确定风廓线雷达噪声功率的方法
CN103983975A (zh) * 2014-06-05 2014-08-13 中国气象局气象探测中心 基于两种雷达的大气运动垂直速度检测方法及系统
CN104316929A (zh) * 2014-09-12 2015-01-28 何平 一种功率谱识别方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
龚树凤: "不确定条件下的机会阵雷达信号处理技术研究", 《中国博士学位论文全文数据库》 *

Also Published As

Publication number Publication date
CN109614579B (zh) 2022-05-17

Similar Documents

Publication Publication Date Title
CN108761418B (zh) 一种基于微多普勒特征的雷达多径目标剔除方法
CN111142086B (zh) 一种pd雷达幅度抖动抑制方法、时间抖动检测方法及系统
JP6260482B2 (ja) 物標検出装置
CN107462873A (zh) 一种雷达干扰快速识别方法
CN103760530A (zh) 基于信号锥的雷达压制式干扰和欺骗式干扰识别方法
KR20160054349A (ko) 클러터를 제거하기 위한 레이더 신호 처리 장치 및 방법
CN107843876A (zh) 一种雷达脉冲重频的分选方法及设备
CN109459745B (zh) 一种利用辐射噪声估计运动声源速度的方法
CN111856420A (zh) 一种多普勒雷达多目标检测方法
CN110133632B (zh) 一种基于cwd时频分析的复合调制信号识别方法
CN102636775A (zh) 基于模糊逻辑识别的风廓线雷达回波谱重构方法
CN106093896A (zh) 一种快速检测和估计脉冲多普勒雷达脉冲抖动的方法
CN109782249B (zh) 一种两目标相关时延估计算法
CN109614579A (zh) 相控阵天气雷达最大值提取信号方法
JP6523790B2 (ja) 目標検出装置
CN105553619B (zh) 基于调制识别辅助的突发信号检测方法
CN109085568B (zh) 一种基于二次混频的调频连续波多目标检测方法
CN107589419A (zh) 一种能量峰值联合宽带高分辨距离像目标提取的方法
CN109541556B (zh) 一种对线性调频信号移频干扰的识别方法
CN106597390B (zh) 一种海杂波抑制方法
CN105891784B (zh) 雷达的风电场杂波抑制中风轮机朝向估计方法
CN107329123A (zh) 一种弱雷达脉冲包络信号检测方法和装置
CN114545384A (zh) 一种应用于跑道异物检测的整合方法及系统
CN103412293A (zh) 一种基于脉冲分选技术的单脉冲雷达水平方向图测试方法
Hyun et al. Two-step pairing algorithm for target range and velocity detection in FMCW automotive radar

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