CN102269810A - 相干激光测风多普勒雷达激光回波信号频谱变换装置 - Google Patents

相干激光测风多普勒雷达激光回波信号频谱变换装置 Download PDF

Info

Publication number
CN102269810A
CN102269810A CN2011101056871A CN201110105687A CN102269810A CN 102269810 A CN102269810 A CN 102269810A CN 2011101056871 A CN2011101056871 A CN 2011101056871A CN 201110105687 A CN201110105687 A CN 201110105687A CN 102269810 A CN102269810 A CN 102269810A
Authority
CN
China
Prior art keywords
module
spectrum
echo signal
data
laser
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.)
Pending
Application number
CN2011101056871A
Other languages
English (en)
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.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
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 Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN2011101056871A priority Critical patent/CN102269810A/zh
Publication of CN102269810A publication Critical patent/CN102269810A/zh
Pending legal-status Critical Current

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

  • Optical Radar Systems And Details Thereof (AREA)

Abstract

一种相干激光测风多普勒雷达激光回波信号频谱变换装置,包括数据分路模块,第一频谱变换模块,第二频谱变换模块,蝶形因子乘法模块,匹配乘法模块,加法模块,功率谱模块。本发明可以实时地得到相干激光多普勒雷达功率谱数据,可以在低时钟频率下获得高数据吞吐率,具有成本低,硬件容易升级,稳定等特点。

Description

相干激光测风多普勒雷达激光回波信号频谱变换装置
技术领域
本发明属于高速信号处理领域,特别是一种相干激光测风多普勒雷达激光回波信号频谱变换装置。
背景技术
在各种气象参数中,如风速、云密度、云顶高、水蒸气浓度,温度和压力,对于许多用户需求,直接风速测量是最有价值的,而且传统的测量手段很难实时获得高分辨的风场垂直分布。经过实验和比较,激光多普勒测风雷达以其探测时空分辨率高等优越性能成为近些年来激光测风遥感重点方向。
多普勒测风激光雷达可以分为相干探测和非相干探测。非相干探测即直接探测,直接探测回波信号,利用发射激光和接收激光的频率差推算多普勒频移,从而获得风场信息。相干探测则是在直接探测的基础上引入一束参考光进行外差相干混频。相干探测相对于直接探测而言,拥有更高的探测灵敏度,显著的提高了接收信号的信噪比,已经有大量的理论和实验研究证明了相干探测的优势。
在多普勒相干测风激光雷达中,对于不同高度的大气层,低高度的大气层由于距离激光雷达近,从激光发射到激光雷达回波信号被激光雷达接收到的时间短,越高的大气层,从激光发射到激光雷达回波信号被激光雷达接收到的时间就越长,因而就可以根据回波信号被接收到的时间信息来得到大气层的距离信息,把回波信号划为一个个距离门,假设在该距离门中风速恒定,通过求取每个距离门中的风速信息就可以得到不同高度的风速信息,而在风速信息的反演算法中,傅里叶变换是整个算法中非常关键的部分。
一般的快速傅里叶变换装置直接采用传统的傅里叶变换方法,需要的系统时钟频率特别高,一般的器件很难满足如此高时钟频率的要求,同时系统时钟频率高时还对电磁兼容性能造成影响,很容易干扰其他设备的工作。而本发明所述的快速傅里叶变换装置由于把输入的信号分成两路,可以把系统时钟频率降为传统方法的一半,同时可以达到非常高的数据吞吐率,并且可以利用成本低廉的器件就可以实现系统的要求,因而在相干激光雷达信号处理中有着很强的应用背景。
发明内容
本发明的目的在于提高一种相干激光测风多普勒雷达激光回波信号频谱变换装置,以实现实时地得到相干激光测风多普勒雷达的功率谱数据,在低时钟频率下获得高数据吞吐率,该装置应具有成本低,硬件容易升级,稳定的特点。
本发明的技术解决方案如下:
一种相干激光测风多普勒雷达激光回波信号频谱变换装置,其特点在于包括数据分路模块,第一频谱变换模块,第二频谱变换模块,蝶形因子乘法模块,匹配乘法模块,加法模块和功率谱模块,上述模块的连接关系为:
输入的相干激光多普勒雷达回波信号x(n)进入数据分路模块,该数据分路模块的输出端分别与所述的第一频谱变换模块,第二频谱变换模块的输入端相连,所述的第一频谱变换模块的输出端经蝶形因子乘法模块后与所述的加法模块的输入端相连,所述的第二频谱变换模块的输出端经所述的匹配乘法模块后与所述的加法模块的输入端相连,该加法模块的输出端与所述的功率谱模块的输入端相连,所述的功率谱模块的输出端为本装置的输出端。
所述的第一频谱变换模块和第二频谱变换模块的点数为
Figure BSA00000482075000021
点,其中的取值范围为:
Figure BSA00000482075000023
所述的匹配乘法模块是一个乘数为常数的乘法模块,该乘数的取值为M,210≤M≤215
所述的蝶形因子乘法模块的两个乘数:一个为第一频谱变换模块输出的傅里叶频谱数据,另一个乘数为蝶形因子
Figure BSA00000482075000024
其中N为第一频谱变换模块的点数的两倍,k表示第一频谱变换模块输出的第k个数据与
Figure BSA00000482075000025
相乘,
Figure BSA00000482075000026
j表示虚数单位。
相干激光多普勒雷达回波信号的功率谱数据为相干激光多普勒雷达回波信号的N点快速傅里叶变换功率谱数据。
所述的相干激光测风多普勒雷达激光回波信号频谱变换装置的数据处理方法,包括下列步骤:
①相干激光测风多普勒雷达激光回波信号x(n),以下简称为回波信号x(n)进入数据分路模块,被分成两路数据:奇数信号x(2b+1)和偶数信号x(2b)分别输入所述的第一频谱变换模块和第二频谱变换模块;
②所述的第一频谱变换模块对所述的奇数信号x(2b+1)的
Figure BSA00000482075000031
点进行傅里叶变换得到傅里叶频谱
Figure BSA00000482075000032
输入到蝶形因子乘法模块,该蝶形因子乘法模块将所述的傅里叶频谱
Figure BSA00000482075000033
与蝶形因子
Figure BSA00000482075000034
相乘,得到蝶形因子频谱数据
Figure BSA00000482075000035
向所述的加法模块输出;
③所述的谱变换模块对输入的偶数信号x(2b)的
Figure BSA00000482075000036
点进行傅里叶变换得到傅里叶频谱信号输入到所述的匹配乘法模块,该匹配乘法模块将所述的频谱信号
Figure BSA00000482075000038
与一个为常数的匹配数M相乘,得到匹配频谱数据
Figure BSA00000482075000039
向所述的加法模块输入;
④所述的加法模块将所述的匹配频谱数据和蝶形因子频谱数据
Figure BSA000004820750000311
相加后,向所述的功率谱模块输入;
⑤所述的功率谱模块对输入的复数数据进行求模运算以及平方运算,得到并输出相干激光多普勒雷达回波信号的N点的功率谱数据。
本发明所述的相干激光测风多普勒雷达激光回波信号频谱变换装置,其工作原理为:
N点傅里叶变换的表达式为: X ( k ) = Σ n = 0 N - 1 x ( n ) w n NK - - - ( 1 )
其中k=0~N-1,n=0~N-1, W N nk = e - j * 2 * π * n * k N
把输入的信号x(n)分解为偶数信号x(2b)和奇数信号x(2b+1),代入傅里叶变换表达式(1)中可以得到
X ( k ) = M * Σ b = 0 N / 2 - 1 x ( 2 b ) W N / 2 bk + M * W N k Σ b = 0 N / 2 - 1 x ( 2 b + 1 ) W N / 2 bk - - - ( 2 )
可以看出前半部分为偶数点的
Figure BSA00000482075000042
点傅里叶变换,后半部分奇数点的
Figure BSA00000482075000043
点傅里叶变换,一个N点傅里叶变换可以分解为两个
Figure BSA00000482075000044
点傅里叶变换,
Figure BSA00000482075000045
即为蝶形因子。由于
Figure BSA00000482075000046
的值都是小数,而硬件不便于小数的表示,因而装置中把
Figure BSA00000482075000047
与常数M相乘(匹配乘法模块),使得扩大为整数,同时也需要在偶数点的
Figure BSA00000482075000049
点傅里叶变换中乘以常数M以便可以得到正确的傅里叶变换数据,由于x(2b)只有
Figure BSA000004820750000410
点,因而式(2)计算出的傅里叶变换只是整个N点傅里叶变换频谱的前
Figure BSA000004820750000411
点,而由于激光雷达回波信号输入的数据都是实数数据,而实数数据的傅里叶变换频谱具有对称性,因而就可以根据前
Figure BSA000004820750000412
点傅里叶变换频谱推算出整个N点的傅里叶变换频谱。
综上所述,本发明的技术效果如下:
1、由于上述各个模块都是采用硬件实现的,因此该激光雷达回波信号处理装置具有速度快,稳定性好等特点。
2、由于采用双路并行的频谱变换处理,因此可以把时钟频率降低到没有并行处理前的一半,同时数据吞吐率提高了一倍,能够满足相干激光雷达回波信号处理的需要。
3、由于时钟频率降低到了没有并行处理前的一半,因此该装置对硬件的要求也相应的降低,从而可以用价格低廉的硬件来实现快速傅里叶变换算法,有效的降低了系统成本。
总之,本发明可以实时地得到相干激光测风多普勒雷达功率谱数据,可以在低时钟频率下获得高数据吞吐率,具有成本低,硬件容易升级,稳定等特点。
附图说明
图1是相干激光测风多普勒雷达激光回波信号频谱变换装置示意图。
图中:1-相干激光多普勒雷达回波信号,2-数据分路模块,3-第一频谱变换模块,4-蝶形因子乘法模块,5-第二频谱变换模块,6-匹配乘法模块,7-加法模块,8-功率谱模块,9-相干激光多普勒雷达回波信号的功率谱数据,10-奇数信号,11-偶数信号,12-蝶形因子频谱数据,13-匹配频谱数据
具体实施方式
下面结合实施例和附图对本发明作进一步说明,但不应以此限制本发明的保护范围。
参考图1,图1是本发明相干激光测风多普勒雷达激光回波信号频谱变换装置,包括数据分路模块2,第一频谱变换模块3,第二频谱变换模块5,蝶形因子乘法模块4,匹配乘法模块6,加法模块7,功率谱模块8。上述模块的连接关系为:
输入的相干激光多普勒雷达回波信号x(n)1进入数据分路模块2,该数据分路模块2的输出端分别与所述的第一频谱变换模块3,第二频谱变换模块5的输入端相连,所述的第一频谱变换模块3的输出端经蝶形因子乘法模块4后与所述的加法模块7的输入端相连,所述的第二频谱变换模块5的输出端经所述的匹配乘法模块6后与所述的加法模块7的输入端相连,该加法模块7的输出端与所述的功率谱模块8的输入端相连,所述的功率谱模块8的输出端为本装置的输出端。
输入的相干激光多普勒雷达回波信号x(n)1进入数据分路模块2,在数据分路模块2中,回波信号x(n)1被分成两路数据输出,一路为输入回波信号x(n)1的奇数信号x(2b+1)10,一路为输入回波信号x(n)1的偶数信号x(2b)11;
奇数信号x(2b+1)10进入第一频谱变换模块3中,经过第一频谱变换模块3变换后输出奇数信号x(2b+1)10的点的傅里叶频谱
Figure BSA00000482075000052
并且输入到蝶形因子乘法模块4中,在蝶形因子乘法模块4中,输入的奇数信号x(2b+1)10的傅里叶频谱
Figure BSA00000482075000053
与蝶形因子
Figure BSA00000482075000054
相乘,相乘后的结果输出得到蝶形因子频谱数据
Figure BSA00000482075000061
数据分路模块2输出的回波信号x(n)1的偶数信号x(2b)11进入第二频谱变换模块5中,经过第二频谱变换模块5变换后输出偶数信号x(2b)11的
Figure BSA00000482075000062
点的傅里叶频谱信号
Figure BSA00000482075000063
并且输入到匹配乘法模块6中,在匹配乘法模块6中,偶数信号x(2b)11的傅里叶频谱信号
Figure BSA00000482075000064
与一个为常数的匹配数M相乘从而输出得到匹配频谱数据
Figure BSA00000482075000065
匹配频谱数据13和蝶形因子频谱数据
Figure BSA00000482075000067
12一起进入复数加法模块7相加后输出的数据进入功率谱模块8中,在功率谱模块8中把复数加法模块7相加后输出的复数数据进行求模运算以及平方运算从而输出相干激光多普勒雷达回波信号的N点的功率谱数据9。
下面是具体实施例
相干激光雷达的回波信号1的数据速率为200M*14bps的数据经过数据分路模块2后被分为奇数信号10(数据速率为100M*14bps)和偶数信号11(数据速率为100M*14bps)。
其中奇数信号10进入第一频谱变换模块3进行
Figure BSA00000482075000068
(不限于该值,取值范围可以为
Figure BSA00000482075000069
)点的傅里叶变换,经过第一频谱变换模块3变换后的频谱信号进入蝶形因子乘法模块4中,经过蝶形因子乘法模块4处理后输出蝶形因子频谱数据12(数据速率为100M*14bps)。
偶数信号11进入第二频谱变换模块5中进行
Figure BSA000004820750000610
(不限于该值,取值范围可以为
Figure BSA000004820750000611
)点的傅里叶变换,经过第二频谱变换模块5变换后的频谱信号进入匹配乘法模块6中,在匹配乘法模块6中与一个常数M=4096(不限于该值,取值范围可以为210≤M≤215)相乘,相乘后的结果为匹配频谱数据13(数据速率为100M*14bps)。
匹配频谱数据13和蝶形因子频谱数据12一起进入复数加法模块7相加后输出的数据进入功率谱模块8中,经过功率谱模块8处理后输出相干激光多普勒雷达回波信号的N=256点傅里叶变换数据(数据速率为100M*14bps)。
实验表明,本发明可满足相干多普勒激光雷达回波信号频谱变换的要求,具有硬件容易升级,集成度高,性能稳定的特点。该装置还可以用于其它需要频谱变换的测量,具有很高的通用性。

Claims (5)

1.一种相干激光测风多普勒雷达激光回波信号频谱变换装置,其特征在于包括数据分路模块(2),第一频谱变换模块(3),第二频谱变换模块(5),蝶形因子乘法模块(4),匹配乘法模块(6),加法模块(7),功率谱模块(8),上述模块的连接关系为:
输入的相干激光多普勒雷达回波信号x(n)(1)进入数据分路模块(2),该数据分路模块(2)的输出端分别与所述的第一频谱变换模块(3),第二频谱变换模块(5)的输入端相连,所述的第一频谱变换模块(3)的输出端经蝶形因子乘法模块(4)后与所述的加法模块(7)的输入端相连,所述的第二频谱变换模块(5)的输出端经所述的匹配乘法模块(6)后与所述的加法模块(7)的输入端相连,该加法模块(7)的输出端与所述的功率谱模块(8)的输入端相连,所述的功率谱模块(8)的输出端为本装置的输出端。
2.根据权利要求1所述的相干激光测风多普勒雷达激光回波信号频谱变换装置,其特征在于所述的第一频谱变换模块(3)和第二频谱变换模块(5)的点数为
Figure FSA00000482074900011
点,其中
Figure FSA00000482074900012
的取值范围为:
3.根据权利要求1所述的相干激光测风多普勒雷达激光回波信号频谱变换装置,其特征在于所述的匹配乘法模块(6)是一个乘数为常数的乘法模块,该乘数的取值为M,210≤M≤215
4.根据权利要求1所述的相干激光测风多普勒雷达激光回波信号频谱变换装置,其特征在于所述的蝶形因子乘法模块(4)的两个乘数:一个为第一频谱变换模块(3)输出的傅里叶频谱数据,另一个乘数为蝶形因子
Figure FSA00000482074900014
其中N为第一频谱变换模块(3)的点数的两倍,k表示第一频谱变换模块(3)输出的第k个数据与
Figure FSA00000482074900015
相乘,
Figure FSA00000482074900016
j表示虚数单位。
5.根据权利要求1所述的相干激光测风多普勒雷达激光回波信号频谱变换装置,其特征在于相干激光多普勒雷达回波信号的功率谱数据(9)为相干激光多普勒雷达回波信号(1)的N点快速傅里叶变换功率谱数据。
CN2011101056871A 2011-04-26 2011-04-26 相干激光测风多普勒雷达激光回波信号频谱变换装置 Pending CN102269810A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101056871A CN102269810A (zh) 2011-04-26 2011-04-26 相干激光测风多普勒雷达激光回波信号频谱变换装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101056871A CN102269810A (zh) 2011-04-26 2011-04-26 相干激光测风多普勒雷达激光回波信号频谱变换装置

Publications (1)

Publication Number Publication Date
CN102269810A true CN102269810A (zh) 2011-12-07

Family

ID=45052167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101056871A Pending CN102269810A (zh) 2011-04-26 2011-04-26 相干激光测风多普勒雷达激光回波信号频谱变换装置

Country Status (1)

Country Link
CN (1) CN102269810A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646422A (zh) * 2016-09-28 2017-05-10 成都凯天电子股份有限公司 增强相干测风雷达多普勒频移信号信噪比的预处理系统
CN107390229A (zh) * 2017-06-14 2017-11-24 南京牧镭激光科技有限公司 一种测风激光雷达信号的处理方法及其处理装置
CN110672875A (zh) * 2019-10-21 2020-01-10 中国科学院电子学研究所 基于Chirp-Z变换的表面水流速度检测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101545977A (zh) * 2009-05-05 2009-09-30 中国科学院上海光学精密机械研究所 相干测风激光雷达信号检测解调装置
CN101825710A (zh) * 2009-12-24 2010-09-08 哈尔滨工业大学 一种2μm全光纤相干激光多普勒测风雷达系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101545977A (zh) * 2009-05-05 2009-09-30 中国科学院上海光学精密机械研究所 相干测风激光雷达信号检测解调装置
CN101825710A (zh) * 2009-12-24 2010-09-08 哈尔滨工业大学 一种2μm全光纤相干激光多普勒测风雷达系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
史永奇: "基于FPGA设计雷达目标检测信号处理器", 《微处理机》, no. 4, 31 August 2007 (2007-08-31) *
张芳沛 等: "相干多普勒测风激光雷达", 《应用光学》, vol. 30, no. 6, 30 November 2009 (2009-11-30) *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646422A (zh) * 2016-09-28 2017-05-10 成都凯天电子股份有限公司 增强相干测风雷达多普勒频移信号信噪比的预处理系统
CN106646422B (zh) * 2016-09-28 2020-06-09 成都凯天电子股份有限公司 增强相干测风雷达多普勒频移信号信噪比的预处理系统
CN107390229A (zh) * 2017-06-14 2017-11-24 南京牧镭激光科技有限公司 一种测风激光雷达信号的处理方法及其处理装置
CN107390229B (zh) * 2017-06-14 2019-02-12 南京牧镭激光科技有限公司 一种测风激光雷达信号的处理方法及其处理装置
CN110672875A (zh) * 2019-10-21 2020-01-10 中国科学院电子学研究所 基于Chirp-Z变换的表面水流速度检测方法

Similar Documents

Publication Publication Date Title
CN103018730B (zh) 分布式子阵波达方向估计方法
CN102628937B (zh) 基于广义keystone变换和非相参积累的雷达检测方法
CN104678372B (zh) 正交频分复用雷达超分辨距离与角度值联合估计方法
CN105229431A (zh) 具有改进的距离确定的料位测量
CN101900761B (zh) 一种高准确度非整周期采样谐波分析测量方法
CN103023586A (zh) 一种天波超视距雷达电离层信道模型
CN111830493A (zh) 水中目标收发分置中高频声目标强度预报系统及方法
CN104833974A (zh) 基于图像谱压缩的sar成像快速后向投影方法
CN105372668A (zh) 一种相位式激光测距方法
CN103869298A (zh) 一种分布式mimo天波超视距雷达海杂波仿真方法
CN103983957A (zh) 一种多普勒偏移测量方法及其装置
CN111487609A (zh) 多频连续波mimo阵列雷达系统及其目标参数估计方法
CN102269810A (zh) 相干激光测风多普勒雷达激光回波信号频谱变换装置
CN101788671A (zh) 应用于外差探测啁啾调幅激光测距装置的多周期调制方法
CN106443614B (zh) 高超声速目标测加速度方法
CN105548959B (zh) 一种基于稀疏重建的多传感器多目标的定位方法
CN102866383A (zh) 一种基于空域自适应滤波的波达方向估计方法
Hysell et al. A multistatic HF beacon network for ionospheric specification in the Peruvian sector
CN103048641B (zh) 具有联合约束优化形式的矢量阵稳健聚焦处理方法
CN108614234B (zh) 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法
CN103852443A (zh) 免标定调制光谱的实现方法
CN113671495A (zh) 基于Zynq平台的太赫兹雷达探测系统以及方法
CN112130140A (zh) 一种伪码调相引信系统提高测速测距精度的方法
CN108020811A (zh) 基于目标源相移差分技术的1维均匀线性阵列测向方法
Zrnic et al. Matched filter criteria and range weighting for weather radar

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20111207