CN105022036A - 风廓线雷达风速测定方法 - Google Patents

风廓线雷达风速测定方法 Download PDF

Info

Publication number
CN105022036A
CN105022036A CN201510530246.4A CN201510530246A CN105022036A CN 105022036 A CN105022036 A CN 105022036A CN 201510530246 A CN201510530246 A CN 201510530246A CN 105022036 A CN105022036 A CN 105022036A
Authority
CN
China
Prior art keywords
signal
frequency
wind speed
wind
profile radar
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
CN201510530246.4A
Other languages
English (en)
Other versions
CN105022036B (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.)
Chengdu University of Information Technology
Original Assignee
Chengdu University of Information Technology
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 Chengdu University of Information Technology filed Critical Chengdu University of Information Technology
Priority to CN201510530246.4A priority Critical patent/CN105022036B/zh
Publication of CN105022036A publication Critical patent/CN105022036A/zh
Application granted granted Critical
Publication of CN105022036B publication Critical patent/CN105022036B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及风廓线雷达。本发明公开了一种风廓线雷达风速测定方法。本发明的风廓线雷达风速测定方法，包括步骤：a、接收回波信号进行处理，得到中频信号；b、对中频信号进行数据采样；c、对采样信号进行汉明窗加权处理；d、对加权信号进行快速傅里叶变换；e、对傅里叶变换后的信号进行样条函数插值，找到频谱幅度最大值f2；f、根据下式计算平均风速Vr：其中，f1为发射信号与本振频率的差频，λ为发射信号波长，f2-f1＝Δf为多普勒频移。本发明通过基于样条函数的频谱插值，使得信号的频谱信息更加精细，应用样条插值后的频谱进行中心频率估计，比传统的中心频率测量方法精度提高很多，从而可以大大提高风廓线雷达风速测定的精确度。

Description

风廓线雷达风速测定方法
技术领域
本发明涉及风廓线雷达，特别涉及风廓线雷达风速测定方法。
背景技术
风廓线雷达靠发射单频脉冲信号测量风场信息，其中平均风速信息是通过测量发射脉冲的中心频率的多普勒频移来实现，如测得的发射脉冲和接收脉冲的频谱中心移动(多普勒频移)为△f＝f2-f1，其中f1和f2分别是发射信号和接收信号中心频率，则平均风速Vr可以根据风速公式：
V r = Δ f λ 2 = ( f 2 - f 1 ) λ 2 ,
计算得到。式中，λ为发射信号波长。
通常，发射信号的中心频率f1就是发射脉冲与本振频率的差，即固定中频，为风廓线雷达内部设定参数。而接收信号中心频率f2是根据接收的回波脉冲进行处理后，通过一定的运算程序进行估计得到的。
风廓线雷达风速测定流程如图1所示，可以看出，回波脉冲的中心频率f2估计精确度，是决定平均风速测量精确度的关键因素。传统方法是对采样后的信号直接进行傅里叶变换，然后寻找频谱峰值频率作为中心频率f2的估计值，中心频率f2计算流程如图2所示。由于发射脉冲宽度有限，所以傅里叶变换后的信号频谱的分辨率有限，这就限制了信号中心频率估计的精确度。
发明内容
本发明所要解决的技术问题，就是提供一种风廓线雷达风速测定方法，通过基于样条函数的频谱插值，使得信号的频谱信息更加精细，提高回波脉冲中心频率f2估计精确度。
本发明解决所述技术问题采用的技术方案是，风廓线雷达风速测定方法，其特征在于，包括以下步骤：
a、接收回波信号进行处理，得到中频信号；
b、对中频信号进行数据采样；
c、对采样信号进行汉明窗加权处理；
d、对加权信号进行快速傅里叶变换；
e、对傅里叶变换后的信号进行样条函数插值，找到频谱幅度最大值f2；
f、根据下式计算平均风速Vr：
V r = ( f 2 - f 1 ) λ 2
其中，f1为发射信号与本振频率的差频，λ为发射信号波长，f2-f1＝Δf为多普勒频移。
具体的，所述步骤a具体为对接收信号进行下变频处理，得到中频信号。
具体的，所述风廓线雷达发射信号波长为1.55μm，且f1＝81.5MHz。
具体的，步骤b中采样频率为400MHz。
具体的，步骤e所述样条函数插值具体为10000倍样条函数插值。
本发明的有益效果是，通过基于样条函数的频谱插值，使得信号的频谱信息更加精细，应用样条插值后的频谱进行中心频率估计，比传统的中心频率测量方法精度提高很多，从而可以大大提高风廓线雷达风速测定的精确度。
附图说明
图1是现有技术风廓线雷达风速测定流程图；
图2是现有技术风廓线雷达回波信号中频计算流程图；
图3是本发明风廓线雷达风速测定流程图。
具体实施方式
下面结合附图及实施例，详细描述本发明的技术方案。
风廓线雷达风速测定方法，首先向探测区发射中心频率已知的电磁脉冲，通过大气散射后回波信号返回风廓线雷达接收单元。该回波信号由于风速的影响产生了多普勒频移，根据多普勒频移可以计算平均风速。本发明的风廓线雷达风速测定方法，主要流程如图3所示，包括步骤：
第一步，风廓线雷达接收回波信号进行处理，获取中频信号。
该步骤主要是对回波信号进行整形滤波等常规处理，以降低各种噪声的影响。然后将回波信号与本振频率信号混频，通过下变频处理得到中频信号。
第二步，对中频信号进行数据采样。
该步骤根据设定的采样频率对连续的中频信号进行数据采样，得到离散的采样信号。
第三步，对采样后的信号进行汉明窗加权处理。
对采样后的信号进行汉明窗加权后，信号的边沿将会更平缓，从而减小信号在快速傅立叶变换后的频谱泄漏，使频谱分辨效果更好。
第四步，对加权信号进行快速傅里叶变换。
该步骤对上述加权信号在频域进行离散变换，得到频谱分布。
第五步，对傅里叶变换后的信号进行样条函数插值，找到频谱幅度最大值f2。
第六步，根据下式计算平均风速Vr：
V r = ( f 2 - f 1 ) λ 2
其中，f1为发射信号与本振频率的差频，λ为发射信号波长，f2-f1＝Δf为多普勒频移。
与现有技术风速测定流程图1比较，本发明的风速测定方法的区别就在于回波脉冲的中频计算方法不同。根据图2现有技术风廓线雷达回波信号中频计算流程图可知，现有技术是对采样后的信号直接进行傅里叶变换，然后在频率-振幅图中寻找振幅最大的频率作为中心频率f2的值。而本发明回波信号中频计算方法是，对采样后的信号先进行汉明窗加权，然后进行FFT(快速傅里叶变换)运算，并通过样条函数插值展宽频谱图，最后根据频谱峰值得到回波信号中频f2。
实施例
设定发射脉冲为固定频率的高频激光脉冲，其波长λ＝1.55μm，发射脉冲与本振频率的差，即f1为81.5MHz，假设发射脉冲发射出去后由于风速的影响产生了1MHz的频移，因此接收脉冲下变频后的中心频率f2为82.5MHz，即△f＝f2-f1＝82.5-81.5＝1MHz，风速为 V r = Δ f . λ 2 = 1 × 10 6 × 1.55 × 10 - 6 2 = 0.775 m / s . 为了验证本发明提出的算法的有效性，仿真产生一个频率为82.5MHz的正弦信号，代表经过下变频后的回波脉冲信号，为了模拟真实的探测环境，在脉冲信号中加入了高斯噪声，信号采样频率为400MHz，采样后分别用传统方法和本发明提出的方法进行中心频率求取，求取过程与结果如下：
传统方法：
(1)对82.5MHz的接收信号进行快速傅里叶变换
(2)找到频谱幅度最大的位置，即谱峰位置，就是中心频率，为81.25MHz
(3)计算发射脉冲信号与回波脉冲信号的频率差，即多普勒频移81.25-81.5＝-0.25MHz
(4)计算风速 V r = Δ f . λ 2 = - 0.25 × 10 6 × 1.55 × 10 - 6 2 = - 0.19375 m / s .
本发明提出的方法：
(1)对82.5MHz的接收信号进行汉明窗加权。汉明窗函数为：
w ( n ) = 0.54 - 0.46 c o s ( 2 π n N ) , 0 ≤ n ≤ N
其中N是窗函数的样点数目，称为窗函数的长度。
(2)对加权后的信号进行快速傅里叶变换
(3)对信号频谱进行10000倍样条函数插值
(4)找到频谱幅度最大的位置，即谱峰位置，就是中心频率，为82.41MHz
(5)计算发射脉冲信号与接收脉冲信号的频率差，即多普勒频移82.41-81.5＝0.91MHz
(6)计算风速 V r = Δ f . λ 2 = 0.91 × 10 6 × 1.55 × 10 - 6 2 = 0.70525 m / s .
从上述中心频率和风速估计结果可以看出，传统方法对脉冲中心频率的估计结果误差较大，所以导致风速估计的误差较大，而本发明提出的方法较为准确的估计出了回波脉冲的中心频率和风速。
为了进一步验证算法的有效性，参照实施例，仿真产生一组接收脉冲下变频后的脉冲信号，加入高斯噪声，对其中心频率进行估计，并由此计算出多普勒频移和计算风速。
仿真时，发射脉冲频率固定为波长λ＝1.55μm的激光脉冲，它与本振的频率差即固定中频为81.5MHz，设定一些与实际风速对应的多普勒频移，得到接收脉冲下变频后的中频脉冲信号，然后对中频脉冲中心频率进行估计。真实中心频率、估计中心频率以及对应的真实风速和估计风速如表1所示。
表1
从表1可以看出，对于两个中心频率间隔较近的接收脉冲，下变频后传统的中心频率估计方法根本无法将它们区分开来，测得的风速不准确，而且分辨率低，存在模糊。而本发明提出的方法可以较为准确的测出风速，而且对于中心频率间隔比较小的脉冲也能区分出来，所以测风的分辨率也较高。

Claims (5)

1.风廓线雷达风速测定方法，其特征在于，包括以下步骤：
a、接收回波信号进行处理，得到中频信号；
b、对中频信号进行数据采样；
c、对采样信号进行汉明窗加权处理；
d、对加权信号进行快速傅里叶变换；
e、对傅里叶变换后的信号进行样条函数插值，找到频谱幅度最大值f2；
f、根据下式计算平均风速Vr：
V r = ( f 2 - f 1 ) λ 2
其中，f1为发射信号与本振频率的差频，λ为发射信号波长，f2-f1＝Δf为多普勒频移。
2.根据权利要求1所述的风廓线雷达风速测定方法，其特征在于，所述步骤a具体为对接收信号进行下变频处理，得到中频信号。
3.根据权利要求1或2所述的风廓线雷达风速测定方法，其特征在于，所述风廓线雷达发射信号波长为1.55μm，且f1＝81.5MHz。
4.根据权利要求3所述的风廓线雷达风速测定方法，其特征在于，步骤b中采样频率为400MHz。
5.根据权利要求3所述的风廓线雷达风速测定方法，其特征在于，步骤e所述样条函数插值具体为10000倍样条函数插值。
CN201510530246.4A 2015-08-26 2015-08-26 风廓线雷达风速测定方法 Expired - Fee Related CN105022036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510530246.4A CN105022036B (zh) 2015-08-26 2015-08-26 风廓线雷达风速测定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510530246.4A CN105022036B (zh) 2015-08-26 2015-08-26 风廓线雷达风速测定方法

Publications (2)

Publication Number Publication Date
CN105022036A true CN105022036A (zh) 2015-11-04
CN105022036B CN105022036B (zh) 2018-02-23

Family

ID=54412141

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510530246.4A Expired - Fee Related CN105022036B (zh) 2015-08-26 2015-08-26 风廓线雷达风速测定方法

Country Status (1)

Country Link
CN (1) CN105022036B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629237A (zh) * 2015-12-25 2016-06-01 中国航天空气动力技术研究院 利用机载雷达对风场进行危险性评估的方法
CN105629238A (zh) * 2015-12-25 2016-06-01 中国航天空气动力技术研究院 利用机载雷达对风场进行数字重建的方法
CN106772440A (zh) * 2017-01-12 2017-05-31 杭州赛尤企业管理咨询有限公司 采用变频激光测风雷达的测风***及变频控制方法
CN109283535A (zh) * 2018-11-16 2019-01-29 北京无线电测量研究所 非均匀风场的探测方法
CN111198361A (zh) * 2018-11-19 2020-05-26 深圳市速腾聚创科技有限公司 数据处理***和方法
RU2786132C1 (ru) * 2022-02-07 2022-12-19 Акционерное общество Центральное конструкторское бюро аппаратостроения Способ формирования и обработки радиолокационных сигналов в импульсно-доплеровской метеорологической рлс
CN117741685A (zh) * 2024-02-21 2024-03-22 南京信息工程大学 一种风速反演方法、***、存储介质及设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1262793A1 (en) * 2001-05-30 2002-12-04 Murata Manufacturing Co., Ltd. Method and apparatus for removing a DC-offset in the frequency spectrum before performing Fourier transform in a radar
CN102323574A (zh) * 2011-06-25 2012-01-18 中国航天科工集团第二研究院二十三所 一种风廓线雷达信号处理方法
CN202433521U (zh) * 2011-10-17 2012-09-12 中国人民解放军理工大学气象学院 风廓线雷达湍流目标检测处理板
CN103630905A (zh) * 2013-08-29 2014-03-12 中国科学院电子学研究所 阵列天线sar极坐标交叠子孔径成像方法
CN104035096A (zh) * 2014-06-06 2014-09-10 南京大学 一种基于多普勒天气雷达的垂直风廓线非线性反演方法
CN104122538A (zh) * 2013-04-24 2014-10-29 何平 确定风廓线雷达噪声功率的方法
CN104133216A (zh) * 2014-07-17 2014-11-05 北京无线电测量研究所 一种获取低空风廓线的雷达探测方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1262793A1 (en) * 2001-05-30 2002-12-04 Murata Manufacturing Co., Ltd. Method and apparatus for removing a DC-offset in the frequency spectrum before performing Fourier transform in a radar
CN102323574A (zh) * 2011-06-25 2012-01-18 中国航天科工集团第二研究院二十三所 一种风廓线雷达信号处理方法
CN202433521U (zh) * 2011-10-17 2012-09-12 中国人民解放军理工大学气象学院 风廓线雷达湍流目标检测处理板
CN104122538A (zh) * 2013-04-24 2014-10-29 何平 确定风廓线雷达噪声功率的方法
CN103630905A (zh) * 2013-08-29 2014-03-12 中国科学院电子学研究所 阵列天线sar极坐标交叠子孔径成像方法
CN104035096A (zh) * 2014-06-06 2014-09-10 南京大学 一种基于多普勒天气雷达的垂直风廓线非线性反演方法
CN104133216A (zh) * 2014-07-17 2014-11-05 北京无线电测量研究所 一种获取低空风廓线的雷达探测方法及装置

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629237A (zh) * 2015-12-25 2016-06-01 中国航天空气动力技术研究院 利用机载雷达对风场进行危险性评估的方法
CN105629238A (zh) * 2015-12-25 2016-06-01 中国航天空气动力技术研究院 利用机载雷达对风场进行数字重建的方法
CN106772440A (zh) * 2017-01-12 2017-05-31 杭州赛尤企业管理咨询有限公司 采用变频激光测风雷达的测风***及变频控制方法
CN106772440B (zh) * 2017-01-12 2023-09-19 杭州赛尤新能源科技有限公司 采用变频激光测风雷达的测风***及变频控制方法
CN109283535A (zh) * 2018-11-16 2019-01-29 北京无线电测量研究所 非均匀风场的探测方法
CN109283535B (zh) * 2018-11-16 2022-12-20 北京无线电测量研究所 非均匀风场的探测方法
CN111198361A (zh) * 2018-11-19 2020-05-26 深圳市速腾聚创科技有限公司 数据处理***和方法
RU2786132C1 (ru) * 2022-02-07 2022-12-19 Акционерное общество Центральное конструкторское бюро аппаратостроения Способ формирования и обработки радиолокационных сигналов в импульсно-доплеровской метеорологической рлс
CN117741685A (zh) * 2024-02-21 2024-03-22 南京信息工程大学 一种风速反演方法、***、存储介质及设备
CN117741685B (zh) * 2024-02-21 2024-05-14 南京信息工程大学 一种风速反演方法、***、存储介质及设备

Also Published As

Publication number Publication date
CN105022036B (zh) 2018-02-23

Similar Documents

Publication Publication Date Title
CN105022036A (zh) 风廓线雷达风速测定方法
CN107861117B (zh) 一种适用于连续波周界监视雷达的多目标参数测量方法
EP2884299B1 (en) Speed determination of a target
CN108469608B (zh) 一种运动平台雷达多普勒质心精确估计方法
US6396436B1 (en) Radar method for measuring distances between and relative speeds of a vehicle and one or more obstacles
JP2020067455A (ja) 妨害信号抑圧を行うｆｍｃｗレーダー
CN109375202B (zh) 一种基于车载毫米波雷达的车辆测距测速方法
CN110850400B (zh) 基于干涉仪测向的lfmcw雷达多目标检测方法
CN104919331A (zh) 雷达装置
CN107843892A (zh) 一种基于最小二乘法的高速目标多普勒测速方法
CN106199582A (zh) 应用于室内的双频连续波雷达人体定位跟踪的方法
CN103278807A (zh) 双通道欠采样线扫频脉冲信号的时延估计方法
CN104297740A (zh) 基于相位分析的雷达目标多普勒谱估计方法
CN106872969A (zh) 基于mtd脉冲积累及滑动处理的雷达目标角度估计方法
CN104007435A (zh) 一种基于中频相邻回波相位差的精确测速方法
CN109164441B (zh) 一种雷达测距的方法
CN110109089B (zh) 一种线性调频连续波探测***距离测量精度的改进方法
CN111007473A (zh) 基于距离频域自相关函数的高速微弱目标检测方法
US20120309324A1 (en) Measurement method and apparatus
CN108120975A (zh) 基于梯形连续波的雷达测速测距方法
CN117008130A (zh) 一种基于毫米波雷达的车辆目标速度变化测量方法及***
CN115308434A (zh) 一种超声波测速方法及***
CN112346045B (zh) 一种多频率点雷达测速方法
CN109521222B (zh) 一种提高激光测速精度的方法
JP2007127529A (ja) 距離測定装置および距離測定方法

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180223

Termination date: 20180826