CN105891821B - 一种机载下视测量目标的自动跟踪方法 - Google Patents
一种机载下视测量目标的自动跟踪方法 Download PDFInfo
- Publication number
- CN105891821B CN105891821B CN201610350347.8A CN201610350347A CN105891821B CN 105891821 B CN105891821 B CN 105891821B CN 201610350347 A CN201610350347 A CN 201610350347A CN 105891821 B CN105891821 B CN 105891821B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- mtr
- msub
- axis
- 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
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
- 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/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
Abstract
公开了一种机载下视测量目标的自动跟踪方法,包括如下步骤:S1、获取飞机和目标的位置信息,根据位置信息确定将地心坐标系的坐标原点移动至飞机处后的目标的第一坐标;S2、获取飞机的姿态信息,根据姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标;S3、基于目标的第四坐标确定雷达天线的方位角和俯仰角,依据方位角和俯仰角调整雷达天线的指向。本发明能够实时自动调整天线波束指向,保证被测量目标始终处于天线主波束中心,无需人工手动操作伺服控制杆来对准目标,显著提高机载下视测量的效率和精度。
Description
技术领域
本发明涉及信号特征控制技术领域,特别涉及一种机载下视测量目标的自动跟踪方法。
背景技术
以下对本发明的相关技术背景进行说明,但这些说明并不一定构成本发明的现有技术。
近年来,由于机载下视测量系统能够获取被测目标全方位、大入射余角条件下的电磁散射特性数据,受到广泛的关注。在机载下视测量过程中,载机与目标是相对运动的,为确保测量的准确性,必须保证机载测量雷达天线始终对准目标。
机载下视测量雷达在工作中需要不断跟踪目标,以得到准确的测量结果。通常这一工作主要由人工手动完成。但由于飞机平台在空中飞行时容易受到气流扰动等影响,特别是在恶劣气候条件下,飞机姿态不稳定,再加上飞行过程中载机机舱环境的复杂性,对人工手动操作提出了较大挑战。
因此,现有技术中需要一种能够解决由于飞机平台在空中飞行时容易受到气流扰动而导致人工手动跟踪目标测量精度差的问题的解决方案。
发明内容
本发明的目的在于提出一种机载下视测量目标的自动跟踪方法,能够实时自动调整天线波束指向,保证被测量目标始终处于天线主波束中心,提高机载下视测量的效率和精度。
根据本发明的机载下视测量目标的自动跟踪方法,包括如下步骤:
S1、获取飞机和目标的位置信息,根据所述位置信息确定将地心坐标系的坐标原点移动至飞机处后的目标的第一坐标;
S2、获取飞机的姿态信息,根据所述姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标;
S3、基于目标的第四坐标确定雷达天线的方位角和俯仰角,依据所述方位角和俯仰角调整雷达天线的指向;
其中,
所述位置信息包括:经度、纬度和高度,所述高度是指当前位置到地表的距离;所述姿态信息包括偏航角、俯仰角和横滚角;
地心坐标系中,Z轴指向地球北极,X轴指向格林威治子午面与地球赤道的交点,Y轴垂直XOZ面;
机体向下坐标系中,X轴指向飞机的机头方向,Y轴指向飞机的机身右侧,Z轴垂直机体向下。
优选地,步骤S2中根据所述姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标包括:
S21、将目标的第一坐标转换成东北天坐标系中的第二坐标;
S22、将目标在东北天坐标系中的第二坐标转换成东南地坐标系中的第三坐标;
S23、根据所述姿态信息,将目标的第三坐标转换成机体向下坐标系中的第四坐标;
其中,东南地坐标系中,X轴指向正东方向,Y轴指向正南方向,Z轴指向地心;东北天坐标系中,X轴指向正东方向,Y轴指向正北方向,Z轴指向地心反向。
优选地,目标的第二坐标为:
[x2,y2,z2]'=M1M0[x1,y1,z1]'
其中,
lon表示飞机的经度,lat表示飞机的纬度,alt表示飞机的高度;x1、y1、z1分别表示将地心坐标系的坐标原点移动至飞机处后目标在X轴、Y轴和Z轴的坐标值;x2、y2、z2分别表示目标在东北天坐标系中的X轴、Y轴和Z轴的坐标值。
优选地,目标的第三坐标为:
[x3,y3,z3]'=M180[x2,y2,z2]'
其中,
x3、y3、z3分别表示目标在东南地坐标系中的X轴、Y轴和Z轴的坐标值。
优选地,目标的第四坐标为:
[x4,y4,z4]'=MROLLMPITCHMYAQ[x3,y3,z3]'
其中,
x4、y4、z4分别表示目标在机体向下坐标系中的X轴、Y轴和Z轴的坐标值;YAW表示飞机的偏航角,PITCH表示飞机的俯仰角,ROLL表示飞机的横滚角。
优选地,步骤S3中基于目标的第四坐标确定雷达天线的方位角和俯仰角包括:
S31、根据雷达在飞机上的安装误差对目标的第四坐标进行校正,得到目标的第五坐标;
S32、根据第五坐标确定雷达天线的方位角和俯仰角。
优选地,目标的第五坐标为:
[x5,y5,z5]'=[x4+x0,y4+y0,z4+z0]'
其中,x5、y5、z5分别表示校正安装误差后的目标在机体向下坐标系中的X轴、Y轴和Z轴的坐标值;x0、y0、z0分别表示在机体向下坐标系的X轴、Y轴和Z轴方向上的飞机惯导设备在天线平台的安装误差值。
优选地,雷达天线的方位角azT为:
雷达天线的俯仰角pitchT为:
其中,a表示在WGS-84地球模型中椭球的长半轴。
优选地,飞机的姿态信息为:
YAW=YAW0+eyaw
PITCH=PITCH0+epitch
ROLL=ROLL0+eroll
其中,YAW0、PITCH0、ROLL0分别表示飞机惯导设备测得的飞机的偏航角、俯仰角和横滚角;eyaw、epitch、eroll分别表示飞机惯导设备在天线平台的偏航角安装误差值、俯仰角安装误差值和横滚角安装误差值。
优选地,飞机在地心坐标系的坐标为:
目标在地心坐标系的坐标为:
其中,
a表示在WGS-84地球模型中椭球的长半轴,a=6378137m;偏心率平方e2=0.0066944;x、y、z分别表示飞机在地心坐标系中的X轴、Y轴和Z轴的坐标值;xT、yT、zT分别表示目标在地心坐标系中的X轴、Y轴和Z轴的坐标值;lon表示飞机的经度,lat表示飞机的纬度,alt表示飞机的高度;lonT表示目标的经度,latT表示目标的纬度,altT表示目标的高度。
根据本发明的机载下视测量目标的自动跟踪方法,根据飞机和目标的位置信息和姿态信息确定目标相对于天线平台的位置关系,从而实时自动调整天线波束指向,保证被测量目标始终处于天线主波束中心。根据本发明的自动跟踪方法在机载下视测量目标时,无需人工手动操作伺服控制杆来对准目标,能显著提高机载下视测量的效率和精度。
附图说明
通过以下参照附图而提供的具体实施方式部分,本发明的特征和优点将变得更加容易理解,在附图中:
图1是根据本发明的机载下视测量目标的自动跟踪方法的流程图;
图2是示出根据本发明的地心坐标系的示意图;
图3是根据本发明的东北天坐标系的示意图;
图4是根据本发明的机体向下坐标系的示意图;
图5是根据本发明的飞机偏航角的示意图;
图6是根据本发明的飞机俯仰角的示意图;
图7是根据本发明的飞机横滚角的示意图。
具体实施方式
下面参照附图对本发明的示例性实施方式进行详细描述。对示例性实施方式的描述仅仅是出于示范目的,而绝不是对本发明及其应用或用法的限制。
人工手动跟踪目标主要通过人眼进行观察,准确性较差。受气流扰动等因素的影响,飞机飞行时的姿态角不稳定,进一步增加了人工手动跟踪目标的操作难度,使得测量精度无法得到保证。为了解决现有技术中的问题,本发明根据飞机和目标的位置信息和姿态信息自动确定和调整雷达天线的波束指向,图1示出了根据本发明的机载下视测量目标的自动跟踪方法的流程图。
根据本发明的自动跟踪方法,在步骤S1中首先获取飞机和目标的在地心坐标系中的位置信息,包括:经度、纬度和高度。图2示出了根据本发明的地心坐标系的示意图,其中,地心坐标系的Z轴指向地球北极、X轴指向格林威治子午面与地球赤道的交点、Y轴垂直XOZ面。图中,点O代表地心坐标系的原点,点P代表飞机的位置,lon表示飞机的经度,lat表示飞机的纬度,alt表示飞机的高度,本发明中的高度是指当前位置到地表的距离。
可以按照如下公式确定飞机和目标在地心坐标系的坐标:
飞机在地心坐标系的坐标为:
目标在地心坐标系的坐标为:
其中,
a表示在WGS-84地球模型中椭球的长半轴,a=6378137m;偏心率平方e2=0.0066944;x、y、z分别表示飞机在地心坐标系中的X轴、Y轴和Z轴的坐标值;xT、yT、zT分别表示目标在地心坐标系中的X轴、Y轴和Z轴的坐标值;lon表示飞机的经度,lat表示飞机的纬度,alt表示飞机的高度;lonT表示目标的经度,latT表示目标的纬度,altT表示目标的高度。
本发明通过飞机和目标的信息确定目标相对于飞机的位置,进行调整雷达天线的波束指向。例如,可以根据飞机和目标在地心坐标系中的位置分别获取飞机和目标在机体向下坐标系中的坐标信息,然后再根据飞机和目标在机体向下坐标系中的坐标信息确定目标与飞机的相对位置。但是这种分别计划飞机和目标在机体向下坐标系中的坐标信息的方式步骤繁琐,工作量大。本发明直接根据飞机和目标在地心坐标系中的位置信息确定目标相对于飞机的相对位置坐标,然后再确定将地心坐标系的坐标原点移动至飞机处后的目标的第一坐标。根据本发明,不仅便于直观了解确定目标相对于飞机的位置,还能简化数据处理过程,降低后续步骤进行坐标转化时的工作量,提高工作效率。
S2、获取飞机的姿态信息,包括偏航角、俯仰角和横滚角,根据飞机的姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标。图4示出了根据本发明的机体向下坐标系的示意图,其中,机体向下坐标系的X轴指向飞机的机头方向、Y轴指向飞机的机身右侧、Z轴垂直机体向下。偏航角指机头方向与正北向的夹角,顺时针为正,旋转轴为Z;俯仰角指机身前后向相对水平面的转角,以Y为旋转轴,右手螺旋为正,抬头为正,低头为负;横滚角指机身左右向相对水平面的转角,以X为旋转轴,右手螺旋为正,左高为正,右高为负。
飞机沿着预设的航行轨迹飞行,为了便于统计分析飞机的速度和航向等信息,在本发明的一些实施例中,将地心坐标系的坐标原点移动至飞机处后的目标的第一坐标后,将该第一坐标转化成以飞机为坐标原点的东北天坐标系中的第二坐标。由于飞机惯导设备一般标称东北天坐标系的Z轴以下,因此在获得第二坐标以后进一步将其转化成东南地坐标系中的第三坐标。具体地址,步骤S2中根据飞机的姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标包括:
S21、将目标的第一坐标转换成东北天坐标系中的第二坐标;
S22、将目标在东北天坐标系中的第二坐标转换成东南地坐标系中的第三坐标;
S23、根据飞机的姿态信息,将目标的第三坐标转换成机体向下坐标系中的第四坐标;
其中,东南地坐标系中,X轴指向正东方向,Y轴指向正南方向,Z轴指向地心;东北天坐标系中,X轴指向正东方向,Y轴指向正北方向,Z轴指向地心反向,参见图3。
优选地,将地心坐标系的坐标原点移动至飞机处后,先以地心坐标系的Z轴为旋转轴,把移动后的地心坐标系按右手螺旋方向旋转(90+lon)度;再以地心坐标系的X轴为旋转轴,把地心坐标系按右手螺旋方向旋转(90-lat)度,得到目标的第二坐标为:
[x2,y2,z2]'=M1M0[x1,y1,z1]'
其中,
lon表示飞机的经度,lat表示飞机的纬度,alt表示飞机的高度;x1、y1、z1分别表示将地心坐标系的坐标原点移动至飞机处后目标在X轴、Y轴和Z轴的坐标值;x2、y2、z2分别表示目标在东北天坐标系中的X轴、Y轴和Z轴的坐标值。
优选地,步骤S22中以东北天坐标系的X轴为旋转轴,将东北天坐标系顺时针旋转180度得到东南地坐标系,目标的第三坐标为:
[x3,y3,z3]'=M180[x2,y2,z2]'
其中,
x3、y3、z3分别表示目标在东南地坐标系中的X轴、Y轴和Z轴的坐标值。
优选地,将目标在东南地坐标系中的第三坐标转换成在机体向下坐标系中的第四坐标的过程中,先以东南地坐标系的Z轴为旋转轴,逆右手螺旋转动(90-YAW),YAW表示飞机的偏航角,参见图5。旋转矩阵为:
然后以东南地坐标系的Y轴为旋转轴,右手螺旋转动PITCH角度,PITCH表示飞机的俯仰角,参见图6。旋转矩阵为:
最后以东南地坐标系的X轴为旋转轴,右手螺旋转动ROLL角度,ROLL表示飞机的横滚角,参见图7。旋转矩阵为:
按照上述方法得到的目标在机体向下坐标系中的第四坐标为:
[x4,y4,z4]'=MROLLMPITCHMYAQ[x3,y3,z3]'。
其中,x4、y4、z4分别表示目标在机体向下坐标系中的X轴、Y轴和Z轴的坐标值。
实际安装过程中,天线平台中心和惯导设备中心有可能不完全重合,存在安装误差,使得获取的飞机的姿态信息与实际值有偏差。为了尽量减小测量的姿态信息与实际姿态信息的偏差,可以按照如下公式对飞机的姿态信息进行校正:
YAW=YAW0+eyaw
PITCH=PITCH0+epitch
ROLL=ROLL0+eroll
其中,YAW0、PITCH0、ROLL0分别表示飞机惯导设备测得的飞机的偏航角、俯仰角和横滚角;eyaw、epitch、eroll分别表示飞机惯导设备在天线平台的偏航角安装误差值、俯仰角安装误差值和横滚角安装误差值。
S3、基于目标的第四坐标确定雷达天线的方位角和俯仰角,依据所述方位角和俯仰角调整雷达天线的指向。
实际安装过程中,天线平台中心和惯导设备中心有可能不完全重合、存在安装误差。为了尽量减小该安装误差对测量结果准确性的影响,在本发明的一些实施例子中,步骤S3中基于目标的第四坐标确定雷达天线的方位角和俯仰角包括:
S31、根据雷达在飞机上的安装误差对目标的第四坐标进行校正,得到目标的第五坐标;
S32、根据第五坐标确定雷达天线的方位角和俯仰角。
优选地,目标的第五坐标为:
[x5,y5,z5]'=[x4+x0,y4+y0,z4+z0]'
其中,x5、y5、z5分别表示校正安装误差后的目标在机体向下坐标系中的X轴、Y轴和Z轴的坐标值;x0、y0、z0分别表示在机体向下坐标系的X轴、Y轴和Z轴方向上的飞机惯导设备在天线平台的安装误差值。
进一步优选地,雷达天线的方位角azT为:
雷达天线的俯仰角pitchT为:
其中,a表示在WGS-84地球模型中椭球的长半轴。
与现有技术相比,本发明飞机和目标的位置信息和姿态信息获得目标相对于天线平台的位置关系,从而实时调整天线波束指向,保证被测量目标始终处于天线主波束中心。这样在机载下视测量目标时,无需人工手动操作伺服控制杆来对准目标,可由程序自动控制完成,提高机载下视测量的效率和精度。
虽然参照示例性实施方式对本发明进行了描述,但是应当理解,本发明并不局限于文中详细描述和示出的具体实施方式,在不偏离权利要求书所限定的范围的情况下,本领域技术人员可以对所述示例性实施方式做出各种改变。
Claims (8)
1.一种机载下视测量目标的自动跟踪方法,其特征在于包括如下步骤:
S1、获取飞机和目标的位置信息,根据所述位置信息确定将地心坐标系的坐标原点移动至飞机处后的目标的第一坐标;
S2、获取飞机的姿态信息,根据所述姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标;
S3、基于目标的第四坐标确定雷达天线的方位角和俯仰角,依据所述方位角和俯仰角调整雷达天线的指向;
其中,
所述位置信息包括:经度、纬度和高度,所述高度是指当前位置到地表的距离;所述姿态信息包括偏航角、俯仰角和横滚角;
地心坐标系中,Z轴指向地球北极,X轴指向格林威治子午面与地球赤道的交点,Y轴垂直XOZ面;
机体向下坐标系中,X轴指向飞机的机头方向,Y轴指向飞机的机身右侧,Z轴垂直机体向下;
其中,步骤S2中根据所述姿态信息将目标的第一坐标转换成机体向下坐标系中的第四坐标包括:
S21、将目标的第一坐标转换成东北天坐标系中的第二坐标;
S22、将目标在东北天坐标系中的第二坐标转换成东南地坐标系中的第三坐标;
S23、根据所述姿态信息,将目标的第三坐标转换成机体向下坐标系中的第四坐标;
其中,东南地坐标系中,X轴指向正东方向,Y轴指向正南方向,Z轴指向地心;东北天坐标系中,X轴指向正东方向,Y轴指向正北方向,Z轴指向地心反向;
其中,目标的第二坐标为:
[x2,y2,z2]'=M1M0[x1,y1,z1]'
其中,
<mrow>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
lon表示飞机的经度,lat表示飞机的纬度;x1、y1、z1分别表示将地心坐标系的坐标原点移动至飞机处后目标在X轴、Y轴和Z轴的坐标值;
x2、y2、z2分别表示目标在东北天坐标系中的X轴、Y轴和Z轴的坐标值。
2.如权利要求1所述的自动跟踪方法,其中,目标的第三坐标为:
[x3,y3,z3]'=M180[x2,y2,z2]'
其中,
<mrow>
<msub>
<mi>M</mi>
<mn>180</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
x3、y3、z3分别表示目标在东南地坐标系中的X轴、Y轴和Z轴的坐标值。
3.如权利要求2所述的自动跟踪方法,其中,目标的第四坐标为:
[x4,y4,z4]′=MROLLMPITCHMYAW[x3,y3,z3]′
其中,
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>Y</mi>
<mi>A</mi>
<mi>W</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>90</mn>
<mo>-</mo>
<mi>Y</mi>
<mi>A</mi>
<mi>W</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>90</mn>
<mo>-</mo>
<mi>Y</mi>
<mi>A</mi>
<mi>W</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>90</mn>
<mo>-</mo>
<mi>Y</mi>
<mi>A</mi>
<mi>W</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>90</mn>
<mo>-</mo>
<mi>Y</mi>
<mi>A</mi>
<mi>W</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>P</mi>
<mi>I</mi>
<mi>T</mi>
<mi>C</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>P</mi>
<mi>I</mi>
<mi>T</mi>
<mi>C</mi>
<mi>H</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>P</mi>
<mi>I</mi>
<mi>T</mi>
<mi>C</mi>
<mi>H</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>P</mi>
<mi>I</mi>
<mi>T</mi>
<mi>C</mi>
<mi>H</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>P</mi>
<mi>I</mi>
<mi>T</mi>
<mi>C</mi>
<mi>H</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>R</mi>
<mi>O</mi>
<mi>L</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>R</mi>
<mi>O</mi>
<mi>L</mi>
<mi>L</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>R</mi>
<mi>O</mi>
<mi>L</mi>
<mi>L</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>R</mi>
<mi>O</mi>
<mi>L</mi>
<mi>L</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>R</mi>
<mi>O</mi>
<mi>L</mi>
<mi>L</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
x4、y4、z4分别表示目标在机体向下坐标系中的X轴、Y轴和Z轴的坐标值;YAW表示飞机的偏航角,PITCH表示飞机的俯仰角,ROLL表示飞机的横滚角。
4.如权利要求3所述的自动跟踪方法,其中,步骤S3中基于目标的第四坐标确定雷达天线的方位角和俯仰角包括:
S31、根据雷达在飞机上的安装误差对目标的第四坐标进行校正,得到目标的第五坐标;
S32、根据第五坐标确定雷达天线的方位角和俯仰角。
5.如权利要求4所述的自动跟踪方法,其中,目标的第五坐标为:
[x5,y5,z5]'=[x4+x0,y4+y0,z4+z0]'
其中,x5、y5、z5分别表示校正安装误差后的目标在机体向下坐标系中的X轴、Y轴和Z轴的坐标值;x0、y0、z0分别表示在机体向下坐标系的X轴、Y轴和Z轴方向上的飞机惯导设备在天线平台的安装误差值。
6.如权利要求5所述的自动跟踪方法,其中,雷达天线的方位角azT为:
雷达天线的俯仰角pitchT为:
<mrow>
<mi>pitchT</mi>
<mo>=</mo>
<mi>a</mi>
<mi>tan</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>x</mi>
<mn>5</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mn>5</mn>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>z</mi>
<mn>5</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,a表示在WGS-84地球模型中椭球的长半轴。
7.如权利要求6所述的自动跟踪方法,其中,飞机的姿态信息为:
YAW=YAW0+eyaw
PITCH=PITCH0+epitch
ROLL=ROLL0+eroll
其中,YAW0、PITCH0、ROLL0分别表示飞机惯导设备测得的飞机的偏航角、俯仰角和横滚角;eyaw、epitch、eroll分别表示飞机惯导设备在天线平台的偏航角安装误差值、俯仰角安装误差值和横滚角安装误差值。
8.如权利要求1所述的自动跟踪方法,其中,飞机在地心坐标系的坐标为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mi>a</mi>
<mi>l</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mo>+</mo>
<mi>a</mi>
<mi>l</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>z</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>a</mi>
<mi>l</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
目标在地心坐标系的坐标为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>T</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
<mo>+</mo>
<msub>
<mi>alt</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>lat</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>lon</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>T</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
<mo>+</mo>
<msub>
<mi>alt</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>lat</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>lon</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>z</mi>
<mi>T</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>alt</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>lat</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,
<mrow>
<mi>N</mi>
<mo>=</mo>
<mi>a</mi>
<mo>/</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mi>a</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
<mo>=</mo>
<mi>a</mi>
<mo>/</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>lat</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
a表示在WGS-84地球模型中椭球的长半轴,a=6378137m;偏心率平方e2=0.0066944;x、y、z分别表示飞机在地心坐标系中的X轴、Y轴和Z轴的坐标值;xT、yT、zT分别表示目标在地心坐标系中的X轴、Y轴和Z轴的坐标值;lon表示飞机的经度,lat表示飞机的纬度,alt表示飞机的高度;lonT表示目标的经度,latT表示目标的纬度,altT表示目标的高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610350347.8A CN105891821B (zh) | 2016-05-24 | 2016-05-24 | 一种机载下视测量目标的自动跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610350347.8A CN105891821B (zh) | 2016-05-24 | 2016-05-24 | 一种机载下视测量目标的自动跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105891821A CN105891821A (zh) | 2016-08-24 |
CN105891821B true CN105891821B (zh) | 2018-04-17 |
Family
ID=56717894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610350347.8A Active CN105891821B (zh) | 2016-05-24 | 2016-05-24 | 一种机载下视测量目标的自动跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105891821B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108061477B (zh) * | 2016-11-08 | 2019-08-13 | 北京机电工程研究所 | 一种导引头与惯组系统间相对安装误差校正方法 |
CN108037523A (zh) * | 2017-10-26 | 2018-05-15 | 北京航空航天大学 | 一种应用于无人机的电子辅助波束对准方法 |
CN108896038B (zh) * | 2018-05-15 | 2022-03-25 | 中国科学院遥感与数字地球研究所 | 月基光学传感器成像方法 |
CN109725306B (zh) * | 2018-11-16 | 2023-02-10 | 中国西安卫星测控中心 | 一种飞行目标弹道确定方法 |
CN110361706B (zh) * | 2019-07-02 | 2021-03-12 | 中国人民解放军陆军炮兵防空兵学院郑州校区 | 一种雷达天线阵面法线方位角测量装置及方法 |
CN110440791A (zh) * | 2019-08-14 | 2019-11-12 | 北京环境特性研究所 | 一种机载雷达rcs测量航线设计方法和装置 |
CN110808447A (zh) * | 2019-10-24 | 2020-02-18 | 迪泰(浙江)通信技术有限公司 | 基于三轴动态跟踪技术的船载卫星天线系统 |
CN111238439B (zh) * | 2020-02-14 | 2022-06-10 | 天津时空经纬测控技术有限公司 | 角度偏差测量系统 |
CN112504277A (zh) * | 2020-11-20 | 2021-03-16 | 中国电子科技集团公司第五十四研究所 | 一种基于无人机数据链的应急导航方法 |
CN113156418B (zh) * | 2021-05-08 | 2023-02-24 | 中国科学院空天信息创新研究院 | 基于蒙特卡洛仿真的雷达目标跟踪精度预测方法 |
CN115017747B (zh) * | 2022-08-09 | 2022-10-21 | 安胜(天津)飞行模拟系统有限公司 | 一种eVTOL飞机的多倍速仿真重定位运行方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1906206A1 (en) * | 2006-09-28 | 2008-04-02 | Mazda Motor Corporation | Obstacle detecting apparatus, method and computer program product for vehicle |
CN102967851A (zh) * | 2012-11-14 | 2013-03-13 | 电子科技大学 | 一种双基地sar的空间同步方法 |
CN103438887A (zh) * | 2013-09-18 | 2013-12-11 | 上海海事大学 | 用于移动机器人定位与环境重构的绝对坐标获取方法 |
CN103471563A (zh) * | 2013-09-27 | 2013-12-25 | 重庆大学 | 分布式相控阵天线的子阵波束指向角度校正方法 |
CN103777201A (zh) * | 2012-10-19 | 2014-05-07 | 中国航天科工集团第二研究院二〇七所 | 基于gps数据的机载sar运动补偿方法 |
CN104459608A (zh) * | 2014-11-27 | 2015-03-25 | 中国船舶重工集团公司第七二四研究所 | 一种舰船纵横摇综合矢量坐标变换天线波束指向修正方法 |
CN104459750A (zh) * | 2014-12-11 | 2015-03-25 | 山东大学 | 一种基于gps/ins的动态指向方法 |
-
2016
- 2016-05-24 CN CN201610350347.8A patent/CN105891821B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1906206A1 (en) * | 2006-09-28 | 2008-04-02 | Mazda Motor Corporation | Obstacle detecting apparatus, method and computer program product for vehicle |
CN103777201A (zh) * | 2012-10-19 | 2014-05-07 | 中国航天科工集团第二研究院二〇七所 | 基于gps数据的机载sar运动补偿方法 |
CN102967851A (zh) * | 2012-11-14 | 2013-03-13 | 电子科技大学 | 一种双基地sar的空间同步方法 |
CN103438887A (zh) * | 2013-09-18 | 2013-12-11 | 上海海事大学 | 用于移动机器人定位与环境重构的绝对坐标获取方法 |
CN103471563A (zh) * | 2013-09-27 | 2013-12-25 | 重庆大学 | 分布式相控阵天线的子阵波束指向角度校正方法 |
CN104459608A (zh) * | 2014-11-27 | 2015-03-25 | 中国船舶重工集团公司第七二四研究所 | 一种舰船纵横摇综合矢量坐标变换天线波束指向修正方法 |
CN104459750A (zh) * | 2014-12-11 | 2015-03-25 | 山东大学 | 一种基于gps/ins的动态指向方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105891821A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105891821B (zh) | 一种机载下视测量目标的自动跟踪方法 | |
CN105698762B (zh) | 一种单机航迹上基于不同时刻观测点的目标快速定位方法 | |
CN105184776B (zh) | 目标跟踪方法 | |
CN106094877B (zh) | 一种无人机着陆导航系统及控制方法 | |
CN109753076A (zh) | 一种无人机视觉追踪实现方法 | |
CN104296751B (zh) | 一种多星敏感器构型布局设计方法 | |
CN106871927A (zh) | 一种无人机光电吊舱安装误差标校方法 | |
CN106595668A (zh) | 一种用于光电吊舱的无源定位算法 | |
CN110470304B (zh) | 一种基于无人机光电平台的高精度目标定位测速方法 | |
CN107505615A (zh) | 星载sar滑动聚束成像模式星体姿态设计方法 | |
CN108427427A (zh) | 一种航天器对地表定向目标姿态角计算方法 | |
CN108061477B (zh) | 一种导引头与惯组系统间相对安装误差校正方法 | |
DE102020113154A1 (de) | Virtueller sensor und kompasskalibrierung | |
CN103955224B (zh) | 一种用于相对运动视线跟踪的姿态控制方法 | |
CN106777489A (zh) | 无人机载光电稳定转台跟踪态建模仿真方法 | |
CN106155081A (zh) | 一种旋翼无人机大范围目标监测和精确定位方法 | |
CN110134134B (zh) | 一种无人机悬停状态下的测风方法 | |
CN107831515A (zh) | 水下定位方法和系统 | |
CN107168355A (zh) | 无人飞行器航线确定方法和装置 | |
CN110161578A (zh) | 一种海洋深部流动地磁测量系统及方法 | |
CN108657467B (zh) | 一种利用虚拟太阳矢量的航天器偏航机动控制方法及系统 | |
CN112945242B (zh) | 一种卫星在轨自主规划任务最佳时间与姿态的方法 | |
CN103760562A (zh) | 一种机载圆迹合成孔径雷达航线的获取方法 | |
CN105739542B (zh) | 一种基于两轴机载雷达伺服平台的兴趣点跟踪控制方法 | |
CN108398956A (zh) | 一种用于船载平台的无人机平衡系统及其控制方法 |
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 |