CN115951332B - 二元脉冲幅度编码测风激光雷达风速测量方法及系统 - Google Patents
二元脉冲幅度编码测风激光雷达风速测量方法及系统 Download PDFInfo
- Publication number
- CN115951332B CN115951332B CN202310246035.2A CN202310246035A CN115951332B CN 115951332 B CN115951332 B CN 115951332B CN 202310246035 A CN202310246035 A CN 202310246035A CN 115951332 B CN115951332 B CN 115951332B
- Authority
- CN
- China
- Prior art keywords
- echo
- pulse amplitude
- signal
- decoding
- binary pulse
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 74
- 238000001228 spectrum Methods 0.000 claims abstract description 65
- 230000003287 optical effect Effects 0.000 claims abstract description 52
- 238000001514 detection method Methods 0.000 claims abstract description 29
- 230000001934 delay Effects 0.000 claims abstract description 19
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims abstract description 18
- 238000012935 Averaging Methods 0.000 claims abstract description 10
- 238000000691 measurement method Methods 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 10
- 239000013307 optical fiber Substances 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 12
- 238000010586 diagram Methods 0.000 description 10
- 238000002592 echocardiography Methods 0.000 description 8
- 239000000835 fiber Substances 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000010183 spectrum analysis Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 108091026890 Coding region Proteins 0.000 description 1
- 108091029480 NONCODE Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本申请涉及二元脉冲幅度编码测风激光雷达风速测量方法及系统,该方法包括:根据二元脉冲幅度编码矩阵对单频连续激光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号;将二元脉冲幅度编码信号经光放大器放大后,通过环形器传输至光学天线并向待测大气区域发射;接收待测大气区域的回波信号后,确定回波信号中回波脉冲序列的解码矩阵;根据解码矩阵对回波脉冲序列进行解码处理,得到不同延时的功率谱;对不同延时的功率谱进行对齐并累加平均处理,得到回波功率谱;根据回波功率谱和单频连续激光的频率计算得到待测大气区域的大气风速。能够满足远距离复杂风场精细探测的现实需求。
Description
技术领域
本发明属于激光雷达测风技术领域,涉及一种二元脉冲幅度编码测风激光雷达风速测量方法及系统。
背景技术
复杂风场是指时间和空间上变化较为剧烈的风场,包括风切变和飞行器扰动场等中小尺度气流场,其全天候探测是研究复杂目标动力学规律的关键问题,也是空间信息获取与气象探测的交叉研究方向,在军民用航空安全、环境监测和重要武器发射保障等方面都具有重要应用需求。而当前测风的主要手段包括微波气象雷达、测风激光雷达以及探空气球等,其中测风激光雷达受益于激光的高单色性、高相干性及高方向性有望实现最高的风场距离分辨率,解决复杂风场的精细探测问题。
测风激光雷达的主要原理是发射的激光载波与待测大气区域的空气中气溶胶粒子相互作用后发生多普勒频移,通过频谱分析检测散射回波的多普勒频移得到风速,间接实现风场信息的感知。然而,在实现本发明的过程中,发明人发现目前传统的测风激光雷达所使用的测量方法,仍难以满足远距离复杂风场精细探测的现实需求。
发明内容
针对上述传统方法中存在的问题,本发明提出了一种二元脉冲幅度编码测风激光雷达风速测量方法以及一种测风激光雷达系统,能够满足远距离复杂风场精细探测的现实需求。
为了实现上述目的,本发明实施例采用以下技术方案:
一方面,提供一种二元脉冲幅度编码测风激光雷达风速测量方法,包括步骤:
根据二元脉冲幅度编码矩阵对单频连续激光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号;二元脉冲幅度编码矩阵为对角元素为0、其余元素为1、行数和列数相等的矩阵;
将二元脉冲幅度编码信号经光放大器放大后,通过环形器传输至光学天线并向待测大气区域发射;
接收待测大气区域的回波信号后,确定回波信号中回波脉冲序列的解码矩阵;
根据解码矩阵对回波脉冲序列进行解码处理,得到不同延时的功率谱;
对不同延时的功率谱进行对齐并累加平均处理,得到回波功率谱;
根据回波功率谱和单频连续激光的频率计算得到待测大气区域的大气风速。
在其中一个实施例中,确定回波信号中回波脉冲序列的解码矩阵的过程,包括:
根据回波信号中回波脉冲的码元强度确定解码矩阵;若码元强度为0,则解码矩阵中对应的元素记为-1,若码元强度为1,则解码矩阵中对应的元素记为1。
在其中一个实施例中,解码处理的解码公式为:
另一方面,还提供一种测风激光雷达系统,包括:
单频连续激光器,用于产生单频连续激光;
第一光纤分束器,用于将单频连续激光分束为本征参考光和信号光;
光脉冲调制器,用于对信号光进行调制输出;
信号发生器,用于产生编码脉冲信号至光脉冲调制器,控制光脉冲调制器根据二元脉冲幅度编码矩阵对信号光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号;
环形器,用于将二元脉冲幅度编码信号传输至光学天线以及用于回传光学天线接收的回波信号;
光学天线,用于将二元脉冲幅度编码信号向待测大气区域发射并接收待测大气区域的回波信号;
第二光纤分束器,用于将本征参考光及回波信号合束;
平衡探测器,用于合束并同时输入的本征参考光及回波信号以平方律检波的方式进行光电转换处理;
数字采集处理器,用于对平衡探测器输出的电信号进行采集,确定回波信号中回波脉冲序列的解码矩阵,根据解码矩阵对回波脉冲序列进行解码处理,得到不同延时的功率谱后,对不同延时的功率谱进行对齐并累加平均处理,得到回波功率谱,根据回波功率谱和单频连续激光的频率计算得到待测大气区域的大气风速。
在其中一个实施例中,数字采集处理器根据回波信号中回波脉冲的码元强度确定解码矩阵;若码元强度为0,则解码矩阵中对应的元素记为-1,若码元强度为1,则解码矩阵中对应的元素记为1。
在其中一个实施例中,解码处理的解码公式为:
在其中一个实施例中,光脉冲调制器包括声光调制器和/或电光调制器。
在其中一个实施例中,信号发生器包括任意波形发生器、晶体管逻辑电平信号发生器或射频信号发生器。
在其中一个实施例中,测风激光雷达系统还包括光放大器,光放大器用于将二元脉冲幅度编码信号放大后输入环形器。
上述技术方案中的一个技术方案具有如下优点和有益效果:
上述二元脉冲幅度编码测风激光雷达风速测量方法及系统,通过利用二元脉冲幅度编码实现单次脉冲发射至待测大气区域后获得多个位置的回波后,其中每个回波的距离分辨率对应单个码元宽度,再通过解码算出单个位置的回波,通过对齐回波位置或延迟回波时间后再累加,即可实现单个位置的多次测量,增强单个位置的测量次数,也即增强了信噪比并且在不影响距离分辨率的基础上实现了输出功率的提高。减小码元宽度可以提高距离分辨率,而编码序列中码元数量增加又可以增大发射能量,增强信噪比,扩大探测范围,能够很好地满足远距离复杂风场精细探测的现实需求,为复杂风场的精细探测提供了新的有效技术手段。
附图说明
为了更清楚地说明本申请实施例或传统技术中的技术方案,下面将对实施例或传统技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为一个实施例中二元脉冲幅度编码测风激光雷达风速测量方法的流程示意图;
图2为一个实施例中基于编码脉冲进行大气探测的原理示意图;
图3为一个实施例中基于二元脉冲幅度编码的编码脉冲示意图;
图4为一个实施例中测风激光雷达系统的结构示意图;
图5为另一个实施例中测风激光雷达系统的结构示意图;
图6为一个实施例中二元脉冲幅度编码测风激光雷达风速测量方法的应用流程示意图;
图7为一个实施例中基于二元脉冲幅度编码的测风激光雷达系统与非编码脉冲测风激光雷达得到的速度-距离探测结果对比示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
除非另有定义,本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。在本申请的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本申请。
需要说明的是,在本文中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本发明的至少一个实施例中。在说明书中的各个位置展示该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。
本领域技术人员可以理解,本文所描述的实施例可以与其它实施例相结合。在本发明说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。
当前的测风激光雷达一般采用脉冲体制,结合飞行时间法判断探测距离,其距离分辨率为:
其中,为发射脉冲的宽度,/>为采样间隔,/>为采样点数,/>为光速。根据信号处理理论,进行频谱分析的采样点数/>越多,频谱分析的速度分辨率越高,但是距离分辨率就越差。在采样点数和采样间隔固定的情况下,要提高距离分辨率就要降低发射脉冲宽度。但是受限于非线性效应,激光脉冲的峰值功率难以提高,当峰值功率大于受激布里渊等非线性效应的阈值时,将会使得输出的光束展宽、严重抖动甚至损坏激光器。
因此脉宽越窄,可发射的单脉冲能量越小,回波的能量也越低,信噪比越低,探测范围越小。而且由于测风激光雷达的最大探测距离为:,如果通过增大发射脉冲的重复频率/>或者减小发射脉冲的重复时间/>来积累更多脉冲提升出射能量,则最大探测距离/>减小。
可见,传统的测风激光雷达受限于脉冲峰值功率的限制,无法兼顾高距离分辨率和远探测距离,也就无法满足远距离复杂风场精细探测的需求。为此,本申请提供了一种新的技术方案,设计原理是发射一组二元编码脉冲序列,以实现单次脉冲发射获得多个位置的回波,再对回波解码提取出单个位置的回波,兼顾了激光雷达的高距离分辨率与大功率输出,能够解决传统单频激光雷达距离分辨和输出功率无法兼顾的问题,在增加脉宽的同时提高输出功率但不降低距离分辨率。
下面将结合本发明实施例图中的附图,对本发明实施方式进行详细说明。
请参阅图1,在一个实施例中,本申请实施例提供了一种二元脉冲幅度编码测风激光雷达风速测量方法,包括如下处理步骤S12至S22:
S12,根据二元脉冲幅度编码矩阵对单频连续激光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号;二元脉冲幅度编码矩阵为对角元素为0、其余元素为1、行数和列数相等的矩阵;
S14,将二元脉冲幅度编码信号经光放大器放大后,通过环形器传输至光学天线并向待测大气区域发射;
S16,接收待测大气区域的回波信号后,确定回波信号中回波脉冲序列的解码矩阵;
S18,根据解码矩阵对回波脉冲序列进行解码处理,得到不同延时的功率谱;
S20,对不同延时的功率谱进行对齐并累加平均处理,得到回波功率谱;
S22,根据回波功率谱和单频连续激光的频率计算得到待测大气区域的大气风速。
可以理解,如图2和图3所示,分别为基于二元脉冲幅度编码的编码脉冲示意图和该编码脉冲进行大气探测的原理图,在实际应用中,可以通过激光调制器对单频连续激光进行二元脉冲幅度编码调制,调制过程中使用的二元脉冲幅度编码矩阵可以表示为:
二元脉冲幅度编码矩阵是一个对角元素为0,其余元素为1,行数和列数相等的矩阵。其中,码元0和1表示信号的有无,即码元0表示信号(单频连续激光)被斩断强度为0,码元1表示有信号强度不为0。其中,的每一行控制一个发射脉冲内每个码元宽度所对应信号的有无。
具体的,假设为m行m列矩阵,m为码元长度,码长代表编码的位数,则每个脉冲可被等分为m段,对于第一行序列/>控制的脉冲幅度即为:第1段无信号,其后m-1段都有信号。相邻的下一个脉冲幅度则由第二行序列控制,即除第2段没有信号外,其余m-1段都有信号。依此类推,每一个脉冲幅度由一行编码序列控制,直至第m行序列结束,再回到/>的第一行序列继续控制接下来的脉冲幅度。
将调制后输出的二元脉冲幅度编码信号经过光放大器放大、环形器传输至光学天线发射到大气中的待测大气区域,并对待测大气区域的回波信号进行接收,对所接收的回波信号进行解码处理。
解码处理过程可以为:先将采集到的回波信号中的回波脉冲按照编码序列分组,即属于二元脉冲幅度编码矩阵中一行序列的回波脉冲分为一组;然后再对每一组回波脉冲,按照系统采用的距离门分段求出累加平均功率谱/>,其分别对应第/>组回波的累加平均功率谱。其中,每个累加平均功率谱都是多个距离分辨率为码元宽度的功率谱的累加结果。按照解码矩阵解码获得不同延时的功率谱,然后对齐并累加平均最终得到距离分辨率为单个码元宽度的回波功率谱。
上述二元脉冲幅度编码测风激光雷达风速测量方法,通过利用二元脉冲幅度编码实现单次脉冲发射至待测大气区域后获得多个位置的回波后,其中每个回波的距离分辨率对应单个码元宽度,再通过解码算出单个位置的回波,通过对齐回波位置或延迟回波时间后再累加,即可实现单个位置的多次测量,增强单个位置的测量次数,也即增强了信噪比并且在不影响距离分辨率的基础上实现了输出功率的提高。减小码元宽度可以提高距离分辨率,而编码序列中码元数量增加又可以增大发射能量,增强信噪比,扩大探测范围,能够很好地满足远距离复杂风场精细探测的现实需求,为复杂风场的精细探测提供了新的有效技术手段。
在一个实施例中,进一步的,关于上述的步骤S16中确定回波信号中回波脉冲序列的解码矩阵的过程,具体可以包括如下处理步骤:
根据回波信号中回波脉冲的码元强度确定解码矩阵;若码元强度为0,则解码矩阵中对应的元素记为-1,若码元强度为1,则解码矩阵中对应的元素记为1。
可以理解,解码矩阵的确定只需要根据回波脉冲的形状即可,这是因为编码矩阵已知且固定,只要通过回波脉冲的形状直接判断出收集的第一个回波脉冲对应编码序列的第几个,即可直接写出解码矩阵,实现解码矩阵的高效且准确获取。如图3所示,以4个码元的编码序列为例(其他码元长度的编码序列同理理解),回波脉冲中(解码前的编码脉冲)的4个码元,第一个码元的码元强度为零、其后码元的码元强度为1,则其编码序列就写为0111。在解码时则是类似的,回波脉冲中哪一个码元所在位置处其码元强度为0,则解码序列中相应位置写成-1,其余码元强度为1的码元则均写成1,因此编码序列为0111的,在解码时就可以写成-1111,如此即可确定解码矩阵中各相应元素对应的取值,从而确定解码矩阵。前述过程也称根据回波脉冲的形状写出解码矩阵的过程,其他码元长度的编码序列情形可同理理解。
在一个实施例中,进一步的,解码处理的解码公式为:
其中,对应延时分别为/>的距离分辨率为单个码元宽度的功率谱,/>为单个码元宽度,/>为解码矩阵,/>分别对应第1,2,…,m组回波的累加平均功率谱。可以理解,码元宽度相当于是一个脉冲的宽度,可以记为τ,其也决定了距离分辨率。解码就是为了获得距离分辨率为单个码元宽度的功率谱,没解码前是每个分组后的序列有一个平均功率谱,是多个不同延时的距离分辨率为码元宽度的功率谱的线性叠加。
具体的,假设获得的解码矩阵为:
需要说明的是,此处列举的解码矩阵示例中,起始第一行为,表示采集分组后的第一组回波脉冲对应着二元脉冲幅度编码矩阵/>中的第二个序列,如此,按前述解码公式解码出不同延时的功率谱且每个功率谱的距离分辨率对应一个码元宽度,实现对接收的回波高效解码输出。
在解码后,可以将解码后不同延时的功率谱进行对齐并累加平均。具体的,将解码后获得的不同延时的功率谱:
分别延时后,进行累加并求平均即可获得估计后的回波功率谱,且回波功率谱的距离分辨率为一个码元宽度。获得回波功率谱后,即可根据本领域既有的大气风速计算公式,利用回波功率谱和单频连续激光的频率来计算得到该待测大气区域的大气风速,进而确定该待测大气区域的风场。
应该理解的是,虽然图1流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且图1的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
请参阅图4,在一个实施例中,提供一种测风激光雷达系统100,包括单频连续激光器11、第一光纤分束器12、光脉冲调制器13、信号发生器14、环形器15、光学天线16、第二光纤分束器17、平衡探测器18和数字采集处理器19。其中,单频连续激光器11用于产生单频连续激光。第一光纤分束器12用于将单频连续激光分束为本征参考光和信号光。光脉冲调制器13用于对信号光进行调制输出。信号发生器14用于产生编码脉冲信号至光脉冲调制器,控制光脉冲调制器根据二元脉冲幅度编码矩阵对信号光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号。环形器15用于将二元脉冲幅度编码信号传输至光学天线以及用于回传光学天线接收的回波信号。光学天线16用于将二元脉冲幅度编码信号向待测大气区域发射并接收待测大气区域的回波信号。
第二光纤分束器17用于将本征参考光及回波信号合束。平衡探测器18用于合束并同时输入的本征参考光及回波信号以平方律检波的方式进行光电转换处理。数字采集处理器19用于对平衡探测器输出的电信号进行采集,确定回波信号中回波脉冲序列的解码矩阵,根据解码矩阵对回波脉冲序列进行解码处理,得到不同延时的功率谱后,对不同延时的功率谱进行对齐并累加平均处理,得到回波功率谱,根据回波功率谱和单频连续激光的频率计算得到待测大气区域的大气风速。
可以理解,上述测风激光雷达系统100中的各部件均可以采用本领域既有的相应部件来搭建该系统,其中,光脉冲调制器13和信号发生器14改为使用上述二元脉冲幅度编码测风激光雷达风速测量方法中的二元脉冲幅度编码调制方式对单频连续激光进行调制输出。而数字采集处理器19则相应地采用上述二元脉冲幅度编码测风激光雷达风速测量方法中的解码处理方式进行解码和计算处理。
上述测风激光雷达系统100,通过利用二元脉冲幅度编码实现单次脉冲发射至待测大气区域后获得多个位置的回波后,其中每个回波的距离分辨率对应单个码元宽度,再通过解码算出单个位置的回波,通过对齐回波位置或延迟回波时间后再累加,即可实现单个位置的多次测量,增强单个位置的测量次数,也即增强了信噪比并且在不影响距离分辨率的基础上实现了输出功率的提高。减小码元宽度可以提高距离分辨率,而编码序列中码元数量增加又可以增大发射能量,增强信噪比,扩大探测范围,能够很好地满足远距离复杂风场精细探测的现实需求,为复杂风场的精细探测提供了新的有效技术手段。
在一个实施例中,进一步的,数字采集处理器19根据回波信号中回波脉冲的码元强度确定解码矩阵;若码元强度为0,则解码矩阵中对应的元素记为-1,若码元强度为1,则解码矩阵中对应的元素记为1。
在一个实施例中,进一步的,解码处理的解码公式为:
在一个实施例中,进一步的,测风激光雷达系统100还包括光放大器20,光放大器20用于将二元脉冲幅度编码信号放大后输入环形器15。
可以理解,在本实施例中,还可以采用光放大器20对二元脉冲幅度编码信号放大后再输出下级单元,从而提高探测能力。具体的,如图5所示,单频连续激光器11的输出端与第一光纤分束器12的输入端连接,第一光纤分束器12将单频连续激光分成两路,其中第一光纤分束器12的第一输出端A1输出信号光,第一光纤分束器12的第二输出端B1输出本征参考光。
第一光纤分束器12的第一输出端A1与光脉冲调制器13的输入端连接,第一光纤分束器12的第二输出端B1和第二光纤分束器17的第一输入端B3连接。信号发生器14的输出端与光脉冲调制器13的驱动信号输入端相连,光脉冲调制器13的输出端和光放大器20的输入端连接,光放大器20的输出端和环形器15的第一输出端A2连接,环形器15的第一输出端A3和光学天线18的输入端连接,光学天线16出射的光束照在待测大气区域中的目标上,目标的后向散射回波经过光学天线16收集,再回到环形器15的第一输出端A3,在环形器15内部传输到第二输出端A4和第二光纤分束器17的第二输入端B2连接。第二光纤分束器17的输出端和平衡探测器18的输入端连接,平衡探测器18的输出端和数字采集处理器19的输入端连接。待测大气区域中的目标为流体粒子。
在一个实施例中,进一步的,光脉冲调制器13包括声光调制器和/或电光调制器。可以理解,本实施例中,光脉冲调制器13可以采用本领域中既有的声光调制器或者电光调制器,也可以采用声光调制器和电光调制器共同组成的组合调制器,以实现所需的调制功能。
在一个实施例中,进一步的,信号发生器14包括任意波形发生器、晶体管逻辑电平信号发生器或射频信号发生器。可以理解,本实施例中,信号发生器14可以是本领域中既有的任意波形发生器,也可以是本领域中既有的晶体管逻辑电平信号发生器,或者还可以是本领域中既有的射频信号发生器,以产生所需的编码脉冲信号。
在一个实施例中,进一步的,编码脉冲信号的码元长度为4 。可以理解,码长长度为4,码元宽度为一个单个码(0或者1)的宽度。本实施例中选用的码元长度为4 ,可以实现较为简便的风速测量效果。
关于测风激光雷达系统100信号处理过程的具体限定,可以参见上文中二元脉冲幅度编码测风激光雷达风速测量方法的相应限定,在此不再赘述。可以理解,上述测风激光雷达系统100的光路中,还可以包括上述实施例未列出的其他组成部件,具体可以根据不同应用场景中具体应用的测风激光雷达系统100的结构组成确定。
在一个实施例中,为了更好地理解本申请的技术方案,本实施例中提供了其中一种详细的实施示例。需要说明的是,该实施示例并非是对本申请的唯一限定,本领域技术人员可以根据本申请的技术方案的构思实现不同场景下的风速测量。
如图6所示,其中表示对采集到的回波,按照先傅里叶变换再平方的方式求出功率谱,以码元长度为4 举例,二元脉冲幅度编码后的激光脉冲发射到大气中,再接收其回波,首先对回波按照编码序列进行分组;然后对每组回波脉冲按照距离门利用傅里叶变换后平方的方法计算其累加平均功率谱;再通过观察回波脉冲部分的形状写出解码矩阵,例如:
若第一组回波的脉冲部分在第二个码元处的强度为0,则解码矩阵为:
即第二个码元为-1,其余都为1,再通过解码公式:
计算出对应不同延时的单个码元宽度的功率谱;最后通过延时处理,对齐每个解码出的功率谱并累加求平均。经过解码,可以获得距离分辨率为单个码元宽度的回波功率谱,但是累加次数大幅提高,发射能量增多,相比于传统发射M次非编码脉冲再累加平均求出的回波功率谱,经过四个码元的二元脉冲幅度编码调制,可将信噪比提升为原来的3/2倍。
本示例中提供的采用基于二元脉冲幅度编码的测风激光雷达系统,与非编码脉冲测风激光雷达得到的速度-距离探测结果对比示意图如图7所示:
使用本申请提出的技术方案测量得到的风速具有更高的信噪比,比以码元宽度发射的使用同样数量的脉冲进行累加的非编码脉冲测风激光雷达能探测更远的距离。从示意图7中可以看到,虚线对应的非编码脉冲在探测距离达到/>之后,由于信噪比更低会受到噪声的干扰导致测量的风速错误。而本申请提出的技术方案可以在更远的测量范围处维持稳定的风速估计。编码和非编码的测量范围差距受限于编码的码元数量M,M越大,编码测风激光雷达的测量范围就提升越大。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可做出若干变形和改进,都属于本申请保护范围。因此本申请专利的保护范围应以所附权利要求为准。
Claims (7)
1.一种二元脉冲幅度编码测风激光雷达风速测量方法,其特征在于,包括步骤:
根据二元脉冲幅度编码矩阵对单频连续激光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号;所述二元脉冲幅度编码矩阵为对角元素为0、其余元素为1、行数和列数相等的矩阵;
将所述二元脉冲幅度编码信号经光放大器放大后,通过环形器传输至光学天线并向待测大气区域发射;
接收所述待测大气区域的回波信号后,确定所述回波信号中回波脉冲序列的解码矩阵;
根据所述解码矩阵对所述回波脉冲序列进行解码处理,得到不同延时的功率谱;
对不同延时的所述功率谱进行对齐并累加平均处理,得到回波功率谱;
根据所述回波功率谱和所述单频连续激光的频率计算得到所述待测大气区域的大气风速;
其中,确定所述回波信号中回波脉冲序列的解码矩阵的过程,包括:
根据所述回波信号中回波脉冲的码元强度确定所述解码矩阵;若所述码元强度为0,则所述解码矩阵中对应的元素记为-1,若所述码元强度为1,则所述解码矩阵中对应的元素记为1。
3.一种测风激光雷达系统,其特征在于,包括:
单频连续激光器,用于产生单频连续激光;
第一光纤分束器,用于将所述单频连续激光分束为本征参考光和信号光;
光脉冲调制器,用于对所述信号光进行调制输出;
信号发生器,用于产生编码脉冲信号至所述光脉冲调制器,控制所述光脉冲调制器根据二元脉冲幅度编码矩阵对所述信号光进行二元脉冲幅度编码调制,得到二元脉冲幅度编码信号;所述二元脉冲幅度编码矩阵为对角元素为0、其余元素为1、行数和列数相等的矩阵;
环形器,用于将所述二元脉冲幅度编码信号传输至光学天线以及用于回传所述光学天线接收的回波信号;
光学天线,用于将所述二元脉冲幅度编码信号向待测大气区域发射并接收所述待测大气区域的回波信号;
第二光纤分束器,用于将所述本征参考光及所述回波信号合束;
平衡探测器,用于合束并同时输入的所述本征参考光及所述回波信号以平方律检波的方式进行光电转换处理;
数字采集处理器,用于对所述平衡探测器输出的电信号进行采集,确定所述回波信号中回波脉冲序列的解码矩阵,根据所述解码矩阵对所述回波脉冲序列进行解码处理,得到不同延时的功率谱后,对不同延时的所述功率谱进行对齐并累加平均处理,得到回波功率谱,根据所述回波功率谱和所述单频连续激光的频率计算得到所述待测大气区域的大气风速;其中,
所述数字采集处理器根据所述回波信号中回波脉冲的码元强度确定所述解码矩阵;若所述码元强度为0,则所述解码矩阵中对应的元素记为-1,若所述码元强度为1,则所述解码矩阵中对应的元素记为1。
5.根据权利要求4所述的测风激光雷达系统,其特征在于,所述光脉冲调制器包括声光调制器和/或电光调制器。
6.根据权利要求4所述的测风激光雷达系统,其特征在于,所述信号发生器包括任意波形发生器、晶体管逻辑电平信号发生器或射频信号发生器。
7.根据权利要求4所述的测风激光雷达系统,其特征在于,还包括光放大器,所述光放大器用于将所述二元脉冲幅度编码信号放大后输入所述环形器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310246035.2A CN115951332B (zh) | 2023-03-15 | 2023-03-15 | 二元脉冲幅度编码测风激光雷达风速测量方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310246035.2A CN115951332B (zh) | 2023-03-15 | 2023-03-15 | 二元脉冲幅度编码测风激光雷达风速测量方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115951332A CN115951332A (zh) | 2023-04-11 |
CN115951332B true CN115951332B (zh) | 2023-05-23 |
Family
ID=85903372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310246035.2A Active CN115951332B (zh) | 2023-03-15 | 2023-03-15 | 二元脉冲幅度编码测风激光雷达风速测量方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115951332B (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003100458A1 (fr) * | 2002-05-29 | 2003-12-04 | Mitsubishi Denki Kabushiki Kaisha | Dispositif radar doppler a laser |
DE102017115710A1 (de) * | 2017-07-12 | 2019-02-07 | Airbus Defence and Space GmbH | LIDAR-Anordnung und LIDAR-Verfahren |
CN108594256B (zh) * | 2018-04-16 | 2021-10-12 | 夏和娣 | 一种基于脉冲编码技术的相干激光雷达 |
CN109323751A (zh) * | 2018-11-14 | 2019-02-12 | 四川鸿禾阳科技有限公司 | 一种脉冲编码的分布式光纤振动传感方法及装置 |
CN110749878B (zh) * | 2019-09-24 | 2021-08-13 | 北京理工大学 | 一种基于数字信道化的脉内频率编码信号参数测量方法 |
CN114325726B (zh) * | 2020-09-27 | 2022-12-13 | 北京一径科技有限公司 | 激光雷达的抗串扰处理方法、存储介质及激光雷达 |
CN112034447B (zh) * | 2020-11-06 | 2021-01-15 | 中国人民解放军国防科技大学 | 一种提高脉冲多普勒雷达检测性能的二相调制方法及装置 |
CN113253301B (zh) * | 2021-07-05 | 2021-10-01 | 中国人民解放军国防科技大学 | 多频脉冲激光雷达信号处理方法及测风雷达系统 |
-
2023
- 2023-03-15 CN CN202310246035.2A patent/CN115951332B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115951332A (zh) | 2023-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108594256B (zh) | 一种基于脉冲编码技术的相干激光雷达 | |
CN104931973B (zh) | 用于激光雷达系统的非对称相位编码测距方法 | |
CN108415010B (zh) | 一种基于梯形lfmcw调制的雷达多目标检测方法 | |
US20130278455A1 (en) | Efficient pulse doppler radar with no blind ranges, range ambiguities, blind speeds, or doppler ambiguities | |
Kurdzo et al. | Optimized NLFM pulse compression waveforms for high-sensitivity radar observations | |
CN110361715B (zh) | 一种脉冲编码装置、编码方法及激光雷达系统 | |
Scotti et al. | In-field experiments of the first photonics-based software-defined coherent radar | |
Ipanov et al. | Radar signals with ZACZ based on pairs of D-code sequences and their compression algorithm | |
CN110646779A (zh) | 相位编码不饱和调制方法、装置、激光雷达测距测速方法及激光雷达系统 | |
JP2019105601A (ja) | レーダ装置及びそのレーダ信号処理方法 | |
CN101738602B (zh) | 伪随机序列在探地雷达应用中回波数据预处理方法 | |
CN103064065A (zh) | 双周期夹断式伪随机码的波形设计和对回波的信号处理方法 | |
Wan et al. | Ground moving target focusing and motion parameter estimation method via MSOKT for synthetic aperture radar | |
CN115951332B (zh) | 二元脉冲幅度编码测风激光雷达风速测量方法及系统 | |
CN113253301B (zh) | 多频脉冲激光雷达信号处理方法及测风雷达系统 | |
Akhtar et al. | Formation of range-doppler maps based on sparse reconstruction | |
CN115236697B (zh) | 分时多频脉冲测风激光雷达系统和风速测量方法 | |
CN211577419U (zh) | 相位编码不饱和调制装置和激光雷达系统 | |
Beltrao et al. | Subpulse processing for long range surveillance noise radars | |
CN117008147A (zh) | 一种基于组合编码的测距测速及成像方法和装置 | |
Kang et al. | A study on pulsed-LFM and pulsed-NLFM waveforms for radar systems | |
CN116165664A (zh) | 一种基于二维波形编码的星载sar距离模糊抑制方法 | |
US20220206127A1 (en) | Phase-coded unsaturated modulation method, apparatus, lidar distance and velocity measurement method and lidar system | |
Xu et al. | A time-division multifrequency wind LiDAR with high spatial resolution and extended detection range | |
CN109581363B (zh) | 一种基于非相干散射雷达的小尺寸空间碎片检测和参数提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |