CN102323574B - 一种风廓线雷达信号处理方法 - Google Patents
一种风廓线雷达信号处理方法 Download PDFInfo
- 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
Links
Images
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
本发明公开了一种风廓线雷达信号处理方法,其具体步骤为:第一步,搭建风廓线雷达信号处理系统,包括:工业控制计算机(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则保存到工业控制计算机的硬盘上。
第四步 系统控制模块进行谱数据的处理、显示和传送
系统控制模块对信号处理模块生成的谱数据进行谱矩、信号功率和信噪比参数的计算,在本地显示器显示谱数据图形和处理的信息,用于对风场回波信号进行监测和评估,同时通过以太网卡将谱数据传送给数据终端机,作进一步的风谱信号识别、提取处理。
谱数据的处理过程如下:
信号功率谱的零阶矩为:
式中:m0为信号功率谱的零阶矩,pi为第i点功率值,求和范围为信号谱带宽。
信号功率谱的一阶矩为:
信号功率谱的二阶矩为:
式中:m1为信号功率谱的一阶矩,m2为信号功率谱的二阶矩,vi为第i点速度值,求和范围为信号谱带宽。
信号功率估计值为:
Pr=m0 (4)
速度谱宽估计值为:
信噪比为:
式中: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,作进一步的风谱信号识别、提取处理。
谱数据的处理过程如下:
信号功率谱的零阶矩为:
式中:m0为信号功率谱的零阶矩,pi为第i点功率值,求和范围为信号谱带宽。
信号功率谱的一阶矩为:
信号功率谱的二阶矩为:
式中:m1为信号功率谱的一阶矩,m2为信号功率谱的二阶矩,vi为第i点速度值,求和范围为信号谱带宽。
信号功率估计值为:
Pr=m0 (4)
速度谱宽估计值为:
信噪比为:
式中: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),作进一步的风谱信号识别、提取处理;
谱数据的处理过程如下:
信号功率谱的零阶矩为:
式中:m0为信号功率谱的零阶矩,pi为第i点功率值,求和范围为信号谱带宽;
信号功率谱的一阶矩为:
信号功率谱的二阶矩为:
式中:m1为信号功率谱的一阶矩,m2为信号功率谱的二阶矩,vi为第i点速度值,求和范围为信号谱带宽;
信号功率估计值为:
Pr=m0 (4)
速度谱宽估计值为:
信噪比为:
式中:PN为噪声功率,求和范围为信号谱所占的频率范围,单位:dB;经过上述处理,得到所处理高度层上的风场回波信号功率和信噪比估计值;至此,实现了风廓线雷达信号处理。
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)
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9599714B2 (en) * | 2012-06-25 | 2017-03-21 | Mitsubishi Electric Corporation | Wind measurement coherent lidar |
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)
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 |
-
2011
- 2011-06-25 CN CN2011101732501A patent/CN102323574B/zh active Active
Cited By (2)
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) | 基于声波实时数据二级筛选的管道泄漏检测方法 | |
CN109490414B (zh) | 一种盆式绝缘子缺陷故障检测方法 | |
CN113419222B (zh) | 基于雷达信号提取桥梁振动频率的方法及系统 | |
CN113008361B (zh) | 一种变电站站界噪声抗环境干扰检测方法及装置 | |
CN102495140A (zh) | 一种面向超声波信号的小波去噪方法及系统 | |
WO2023005821A1 (zh) | 生命体检测方法、终端及存储介质 | |
CN110167114A (zh) | 一种基于帧同步信号的水下声通信节点唤醒信号检测方法 | |
CN111580038B (zh) | 一种基于fpga的声学水下信标信号处理系统 | |
CN102176121A (zh) | 数字超声经颅多普勒数字解调和信号处理方法及装置 | |
CN104237892B (zh) | 一种风廓线雷达大气温度探测系统及方法 | |
CN103063169A (zh) | 一种利用超声波监测危岩体崩塌的装置 | |
CN116796571A (zh) | 基于声发射特征的岩石构造应力预测系统 | |
CN104950292A (zh) | 一种气象雷达的二次回波识别方法及装置 | |
KR101295061B1 (ko) | Fmcw 레이더 신호처리모듈 및 그의 메모리 관리 방법 | |
CN112034443B (zh) | 一种雷达探测盲区计算方法及装置、存储介质和电子设备 | |
CN103185881B (zh) | 基于波导不变量的混响抑制方法 | |
CN201555950U (zh) | 基于dsp的雷达信号处理器 | |
CN102901953A (zh) | 一种相关峰锐化方法及装置 | |
CN202735511U (zh) | 一种风廓线雷达信号处理系统 | |
CN202334539U (zh) | 一种用于数字信道化接收机的瞬时频域测量系统 | |
RU180088U1 (ru) | Устройство построения радиолокационного изображения с помощью радиолокационной станции с синтезированной апертурой | |
CN109540515A (zh) | 一种振动信号处理方法 | |
CN109324313A (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 | ||
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 |