CN104596515B - 一种惯导/激光测速仪组合导航的半实物仿真方法 - Google Patents
一种惯导/激光测速仪组合导航的半实物仿真方法 Download PDFInfo
- Publication number
- CN104596515B CN104596515B CN201410539221.6A CN201410539221A CN104596515B CN 104596515 B CN104596515 B CN 104596515B CN 201410539221 A CN201410539221 A CN 201410539221A CN 104596515 B CN104596515 B CN 104596515B
- Authority
- CN
- China
- Prior art keywords
- laser velocimeter
- inertial navigation
- carrier
- velocity
- gyro
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/50—Systems of measurement based on relative movement of target
- G01S17/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/86—Combinations of lidar systems with systems other than lidar, radar or sonar, e.g. with direction finders
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Manufacturing & Machinery (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
Abstract
一种惯导/激光测速仪组合导航的半实物仿真方法,有七大步骤:一、将惯导/激光测速仪安装载体上并上电启动;二、装订初始参数至导航计算机;三、载体静止,采集保存1小时陀螺和加速度计输出数据;求其平均值并分别与平均值作差,得到陀螺和加速度计实际噪声;四、载体静止,惯导进行5分钟静态初始对准;五、载体运动,惯导计算,采集并保存运动后50秒内惯导输出速度值、姿态矩阵等数据;六、把速度投影到激光测速仪速度方向上得参考速度;将其速度与参考速度作差,得其输出速度的噪声值;七、设定仿真轨迹,得姿态矩阵、速度分量及速度标准值;将各相关数据与步骤3、6中噪声数据融合,得到器件的模拟数据进行惯导/激光测速仪组合导航仿真。
Description
技术领域:
本发明涉及一种惯导/激光测速仪组合导航的半实物仿真方法,属于惯性导航技术领域。
背景技术:
激光测速仪根据多普勒频移可以实时测量运动载体相对地面移动的速度,且其测量误差不随时间积累,与惯导***组合能优势互补,能实现全自主、高精度导航定位。
在进行惯导/激光测速仪组合导航研究时,一般是设定仿真轨迹(包含载体的运动速度、位置、姿态),利用轨迹发生器生成加速度计、陀螺的标准值,由设定仿真轨迹中载体的运动速度等信息计算得到激光测速仪的速度标准值,并给定激光测速仪的标度误差、安装角及安装角误差以及噪声,给定加速度计和陀螺的标度误差、零偏、失准角误差以及噪声,给定惯导的安装角及安装角误差,将以上数据融合得到用于组合导航仿真的器件模拟数据,然后采用这些模拟数据进行组合导航的仿真计算与分析。其中,给定的仿真噪声不能全面反映激光测速仪、加速度计、陀螺的噪声特性,导致仿真结果不能较全面地反映实际情况。为了更好地进行惯导/激光测速仪组合导航研究,本文提出了一种惯导/激光测速仪组合导航的半实物仿真方法。
发明内容
1、发明目的:本发明的目的是提供了一种惯导/激光测速仪组合导航的半实物仿真方法,它能够更好地模拟惯导/激光测速仪组合导航过程。
2、技术方案:
本发明一种惯导/激光测速仪组合导航的半实物仿真方法,该方法具体步骤如下:
步骤1、将惯导/激光测速仪组合***安装到载体上,并上电启动。
步骤2、装订初始参数(包括初始的经度、纬度、高度、激光测速仪标度值、激光测速仪安装角、惯导安装角)至导航计算机。
步骤3、保持载体静止,采集保存1小时陀螺和加速度计的输出数据。求取陀螺和加速度计的输出数据的平均值,将采集的陀螺和加速度计数据分别与其平均值作差,得到陀螺和加速度计的实际噪声。
步骤4、载体保持静止,惯导进行5分钟静态初始对准。
步骤5、完成对准后载体开始运动,运动过程中惯导进行惯性导航计算,采集并保存开始运动后50秒内惯导输出的速度值、惯导输出的姿态矩阵、激光测速仪输出的速度。
步骤6、惯导输出的速度为东北天坐标系下的速度,由惯导输出姿态矩阵将东北天坐标系下的速度转换到惯组坐标系下,再根据惯导与载体、激光测速仪与载体之间的位置关系,把惯组坐标系下的速度投影到激光测速仪的速度方向上得到参考速度。将激光测速仪的速度与参考速度作差,得到激光测速仪输出速度的噪声值。
步骤7、设定仿真轨迹(包含载体的运动速度、位置、姿态)。由设定仿真轨迹中载体的姿态得到仿真载体的姿态矩阵,根据得到的姿态矩阵和仿真轨迹中载体的运动速度得到载体在载体系下的速度分量,再由载体的速度分量和给定的激光测速仪的安装角、惯导安装角计算得到激光测速仪的速度标准值;利用轨迹发生器生成加速度计的标准值、陀螺的标准值;并给定加速度计和陀螺的标度误差、零偏、失准角误差;给定激光测速仪的标度误差、安装角误差;给定惯导的安装角误差。将这些数据与步骤3、步骤6中得到的噪声数据融合,得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真。
其中,步骤3中所述的“得到陀螺、加速度计的实际噪声”,其实现过程说明如下:
定义wg为陀螺噪声,wa为加速度计噪声,采集陀螺x、y、z三个轴输出的数据分别为ωx(1)、ωx(2)…ωx(n),ωy(1)、ωy(2)…ωy(n),ωz(1)、ωz(2)…ωz(n),平均值分别为采集加速度计x、y、z三个轴输出的数据分别为ax(1)、ax(2)…ax(n),ay(1)、ay(2)…ay(n),az(1)、az(2)…az(n),平均值分别为
陀螺噪声计算公式如下:
加速度计噪声计算公式如下:
其中,步骤6中所述的“得到激光测速仪输出速度的噪声值”,其实现过程说明如下:
载体开始运动,采集得到了50秒内激光测速仪的速度输出惯导输出东北天坐标系下速度惯导输出的姿态矩阵导航计算机装订的激光测速仪的安装角为θ,惯导的安装角为α、β、γ;惯导在惯组坐标系下的速度输出分别为Vx(1)、Vx(2)…Vx(n),Vy(1)、Vy(2)…Vy(n),Vz(1)、Vz(2)…Vz(n)。惯导在惯组坐标系下的速度投影到激光测速仪的速度方向上的参考速度为惯导在东北天坐标系下输出的速度转到惯组坐标系下,计算公式如下:
惯导在惯组坐标系下的速度投影到激光测速仪的速度方向上的参考速度为Vm(i),计算公式如下:
Vm(i)=[Vx(i)cosα+Vy(i)cosβ+Vz(i)cosγ]cosθ(i=1...n)
激光测速仪输出的速度噪声值为wL(i),计算公式如下:
其中,步骤7中所述的“得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真”,其实现过程说明如下:
设定仿真轨迹(包含载体的运动速度、位置、姿态)。由设定仿真轨迹中载体的姿态得到仿真载体的姿态矩阵,根据得到的姿态矩阵和仿真轨迹中载体的运动速度得到载体在载体系下的速度分量,再由载体的速度分量和给定的激光测速仪的安装角、惯导安装角计算得到激光测速仪的速度标准值;利用轨迹发生器生成加速度计的标准值、陀螺的标准值;并给定加速度计和陀螺的标度误差、零偏、失准角误差;给定激光测速仪的标度误差、安装角误差;给定惯导的安装角误差。将这些数据与步骤3、步骤6中得到的噪声数据融合,得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真。
定义:仿真轨迹中导航坐标系下载体运动速度向量为Vn,在导航坐标系下三个方向的分量分别为仿真轨迹中载体的运动速度在载体系下x、y、z三个轴的分量分别为由仿真轨迹中载体的姿态得到的载体姿态矩阵为给定激光测速仪的安装角θf,激光测速仪安装角误差为δθf,惯导的安装角分别为αb、βb、γb以及惯导的安装角误差分别为δαb、δβb、δγb,激光测速仪的标度误差值δKL;加速度计x、y、z三个轴的标准值分别为atx(1)、atx(2)…atx(n),aty(1)、aty(2)…aty(n),atz(1)、atz(2)…atz(n);加速度计x、y、z三个轴的标度误差分别为δKax、δKay、δKaz;加速度计的失准角误差为δMaxy、δMaxz、δMayx、δMayz、δMazx、δMazy;加速度计x、y、z三个轴的零偏分别为δBax、δBay、δBaz;陀螺x、y、z三个轴的标准值分别为ωtx(1)、ωtx(2)…ωtx(n),ωty(1)、ωty(2)…ωty(n),ωtz(1)、ωtz(2)…ωtz(n);陀螺x、y、z三个轴的标度误差分别为δKgx、δKgy、δKgz;陀螺的失准角误差为δMgxy、δMgxz、δMgyx、δMgyz、δMgzx、δMgzy;陀螺x、y、z三个轴的零偏分别为δBgx、δBgy、δBgz。陀螺x、y、z三个轴的噪声分别为wgx(1)、wgx(2)…wgx(n),wgy(1)、wgy(2)…wgy(n),wgz(1)、wgz(2)…wgz(n);加速度计x、y、z三个轴的噪声分别为wax(1)、wax(2)…wax(n),way(1)、way(2)…way(n),waz(1)、waz(2)…waz(n);激光测速仪的速度噪声为wL(1)、wL(2)…wL(n)。
仿真轨迹中载体在导航系下的速度转换到载体系下x、y、z三个轴的速度计算如下:
计算得到激光测速仪的速度标准值为计算公式如下:
由计算得到激光测速仪的速度标准值、给定的激光测速仪的标度误差值、安装角误差、惯导安装角误差以及提取的激光测速仪的速度噪声,综合得到用于组合导航半实物仿真的激光测速仪的速度计算公式如下:
由计算得到的陀螺、加速度计标准值,给定的陀螺、加速度计零偏、标度误差、失准角误差以及提取的陀螺、加速度计的噪声,综合得到用于组合导航半实物仿真的陀螺、加速度计数据计算公式如下:
采用上述得到的激光测速仪模拟数据加速度计模拟数据af、陀螺模拟数据ωf,结合由仿真轨迹中给定的载体初始位置、初始姿态、初始速度,进行组合导航半实物仿真。
3、优点及功效:该方法的优点是用激光测速仪、陀螺、加速度计的实际噪声替代仿真噪声,然后进行惯导/激光测速仪组合导航的半实物仿真,为研究惯导/激光测速仪组合导航提供了更好的支持。
附图说明
图1为本发明噪声提取框图
图2为本发明惯导/激光测速仪组合导航的半实物仿真方法框图;
图3为本发明惯导/激光测速仪组合导航的半实物仿真方法流程图;
图1中符号说明如下:
惯导输出的东、北、天方向的速度
Vx、Vy、Vz:惯导在惯组坐标系下x、y、z三个方向的速度分量
激光测速仪输出的速度
Vm:激光测速仪测速方向上的参考速度
a:加速度计输出值
加速度计输出值的平均值
ω:陀螺输出值
陀螺输出值的平均值
wL:激光测速仪的噪声
wa:加速度计的噪声
wg:陀螺的噪声
图2中符号说明如下:
计算得到的激光测速仪的速度标准值
δKL:给定的激光测速仪的标度误差
ωt:轨迹发生器生成的陀螺数据标准值
at:轨迹发生器生成的加速度计数据标准值
δBg:给定的陀螺零偏
δBa:给定的加速度计零偏
ωc:包含标度误差、安装角误差的陀螺数据
ac:包含标度误差、安装角误差的加速度计数据
用于组合导航半实物仿真的激光测速仪的速度值
ωf:用于组合导航半实物仿真的陀螺数据
af:用于组合导航半实物仿真的加速度计数据
P:仿真轨迹中设定的载***置
A:仿真轨迹中设定的载体姿态
V:仿真轨迹中设定的载体速度
P0:仿真轨迹中设定的载体初始位置
A0:仿真轨迹中设定的载体初始姿态
V0:仿真轨迹中设定的载体初始速度
wL:激光测速仪的速度噪声
wa:加速度计噪声
wg:陀螺噪声
具体实施方式
见图1—图3,本发明一种惯导/激光测速仪组合导航的半实物仿真方法,该方法具体步骤如下:
步骤1、将惯导/激光测速仪组合***安装到载体上,并上电启动。
步骤2、装订初始参数(包括初始的经度、纬度、高度、激光测速仪标度值、激光测速仪安装角、惯导安装角)至导航计算机。
步骤3、保持载体静止,采集保存1小时陀螺和加速度计的输出数据。求取陀螺和加速度计的输出数据的平均值,将采集的陀螺和加速度计数据分别与其平均值作差,得到陀螺和加速度计的实际噪声。
步骤4、载体保持静止,惯导进行5分钟静态初始对准。
步骤5、完成对准后载体开始运动,运动过程中惯导进行惯性导航计算,采集并保存开始运动后50秒内惯导输出的速度值、惯导输出的姿态矩阵、激光测速仪输出的速度。
步骤6、惯导输出的速度为东北天坐标系下的速度,由惯导输出姿态矩阵将东北天坐标系下的速度转换到惯组坐标系下,再根据惯导与载体、激光测速仪与载体之间的位置关系,把惯组坐标系下的速度投影到激光测速仪的速度方向上得到参考速度。将激光测速仪的速度与参考速度作差,得到激光测速仪输出速度的噪声值。
步骤7、设定仿真轨迹(包含载体的运动速度、位置、姿态)。由设定仿真轨迹中载体的姿态得到仿真载体的姿态矩阵,根据得到的姿态矩阵和仿真轨迹中载体的运动速度得到载体在载体系下的速度分量,再由载体的速度分量和给定的激光测速仪的安装角、惯导安装角计算得到激光测速仪的速度标准值;利用轨迹发生器生成加速度计的标准值、陀螺的标准值;并给定加速度计和陀螺的标度误差、零偏、失准角误差;给定激光测速仪的标度误差、安装角误差;给定惯导的安装角误差。将这些数据与步骤3、步骤6中得到的噪声数据融合,得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真。
其中,步骤3中所述的“得到陀螺、加速度计的实际噪声”,其实现过程说明如下:
定义wg为陀螺噪声,wa为加速度计噪声,采集陀螺x、y、z三个轴输出的数据分别为ωx(1)、ωx(2)…ωx(n),ωy(1)、ωy(2)…ωy(n),ωz(1)、ωz(2)…ωz(n),平均值分别为采集加速度计x、y、z三个轴输出的数据分别为ax(1)、ax(2)…ax(n),ay(1)、ay(2)…ay(n),az(1)、az(2)…az(n),平均值分别为
陀螺噪声计算公式如下:
加速度计噪声计算公式如下:
其中,步骤6中所述的“得到激光测速仪输出速度的噪声值”,其实现过程说明如下:
载体开始运动,采集得到了50秒内激光测速仪的速度输出惯导输出东北天坐标系下速度惯导输出的姿态矩阵导航计算机装订的激光测速仪的安装角为θ,惯导的安装角为α、β、γ;惯导在惯组坐标系下的速度输出分别为Vx(1)、Vx(2)…Vx(n),Vy(1)、Vy(2)…Vy(n),Vz(1)、Vz(2)…Vz(n)。惯导在惯组坐标系下的速度投影到激光测速仪的速度方向上的参考速度为惯导在东北天坐标系下输出的速度转到惯组坐标系下,计算公式如下:
惯导在惯组坐标系下的速度投影到激光测速仪的速度方向上的参考速度为Vm(i),计算公式如下:
Vm(i)=[Vx(i)cosα+Vy(i)cosβ+Vz(i)cosγ]cosθ(i=1…n)
激光测速仪输出的速度噪声值为wL(i),计算公式如下:
其中,步骤7中所述的“得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真”,其实现过程说明如下:
设定仿真轨迹(包含载体的运动速度、位置、姿态)。由设定仿真轨迹中载体的姿态得到仿真载体的姿态矩阵,根据得到的姿态矩阵和仿真轨迹中载体的运动速度得到载体在载体系下的速度分量,再由载体的速度分量和给定的激光测速仪的安装角、惯导安装角计算得到激光测速仪的速度标准值;利用轨迹发生器生成加速度计的标准值、陀螺的标准值;并给定加速度计和陀螺的标度误差、零偏、失准角误差;给定激光测速仪的标度误差、安装角误差;给定惯导的安装角误差。将这些数据与步骤3、步骤6中得到的噪声数据融合,得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真。
定义:仿真轨迹中导航坐标系下载体运动速度向量为Vn,在导航坐标系下三个方向的分量分别为仿真轨迹中载体的运动速度在载体系下x、y、z三个轴的分量分别为由仿真轨迹中载体的姿态得到的载体姿态矩阵为给定激光测速仪的安装角θf,激光测速仪安装角误差为δθf,惯导的安装角分别为αb、βb、γb以及惯导的安装角误差分别为δαb、δβb、δγb,激光测速仪的标度误差值δKL;加速度计x、y、z三个轴的标准值分别为atx(1)、atx(2)…atx(n),aty(1)、aty(2)…aty(n),atz(1)、atz(2)…atz(n);加速度计x、y、z三个轴的标度误差分别为δKax、δKay、δKaz;加速度计的失准角误差为δMaxy、δMaxz、δMayx、δMayz、δMazx、δMazy;加速度计x、y、z三个轴的零偏分别为δBax、δBay、δBaz;陀螺x、y、z三个轴的标准值分别为ωtx(1)、ωtx(2)…ωtx(n),ωty(1)、ωty(2)…ωty(n),ωtz(1)、ωtz(2)…ωtz(n);陀螺x、y、z三个轴的标度误差分别为δKgx、δKgy、δKgz;陀螺的失准角误差为δMgxy、δMgxz、δMgyx、δMgyz、δMgzx、δMgzy;陀螺x、y、z三个轴的零偏分别为δBgx、δBgy、δBgz。陀螺x、y、z三个轴的噪声分别为wgx(1)、wgx(2)…wgx(n),wgy(1)、wgy(2)…wgy(n),wgz(1)、wgz(2)…wgz(n);加速度计x、y、z三个轴的噪声分别为wax(1)、wax(2)…wax(n),way(1)、way(2)…way(n),waz(1)、waz(2)…waz(n);激光测速仪的速度噪声为wL(1)、wL(2)…wL(n)。
仿真轨迹中载体在导航系下的速度转换到载体系下x、y、z三个轴的速度计算如下:
计算得到激光测速仪的速度标准值为计算公式如下:
由计算得到激光测速仪的速度标准值、给定的激光测速仪的标度误差值、安装角误差、惯导安装角误差以及提取的激光测速仪的速度噪声,综合得到用于组合导航半实物仿真的激光测速仪的速度计算公式如下:
由计算得到的陀螺、加速度计标准值,给定的陀螺、加速度计零偏、标度误差、失准角误差以及提取的陀螺、加速度计的噪声,综合得到用于组合导航半实物仿真的陀螺、加速度计数据计算公式如下:
采用上述得到的激光测速仪模拟数据加速度计模拟数据af、陀螺模拟数据ωf,结合由仿真轨迹中给定的载体初始位置、初始姿态、初始速度,进行组合导航半实物仿真。
Claims (4)
1.一种惯导/激光测速仪组合导航的半实物仿真方法,其特征在于:该方法具体步骤如下:
步骤1、将惯导/激光测速仪组合***安装到载体上,并上电启动;
步骤2、装订初始参数,包括初始的经度、纬度、高度、激光测速仪标度值、激光测速仪安装角、惯导安装角至导航计算机;
步骤3、保持载体静止,采集保存1小时陀螺和加速度计的输出数据;求取陀螺和加速度计的输出数据的平均值,将采集的陀螺和加速度计数据分别与其平均值作差,得到陀螺和加速度计的实际噪声;
步骤4、载体保持静止,惯导进行5分钟静态初始对准;
步骤5、完成对准后载体开始运动,运动过程中惯导进行惯性导航计算,采集并保存开始运动后50秒内惯导输出的速度值、惯导输出的姿态矩阵、激光测速仪输出的速度;
步骤6、惯导输出的速度为东北天坐标系下的速度,由惯导输出姿态矩阵将东北天坐标系下的速度转换到惯组坐标系下,再根据惯导与载体、激光测速仪与载体之间的位置关系,把惯组坐标系下的速度投影到激光测速仪的速度方向上得到参考速度;将激光测速仪的速度与参考速度作差,得到激光测速仪输出速度的噪声值;
步骤7、设定仿真轨迹,包含载体的运动速度、位置、姿态,由设定仿真轨迹中载体的姿态得到仿真载体的姿态矩阵,根据得到的姿态矩阵和仿真轨迹中载体的运动速度得到载体在载体系下的速度分量,再由载体的速度分量和给定的激光测速仪的安装角、惯导安装角计算得到激光测速仪的速度标准值;利用轨迹发生器生成加速度计的标准值、陀螺的标准值;并给定加速度计和陀螺的标度误差、零偏、失准角误差;给定激光测速仪的标度误差、安装角误差;给定惯导的安装角误差;将这些数据与步骤3、步骤6中得到的噪声数据融合,得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真。
2.根据权利要求1所述的一种惯导/激光测速仪组合导航的半实物仿真方法,其特征在于:步骤3中所述的“得到陀螺、加速度计的实际噪声”,其实现过程说明如下:
定义wg为陀螺噪声,wa为加速度计噪声,采集陀螺x、y、z三个轴输出的数据分别为ωx(1)、ωx(2)…ωx(n),ωy(1)、ωy(2)…ωy(n),ωz(1)、ωz(2)…ωz(n),平均值分别为采集加速度计x、y、z三个轴输出的数据分别为ax(1)、ax(2)…ax(n),ay(1)、ay(2)…ay(n),az(1)、az(2)…az(n),平均值分别为
陀螺噪声计算公式如下:
其中,i=1...n;
加速度计噪声计算公式如下:
其中,i=1...n。
3.根据权利要求1所述的一种惯导/激光测速仪组合导航的半实物仿真方法,其特征在于:步骤6中所述的“得到激光测速仪输出速度的噪声值”,其实现过程说明如下:
载体开始运动,采集得到了50秒内激光测速仪的速度输出惯导输出东北天坐标系下速度惯导输出的姿态矩阵导航计算机装订的激光测速仪的安装角为θ,惯导的安装角为α、β、γ;惯导在惯组坐标系下的速度输出分别为Vx(1)、Vx(2)…Vx(n),Vy(1)、Vy(2)…Vy(n),Vz(1)、Vz(2)…Vz(n);惯导在惯组坐标系下的速度投影到激光测速仪的速度方向上的参考速度为惯导在东北天坐标系下输出的速度转到惯组坐标系下,计算公式如下:
其中,i=1...n;
惯导在惯组坐标系下的速度投影到激光测速仪的速度方向上的参考速度为Vm(i),计算公式如下:
Vm(i)=[Vx(i)cosα+Vy(i)cosβ+Vz(i)cosγ]cosθ;其中,i=1...n;
激光测速仪输出的速度噪声值为wL(i),计算公式如下:
其中,i=1...n。
4.根据权利要求1所述的一种惯导/激光测速仪组合导航的半实物仿真方法,其特征在于:步骤7中所述的“得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真”,其实现过程说明如下:
设定仿真轨迹,包含载体的运动速度、位置、姿态,由设定仿真轨迹中载体的姿态得到仿真载体的姿态矩阵,根据得到的姿态矩阵和仿真轨迹中载体的运动速度得到载体在载体系下的速度分量,再由载体的速度分量和给定的激光测速仪的安装角、惯导安装角计算得到激光测速仪的速度标准值;利用轨迹发生器生成加速度计的标准值、陀螺的标准值;并给定加速度计和陀螺的标度误差、零偏、失准角误差;给定激光测速仪的标度误差、安装角误差;给定惯导的安装角误差;将这些数据与步骤3、步骤6中得到的噪声数据融合,得到用于惯导/激光测速仪组合导航半实物仿真的器件模拟数据,然后采用这些模拟数据进行惯导/激光测速仪组合导航仿真;
定义:仿真轨迹中导航坐标系下载体运动速度向量为Vn,在导航坐标系下三个方向的分量分别为仿真轨迹中载体的运动速度在载体系下x、y、z三个轴的分量分别为由仿真轨迹中载体的姿态得到的载体姿态矩阵为给定激光测速仪的安装角θf,激光测速仪安装角误差为δθf,惯导的安装角分别为αb、βb、γb以及惯导的安装角误差分别为δαb、δβb、δγb,激光测速仪的标度误差值δKL;加速度计x、y、z三个轴的标准值分别为atx(1)、atx(2)…atx(n),aty(1)、aty(2)…aty(n),atz(1)、atz(2)…atz(n);加速度计x、y、z三个轴的标度误差分别为δKax、δKay、δKaz;加速度计的失准角误差为δMaxy、δMaxz、δMayx、δMayz、δMazx、δMazy;加速度计x、y、z三个轴的零偏分别为δBax、δBay、δBaz;陀螺x、y、z三个轴的标准值分别为ωtx(1)、ωtx(2)…ωtx(n),ωty(1)、ωty(2)…ωty(n),ωtz(1)、ωtz(2)…ωtz(n);陀螺x、y、z三个轴的标度误差分别为δKgx、δKgy、δKgz;陀螺的失准角误差为δMgxy、δMgxz、δMgyx、δMgyz、δMgzx、δMgzy;陀螺x、y、z三个轴的零偏分别为δBgx、δBgy、δBgz;陀螺x、y、z三个轴的噪声分别为wgx(1)、wgx(2)…wgx(n),wgy(1)、wgy(2)…wgy(n),wgz(1)、wgz(2)…wgz(n);加速度计x、y、z三个轴的噪声分别为wax(1)、wax(2)…wax(n),way(1)、way(2)…way(n),waz(1)、waz(2)…waz(n);激光测速仪的速度噪声为wL(1)、wL(2)…wL(n);
仿真轨迹中载体在导航系下的速度转换到载体系下x、y、z三个轴的速度计算如下:
其中,i=1...n;
计算得到激光测速仪的速度标准值为计算公式如下:
其中,i=1...n;
由计算得到激光测速仪的速度标准值、给定的激光测速仪的标度误差值、安装角误差、惯导安装角误差以及提取的激光测速仪的速度噪声,综合得到用于组合导航半实物仿真的激光测速仪的速度计算公式如下:
其中,i=1...n;
由计算得到的陀螺、加速度计标准值,给定的陀螺、加速度计零偏、标度误差、失准角误差以及提取的陀螺、加速度计的噪声,综合得到用于组合导航半实物仿真的陀螺、加速度计数据计算公式如下:
其中,i=1...n;
其中,i=1...n;
采用上述得到的激光测速仪模拟数据加速度计模拟数据af、陀螺模拟数据ωf,结合由仿真轨迹中给定的载体初始位置、初始姿态、初始速度,进行组合导航半实物仿真。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410539221.6A CN104596515B (zh) | 2014-10-13 | 2014-10-13 | 一种惯导/激光测速仪组合导航的半实物仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410539221.6A CN104596515B (zh) | 2014-10-13 | 2014-10-13 | 一种惯导/激光测速仪组合导航的半实物仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104596515A CN104596515A (zh) | 2015-05-06 |
CN104596515B true CN104596515B (zh) | 2017-06-06 |
Family
ID=53122451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410539221.6A Expired - Fee Related CN104596515B (zh) | 2014-10-13 | 2014-10-13 | 一种惯导/激光测速仪组合导航的半实物仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104596515B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101122780A (zh) * | 2007-09-04 | 2008-02-13 | 北京控制工程研究所 | 月球软着陆制导、导航与控制半物理仿真试验*** |
CN102707092A (zh) * | 2012-06-19 | 2012-10-03 | 北京航空航天大学 | 一种基于角速率台的单波束激光测速仪标定方法 |
CN102706367A (zh) * | 2012-06-19 | 2012-10-03 | 北京航空航天大学 | 一种用于组合导航的单波束激光测速仪精度测试与计算方法 |
CN103278163A (zh) * | 2013-05-24 | 2013-09-04 | 哈尔滨工程大学 | 一种基于非线性模型的sins/dvl组合导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2965928B1 (fr) * | 2010-10-08 | 2013-05-31 | Thales Sa | Systeme pour la determination de la vitesse air d'un aeronef |
-
2014
- 2014-10-13 CN CN201410539221.6A patent/CN104596515B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101122780A (zh) * | 2007-09-04 | 2008-02-13 | 北京控制工程研究所 | 月球软着陆制导、导航与控制半物理仿真试验*** |
CN102707092A (zh) * | 2012-06-19 | 2012-10-03 | 北京航空航天大学 | 一种基于角速率台的单波束激光测速仪标定方法 |
CN102706367A (zh) * | 2012-06-19 | 2012-10-03 | 北京航空航天大学 | 一种用于组合导航的单波束激光测速仪精度测试与计算方法 |
CN103278163A (zh) * | 2013-05-24 | 2013-09-04 | 哈尔滨工程大学 | 一种基于非线性模型的sins/dvl组合导航方法 |
Non-Patent Citations (1)
Title |
---|
基于激光测速仪的高精度定位定向技术;张小跃等;《强激光与粒子束》;20140331;第26卷(第3期);031012-1-031012-5 * |
Also Published As
Publication number | Publication date |
---|---|
CN104596515A (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105259787B (zh) | 一种组合导航半物理仿真测试同步控制方法 | |
CN105371844B (zh) | 一种基于惯性/天文互助的惯性导航***初始化方法 | |
CN101979277B (zh) | 卫星磁测磁控***的全实物验证平台与工作方法 | |
CN104374388B (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN101424534B (zh) | 惯性/重力组合导航半实物模拟装置 | |
CN106840196A (zh) | 一种捷联惯性导航计算机测试***及实现方法 | |
CN106885570A (zh) | 一种基于鲁棒sckf滤波的紧组合导航方法 | |
CN106643737A (zh) | 风力干扰环境下四旋翼飞行器姿态解算方法 | |
CN105136166B (zh) | 一种指定惯导位置精度的捷联惯导***误差模型仿真方法 | |
CN103017765B (zh) | 应用于微机械组合导航***的偏航角修正方法和修正装置 | |
CN106019323A (zh) | 卫星导航接收机仿真测试*** | |
CN103414451B (zh) | 一种应用于飞行器姿态估计的扩展卡尔曼滤波方法 | |
CN102997935B (zh) | 一种基于光学和惯性组合测量的自主gnc仿真试验*** | |
CN106441301A (zh) | 一种飞行器发射初始参数获取方法及*** | |
CN105242679A (zh) | 一种四旋翼飞行器的控制***设计方法 | |
CN103604430A (zh) | 一种基于边缘化ckf重力辅助导航的方法 | |
CN103454662B (zh) | 一种基于ckf的sins/北斗/dvl组合对准方法 | |
CN106767931A (zh) | 一种基于角加速度给定的无陀螺惯性导航***的验证方法 | |
CN108827345A (zh) | 一种基于杆臂挠曲变形补偿的机载武器传递对准方法 | |
CN105865432B (zh) | 一种针对陀螺仪多源噪声的混合滤波方法与测试平台 | |
CN104596540B (zh) | 一种惯导/里程计组合导航的半实物仿真方法 | |
CN105069237B (zh) | 硬件接口激励的惯性/卫星深组合导航***动态测试方法 | |
CN104729530B (zh) | 一种惯导/北斗组合导航的半实物仿真方法 | |
CN109029499A (zh) | 一种基于重力视运动模型的加速度计零偏迭代寻优估计方法 | |
CN110940336B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170606 Termination date: 20181013 |