CN102323574B - 一种风廓线雷达信号处理方法 - Google Patents

一种风廓线雷达信号处理方法 Download PDF

Info

Publication number
CN102323574B
CN102323574B CN2011101732501A CN201110173250A CN102323574B CN 102323574 B CN102323574 B CN 102323574B CN 2011101732501 A CN2011101732501 A CN 2011101732501A CN 201110173250 A CN201110173250 A CN 201110173250A CN 102323574 B CN102323574 B CN 102323574B
Authority
CN
China
Prior art keywords
signal
signal processing
control module
system control
spectrum
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
CN2011101732501A
Other languages
English (en)
Other versions
CN102323574A (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.)
Aerospace new weather Technology Co., Ltd
Original Assignee
No23 Institute Of No20 Academy Casic
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 No23 Institute Of No20 Academy Casic filed Critical No23 Institute Of No20 Academy Casic
Priority to CN2011101732501A priority Critical patent/CN102323574B/zh
Publication of CN102323574A publication Critical patent/CN102323574A/zh
Application granted granted Critical
Publication of CN102323574B publication Critical patent/CN102323574B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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

本发明公开了一种风廓线雷达信号处理方法,其具体步骤为:第一步,搭建风廓线雷达信号处理***,包括:工业控制计算机(1)、数字中频接收机板(2)、***控制模块(3)和信号处理模块(4);第二步,数字中频接收机板(2)采集雷达回波信号;第三步,信号处理模块(4)进行信号处理;第四步,***控制模块(3)进行谱数据的处理、显示和传送;经过上述处理,得到所处理高度层上的风场回波信号功率和信噪比估计值;至此,实现了风廓线雷达信号处理。本方法数据处理精度高、存储容量大,减小了硬件复杂度和成本,提高了***性能和可靠性,便于信号处理算法的扩展升级,通用性强,产品性价比高,保证了产品批量化的需要。

Description

一种风廓线雷达信号处理方法
技术领域
本发明涉及一种雷达信号处理方法,特别是一种风廓线雷达信号处理方法。
背景技术
风廓线雷达(Wind Profiling Radar,WPR)是通过探测大气湍流对电磁波的反射和散射来连续探测大气风场等物理量的一种全天候大气遥感探测设备,具有较高的测风实时性和精度,是中小尺度大气探测的关键设备之一,在我国得到了快速发展和应用。
原有的风廓线雷达信号处理方法是使用专用的进口信号处理器,进口信号处理器包括:AD转换部件、DSP信号处理部件,硬件复杂成本较高。由于硬件资源有限,在算法扩展、计算精度、数据存储等方面都受到限制。在产品设计时需要专用的开发工具,程序固化在器件中,数据处理位数只有23位,不能录取原始数据和中间处理数据,在风廓线雷达中的应用有其局限性。
信号处理器是风廓线雷达***的关键设备之一,其性能好坏直接影响到测风的准确性和精度。原有信号处理器的硬件电路会带来额外的干扰信号,增加了故障几率,影响到***的性能和可靠性,并且批量化受到国外生产厂家的限制,***升级和扩展困难,这些问题制约了风廓线雷达产品的大批量生产和推广使用。
发明内容
本发明的目的在于提供一种风廓线雷达信号处理方法,解决原有风廓线雷达信号处理方法依赖进口器件、成本高、故障几率大以及升级困难的问题。
一种风廓线雷达信号处理方法,其具体步骤为:
第一步 搭建风廓线雷达信号处理***
风廓线雷达信号处理***包括:工业控制计算机、数字中频接收机板、***控制模块和信号处理模块,其中工业控制计算机包括:PCI总线、共享内存和以太网卡。
数字中频接收机板与工业控制计算机通过PCI总线连接,数字中频接收机板与风廓线雷达***中的模拟接收机连接;工业控制计算机与数据终端机通过以太网卡连接;***控制模块与工业控制计算机通过PCI总线连接,信号处理模块与工业控制计算机通过PCI总线连接。工业控制计算机为***的运行平台,数字中频接收机板用于采集雷达回波信号,***控制模块用于控制和调度***的工作流程,信号处理模块用于对回波采样信号进行信号处理,***控制模块和信号处理模块保存在工业控制计算机上,***运行时加载到共享内存中。
第二步 数字中频接收机板采集雷达回波信号
***控制模块首先进行数字中频接收机板的初始化和参数加载,然后启动工作时序,数字中频接收机板按照加载参数中设置的采样率、抽取因子、滤波系数,直接对雷达回波的中频信号进行A/D变换和采集,将采集的时域数据先存入数字中频接收机板内的FIFO缓存中,通过DMA中断方式传输到工业控制计算机内存中,由***控制模块接收,再保存到工业控制计算机的共享内存中。***控制模块边接收时域数据边进行计数,当计数值达到预设的数值N时,N=FFT点数×时域积累数×谱平均数,停止信号采集的工作时序,***控制模块完成信号接收,然后对预设的时域数据保存标志作判断,如果标志为1,则保存到工业控制计算机的硬盘上。
第三步 信号处理模块进行信号处理
***控制模块与信号处理模块使用共享内存交换接收到的时域数据,***控制模块调用信号处理模块,向信号处理模块传送信号处理参数和共享内存地址,由信号处理模块利用信号处理参数对共享内存中的时域数据进行时域积累、脉冲压缩、去直流、加窗、FFT变换、谱平均的一系列信号处理过程,处理完毕后将生成的谱数据保存到共享内存中,向***控制模块返回处理完毕的标志信息,***控制模块接收处理后生成的谱数据,并对预设的谱数据保存标志作判断,如果标志为1则保存到工业控制计算机的硬盘上。
第四步 ***控制模块进行谱数据的处理、显示和传送
***控制模块对信号处理模块生成的谱数据进行谱矩、信号功率和信噪比参数的计算,在本地显示器显示谱数据图形和处理的信息,用于对风场回波信号进行监测和评估,同时通过以太网卡将谱数据传送给数据终端机,作进一步的风谱信号识别、提取处理。
谱数据的处理过程如下:
信号功率谱的零阶矩为:
m 0 = Σ i = 1 k p i - - - ( 1 )
式中:m0为信号功率谱的零阶矩,pi为第i点功率值,求和范围为信号谱带宽。
信号功率谱的一阶矩为:
m 1 = Σ i = 1 k v i p i - - - ( 2 )
信号功率谱的二阶矩为:
m 2 = Σ i = 1 k v i 2 p i - - - ( 3 )
式中:m1为信号功率谱的一阶矩,m2为信号功率谱的二阶矩,vi为第i点速度值,求和范围为信号谱带宽。
信号功率估计值为:
Pr=m0          (4)
速度谱宽估计值为:
w ~ = 2 * [ m 2 m 0 - ( m 1 m 0 ) 2 ] 1 / 2 - - - ( 5 )
信噪比为:
SNR = 10 log ( Σ i = 1 k p i / P N ) - - - ( 6 )
式中:PN为噪声功率,求和范围为信号谱所占的频率范围,单位:dB。
经过上述处理,得到所处理高度层上的风场回波信号功率和信噪比估计值。
至此,实现了风廓线雷达信号处理。
本发明提供的风廓线雷达信号处理方法采用工业控制计算机、数字中频接收机板和处理模块实现,数据位数达到64位,数据处理精度高、存储容量大,减小了硬件复杂度和成本,提高了***性能和可靠性,便于信号处理算法的扩展升级,通用性强,产品性价比高,保证了产品批量化的需要。
附图说明
图1一种风廓线雷达信号处理方法的风廓线雷达信号处理***示意图。
1.工业控制计算机 2.数字中频接收机板 3.***控制模块 4.信号处理模块5.PCI总线 6.共享内存 7.以太网卡 8.模拟接收机 9.数据终端机
具体实施方式
一种风廓线雷达信号处理方法,其具体步骤为:
第一步 搭建风廓线雷达信号处理***
风廓线雷达信号处理***包括:工业控制计算机1、数字中频接收机板2、***控制模块3和信号处理模块4,其中工业控制计算机1包括:PCI总线5、共享内存6和以太网卡7。
数字中频接收机板2与工业控制计算机1通过PCI总线5连接,数字中频接收机板2与风廓线雷达***中的模拟接收机8连接;工业控制计算机1与数据终端机9通过以太网卡7连接;***控制模块3与工业控制计算机1通过PCI总线5连接,信号处理模块4与工业控制计算机1通过PCI总线5连接。工业控制计算机1为***的运行平台,数字中频接收机板2用于采集雷达回波信号,***控制模块3用于控制和调度***的工作流程,信号处理模块4用于对回波采样信号进行信号处理,***控制模块3和信号处理模块4保存在工业控制计算机1上,***运行时加载到共享内存6中。
第二步 数字中频接收机板2采集雷达回波信号
***控制模块3首先进行数字中频接收机板2的初始化和参数加载,然后启动工作时序,数字中频接收机板2按照加载参数中设置的采样率、抽取因子、滤波系数,直接对雷达回波的中频信号进行A/D变换和采集,将采集的时域数据先存入数字中频接收机板2内的FIFO缓存中,通过DMA中断方式传输到工业控制计算机1内存中,由***控制模块3接收,再保存到工业控制计算机1的共享内存6中。***控制模块3边接收时域数据边进行计数,当计数值达到预设的数值N时,N=FFT点数×时域积累数×谱平均数,停止信号采集的工作时序,***控制模块3完成信号接收,然后对预设的时域数据保存标志作判断,如果标志为1,则保存到工业控制计算机1的硬盘上。
第三步 信号处理模块4进行信号处理
***控制模块3与信号处理模块4使用共享内存6交换接收到的时域数据,***控制模块3调用信号处理模块4,向信号处理模块4传送信号处理参数和共享内存地址,由信号处理模块4利用信号处理参数对共享内存中的时域数据进行时域积累、脉冲压缩、去直流、加窗、FFT变换、谱平均的一系列信号处理过程,处理完毕后将生成的谱数据保存到共享内存中,向***控制模块返回处理完毕的标志信息,***控制模块3接收处理后生成的谱数据,并对预设的谱数据保存标志作判断,如果标志为1则保存到工业控制计算机1的硬盘上。
第四步 ***控制模块3进行谱数据的处理、显示和传送
***控制模块3对信号处理模块4生成的谱数据进行谱矩、信号功率和信噪比参数的计算,在本地显示器显示谱数据图形和处理的信息,用于对风场回波信号进行监测和评估,同时通过以太网卡将谱数据传送给数据终端机9,作进一步的风谱信号识别、提取处理。
谱数据的处理过程如下:
信号功率谱的零阶矩为:
m 0 = Σ i = 1 k p i - - - ( 1 )
式中:m0为信号功率谱的零阶矩,pi为第i点功率值,求和范围为信号谱带宽。
信号功率谱的一阶矩为:
m 1 = Σ i = 1 k v i p i - - - ( 2 )
信号功率谱的二阶矩为:
m 2 = Σ i = 1 k v i 2 p i - - - ( 3 )
式中:m1为信号功率谱的一阶矩,m2为信号功率谱的二阶矩,vi为第i点速度值,求和范围为信号谱带宽。
信号功率估计值为:
Pr=m0         (4)
速度谱宽估计值为:
w ~ = 2 * [ m 2 m 0 - ( m 1 m 0 ) 2 ] 1 / 2 - - - ( 5 )
信噪比为:
SNR = 10 log ( Σ i = 1 k p i / P N ) - - - ( 6 )
式中:PN为噪声功率,求和范围为信号谱所占的频率范围,单位:dB。
经过上述处理,得到所处理高度层上的风场回波信号功率和信噪比估计值。
至此,实现了风廓线雷达信号处理。

Claims (1)

1.一种风廓线雷达信号处理方法,其特征在于具体步骤为:
第一步 搭建风廓线雷达信号处理***
风廓线雷达信号处理***包括:工业控制计算机(1)、数字中频接收机板(2)、***控制模块(3)和信号处理模块(4),其中工业控制计算机(1)包括:PCI总线(5)、共享内存(6)和以太网卡(7);
数字中频接收机板(2)与工业控制计算机(1)通过PCI总线(5)连接,数字中频接收机板(2)与风廓线雷达***中的模拟接收机(8)连接;工业控制计算机(1)与数据终端机(9)通过以太网卡(7)连接;***控制模块(3)与工业控制计算机(1)通过PCI总线(5)连接,信号处理模块(4)与工业控制计算机(1)通过PCI总线(5)连接;工业控制计算机(1)为***的运行平台,数字中频接收机板(2)用于采集雷达回波信号,***控制模块(3)用于控制和调度***的工作流程,信号处理模块(4)用于对回波采样信号进行信号处理;***控制模块(3)和信号处理模块(4)保存在工业控制计算机(1)上,***运行时加载到共享内存(6)中;
第二步 数字中频接收机板(2)采集雷达回波信号
***控制模块(3)首先进行数字中频接收机板(2)的初始化和参数加载,然后启动工作时序,数字中频接收机板(2)按照加载参数中设置的采样率、抽取因子、滤波系数,直接对雷达回波的中频信号进行A/D变换和采集,将采集的时域数据先存入数字中频接收机板(2)内的FIFO缓存中,通过DMA中断方式传输到工业控制计算机(1)内存中,由***控制模块(3)接收,再保存到工业控制计算机(1)的共享内存(6)中;***控制模块(3)边接收时域数据边进行计数,当计数值达到预设的数值N时,N=FFT点数×时域积累数×谱平均数,停止信号采集的工作时序,***控制模块(3)完成信号接收,然后对预设的时域数据保存标志作判断,如果标志为1,则保存到工业控制计算机(1)的硬盘上;
第三步 信号处理模块(4)进行信号处理
***控制模块(3)与信号处理模块(4)使用共享内存(6)交换接收到的时域数据,***控制模块(3)调用信号处理模块(4),向信号处理模块(4)传送信号处理参数和共享内存地址,由信号处理模块(4)利用信号处理参数对共享内存中的时域数据进行时域积累、脉冲压缩、去直流、加窗、FFT变换、谱平均的一系列信号处理过程,处理完毕后将生成的谱数据保存到共享内存中,向***控制模块返回处理完毕的标志信息,***控制模块(3)接收处理后生成的谱数据,并对预设的谱数据保存标志作判断,如果标志为1则保存到工业控制计算机(1)的硬盘上;
第四步 ***控制模块(3)进行谱数据的处理、显示和传送
***控制模块(3)对信号处理模块(4)生成的谱数据进行谱矩、信号功率和信噪比参数的计算,在本地显示器显示谱数据图形和处理的信息,用于对风场回波信号进行监测和评估,同时通过以太网卡将谱数据传送给数据终端机(9),作进一步的风谱信号识别、提取处理;
谱数据的处理过程如下:
信号功率谱的零阶矩为:
m 0 = Σ i = 1 k p i - - - ( 1 )
式中:m0为信号功率谱的零阶矩,pi为第i点功率值,求和范围为信号谱带宽;
信号功率谱的一阶矩为:
m 1 = Σ i = 1 k v i p i - - - ( 2 )
信号功率谱的二阶矩为:
m 2 = Σ i = 1 k v i 2 p i - - - ( 3 )
式中:m1为信号功率谱的一阶矩,m2为信号功率谱的二阶矩,vi为第i点速度值,求和范围为信号谱带宽;
信号功率估计值为:
Pr=m0             (4)
速度谱宽估计值为:
w ~ = 2 * [ m 2 m 0 - ( m 1 m 0 ) 2 ] 1 / 2 - - - ( 5 )
信噪比为:
SNR = 10 log ( Σ i = 1 k p i / P N ) - - - ( 6 )
式中:PN为噪声功率,求和范围为信号谱所占的频率范围,单位:dB;经过上述处理,得到所处理高度层上的风场回波信号功率和信噪比估计值;至此,实现了风廓线雷达信号处理。
CN2011101732501A 2011-06-25 2011-06-25 一种风廓线雷达信号处理方法 Active CN102323574B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101732501A CN102323574B (zh) 2011-06-25 2011-06-25 一种风廓线雷达信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101732501A CN102323574B (zh) 2011-06-25 2011-06-25 一种风廓线雷达信号处理方法

Publications (2)

Publication Number Publication Date
CN102323574A CN102323574A (zh) 2012-01-18
CN102323574B true CN102323574B (zh) 2013-03-06

Family

ID=45451354

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101732501A Active CN102323574B (zh) 2011-06-25 2011-06-25 一种风廓线雷达信号处理方法

Country Status (1)

Country Link
CN (1) CN102323574B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105527624A (zh) * 2014-09-29 2016-04-27 北京敏视达雷达有限公司 一种雷达回波动态估算噪声的方法和气象雷达***
CN107870332A (zh) * 2016-09-26 2018-04-03 北京无线电测量研究所 一种风廓线雷达大气湍流能量探测***及方法

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014002564A1 (ja) * 2012-06-25 2014-01-03 三菱電機株式会社 風計測コヒーレントライダ装置
CN102890272B (zh) * 2012-11-05 2013-12-04 中国航天科工集团第二研究院二十三所 一种毫米波云雷达信号处理方法
CN104122538B (zh) * 2013-04-24 2017-08-04 何平 确定风廓线雷达噪声功率的方法
CN104237892B (zh) * 2014-09-17 2017-02-15 北京无线电测量研究所 一种风廓线雷达大气温度探测***及方法
CN104345312B (zh) * 2014-10-27 2017-01-25 北京无线电测量研究所 一种毫米波云雷达的数据融合方法及***
CN104459650B (zh) * 2014-11-26 2017-02-22 北京无线电测量研究所 一种毫米波测云雷达实时标校***及方法
CN105022036B (zh) * 2015-08-26 2018-02-23 成都信息工程大学 风廓线雷达风速测定方法
CN106054195B (zh) * 2016-05-31 2018-07-06 中国民航大学 基于空时最优处理器的湍流谱宽估计方法
CN106556836B (zh) * 2016-10-18 2019-04-16 北京无线电测量研究所 一种风廓线雷达有降水无降水的自适应探测***及方法
CN111868565A (zh) * 2018-03-20 2020-10-30 三菱电机株式会社 雷达信号处理装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4819343B2 (ja) * 2004-11-09 2011-11-24 Necエンジニアリング株式会社 ウィンドプロファイラ・システム
CN100422708C (zh) * 2004-12-14 2008-10-01 中国科学院安徽光学精密机械研究所 近地面湍流廓线仪
JP2008014837A (ja) * 2006-07-06 2008-01-24 Toshiba Corp レーダ装置とその信号処理方法
US20080169975A1 (en) * 2007-01-12 2008-07-17 Young Paul Yee Process for generating spatially continuous wind profiles from wind profiler measurements

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105527624A (zh) * 2014-09-29 2016-04-27 北京敏视达雷达有限公司 一种雷达回波动态估算噪声的方法和气象雷达***
CN107870332A (zh) * 2016-09-26 2018-04-03 北京无线电测量研究所 一种风廓线雷达大气湍流能量探测***及方法

Also Published As

Publication number Publication date
CN102323574A (zh) 2012-01-18

Similar Documents

Publication Publication Date Title
CN102323574B (zh) 一种风廓线雷达信号处理方法
CN104483011A (zh) 一种旋转机械多通道振动信号的在线检测与分析***及方法
CN106523928B (zh) 基于声波实时数据二级筛选的管道泄漏检测方法
CN102830409B (zh) 一种基于压缩感知的导航信号采集方法
CN109490414B (zh) 一种盆式绝缘子缺陷故障检测方法
CN113419222B (zh) 基于雷达信号提取桥梁振动频率的方法及***
CN113008361B (zh) 一种变电站站界噪声抗环境干扰检测方法及装置
CN102495140A (zh) 一种面向超声波信号的小波去噪方法及***
WO2023005821A1 (zh) 生命体检测方法、终端及存储介质
CN110167114A (zh) 一种基于帧同步信号的水下声通信节点唤醒信号检测方法
CN111580038B (zh) 一种基于fpga的声学水下信标信号处理***
CN203133273U (zh) 基于cpci总线的高频地波雷达数据采集和处理装置
CN102176121A (zh) 数字超声经颅多普勒数字解调和信号处理方法及装置
CN103063169A (zh) 一种利用超声波监测危岩体崩塌的装置
CN108680843A (zh) 一种基于无线传输技术的开关柜局部放电检测方法
CN210982721U (zh) 基于fpga的低小慢目标探测雷达静态杂波抑制装置
CN116796571A (zh) 基于声发射特征的岩石构造应力预测***
CN104237892A (zh) 一种风廓线雷达大气温度探测***及方法
CN112034443B (zh) 一种雷达探测盲区计算方法及装置、存储介质和电子设备
CN103185881B (zh) 基于波导不变量的混响抑制方法
CN103166887B (zh) 一种边带调制信号分类的方法及装置
CN102901953A (zh) 一种相关峰锐化方法及装置
CN202334539U (zh) 一种用于数字信道化接收机的瞬时频域测量***
CN201555950U (zh) 基于dsp的雷达信号处理器
RU180088U1 (ru) Устройство построения радиолокационного изображения с помощью радиолокационной станции с синтезированной апертурой

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20191216

Address after: 214000 28 Weiming Road, Binhu District, Wuxi City, Jiangsu Province

Patentee after: Aerospace new weather Technology Co., Ltd

Address before: 100854 box 203, box 45, No. 142, Beijing, Beijing, Haidian District

Patentee before: No.23 Institute of No.20 Academy, CASIC

TR01 Transfer of patent right