CN105527624B - 一种雷达回波动态估算噪声的方法和气象雷达*** - Google Patents

一种雷达回波动态估算噪声的方法和气象雷达*** Download PDF

Info

Publication number
CN105527624B
CN105527624B CN201410515512.1A CN201410515512A CN105527624B CN 105527624 B CN105527624 B CN 105527624B CN 201410515512 A CN201410515512 A CN 201410515512A CN 105527624 B CN105527624 B CN 105527624B
Authority
CN
China
Prior art keywords
power spectrum
clock signal
range bin
noise
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.)
Active
Application number
CN201410515512.1A
Other languages
English (en)
Other versions
CN105527624A (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.)
BEIJING METSTAR RADAR Co Ltd
Original Assignee
BEIJING METSTAR RADAR 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 BEIJING METSTAR RADAR Co Ltd filed Critical BEIJING METSTAR RADAR Co Ltd
Priority to CN201410515512.1A priority Critical patent/CN105527624B/zh
Publication of CN105527624A publication Critical patent/CN105527624A/zh
Application granted granted Critical
Publication of CN105527624B publication Critical patent/CN105527624B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种雷达回波动态估算噪声的方法和气象雷达***，以解决现有的气象雷达***不能实时采集动态变化的***噪声，因而气象预测的准确度不够高的技术问题。上述方法可包括：将雷达回波信号转换为IQ时序信号，对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱spM，M为正整数；按升序对功率谱spM进行排序，记为spS；设变量mc和b0，初始化时令mc＝N，b0＝0；判断mc>2，且b0<1时，执行以下循环计算：令mc＝mc‑1，计算b1和b2及新的b0，其中： b0＝b2/b1；每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为

Description

一种雷达回波动态估算噪声的方法和气象雷达***
技术领域
本发明涉及雷达技术领域，具体涉及一种雷达回波动态估算噪声的方法和气象雷达***。
背景技术
在气象雷达信号处理中，需要将气象雷达***的***噪声从雷达回波信号中去除，以保障气象预测的准确度。
气象雷达***一般有两种扫描模式，即平面位置扫描和体扫。平面位置扫描是指：雷达固定仰角，进行一次360度扫描。体扫是指：一组从低到高的不同仰角的平面位置扫描。
目前，通常采用以下方法采集气象雷达***的***噪声，采集方法包括：在***相邻体扫中间，于天线仰角比较高的位置(如19.5度)，控制***的发射机不发射信号，此时采集***内的信号，作为***噪声，接下来的体扫就以此噪声为***噪声。
气象雷达信号处理中，将雷达回波信号沿射线方向按距离分成的小的距离单元，被称为距离库(range bin)。实际的气象雷达***中，不同仰角的***噪声是不同的，不同距离库的***噪声也有可能是不同的，因此，气象雷达***的***噪声是动态变化的。
但是，上述的现有***噪声采集方法只采集一个固定的***噪声，在信号处理中只使用这一个固定的***噪声，不能反映实际上动态变化的***噪声，因此，会在后续的信号处理中来带偏差，导致信号处理结果出现误差，不能进行准确的气象预测。
发明内容
有鉴于此，本发明实施例提供一种雷达回波动态估算噪声的方法和气象雷达***，以解决现有的气象雷达***不能实时采集动态变化的***噪声，因而气象预测的准确度不够高的技术问题。
本发明第一方面提供一种雷达回波动态估算噪声的方法，包括：
将雷达回波信号转换为IQ时序信号，对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱spM ，M为正整数；
按升序对功率谱spM 进行排序，记排序后的功率谱为spS；
设变量mc和b0，初始化时令mc＝N，b0＝0；
判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：b0＝b2/b1；
每一次循环计算后，判断是否有b1≤0，并记录当前的mc；
如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为
本发明第二方面提供一种气象雷达***，包括：
信号转换模块，用于将雷达回波信号转换为IQ时序信号；
功率谱模块，用于对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱spM ，M为正整数；按升序对功率谱spM 进行排序，记排序后的功率谱为spS；
循环计算模块，用于设变量mc和b0，初始化时令mc＝N，b0＝0；判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：
b0＝b2/b1；
噪声计算模块，用于每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为
本发明实施例通过采用上述技术方案，可以实时精确的估算雷达回波信号每个距离库的***噪声，能够精确计算每个距离库的噪声，得到合适的信噪比，避免使用固定的***噪声带来的计算偏差，可以提高气象预测的准确度。
附图说明
图1是本发明实施例提供的雷达回波动态估算噪声的方法的流程图；
图2是本发明实施例提供的气象雷达***的示意图。
具体实施方式
本发明实施例提供一种雷达回波动态估算噪声的方法和气象雷达***，以解决现有的气象雷达***不能实时采集动态变化的***噪声，因而气象预测的准确度不够高的技术问题。下面结合附图分别进行详细说明。
实施例一、
请参考图1，本发明实施例提供一种地物杂波残留滤除方法，包括：
101、计算功率谱：将雷达回波信号转换为IQ时序信号，对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换(Fast Fourier Transform，FFT)，并计算功率谱系数，得到功率谱spM，M为正整数。
进行气象预测时，气象雷达***采用平面位置扫描和/或体扫模式进行扫描，得到雷达回波信号。后续，根据对雷达回波信号的分析计算，预测气象。实际的信号处理中，也将雷达回波信号转换为IO时序信号，IQ时序信号是指将雷达回波信号相移90度得到的信号，后续基于IO时序信号进行处理。
雷达回波信号处理中，将雷达回波信号沿射线方向按距离分成的小的距离单元，被称为距离库(range bin)，后续将对每个距离库对应的IQ时序信号分别进行处理。本发明实施例中以第M个距离库为例，并假定第M个距离库的IQ 时序信号的采样数为N，M和N均为正整数。
本实施例中记第M个距离库的IQ时序信号为x(n)，n为不大于N的正整数，表示第n个采样点。通过对x(n)进行TFT处理，可得到信号其中，k＝0,1,…,N-1，然后根据公式spM(k) ＝ |yM(k)|2 计算功率谱系数，即可得到功率谱spM。
102、功率谱排序：按升序对功率谱spM进行排序，记排序后的功率谱为spS。
本发明实施例中，基于功率谱spM估算***噪声，需要将功率谱spM按照升序进行排序，可记排序后的功率谱为spS。功率谱为spS中包括N个功率谱系数，分别对应于N个采样点。可定义功率谱spS中，功率最低的功率谱系数为第1个，功率最高的功率谱系数为第N个。
103、循环计算：设变量mc和b0，初始化时令mc＝N，b0＝0；判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：
b0＝b2/b1。
本发明实施例中，采用一定的规则从功率谱中选择出功率最低的若干个功率谱系数，来计算***噪声，方法如下：
设变量mc和b0，初始化时令mc＝N，b0＝0；
判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：b0＝b1/b2；
这里所说的循环计算是指，每一次令mc＝mc-1，并按照上述公式计算出b1 和b2及新的b0后，判断是否有mc&gt;2，且b0&lt;1，若是，则再次令mc＝mc-1，并按照上述公式计算出b1和b2及新的b0，再次判断是否有mc&gt;2，且b0&lt;1，若是，按照上述方式继续计算，直到不满足mc&gt;2，或b0&lt;1。
该循环计算过程中，mc不断的变化，具体为从N一次减小，每次减小1。
104、***噪声计算：每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为
本发明实施例中，在每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算***噪声，即，根据功率最低的第1 个至第mc个功率谱系数计算***噪声，计算公式为
至此，即计算得到第M个距离库对应的***噪声。采用上述方法依次对每个距离库进行计算，可得到每个距离库的***噪声。
由上可见，本发明实施例提供了一种雷达回波动态估算噪声的方法，该方法可以实时精确的估算雷达回波信号每个距离库的***噪声，能够精确计算每个距离库的噪声，得到合适的信噪比，避免使用固定的***噪声带来的计算偏差，可以提高气象预测的准确度。
实施例二、
请参考图2，本发明实施例提供一种气象雷达***，包括：
信号转换模块201，用于将雷达回波信号转换为IQ时序信号；
功率谱模块202，用于对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换，并计算功率谱系数，得到功率谱spM，M为正整数；按升序对功率谱spM进行排序，记排序后的功率谱为spS；
循环计算模块203，用于设变量mc和b0，初始化时令mc＝N，b0＝0；判断 mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：
b0＝b2/b1；
噪声计算模块204，用于每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为
可选的，所述信号转换模块201，具体用于将雷达回波信号相移90度得到 IQ时序信号。
可选的，所述功率谱模块202，具体用于对第M个距离库的采样数为N的IQ 时序信号进行傅里叶变换，得到信号其中， k＝0,1,…,N-1；根据公式spM(k)＝|yM(k)|2计算功率谱系数，得到功率谱spM。
可选的，所述循环计算模块203，具体用于：判断mc&gt;2，且b0&lt;1时，令 mc＝mc-1，计算b1和b2及新的b0，其中： b0＝b2/b1；重复上述步骤，直到mc不大于2或者b0不小于1。
由上可见，本发明实施例提供了一种气象雷达***，关于该***的更详细的说明，请参考实施例一中的记载，该***可以实时精确的估算雷达回波信号每个距离库的***噪声，能够精确计算每个距离库的噪声，得到合适的信噪比，避免使用固定的***噪声带来的计算偏差，可以提高气象预测的准确度。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤可以通过硬件来完成，也可以通过程序指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器、随机读取存储器、磁盘或光盘等。
以上对本发明实施例所提供的一种雷达回波动态估算噪声的方法和气象雷达***进行了详细介绍，但以上实施例的说明只是用于帮助理解本发明的方法及其核心思想，不应理解为对本发明的限制。本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (8)

1.一种雷达回波动态估算噪声的方法，其特征在于，包括：
将雷达回波信号转换为IQ时序信号，对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱spM，M为正整数；
按升序对功率谱spM进行排序，记排序后的功率谱为spS；
设变量mc和b0，初始化时令mc＝N，b0＝0；
判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：
b0＝b2/b1，spS(n)表示spS中的第n个功率谱系数，n为不大于N的正整数；
每一次循环计算后，判断是否有b1≤0，并记录当前的mc；
如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为
2.根据权利要求1所述的方法，其特征在于，所述将雷达回波信号转换为IQ时序信号包括：
将雷达回波信号相移90度得到IQ时序信号。
3.根据权利要求1所述的方法，其特征在于，所述对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱spM包括：
对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换，得到信号
其中，k＝0,1,…,N-1；
根据公式spM(k)＝|yM(k)|2计算功率谱系数，得到功率谱spM。
4.根据权利要求1、2或3所述的方法，其特征在于，所述判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及b0，其中：
b0＝b2/b1包括：
判断mc&gt;2，且b0&lt;1时，令mc＝mc-1，计算b1和b2及新的b0，其中：
b0＝b2/b1；
重复上述步骤，直到mc不大于2或者b0不小于1。
5.一种气象雷达***，其特征在于，包括：
信号转换模块，用于将雷达回波信号转换为IQ时序信号；
功率谱模块，用于对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换并计算功率谱系数，得到功率谱spM ，M为正整数；按升序对功率谱spM 进行排序，记排序后的功率谱为spS；
循环计算模块，用于设变量mc和b0，初始化时令mc＝N，b0＝0；判断mc&gt;2，且b0&lt;1时，执行以下循环计算：令mc＝mc-1，计算b1和b2及新的b0，其中：
b0＝b2/b1，spS(n)表示spS中的第n个功率谱系数，n为不大于N的正整数；
噪声计算模块，用于每一次循环计算后，判断是否有b1≤0，并记录当前的mc；如果b1≤0，则采用当前的mc估算***噪声，得到第M个距离库的***噪声为
6.根据权利要求5所述的***，其特征在于，
所述信号转换模块，具体用于将雷达回波信号相移90度得到IQ时序信号。
7.根据权利要求5所述的***，其特征在于，
所述功率谱模块，具体用于对第M个距离库的采样数为N的IQ时序信号进行傅里叶变换，得到信号其中，k＝0,1,…,N-1；根据公式spM(k)＝|yM(k)| 2计算功率谱系数，得到功率谱spM 。
8.根据权利要求5、6或7所述的***，其特征在于，
所述循环计算模块，具体用于判断mc&gt;2，且b0&lt;1时，令mc＝mc-1，计算b1和b2及新的b0，其中：b0＝b2/b1；重复上述步骤，直到mc不大于2或者b0不小于1。
CN201410515512.1A 2014-09-29 2014-09-29 一种雷达回波动态估算噪声的方法和气象雷达*** Active CN105527624B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410515512.1A CN105527624B (zh) 2014-09-29 2014-09-29 一种雷达回波动态估算噪声的方法和气象雷达***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410515512.1A CN105527624B (zh) 2014-09-29 2014-09-29 一种雷达回波动态估算噪声的方法和气象雷达***

Publications (2)

Publication Number Publication Date
CN105527624A CN105527624A (zh) 2016-04-27
CN105527624B true CN105527624B (zh) 2019-02-01

Family

ID=55769961

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410515512.1A Active CN105527624B (zh) 2014-09-29 2014-09-29 一种雷达回波动态估算噪声的方法和气象雷达***

Country Status (1)

Country Link
CN (1) CN105527624B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106249207B (zh) * 2016-07-11 2018-10-12 北京敏视达雷达有限公司 一种估算噪声功率的方法及装置
CN108828542B (zh) * 2018-08-03 2022-07-29 中国航空工业集团公司雷华电子技术研究所 一种机载气象雷达噪声周期动态校准方法
CN109358331B (zh) * 2018-10-15 2023-03-03 成都信息工程大学 气象雷达实时动态噪声功率检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101881824A (zh) * 2009-05-05 2010-11-10 何平 功率谱密度数据噪声门限的客观快速确定方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5508706A (en) * 1991-09-30 1996-04-16 Trw Inc. Radar signal processor
CN101325007B (zh) * 2008-07-24 2011-04-06 中国科学院上海微***与信息技术研究所 基于毫米波交通雷达的公路背景识别与公路背景抑制方法
CN102323574B (zh) * 2011-06-25 2013-03-06 中国航天科工集团第二研究院二十三所 一种风廓线雷达信号处理方法
CN203133274U (zh) * 2013-03-27 2013-08-14 武汉大学 一种基于全数字高频雷达的环境噪声监测***

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101881824A (zh) * 2009-05-05 2010-11-10 何平 功率谱密度数据噪声门限的客观快速确定方法

Also Published As

Publication number Publication date
CN105527624A (zh) 2016-04-27

Similar Documents

Publication Publication Date Title
US9097805B2 (en) Systems and methods for calibrating dual polarization radar systems
CN103630888B (zh) 基于对称三角lfmcw雷达的高精度实时微波测速测距装置
CN103926474B (zh) 相控阵天线单元特性近场测量方法
CN110646773B (zh) 一种基于汽车毫米波雷达的弱目标检测方法、跟踪方法及***
CN109471095B (zh) 一种基于快速迭代插值的fmcw雷达距离估计方法
CN103364770B (zh) 基于矩阵填充的雷达目标检测***及其检测方法
Gong et al. A three-step dealiasing method for Doppler velocity data quality control
CN103837866B (zh) 一种岸对海雷达自动录取能力的测评方法及***
CN105527624B (zh) 一种雷达回波动态估算噪声的方法和气象雷达***
CN103616687A (zh) 分段线性估计的多项式拟合isar包络对齐方法
CN107247193B (zh) 天线近场测试方法及装置
CN109613503A (zh) 雷达回波信号的标校方法和装置
CN108414992B (zh) 一种基于相位信息杂波图的目标检测方法
CN111812636B (zh) 一种基于权重融合选择的粒子滤波检测前跟踪方法
Cosoli et al. A real-time and offline quality control methodology for SeaSonde high-frequency radar currents
CN111352083A (zh) 一种高频地波雷达多接收通道增益自动校准方法及装置
CN102540153A (zh) 一种基于站间直达波干扰的阵列幅相误差校准方法
CN102565755A (zh) 利用所测频谱数据进行宽带无线电测向的方法
CN105652254A (zh) 室外场rcs测量方法及***
CN107171748B (zh) 欠采样的多阵列协同测频与直接定位方法
CN110426676B (zh) 一种基于超声波的定位***及定位方法
CN107329128B (zh) 一种基于快速存储技术的超分辨高精度低空测角方法
CN117075035B (zh) 自旋短基线高精度测向方法、***、设备、介质及终端
CN112285686B (zh) 多基地雷达异步数据的实时在线融合方法
CN117647811B (zh) 运动目标参数估计与重定位方法、***、终端及存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant