CN114325666A - 高速运动下的相干外差测风雷达修正方法 - Google Patents

高速运动下的相干外差测风雷达修正方法 Download PDF

Info

Publication number
CN114325666A
CN114325666A CN202111681292.6A CN202111681292A CN114325666A CN 114325666 A CN114325666 A CN 114325666A CN 202111681292 A CN202111681292 A CN 202111681292A CN 114325666 A CN114325666 A CN 114325666A
Authority
CN
China
Prior art keywords
radar
axis
coordinate system
circulator
wind
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
CN202111681292.6A
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.)
Sichuan Zhizhou Technology Co ltd
Original Assignee
Sichuan Zhizhou Technology Co ltd
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 Sichuan Zhizhou Technology Co ltd filed Critical Sichuan Zhizhou Technology Co ltd
Priority to CN202111681292.6A priority Critical patent/CN114325666A/zh
Publication of CN114325666A publication Critical patent/CN114325666A/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

高速运动下的相干外差测风雷达修正方法
技术领域
本发明涉及气象雷达探测技术领域,尤其涉及高速运动下的相干外差测风雷达修正方法。
背景技术
目前现有的相干外差测风雷达作为一种大气风场的传感器,多为地基式,即相对地面静态环境下探测大气中的风速与风向,若应用在高速运动平台,该传感器测得的风场速度矢量为气溶胶粒子相对雷达光轴的运动参数,是由雷达本身相对地面的高速运动及气溶胶粒子相对地面的运动(即大气风场)耦合的结果,因此,不可直接搭载在高速运动的平台上获得大气风场。
发明内容
针对上述问题,本发明所提供的方法可修正雷达本身相对地面的运动,获得真实的大气风场,即北东天坐标下的风场。
本发明通过以下技术方案实现:
高速运动下的相干外差测风雷达修正方法,雷达在高速运动的载体平台上进行数据采集,采集的数据通过运动修正算法得到大气的三维风场信息;
所述运动修正算法包括以下步骤:
S1、建立3组模型坐标系,并根据3组模型坐标系得到相关数据,其中,相关数据包括2组模型坐标系的旋转角度以及雷达本振光R发射产生的四束光R1、R2、R3、R4的光轴坐标;
S2、根据S1的旋转角度得到每个旋转角度的矢量旋转表达式以及单位矢量的旋转表达式;
S3、根据S1所述的四束光R1、R2、R3、R4的光轴位置坐标得到本振光R的光轴旋转后在直角坐标系的表达式R’
S4、根据S3得到的R’与真实风速矢量以及雷达速度矢量计算出中四束光的各光轴的相对速度标量;
S5、采用最小二乘法解出S4中的相对速度标量的参数值,得到大气的三维风场信息,所述三维风场信息包括垂直风速和水平风速。
进一步的,步骤S1所述的3组模型坐标系为地理坐标系、载体坐标系和雷达坐标系,所述地理坐标系O-xgygzg采用东北天坐标系,其中原点为载体重心,xg轴指向东,yg轴指向正北,zg轴指向天,所述载体坐标系O-xbybzb的原点为载体中心,xb轴沿载体横轴向右,yb轴沿载体横轴向前,zb轴沿载体立轴向上,所述雷达坐标系O-xryrzr为右手坐标系,其中原点为雷达底部安装面的中心,zr轴沿雷达安装面垂直光束出射的方向,雷达安装面与载体安装面重合时,yr轴指向yb轴,xr轴指向-xb轴,zr轴指向-zb轴,所述R1、R2、R3、R4从雷达坐标系原点射出后与雷达坐标系的正z轴夹角均为15°,R1-R4的光轴坐标由公式(1)-(4)得到:
Figure BDA0003451766450000021
Figure BDA0003451766450000022
Figure BDA0003451766450000023
Figure BDA0003451766450000024
由(1)-(4)公式可得到R的矩阵表达式(5)为:
Figure BDA0003451766450000025
进一步的,步骤S1的旋转角度由载体坐标系相对于地理坐标系在xgyg平面绕zg轴旋转得到航向角θ,载体坐标系相对于地理坐标系在ygzg平面绕xg轴旋转得到俯仰角β,载体坐标系相对于地理坐标系在xgzg平面绕yg轴旋转得到横滚角γ,单位矢量在载体坐标系中,分别依次经过航向角θ,俯仰角β,横滚角γ的顺序旋转,得到步骤S2所述的矢量旋转表达式(6)-(8)如下:
Figure BDA0003451766450000031
Figure BDA0003451766450000032
Figure BDA0003451766450000033
单位矢量为r=i+0j+0k(r=[1 0 0]),经过旋转后的矢量旋转表达式为该矢量矩阵右乘旋转矩阵,具体表达式(8)如下:
Figure BDA0003451766450000034
矩阵相乘满足结合律,则有公式(9)如下:
Figure BDA0003451766450000035
H为航向角度且满足三轴旋转矩阵,由H得到单位矢量r的旋转表达式(10)如下:
r'=rH (10)。
进一步的,步骤S3所述本振光R经过旋转后在载体坐标系的矩阵表达式R'(11)如下:
Figure BDA0003451766450000036
Figure BDA0003451766450000041
进一步的,基于地理坐标系的真实风速矢量为Vwind,雷达的速度矢量为Vins,雷达四束光测得的相对速度标量为Vr1、、Vr2、Vr3、Vr4,则有真实风速矢量与雷达的速度矢量的差值投影到各个光轴的值等于雷达四束光光轴的相对速度标量值,表达式(12)-(16)如下:
由vr=(Vwind-Vins)·R' (12)
得到:
vr1=(U-Vx)×mo1+(V-Vy)×mo2+(W-Vz)×mo3 (13)
vr2=(U-Vx)×m11+(V-Vy)×m12+(W-Vz)×m13 (14)
vr3=(U-Vx)×m21+(V-Vy)×m22+(W-Vz)×m23 (15)
vr4=(U-Vx)×m31+(V-Vy)×m32+(W-Vz)×m33 (16)
变换可得表达式(17)-(20)如下:
vr1+Vx×mo1+Vy×mo2+Vz×mo3=U×mo1+V×mo2+W×mo3 (17)
vr2+Vx×m11+Vy×m12+Vz×m13=U×m11+V×m12+W×m13 (18)
vr3+Vx×m21+Vy×m22+Vz×m23=U×m21+V×m22+W×m23 (19)
vr4+Vx×m31+Vy×m32+Vz×m33=U×m31+V×m32+W×m33 (20)
进一步的,采用最小二乘法解出满足表达式(17)-(20)的U、V、W值,其中垂直风速为W值,水平风速公式(21)如下:
Figure BDA0003451766450000042
由公式(21)可得风向为矢量(-U-V)与y轴的风向夹角公式(22)-(23)如下:
当-U≥0时,
Figure BDA0003451766450000043
当-U≤0时,
Figure BDA0003451766450000051
进一步的,所述雷达采用四波束同步收发系统在高速运动的载体平台上进行数据采集,所述四波束同步收发系统包括窄线宽种子光、第一声光调制器、第二声光调制器、多路输出放大器、第一环形器、第二环形器、第三环形器、第四环形器、多轴收/发天线、分束器、第一合束器、第二合束器、第三合束器、第四合束器、第一探测器、第二探测器、第三探测器、第四探测器、多通道放大器和多通道数据采集卡,所述窄线宽种子光分别连接第一声光调制器和分束器,第一声光调制器连接第二声光调制器,第二声光调制器连接多路输出放大器,多路输出放大器分别连接第一环形器、第二环形器、第三环形器和第四环形器,第一环形器、第二环形器、第三环形器和第四环形器均连接多轴收/发天线,所述第一环形器、第一合束器和第一探测器依次连接,所述第二环形器、第二合束器和第二探测器依次连接,所述第三环形器、第三合束器和第三探测器依次连接,所述第四环形器、第四合束器和第四探测器依次连接,所述第一探测器、第二探测器、第三探测器和第四探测器通过不同通道连接多通道放大器,多通道放大器通过不同通道连接多通道数据采集卡,所述分束器分别连接第一合束器、第二合束器、第三合束器和第四合束器。
进一步的,所述四波束同步收发系统同时以一定立体角向大气中发射4束激光,并接受各自光轴路上的后向反射光进行光束扫描得到相干外差信号。
本发明的有益效果:
本发明提出的高速运动下的相干外差测风雷达修正方法,可将相干外差测风技术应用于高速运动平台,并修正位于高速平台上的雷达相对于地面的运动数据,从而获得真实的大气风场的三维信息。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简要介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提出的高速运动下的相干外差测风雷达修正方法的四波束同步收发系统结构框图;
图2为本发明提出的高速运动下的相干外差测风雷达修正方法的高速运动下的相干外差测风雷达原理流程图;
图3为本发明提出的高速运动下的相干外差测风雷达修正方法的坐标图;
图4为本发明提出的高速运动下的相干外差测风雷达修正方法的四束光轴坐标图;
图5为本发明提出的高速运动下的相干外差测风雷达修正方法的四束光轴坐标在xy面的投影坐标图;
图6为本发明提出的高速运动下的相干外差测风雷达修正方法的雷达与载体安装坐标图。
图中,1-窄线宽种子光,2-第一声光调制器,3-第二声光调制器,4-多路输出放大器,5-第一环形器,6-第二环形器,7-第三环形器,8-第四环形器,9-多轴收/发天线,10-分束器,11-第一合束器,12-第二合束器,13-第三合束器,14-第四合束器,15-第一探测器,16-第二探测器,17-第三探测器,18-第四探测器,19-多通道放大器,20-多通道数据采集卡。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
本发明以激光相干外差测风雷达搭载于机载运动平台为例,雷达向下安装探测机下三维风场,其他搭载或安装方式以此方法类推如雷达光束向上安装,可探测上方大气三维风场场,相干外差技术可探测激光径向矢量方向每个物理分辨率距离内的气溶胶粒子与飞机在激光径向方向上的相对运动速度,由于飞机的飞行姿态、航向变化以及飞机抖动,使得激光束的径向矢量在地理坐标系中的位置实时变化,而相干外差多普勒理论探测到的气溶胶与飞机相对运动速度的大小由飞机速度和大气的风场共同作用,因此想要获得大气的三维风场,需要按照以下方法进行:
如图1,首先,雷达光源在一个周期的光束扫描时间内,假设载体平台的姿态和速度不变,考虑高速运动,本发明采用四波束收发同步扫描方案,光学收/发天线为一个四轴天线,与同时输出四束相同激光束的脉冲激光源相连,构成四波束同步光路收/发系统。
其次,雷达选取声光调制器对本振光进行移频时,结合运动载体的运动速度,选取合适的向右移频量,确保载体的运动速度分量不会超出测风雷达探测的速度量程。
再次,增加运动传感器,该运动传感器根据实际情况安装于雷达或载体上探测运动速度、航向以及姿态角等参数,该传感器主要由高精度惯导及GPS或北斗组成。若雷达安装的载体提供运动参数数据,只需在雷达上预留通信接口,接收载体运动参数即可,同时注意雷达与载体的安装方式与角度,结构安装采用硬连接方式安装在载体的轴系上;若雷达自身携带运动传感器,则雷达与载体通过减震器或减振隔离装置安装。
最后,雷达同步采集四束光路的相干外差信号以及运动传感器参数,通过运动修正算法,获得大气三维风场。
实施例2
如图2、图3、图4、图6,以机载为例,因目标载体飞机自身搭载有惯导系统,能够实时地向机载激光测风雷达提供实时飞行姿态、航向、速度、高度等数据。具体运动修正算法内容及步骤如下:
一)坐标系建模理论
首先,建立相关坐标系。定义三个坐标系,其一地理坐标系,地理坐标系O-xg yg zg采用东北天坐标系。原点为载体重心,xg轴指向东,yg轴指向正北,zg轴指向天。该坐标系也是飞机惯导和大气风场的坐标系。其二,飞机载体坐标系,载体坐标系O-xb yb zb,原点位于飞机中心,xb轴沿载体横轴向右(即指向右侧翅膀),yb轴沿载体纵轴向前(指向机头),zb轴沿飞机立轴向上。该系与飞机固连,载体坐标系相对于地理坐标系的方位关系就是飞机的航向和姿态。其三,雷达坐标系,雷达坐标系O-xr yr zr,是一个右手坐标系,原点位于底部安装平面的中心,zr轴沿着安装面垂直顺着光束出射的方向。
如图2,航向角θ为机体坐标系相对于地理坐标系在yg xg平面绕zg轴旋转的角度,范围0到360°,北向东为正,yb轴在yg xg平面的投影与yg平行(即机头在yg xg平面的投影指向正北)为0°。俯仰角β为机体坐标系相对于地理坐标系在yg zg平面绕xg轴旋转的角度,范围0到±90°,向天为正,向地为负,yb轴在yg zg平面上的投影yg轴重合为0°。横滚角γ为机体坐标系相对于地理坐标系在Xg Zg平面上绕yg轴旋转的角度,范围为0到±180°,东向天为正,东向地为负,xb轴在平面zg xg平面的投影,与xg轴重合为0°,飞机载体在飞行过程中,惯导给出的飞行速度Vins(Vx Vy Vz),(Vx向东的飞行速率,Vy向北的飞行速率,Vz向上的飞行速率)、飞机姿态及航向H(θβγ)(航向角θ,俯仰β,横滚γ)。
安装后,雷达安装面与飞机安装面重合,yr轴指向yb轴,xr轴指向-xb轴,zr轴指向-zb轴,即雷达坐标系与飞机坐标系平行如图6所示。
光轴的位置坐标,雷达四束光与正Z轴夹角为15°,均匀分布在XYZ直角坐标系的正Z立体空间内,四束光是由激光源提供,来自同一本振光,通过设计及安装,设定四束光坐标及在XY面投影如图5所示。
四束光轴R1-R4的光轴坐标由以下公式可得:
Figure BDA0003451766450000081
Figure BDA0003451766450000082
Figure BDA0003451766450000083
Figure BDA0003451766450000084
R的矩阵表达式为:
Figure BDA0003451766450000085
二)矢量旋转表达
单位矢量在坐标系某一位置,分别依次经过航向旋转θ,俯仰旋转β,横滚旋转γ的顺序旋转,得到各自的旋转矩阵如下:
Figure BDA0003451766450000091
(原始单位矩阵,每行的含义为各个轴的单位矢量在XYZ坐标系的表达的,每列的含义为每个坐标轴下的三个轴单位矢量的分量值)
Figure BDA0003451766450000092
Figure BDA0003451766450000093
Figure BDA0003451766450000094
假设一单位矢量r=i+0j+0k(r=[1 0 0],经过旋转后的矢量表达式,应为该矢量矩阵右乘旋转矩阵,具体表达式如下:
Figure BDA0003451766450000095
矩阵相乘满足结合律,则有H为:
Figure BDA0003451766450000096
H为三轴旋转矩阵,则有:
r’=rH。
三)光轴经旋转后在直角坐标系的表达式
光轴R经过旋转后在XYZ坐标系的矩阵为:
Figure BDA0003451766450000101
四)风速求取
设自然界的真实风速矢量为Vwind(u v w)(东北天坐标系),惯导测得雷达的速度矢量为Vins(Vx Vy Vz)(东北天坐标系),雷达四束光测得的相对速度标量为Vr1、Vr2、Vr3、Vr4(顺着光矢量风向为正,逆着光轴方向为负)。
则有,真实风矢量与惯导运动矢量的差值投影到各个光轴的数值为雷达各个光轴测得的相对速度标量值,表达式如下:
vr=(Vwind-Vins)·R'
得到:
vr1=(U-Vx)×mo1+(V-Vy)×mo2+(W-Vz)×mo3
vr2=(U-Vx)×m11+(V-Vy)×m12+(W-Vz)×m13
vr3=(U-Vx)×m21+(V-Vy)×m22+(W-Vz)×m23
vr4=(U-Vx)×m31+(V-Vy)×m32+(W-Vz)×m33
变换得:
vr1+Vx×mo1+Vy×mo2+Vz×mo3=U×mo1+V×mo2+W×mo3
vr2+Vx×m11+Vy×m12+Vz×m13=U×m11+V×m12+W×m13
vr3+Vx×m21+Vy×m22+Vz×m23=U×m21+V×m22+W×m23
vr4+Vx×m31+Vy×m32+Vz×m33=U×m31+V×m32+W×m33
五)风场三维信息转换
采用最小二乘法解出满足以上4个等式的U、V、W值。
自然界风场表达规则如下:
1)垂直气流(垂直风)
垂直气流指东北天坐标下,垂直指向天的速率及方向,即W值。
2)水平风
水平风速和风向指水平风向的合风速及合风速的来风方向,风向范围为0到360°(北吹南的风为0°,称为北风;东吹向西的风为90°,称为东风)。
Figure BDA0003451766450000111
风向为矢量(-U-V)与y轴(北)的风向夹角,公式如下:
当-U≥0时,
Figure BDA0003451766450000112
当-U<0时,
Figure BDA0003451766450000113
实施例2
本实施例提出用于发射四束光并进行采集的四波束同步收发系统。
所述四波束同步收发系统包括窄线宽种子光1、第一声光调制器2、第二声光调制器3、多路输出放大器4、第一环形器5、第二环形器6、第三环形器7、第四环形器8、多轴收/发天线9、分束器10、第一合束器11、第二合束器12、第三合束器13、第四合束器14、第一探测器15、第二探测器16、第三探测器17、第四探测器18、多通道放大器19和多通道数据采集卡20,所述窄线宽种子光1分别连接第一声光调制器2和分束器10,第一声光调制器2连接第二声光调制器3,第二声光调制器3连接多路输出放大器4,多路输出放大器4分别连接第一环形器5、第二环形器6、第三环形器7和第四环形器8,第一环形器5、第二环形器6、第三环形器7和第四环形器8均连接多轴收/发天线9,所述第一环形器5、第一合束器11和第一探测器15依次连接,所述第二环形器6、第二合束器12和第二探测器16依次连接,所述第三环形器7、第三合束器13和第三探测器17依次连接,所述第四环形器8、第四合束器14和第四探测器18依次连接,所述第一探测器15、第二探测器16、第三探测器17和第四探测器18通过不同通道连接多通道放大器19,多通道放大器19通过不同通道连接多通道数据采集卡20,所述分束器10分别连接第一合束器11、第二合束器12、第三合束器13和第四合束器14,所述四波束同步收发系统同时以一定立体角向大气中发射4束激光,并接受各自光轴路上的后向反射光进行光束扫描得到相干外差信号。
实施例3
本实施例提出四波束同步收发系统的工作原理。
窄线宽种子光1分出两束光,一束光为信号光,经过第一声光调制器2和第二声光调制器3调制,变成脉冲信号,并完成移频,经移频后的脉冲信号进入具有多级预放大的多路输出放大器4,输出四束相同的激光束,各自经过环形器后链接到多轴收/发天线9,四轴天线将激光束发射到大气中,同时接受来自气溶胶粒子反射回的回波信号,再经环形器后分别连接到第一合束器11、第二合束器12、第三合束器13、第四合束器14,窄线宽种子光1分出的另一束光作为参考光经分束器10分出四束相同的参考光,与经天线和环形器接收的四个回波信号一起进入四个合束器,在合束器中进行合束相干,再经第一探测器15、第二探测器16、第三探测器17、第四探测器18转化为外差电信号,转化后的电信号非常微弱,需经多通道放大器19放大后,被多通道数据采集卡20采集后转换成数字信号完成数据采集。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。
本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (8)

1.高速运动下的相干外差测风雷达修正方法,其特征在于,雷达在高速运动的载体平台上进行相干外差信号采集,采集的相干外差信号通过运动修正算法得到大气的三维风场信息;
所述运动修正算法包括以下步骤:
S1、建立3组模型坐标系,并根据3组模型坐标系得到相关数据,其中,相关数据包括2组模型坐标系的旋转角度以及雷达本振光R发射产生的四束光R1、R2、R3、R4的光轴坐标;
S2、根据S1的旋转角度得到每个旋转角度的矢量旋转表达式以及单位矢量的旋转表达式;
S3、根据S1所述的四束光R1、R2、R3、R4的光轴位置坐标得到本振光R的光轴旋转后在直角坐标系的表达式R’
S4、根据S3得到的R’与真实风速矢量以及雷达速度矢量计算出中四束光的各光轴的相对速度标量;
S5、采用最小二乘法解出S4中的相对速度标量的参数值,得到大气的三维风场信息,所述三维风场信息包括垂直风速和水平风速。
2.根据权利要求1所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,步骤S1所述的3组模型坐标系为地理坐标系、载体坐标系和雷达坐标系,所述地理坐标系O-xgygzg采用东北天坐标系,其中原点为载体重心,xg轴指向东,yg轴指向正北,zg轴指向天,所述载体坐标系O-xbybzb的原点为载体中心,xb轴沿载体横轴向右,yb轴沿载体横轴向前,zb轴沿载体立轴向上,所述雷达坐标系O-xryrzr为右手坐标系,其中原点为雷达底部安装面的中心,zr轴沿雷达安装面垂直光束出射的方向,雷达安装面与载体安装面重合时,yr轴指向yb轴,xr轴指向-xb轴,zr轴指向-zb轴,所述R1、R2、R3、R4从雷达坐标系原点射出后与雷达坐标系的正z轴夹角均为15°,R1-R4的光轴坐标由公式(1)-(4)得到:
Figure FDA0003451766440000011
Figure FDA0003451766440000012
Figure FDA0003451766440000013
Figure FDA0003451766440000014
由(1)-(4)公式可得到R的矩阵表达式(5)为:
Figure FDA0003451766440000021
3.根据权利要求2所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,步骤S1的旋转角度由载体坐标系相对于地理坐标系在xgyg平面绕zg轴旋转得到航向角θ,载体坐标系相对于地理坐标系在ygzg平面绕xg轴旋转得到俯仰角β,载体坐标系相对于地理坐标系在xgzg平面绕yg轴旋转得到横滚角γ,单位矢量在载体坐标系中,分别依次经过航向角θ,俯仰角β,横滚角γ的顺序旋转,得到步骤S2所述的矢量旋转表达式(6)-(8)如下:
Figure FDA0003451766440000022
Figure FDA0003451766440000023
Figure FDA0003451766440000024
单位矢量为r=i+0j+0k(r=[1 0 0]),经过旋转后的矢量旋转表达式为该矢量矩阵右乘旋转矩阵,具体表达式(8)如下:
Figure FDA0003451766440000025
矩阵相乘满足结合律,则有公式(9)如下:
Figure FDA0003451766440000026
H为航向角度且满足三轴旋转矩阵,由H得到单位矢量r的旋转表达式(10)如下:
r'=rH (10)。
4.根据权利要求3所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,步骤S3所述本振光R经过旋转后在载体坐标系的矩阵表达式R'(11)如下:
Figure FDA0003451766440000031
5.根据权利要求4所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,基于地理坐标系的真实风速矢量为Vwind,雷达的速度矢量为Vins,雷达四束光测得的相对速度标量为Vr1、、Vr2、Vr3、Vr4,则有真实风速矢量与雷达的速度矢量的差值投影到各个光轴的值等于雷达四束光光轴的相对速度标量值,表达式(12)-(16)如下:
由vr=(Vwind-Vins)·R' (12)
得到:
vr1=(U-Vx)×mo1+(V-Vy)×mo2+(W-Vz)×mo3 (13)
vr2=(U-Vx)×m11+(V-Vy)×m12+(W-Vz)×m13 (14)
vr3=(U-Vx)×m21+(V-Vy)×m22+(W-Vz)×m23 (15)
vr4=(U-Vx)×m31+(V-Vy)×m32+(W-Vz)×m33 (16)
变换可得表达式(17)-(20)如下:
vr1+Vx×mo1+Vy×mo2+Vz×mo3=U×mo1+V×mo2+W×mo3(17)
vr2+Vx×m11+Vy×m12+Vz×m13=U×m11+V×m12+W×m13(18)
vr3+Vx×m21+Vy×m22+Vz×m23=U×m21+V×m22+W×m23(19)
vr4+Vx×m31+Vy×m32+Vz×m33=U×m31+V×m32+W×m33(20)
6.根据权利要求5所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,采用最小二乘法解出满足表达式(17)-(20)的U、V、W值,其中垂直风速为W值,水平风速公式(21)如下:
Figure FDA0003451766440000041
由公式(21)可得风向为矢量(-U-V)与y轴的夹角公式(22)-(23)如下:
当-U≥0时,
Figure FDA0003451766440000042
当-U≤0时,
Figure FDA0003451766440000043
7.根据权利要求1所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,所述雷达采用四波束同步收发系统在高速运动的载体平台上进行数据采集,所述四波束同步收发系统包括窄线宽种子光、第一声光调制器、第二声光调制器、多路输出放大器、第一环形器、第二环形器、第三环形器、第四环形器、多轴收/发天线、分束器、第一合束器、第二合束器、第三合束器、第四合束器、第一探测器、第二探测器、第三探测器、第四探测器、多通道放大器和多通道数据采集卡,所述窄线宽种子光分别连接第一声光调制器和分束器,第一声光调制器连接第二声光调制器,第二声光调制器连接多路输出放大器,多路输出放大器分别连接第一环形器、第二环形器、第三环形器和第四环形器,第一环形器、第二环形器、第三环形器和第四环形器均连接多轴收/发天线,所述第一环形器、第一合束器和第一探测器依次连接,所述第二环形器、第二合束器和第二探测器依次连接,所述第三环形器、第三合束器和第三探测器依次连接,所述第四环形器、第四合束器和第四探测器依次连接,所述第一探测器、第二探测器、第三探测器和第四探测器通过不同通道连接多通道放大器,多通道放大器通过不同通道连接多通道数据采集卡,所述分束器分别连接第一合束器、第二合束器、第三合束器和第四合束器。
8.根据权利要求7所述的一种高速运动下的相干外差测风雷达修正方法,其特征在于,所述四波束同步收发系统同时以一定立体角向大气中发射4束激光,并接受各自光轴路上的后向反射光进行光束扫描得到相干外差信号。
CN202111681292.6A 2021-12-31 2021-12-31 高速运动下的相干外差测风雷达修正方法 Pending CN114325666A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111681292.6A CN114325666A (zh) 2021-12-31 2021-12-31 高速运动下的相干外差测风雷达修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111681292.6A CN114325666A (zh) 2021-12-31 2021-12-31 高速运动下的相干外差测风雷达修正方法

Publications (1)

Publication Number Publication Date
CN114325666A true CN114325666A (zh) 2022-04-12

Family

ID=81023090

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111681292.6A Pending CN114325666A (zh) 2021-12-31 2021-12-31 高速运动下的相干外差测风雷达修正方法

Country Status (1)

Country Link
CN (1) CN114325666A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117420570A (zh) * 2023-12-19 2024-01-19 南京牧镭激光科技股份有限公司 一种激光雷达姿态补偿测风方法、设备以及存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117420570A (zh) * 2023-12-19 2024-01-19 南京牧镭激光科技股份有限公司 一种激光雷达姿态补偿测风方法、设备以及存储介质
CN117420570B (zh) * 2023-12-19 2024-03-29 南京牧镭激光科技股份有限公司 一种激光雷达姿态补偿测风方法、设备以及存储介质

Similar Documents

Publication Publication Date Title
CN110108984B (zh) 电力巡线激光雷达系统多传感器的空间关系同步方法
CN113311436B (zh) 一种移动平台上激光测风雷达运动姿态测风订正方法
CN112197761B (zh) 一种高精度多旋翼机协同定位方法及系统
CN108375382A (zh) 基于单目视觉的位置姿态测量系统精度检校方法和装置
CN107132542B (zh) 一种基于光学和多普勒雷达的小天体软着陆自主导航方法
CN109911188A (zh) 非卫星导航定位环境的桥梁检测无人机系统
CN107765226B (zh) 一种sar卫星雷达回波模拟方法、系统和介质
CN110988851B (zh) 一种基于星位优化的异轨单星分时测频定位方法
KR100860767B1 (ko) 항공 레이저 측량 데이터를 이용한 수치도화 제작 장치 및방법
CN111829512B (zh) 一种基于多传感器数据融合的auv导航定位方法及系统
CN111380518B (zh) 一种引入径向速度的sins/usbl紧组合导航定位方法
CN108896957A (zh) 一种无人机操控信号源的定位系统及方法
CN109581365A (zh) 一种基于多旋翼无人机的生命探测系统、方法
CN111366913A (zh) 星载激光雷达光轴指向测量误差的标定方法
CN109633724A (zh) 基于单星与多地面站联合测量的无源目标定位方法
CN109471447A (zh) 无人机导航方法、装置、无人机和数据可读存储装置
CN108444476A (zh) 一种考虑水声通信延迟的多水下无人航行器极区协同导航方法
CN111694012A (zh) 一种基于机载激光雷达的三维地形在线生成方法及系统
CN110244308A (zh) 一种适用于无人机测高定姿的激光传感器及其工作方法
CN110220533A (zh) 一种基于传递对准的机载光电吊舱设备失准角标定方法
CN112710303A (zh) 由运动平台运动引起目标在视场中姿态角θ变化的确定方法
CN103823209B (zh) 用于轻小型合成孔径雷达系统中低成本运动误差测量装置
CN113409400A (zh) 一种基于自动跟踪的机载光电系统目标地理定位方法
CN114325666A (zh) 高速运动下的相干外差测风雷达修正方法
CN102262230A (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