基于遗传算法的动平台相干测风激光雷达的反演方法技术

一种基于遗传算法的动平台相干测风激光雷达的反演方法，该激光雷达包括脉冲激光器、分光器、耦合器、声光调制器、放大器、环形器、收发望远镜、透射型扫瞄镜、平衡探测器、A/D转换器、由惯导系统和全球定位系统组成的POS系统和上位机数据处理器。本发明专利技术基于遗传算法的动平台相干测风激光雷达反演方法包括以下步骤：确定遗传算法相关参数，包括风矢量的解域范围和累加的总方位数，然后确定种群规模、交叉概率、变异概率、最大进化代数，确定动平台条件下通过透射型扫描镜扫描所获得的频谱信息与校正后大气风场之间的表达式；对风速进行个体初始化；采用二进制编码方式将风速转化为基因进行编码；进行选择、交叉与变异操作；判断是否达到最大进化代数和是否收敛；若运动平台高度会发生变化，则利用POS系统所得动平台高度信息对高度进行校正，最终得到在地面坐标系下的三维风场。本发明专利技术的特点是不依赖于视向风速的计算，并利用遗传算法通过全局寻优来反演动平台测风雷达所测三维风场，提高了在低信噪比情况下的风场反演精度。下的风场反演精度。下的风场反演精度。

【技术实现步骤摘要】
基于遗传算法的动平台相干测风激光雷达的反演方法


[0001]本专利技术属于激光雷达数据处理领域，特别是一种基于遗传算法的动平台相干测风激光雷达的反演方法。

技术介绍

[0002]大气风场作为大气的一项重要参数，在气象学、天气预报、环境监测、航空领域等方面有重要应用，因此能够精确地实现对于三维风场的反演具有重要意义。相干多普勒测风激光雷达是研究大气风场的重要手段，与其他风场监测手段不同，多普勒测风激光雷达是利用大气环境中的气溶胶粒子和大气分子的后向散射信号来进行探测的。该探测方式具有信噪比高、探测灵敏度高、精度高和对背景光不敏感等优点。将相干测风激光雷达安装于动平台上进行探测，能够扩大雷达系统的探测范围。但目前常用的三维风场反演算法，如最小二乘法等在反演过程中依赖于视向风速的反演精度，当气溶胶浓度较低或风速较小，即信号较弱时，往往会影响到视向风速的反演精度，进一步对三维风场的反演结果准确性产生影响。

技术实现思路

[0003]本专利技术的目的在于克服现有动平台相干多普勒测风激光雷达反演过程中依赖于视向风速反演精度的问题，提供一种基于遗传算法的动平台相干测风激光雷达的反演方法，该方法不通过视向风速的求解来反演风场，而是直接通过雷达在平台运动过程中所获得的多方位频谱信息来获得三维风场。因此，该算法避免了在信噪比较低的情况下所计算出的视向风速不准确，进一步影响三维风场反演结果精确度的情况。直接利用多个不同方位的频谱之间的相关性，对不同扫描角频谱信息中的有用信号进行积累，同时对噪声信号进行平均，使得该算法在信噪比较低时有更好的反演精度。
[0004]本专利技术技术解决方案如下：
[0005]一种基于遗传算法的动平台相干测风激光雷达的反演方法，其特点在于，该方法包括以下步骤：
[0006]1)、建立相干多普勒测风激光雷达动平台，该相干多普勒测风激光雷达动平台包括脉冲激光器、分光器、耦合器、声光调制器、放大器、环形器、收发望远镜、透射型扫瞄镜、平衡探测器、A/D转换器、由惯导系统和全球定位系统组成的POS系统和上位机数据处理器，所述的脉冲激光器发出的光束经所述的分光器分为两路光束：其中一光束直接与所述的耦合器的1端口相连；另一光束经所述的声光调制器进行移频后通过光纤与放大器相连进行放大，之后经所述的环形器的1端口、2端口输出后输入所述的收发望远镜，利用透射型扫瞄镜将光发射到大气中，进行VAD(velocity azimuth display)扫描，即以固定仰角的方式进行不同方位角的扫描，回波信号通过同一收发望远镜(7)进行接收，接收后经过所述的环形器的2端口、3端后进入所述的耦合器的2端口，与进入所述的耦合器1端口的另一光束进行拍频分成两路经所述的平衡探测器)进行信号的探测，利用所述的A/D转换器进行信号采
集，将模拟信号转换为数字信号进入所述的上位机数据处理器；
[0007]2)、对于相干多普勒测风激光雷达动平台而言，考虑到平台运动的影响，再根据激光雷达几何关系我们可以得到每个扫描方向的风速为：
[0008][0009]其中，V＝[V
x V
y V
z
]为三维风场矢量，i为透射型扫瞄镜(8)所扫描的当前方位数，V
LOSi
为不同扫描方向的视向风速；α
i
为系统的视向仰角，β
i
为视向方位角。V
xiA
、V
yiA
、V
ziA
分别为不同方位探测时运动平台的北向、东向和地向速度；α
Gi
和β
Gi
为大地坐标系下的视向仰角和方位角；通过雷达系统中由惯导系统和全球定位系统组成的POS系统(11)，我们可以直接获得平台的运动速度及在动平台坐标系下的视向仰角α
P
和方位角β
P
，需要利用坐标变换将其转化为地面坐标系下的α
G
和β
G
；根据相干多普勒测风原理，我们可以得到视向风速与风场引起的多普勒频移之间的关系为：
[0010][0011]式中，f
i
为不同扫描方向所获信号的多普勒频移值；f
IF
为系统中频，λ为系统所用光源波长；之后，我们便可以按下式得到当前风场下所有方位多普勒频移对应的频谱强度累加和的大小S(V)：
[0012][0013]式中，W
i
(f
i
)为以当前假设风场计算出的单方位多普勒频移所对应的频谱强度；N为参与计算的总方位数；当S(V)最大时所对应的V的解值便为当前风场信息；
[0014]3)、确定遗传算法的相关参数，包括风矢量的求解范围和累加的总方位数，一般情况下水平风速求解范围在
±
20m/s，垂直风速求解范围在
±
3m/s，累加的总方位数则根据扫瞄镜的扫描方位数确定。然后确定种群规模、交叉概率、变异概率、最大进化代数。种群规模不宜过大也不宜过小，一般选取0
‑
100，交叉概率和变异概率的选取一般受种群大小、编码长度等因素的影响，交叉概率取值一般为0.4
‑
0.99，变异概率一般为0.001－0.2。最大进化代数则取决于所希望的循环次数，最大进化代数越大越接近最优解，但相应的计算时间也会延长；
[0015]4)、对风场进行初始化，确定初始风场V的值：第一距离门的风场采用最小二乘法风场反演结果，之后其余距离门以上一距离门求解风速为中心，在小范围内随机选择风速进行初始化；若上一距离门中没有求解出风速信息并且最小二乘法也无法提取风速信息，则采用随机初始化，即在给定风场求解范围随机确定风速；
[0016]5)、然后采用二进制编码方式，将风速转化为基因进行编码；
[0017]6)、将风矢量转换为二进制数的基因序列后进行选择操作、交叉操作与变异操作：选择操作为利用风速与频移量之间的关系计算累加频谱强度S(V)的大小，按照S(V)的大小从小到大对风速进行排序，将排序后的风速平均分为前、中、后3段,分别按照0.6，0.8，1的比例对每段的风速进行选择，采用最后段的风速对被淘汰的风速进行补充：交叉操作为将两个相邻风速间的二进制数的基因序列根据交叉概率进行交叉互换形成新的风速；变异操作为按变异概率随机地改变风速的单个基因值，防止整个算法在求解时陷入局部最优；
[0018]7)、最后判断是否达到步骤2)中所确定的最大进化代数，即算法是否完成循环：
[0019]如达到最大进化代数，则进入下一步；未达到最大进化代数则再返回步骤6)再次进行选择、交叉与变异操作，如达到最大进化代数，则进入下一步；未达到最大进化代数则再返回步骤6)再次进行选择、交叉与变异操作，
[0020]8)根据S(V)所对应最优解V的保持代数判断算法是否收敛：
[0021]若未收敛且执行算法次数少于3次，则返回步骤4)对风场进行初始化；
[0022]若执行算法3次均未收敛，则无法求解本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于遗传算法的动平台相干测风激光雷达的反演方法，其特征在于，该方法包括以下步骤：1)、建立相干多普勒测风激光雷达动平台，该相干多普勒测风激光雷达动平台包括脉冲激光器(1)、分光器(2)、耦合器(3)、声光调制器(4)、放大器(5)、环形器(6)、收发望远镜(7)、透射型扫瞄镜(8)、平衡探测器(9)、A/D转换器(10)、由惯导系统和全球定位系统组成的POS系统(11)和上位机数据处理器(12)，所述的脉冲激光器(1)发出的光束经所述的分光器(2)分为两路光束：其中一光束直接与所述的耦合器(3)的1端口相连；另一光束经所述的声光调制器(4)进行移频后通过光纤与放大器(5)相连进行放大，之后经所述的环形器(6)的1端口、2端口输出后输入所述的收发望远镜(7)，利用透射型扫瞄镜(8)将光发射到大气中，进行VAD(velocity azimuth display)扫描，即以固定仰角的方式进行不同方位角的扫描，回波信号通过同一收发望远镜(7)进行接收，接收后经过所述的环形器(6)的2端口、3端后进入所述的耦合器(3)的2端口，与进入所述的耦合器(3)1端口的另一光束进行拍频分成两路经所述的平衡探测器(9)进行信号的探测，利用所述的A/D转换器(10)进行信号采集，将模拟信号转换为数字信号进入所述的上位机数据处理器(12)；2)、对于相干多普勒测风激光雷达动平台而言，考虑到平台运动的影响，再根据激光雷达几何关系我们可以得到每个扫描方向的风速为：其中，V＝[V
x V
y V
z
]为三维风场矢量，i为透射型扫瞄镜(8)所扫描的当前方位数，V
LOSi
为不同扫描方向的视向风速；α
i
为系统的视向仰角，β
i
为视向方位角。V
xiA
、V
yiA
、V
ziA
分别为不同方位探测时运动平台的北向、东向和地向速度；α
Gi
和β
Gi
为大地坐标系下的视向仰角和方位角；通过雷达系统中由惯导系统和全球定位系统组成的POS系统(11)，我们可以直接获得平台的运动速度及在动平台坐标系下的视向仰角α
P
和方位角β
P
，需要利用坐标变换将其转化为地面坐标系下的α
G
和β
G
；根据相干多普勒测风原理，我们可以得到视向风速与风场引起的多普勒频移之间的关系为：式中，f
i
为不同扫描方向所获信号的多普勒频移值；f
IF
为系统中频，λ为系统所用光源波长；之后，我们便可以按下式得到当前风场下所有方位多普勒频移对应的频谱强度累加和的大小S(V)：式中，W
i
(f
i
)为以当前假设风场计算出的单方位多普勒频移所对应的频谱强度；N...

【专利技术属性】
技术研发人员：赵源，竹孝鹏，刘继桥，朱小磊，陈卫标，
申请(专利权)人：中国科学院上海光学精密机械研究所，
类型：发明
国别省市：