CN101975955B - 一种gnss模拟器中的通用三维载体运动轨迹生成方法 - Google Patents
一种gnss模拟器中的通用三维载体运动轨迹生成方法 Download PDFInfo
- Publication number
- CN101975955B CN101975955B CN 201010296636 CN201010296636A CN101975955B CN 101975955 B CN101975955 B CN 101975955B CN 201010296636 CN201010296636 CN 201010296636 CN 201010296636 A CN201010296636 A CN 201010296636A CN 101975955 B CN101975955 B CN 101975955B
- Authority
- CN
- China
- Prior art keywords
- carrier
- motion
- coordinate system
- centerdot
- delta
- 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
- 230000033001 locomotion Effects 0.000 title claims abstract description 220
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004088 simulation Methods 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 230000001133 acceleration Effects 0.000 claims description 32
- 230000008859 change Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 description 9
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- HEFNNWSXXWATRW-UHFFFAOYSA-N Ibuprofen Chemical compound CC(C)CC1=CC=C(C(C)C(O)=O)C=C1 HEFNNWSXXWATRW-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种GNSS模拟器中的通用三维载体运动轨迹生成方法,它有六大步骤:一、预生成载体运动轨迹用户配置文件,设置载体运动分段与轨迹参数;二、在载体运动轨迹模拟开始时,以设置的初始位置为原点建立地理坐标系;三、载体运动分段为空间直线运动类型的实时轨迹模拟计算;四、载体运动分段为空间圆周运动类型的实时轨迹模拟计算;五、对于运动分段之间的衔接,以上一运动段的末状态为当前运动段的初状态;六、将载体即时位置转换到大地坐标系实时显示输出,并将载体即时位置、速度转换到ECEF坐标系支持模拟器进一步功能实现。本发明能建立通用的三维载体运动轨迹产生模型,它在卫星导航技术领域里具有较好的实用价值和广阔的应用前景。
Description
(一)技术领域
本发明涉及全球卫星导航***(GNSS)卫星信号模拟器中接收机载体的运动轨迹模拟产生方法,尤其涉及一种GNSS模拟器中的通用三维载体运动轨迹生成方法,属于卫星导航技术领域。
(二)背景技术
GNSS卫星信号模拟器通过产生卫星信号接收机载体的运动轨迹,实时计算出接收机的位置,并结合卫星的实时位置判断卫星可见性,计算卫星信号功率、伪距、多普勒、载波相位等接收机实时观测量,以此实现对导航数据的扩频码调制和载波调制,模拟产生实际应用场景中接收机接收到的卫星导航信号,成为卫星信号接收机研制开发与测试验证中的仿真信号源。
相比直接利用真实卫星信号,GNSS卫星信号模拟器能够提供精确可控、可复现的仿真环境以及非正常情况的测试条件,使接收机的研发效率得以保证。随着美国GPS***现代化和俄罗斯GLONASS***改进升级,以及欧洲Galileo***和我国北斗卫星导航***的建设(上述全球卫星导航***可通称GNSS***),GNSS模拟器持续受到工业及军事部门的关注。
国外已研制出多种型号的GNSS卫星信号模拟器,如,美国CAST公司的GPS模拟器、英国Spirent公司的GSS系列模拟器、德国IFEN公司的GNSS模拟器,但大都功能固定,所能模拟的导航卫星***、信号频点与结构、载体动态等严格受限。
模拟用户动态是GNSS模拟器的重要环节之一,精确模拟出接收机载体的运动模式,有利于对接收机的动态性能进行测量评定;而动态导航信号的模拟产生能够支持接收机的捕获跟踪性能测试,以及定位精度等性能的分析评估。接收机载体运动轨迹的模拟可以基于信号模拟器内置的载体运动模型,也可以实测轨迹数据文件的形式载入GNSS信号模拟器,但内置载体运动模型的实现方式更为灵活,在精确性与可重复性方面更为突出。
(三)发明内容
1、目的
在现有的GNSS模拟器产品中,有关接收机载体运动轨迹生成,通常首先划分若干载体类型,在载体类型之下再划分以载体运动特性区分的运动模式,而在运动模式之下支持用户进行与运动模式关联的运动轨迹参数自定义,依据上述流程最终完成载体轨迹的定义,模拟器据此生成载体运动轨迹。以Spirent GSS系列模拟器为例,其载体类型的用户选项分为:简单运动(仅包括平面圆周运动和平面矩形轨迹运动)、车辆、舰船、飞机。在上述载体类型的选项之下,用户需要进行与载体类型关联的运动模式选择,以飞机载体类型为例,依据不同飞行阶段的运动特性,其运动模式分为:爬升、加速直线、匀速直线、匀速转弯、加速转弯、盘旋等。以匀速转弯运动模式为例,用户可定义的运动轨迹参数包括:运动段内飞机航向角的变化量、向心加速度;以加速直线运动模式为例,用户可定义的运动轨迹参数包括:加速段的持续时间或持续距离,加速度或末速度。
扩展性是上述模拟器接收机载体运动轨迹生成技术的主要问题,即载体类型、运动模式、可定义的轨迹参数受到繁琐且特定的限制。由模拟器支持的接收机测试与试验,通常只是利用实际GNSS信号进行测试之前的关键步骤,从而面向特定的应用场景与典型的测试试验成为一般GNSS模拟器的特点之一,这是GNSS模拟器载体动态模式固定并采用概略轨迹模型的主要原因。对于仍在建设中的GNSS***(如,Galileo***,北斗卫星导航***),通常情况下利用实际的卫星信号并不可行,由于***建设周期的影响,实际信号测试可能远滞后于模拟器测试,从而对相应的GNSS模拟器提出了更高的要求,在***建成并具备运行能力之前更加全面的测试试验需要模拟器的载体运动轨迹模拟具备扩展能力、具备反映各种动态的通用性、具备模拟复杂轨迹的能力。
本发明的目的是为了克服现有GNSS模拟器载体运动轨迹生成技术的局限,提供一种GNSS模拟器中的通用三维载体运动轨迹生成方法,它能将任意复杂形式的载体运动轨迹分解为空间直线运动段与空间圆周运动段的组合进行近似,借助姿态描述及多个坐标系的相互转换,建立通用的三维载体运动轨迹产生模型。
2、技术方案
本发明的目的通过以下技术方案实现:一种GNSS模拟器中的通用三维载体运动轨迹生成方法,该方法具体步骤如下:
步骤一:预生成载体运动轨迹用户配置文件,配置文件定义了载体运动分段与轨迹参数设置;
关于载体运动分段的轨迹参数设置,对于直线运动类型,需要设置的轨迹参数包括:直线运动类型标识,运动段持续时间,东向加速度,北向加速度,天向加速度;对于圆周运动类型,需要设置的轨迹参数包括:圆周运动类型标识,运动段持续时间,运动方向(顺时针或逆时针运动),圆周半径,线加速度,圆周运动平面方向相对地理坐标系的方向角。
步骤二:在载体运动轨迹模拟开始时,以设置的初始位置(即经度、纬度和高程)为原点建立地理坐标系,初始位置(x0,y0,z0)即为(0,0,0);
地理坐标系的x轴(又称e轴)沿当地纬线切线指向东,y轴(又称n轴)沿当地经线切线指向北,z轴(又称u轴)沿当地地理垂线(地球参考椭球表面法线)指向天顶。对于动态范围在180km内的载体运动轨迹模拟,这个初始设置的地理坐标系可视为三个轴线方向维持不变,所引起的误差可忽略不计,则在此距离范围内初始的地理坐标系将被保持。在载体运动轨迹模拟结束之前,如果载体动态超出这一距离范围,则需要更新地理坐标系;
步骤三:获取载体运动第一分段的轨迹参数,如果为空间直线运动类型,在地理坐标系中,根据初速度(ve,0,vn,0,vu,0)、加速度(ae,an,au)的设置,可以确定载体的即时位置、速度;
运动段持续时间内,根据tk-1仿真历元时刻载体的位置(xk-1,yk-1,zk-1)、速度(ve,k-1,vn,k-1,vu,k-1)来计算tk时刻载***置(xk,yk,zk)、速度(ve,k,vn,k,vu,k)的方法如下:
(a)计算时间间隔Δt
Δt=tk-tk-1
(b)计算tk-1时刻到tk时刻的距离变化量(Δse,Δsn,Δsu)
(c)计算tk时刻载体速度(ve,k,vn,k,vu,k)
(d)计算tk时刻载***置(xk,yk,zk)
步骤四:如果载体运动第一分段为空间圆周运动类型,根据运动平面相对地理坐标系方向角γ(取运动平面方向与地理坐标系z轴的夹角,以z轴右侧为正角度)的设置,首先确定地理坐标系中圆周运动的二维平面,该平面的x轴即为地理坐标系x轴(e轴),y轴与x轴垂直;
其次,根据初速度(ve,0,vn,0,vu,0)、圆周半径R和线加速度a的设置,可以确定圆周运动平面内载体的即时位置、速度。运动段持续时间内,计算tk仿真历元时刻圆周运动平面内载***置速度(vx,k,vy,k)的方法如下:
(a)计算时间间隔Δt
Δt=tk-tk-1
(b)计算tk-1时刻到tk时刻的转角变化量Δφ
其中,|.|为绝对值算子,
而
(c)计算tk时刻运动平面内载体方位角φk
φk=φk-1+Δφ
而初始方位
(d)计算tk时刻运动平面内载体速度(vx,k,vy,k)
而|vk|=|vk-1|+a·Δt
而
最后,将运动平面内载体的位置和速度转换到地理坐标系中的相应运动状态参数,转换方法如下:
(b)计算tk时刻地理坐标系载体速度(ve,k,vn,k,vu,k)
(c)计算tk时刻地理坐标系载***置(xk,yk,zk)
步骤五:根据直线或圆周运动类型的用户设置,按照步骤三、步骤四,依次处理后续运动分段。对于运动分段之间的衔接,以上一运动段的末状态为当前运动段的初状态,即上一运动段的载体末位置(xm,ym,zm)、末速度(ve,m,vn,m,vu,m)分别成为当前运动段的初始位置(x0,y0,z0)、初始速度(ve,0,vn,0,vu,0),即
步骤六:在模拟器的每个仿真历元,将地理坐标系中的载体即时位置转换到大地坐标系,用于用户界面实时显示输出;并将载体即时位置、速度转换到地球直角坐标系(又称ECEF坐标系)以进一步支持卫星可见性判断,伪距、多普勒等信号状态参数计算等模拟器功能实现;
便于数值计算起见,考虑RN和RM的近似算法,即
而扁率
(b)计算tk时刻ECEF坐标系中相应的载***置坐标(xe,k,ye,k,ze,k)
(c)计算tk时刻ECEF坐标系中相应的载体速度(vx,k,vy,k,vz,k)
3、优点及功效
由上述技术方案可以看出,与国内外已有的GNSS模拟器载体运动轨迹生成方法相比,本发明提出的载体轨迹生成方法具有以下技术优点:
(1)载体运动轨迹模拟具有通用性,不区分载体类型,认为载体轨迹自身的动态特征即反映载体类型;以空间直线运动、空间圆周运动两种运动模式涵盖适用于GNSS模拟器的卫星信号接收机载体的运动模式。
(2)基于运动轨迹分段分解近似并引入姿态描述,具有适用于GNSS模拟器的复杂载体运动轨迹的模拟能力,对复杂轨迹的近似度仅取决于运动段划分的精细程度,支持载体动态细节的小尺度精确描述。
(3)采用载体运动轨迹用户配置文件方式,载体轨迹参数的配置灵活、便捷,支持用户可完全自定义的载体运动轨迹模拟产生,使载体轨迹模拟具有扩展性。
(四)附图说明
图1为本发明的通用三维载体运动轨迹生成方法流程框图
图2为本发明的载体运动轨迹用户配置文件格式示意图
图3为本发明的载体运动轨迹生成示例图
图中符号说明如下:
+e:东向 +n:北向 +u:天向
(五)具体实施方式
本发明一种适用于GNSS模拟器的通用三维载体运动轨迹生成方法,所述方法流程框图如图1所示。该方法通过GNSS模拟器的软件平台实现,在方法具体实施之前,基于GNSS模拟器的人机交互界面,选择所需仿真的导航***及频点,设置仿真起始时间和仿真时间长度,选取星历源文件,设置接收机载体动态,选择是否引入将折算为测距误差的各种误差项,等等。当设置接收机载体动态时,需要选择载体运动轨迹用户配置文件,从而进入本发明所述方法的流程。
该方法具体通过以下步骤实现:
步骤一:预生成载体运动轨迹用户配置文件,配置文件定义了载体运动分段与轨迹参数设置。有关用户配置文件的格式说明,如图2所示。
关于载体运动分段的轨迹参数设置,对于直线运动类型,需要设置的轨迹参数包括:直线运动类型标识,运动段持续时间,东向加速度,北向加速度,天向加速度;对于圆周运动类型,需要设置的轨迹参数包括:圆周运动类型标识,运动段持续时间,运动方向(顺时针或逆时针运动),圆周半径,线加速度,圆周运动平面方向相对地理坐标系的方向角。
步骤二:在载体运动轨迹模拟开始时,以设置的初始位置(即经度、纬度和高程)为原点建立地理坐标系,初始位置(x0,y0,z0)即为(0,0,0)。
地理坐标系的x轴(又称e轴)沿当地纬线切线指向东,y轴(又称n轴)沿当地经线切线指向北,z轴(又称u轴)沿当地地理垂线(地球参考椭球表面法线)指向天顶。对于动态范围在180km内的载体运动轨迹模拟,这个初始设置的地理坐标系可视为三个轴线方向维持不变,所引起的误差可忽略不计,则在此距离范围内初始的地理坐标系将被保持。在载体运动轨迹模拟结束之前,如果载体动态超出这一距离范围,则需要更新地理坐标系。
步骤三:获取载体运动第一分段的轨迹参数,如果为空间直线运动类型,在地理坐标系中,根据初速度(ve,0,vn,0,vu,0)、加速度(ae,an,au)的设置,可以确定载体的即时位置、速度。
运动段持续时间内,根据tk-1仿真历元时刻载体的位置(xk-1,yk-1,zk-1)、速度(ve,k-1,vn,k-1,vu,k-1)来计算tk时刻载***置(xk,yk,zk)、速度(ve,k,vn,k,vu,k)的方法如下:
(a)计算时间间隔Δt
Δt=tk-tk-1
(b)计算tk-1时刻到tk时刻的距离变化量(Δse,Δsn,Δsu)
(c)计算tk时刻载体速度(ve,k,vn,k,vu,k)
(d)计算tk时刻载***置(xk,yk,zk)
步骤四:如果载体运动第一分段为空间圆周运动类型,根据运动平面相对地理坐标系方向角γ(取运动平面方向与地理坐标系z轴的夹角,以z轴右侧为正角度)的设置,首先确定地理坐标系中圆周运动的二维平面,该平面的x轴即为地理坐标系x轴(e轴),y轴与x轴垂直。
其次,根据初速度(ve,0,vn,0,vu,0)、圆周半径R和线加速度a的设置,可以确定圆周运动平面内载体的即时位置、速度。运动段持续时间内,计算tk仿真历元时刻圆周运动平面内载***置速度(vx,k,vy,k)的方法如下:
(a)计算时间间隔Δt
Δt=tk-tk-1
(b)计算tk-1时刻到tk时刻的转角变化量Δφ
而
(c)计算tk时刻运动平面内载体方位角φk
φk=φk-1+Δφ
而初始方位
(d)计算tk时刻运动平面内载体速度(vx,k,vy,k)
而|vk|=|vk-1|+a·Δt
而
最后,将运动平面内载体的位置和速度转换到地理坐标系中的相应运动状态参数,转换方法如下:
(b)计算tk时刻地理坐标系载体速度(ve,k,vn,k,vu,k)
(c)计算tk时刻地理坐标系载***置(xk,yk,zk)
步骤五:根据直线或圆周运动类型的用户设置,按照步骤三、步骤四,依次处理后续运动分段。对于运动分段之间的衔接,以上一运动段的末状态为当前运动段的初状态,即上一运动段的载体末位置(xm,ym,zm)、末速度(ve,m,vn,m,vu,m)分别成为当前运动段的初始位置(x0,y0,z0)、初始速度(ve,0,vn,0,vu,0),即
步骤六:在模拟器的每个仿真历元,将地理坐标系中的载体即时位置转换到大地坐标系,用于用户界面实时显示输出;并将载体即时位置、速度转换到地球直角坐标系(又称ECEF坐标系)以进一步支持卫星可见性判断,伪距、多普勒等信号状态参数计算等模拟器功能实现。
其中, a为地球参考椭球半长轴,e为地球参考椭球偏心率。
便于数值计算起见,考虑RN和RM的近似算法,即
而扁率
(b)计算tk时刻ECEF坐标系中相应的载***置坐标(xe,k,ye,k,ze,k)
(c)计算tk时刻ECEF坐标系中相应的载体速度(vx,k,vy,k,vz,k)
实施例
利用本发明提供的提供通用三维载体运动轨迹生成方法,以模拟产生如图3所示的载体运动轨迹为例,并以WGS-84地球坐标系示例,说明具体实现步骤。
步骤一:GNSS模拟器用户设置载体运动轨迹参数,生成载体运动轨迹用户配置文件。
例如,设置载体初始位置:经度116°,纬度39°,高程100m;设置地理坐标系中载体初始速度:东向100m/s,北向0m/s,天向0m/s;设置载体轨迹分段及参数,该轨迹由10个直线运动分段和圆周运动分段组成,以前5个分段为例说明,分别为:
(a)第一分段为直线运动,东向、北向及天向加速度均为0,历时5s;
(b)第二分段为圆周运动,运动平面相对地理坐标系的方向角为0,圆周半径为200m,顺时针运动方向,线加速度为0,历时20s;
(c)第三分段为直线运动,东向、北向及天向加速度均为0,历时5s;
(d)第四分段为圆周运动,运动平面相对地理坐标系的方向角为45°,圆周半径为200m,顺时针运动方向,线加速度为40m/s2,历时5s;
(e)第五分段为直线运动,东向加速度为65m/s2、北向加速度为34m/s2、天向加速度为-40m/s2,历时5s。
步骤二:读取载体运动轨迹用户配置文件,以设置的初始位置为原点建立地理坐标系。该地理坐标系的x轴(又称e轴)沿39°纬线切线指向东,y轴(又称n轴)沿116°经线切线指向北,z轴(又称u轴)与x和y轴所确定的平面垂直指向天顶。
步骤三:依据步骤一中的运动分段轨迹参数设置,实时模拟载体轨迹、计算载体即时速度及位置坐标。
对于第一分段直线运动,载体从前述地理坐标系原点(0,0,0)起始,以100m/s初速度东向匀速直线运动,历时5s(即Δt=5s),则在该运动段的结束时刻,载体的距离变化量、速度和位置坐标分别为:
步骤四:对于第二分段圆周运动,载体从上一运动段的末位置(500,0,0)起始,以上一运动段的末速率为初速率顺时针匀速圆周运动,运动平面维持在前述地理坐标系的水平面内(即γ=0),历时20s(即Δt=20s),则在该运动段的起始时刻,载体的初始方位、初速率、初始位置分别为:
φ0=0.5πrad,
而在该运动段的结束时刻,载体的转角变化量、方位角、运动平面速度和运动平面坐标分别为:
Δφ=-10rad,φm=-8.4292rad,
进一步实现至地理坐标系的转换,则坐标旋转矩阵、地理坐标系载体速度及位置坐标分别为:
步骤五:按照步骤三、步骤四,依次处理后续运动分段。
(a)对于第三分段,载体从上一运动段的末位置(391.1958,-367.8143,0)起始,以(-83.9027,54.4021,0)初速度匀速直线运动,历时5s(即Δt=5s),则在该运动段的结束时刻,载体的距离变化量、速度和位置坐标分别为:
(b)对于第四分段,载体从上一运动段的末位置(-28.3400,-95.8038,0)起始,以上一运动段的末速率为初速率顺时针加速圆周运动,运动平面与前述地理坐标系的水平面呈45度方向角(即γ=45),历时5s(即Δt=5s),则在该运动段的起始时刻,载体的初始方位、初速率、初始位置分别为:
φ0=4.1372rad,v0=100m/s,
而在该运动段的结束时刻,载体的转角变化量、方位角、运动平面速度和运动平面坐标分别为:
Δφ=-5rad,φm=-0.8628rad,
进一步实现至地理坐标系的转换,则坐标旋转矩阵、地理坐标系载体速度和位置坐标分别为:
(c)对于第五分段,载体从上一运动段的末位置(210.5218,-56.5169,56.5169)起始,以(-227.9064,-137.9469,137.9469)初速度加速直线运动,历时5s(即Δt=5s),则在该运动段的结束时刻,载体的距离变化量、速度和位置坐标分别为:
步骤六:在模拟器的每个仿真历元,将地理坐标系中的载体即时位置转换到大地坐标系,用于用户界面实时显示输出。由于以WGS-84地球坐标系示例,WGS-84中地球椭球的主要参数定义如下:半长轴a=6378137m,扁率f=1/298.257223563,则在上述运动分段的结束时刻,用户界面实时显示的大地坐标系载***置坐标与地理坐标系载体速度,见表1;而用于支持模拟器进一步功能实现的ECEF坐标系中相应的载***置坐标及速度,见表2。
表1各运动段结束时刻的大地坐标系载***置与地理坐标系载体速度
表2各运动段结束时刻的ECEF坐标系载***置及速度
Claims (1)
1.一种GNSS模拟器中的通用三维载体运动轨迹生成方法,其特征在于:该方法具体步骤如下:
步骤一:预生成载体运动轨迹用户配置文件,配置文件定义了载体运动分段与轨迹参数设置;
关于载体运动分段的轨迹参数设置,对于直线运动类型,需要设置的轨迹参数包括:直线运动类型标识,运动段持续时间,东向加速度,北向加速度,天向加速度;对于圆周运动类型,需要设置的轨迹参数包括:圆周运动类型标识,运动段持续时间,运动方向即顺时针或逆时针运动,圆周半径,线加速度,圆周运动平面方向相对地理坐标系的方向角;
步骤二:在载体运动轨迹模拟开始时,以设置的初始位置即经度、纬度和高程为原点建立地理坐标系,初始位置(x0,y0,z0)即为(0,0,0);
地理坐标系的x轴即e轴沿当地纬线切线指向东,y轴即n轴沿当地经线切线指向北,z轴即u轴沿当地地理垂线即地球参考椭球表面法线指向天顶;对于动态范围在180km内的载体运动轨迹模拟,这个初始设置的地理坐标系视为三个轴线方向维持不变;所引起的误差忽略不计,则在此距离范围内初始的地理坐标系将被保持;在载体运动轨迹模拟结束之前,如果载体动态超出这一距离范围,则需要更新地理坐标系;
步骤三:获取载体运动第一分段的轨迹参数,如果为空间直线运动类型,在地理坐标系中,根据初速度(ve,0,vn,0,vu,0)、加速度(ae,an,au)的设置,确定载体的即时位置、速度;
运动段持续时间内,根据tk-1仿真历元时刻载体的位置(xk-1,yk-1,zk-1)、速度(ve ,k-1,vn,k-1,vu,k-1)来计算tk时刻载***置(xk,yk,zk)、速度(ve,k,vn,k,vu,k)的方法如下:
(a)计算时间间隔Δt
Δt=tk-tk-1
(b)计算tk-1时刻到tk时刻的距离变化量(Δse,Δsn,Δsu)
(c)计算tk时刻载体速度(ve,k,vn,k,vu,k)
(d)计算tk时刻载***置(xk,yk,zk)
步骤四:如果载体运动第一分段为空间圆周运动类型,根据运动平面相对地理坐标系方向角γ即取运动平面方向与地理坐标系z轴的夹角,以z轴右侧为正角度的设置,首先确定地理坐标系中圆周运动的二维平面,该平面的x轴即为地理坐标系x轴,y轴与x轴垂直;
其次,根据初速度(ve,0,vn,0,vu,0)、圆周半径R和线加速度a的设置,确定圆周运动平面内载体的即时位置、速度;运动段持续时间内,计算tk仿真历元时刻圆周运动平面内载***置速度(vx,k,vy,k)的方法如下:
(a)计算时间间隔Δt
Δt=tk-tk-1
(b)计算tk-1时刻到tk时刻的转角变化量Δφ
而
(c)计算tk时刻运动平面内载体方位角φk
φk=φk-1+Δφ
而初始方位
(d)计算tk时刻运动平面内载体速度(vx,k,vy,k)
而|vk|=|vk-1|+a·Δt
而
其中,λ为载体坐标对应的经度;最后,将运动平面内载体的位置和速度转换到地理坐标系中的相应运动状态参数,转换方法如下:
(b)计算tk时刻地理坐标系载体速度(ve,k,vn,k,vu,k)
(c)计算tk时刻地理坐标系载***置(xk,yk,zk)
步骤五:根据直线或圆周运动类型的用户设置,按照步骤三、步骤四,依次处理后续运动分段;对于运动分段之间的衔接,以上一运动段的末状态为当前运动段的初状态,即上一运动段的载体末位置(xm,ym,zm)、末速度(ve,m,vn,m,vu,m)分别成为当前运动段的初始位置(x0,y0,z0)、初始速度(ve,0,vn,0,vu,0),即
步骤六:在模拟器的每个仿真历元,将地理坐标系中的载体即时位置转换到大地坐标系,用于用户界面实时显示输出;并将载体即时位置、速度转换到地球直角坐标系即ECEF坐标系以进一步支持卫星可见性判断,伪距、多普勒信号状态参数计算模拟器功能实现;
便于数值计算起见,考虑RN和RM的近似算法,即
而扁率
(b)计算tk时刻ECEF坐标系中相应的载***置坐标(xe,k,ye,k,ze,k)
(c)计算tk时刻ECEF坐标系中相应的载体速度(vx,k,vy,k,vz,k)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010296636 CN101975955B (zh) | 2010-09-29 | 2010-09-29 | 一种gnss模拟器中的通用三维载体运动轨迹生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010296636 CN101975955B (zh) | 2010-09-29 | 2010-09-29 | 一种gnss模拟器中的通用三维载体运动轨迹生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101975955A CN101975955A (zh) | 2011-02-16 |
CN101975955B true CN101975955B (zh) | 2013-04-24 |
Family
ID=43575857
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010296636 Active CN101975955B (zh) | 2010-09-29 | 2010-09-29 | 一种gnss模拟器中的通用三维载体运动轨迹生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101975955B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105787003B (zh) * | 2013-08-26 | 2018-11-20 | 上海安吉四维信息技术有限公司 | 实现逼真的三维城市道路地图的方法 |
CN103645933B (zh) * | 2013-12-04 | 2017-02-01 | 北京经纬恒润科技有限公司 | 一种控制方法及装置 |
CN106054149B (zh) * | 2016-07-22 | 2019-04-05 | 中国船舶重工集团公司第七二四研究所 | 一种雷达机动目标三维航迹模拟方法 |
CN108801238B (zh) * | 2017-04-27 | 2021-08-13 | 电子科技大学中山学院 | 一种基于gpx的路径优化方法 |
CN109227543B (zh) * | 2018-10-19 | 2020-07-28 | 华南理工大学 | 一种位置受限柔性关节机器人基于模式的智能控制方法 |
CN109581444B (zh) * | 2018-11-01 | 2023-01-13 | 西北大学 | 一种gps轨迹分段及语义标注方法 |
CN109559375A (zh) * | 2018-11-21 | 2019-04-02 | 北京像素软件科技股份有限公司 | 三维运动轨迹模拟方法及装置、电子设备及存储介质 |
CN109343089B (zh) * | 2018-12-20 | 2020-07-17 | 广州广电计量检测股份有限公司 | 定位设备的性能测试设备、测试方法及其装置 |
-
2010
- 2010-09-29 CN CN 201010296636 patent/CN101975955B/zh active Active
Non-Patent Citations (4)
Title |
---|
姜瑞健等.导航信号模拟器上位机软件设计与实现.《遥测遥控》.2008,第29卷(第6期),53-58. |
导航信号模拟器上位机软件设计与实现;姜瑞健等;《遥测遥控》;20081130;第29卷(第6期);53-58 * |
赵昀等.软件GNSS信号模拟器架构与中频信号生成.《北京航空航天大学学报》.2009,第35卷(第7期),789-793. |
软件GNSS信号模拟器架构与中频信号生成;赵昀等;《北京航空航天大学学报》;20090731;第35卷(第7期);789-793 * |
Also Published As
Publication number | Publication date |
---|---|
CN101975955A (zh) | 2011-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101975955B (zh) | 一种gnss模拟器中的通用三维载体运动轨迹生成方法 | |
CN113311436B (zh) | 一种移动平台上激光测风雷达运动姿态测风订正方法 | |
CN103675844B (zh) | 一种gnss/ins组合导航同步模拟*** | |
CN101858980B (zh) | 一种车载基于gps软件接收机的智能超紧组合导航方法 | |
CN103487820B (zh) | 一种车载捷联/卫星紧组合无缝导航方法 | |
CN105259787A (zh) | 一种组合导航半物理仿真测试同步控制方法 | |
CN102608642A (zh) | 北斗/惯性组合导航*** | |
CN111123307B (zh) | 一种用于检测bds用户设备动态性能的***和方法 | |
CN106773783B (zh) | 一种连续运行参考站***的半物理仿真测试***及方法 | |
CN102736520A (zh) | 一种卫星导航***原理仿真方法和卫星信号模拟器 | |
CN103777218A (zh) | Gnss/ins超紧组合导航***的性能评估***及方法 | |
CN106842256B (zh) | 一种利用gnss单星信号的导航定位方法 | |
CN102508268B (zh) | 一种地基导航信号观测数据生成方法 | |
CN106405592B (zh) | 车载北斗载波相位周跳检测与修复方法及*** | |
Breuer et al. | High precision localisation in customised GNSS receiver for railway applications | |
Branzanti et al. | Variometric approach for real-time GNSS navigation: First demonstration of Kin-VADASE capabilities | |
CN112596077A (zh) | 一种针对终端载体为低轨卫星的卫星导航信号仿真方法 | |
CN105988129A (zh) | 一种基于标量估计算法的ins/gnss组合导航方法 | |
CN102914781A (zh) | 一种格洛纳斯卫星信号星历电文生成方法及装置 | |
CN103033822A (zh) | 移动信息确定装置、方法以及接收机 | |
CN111487660A (zh) | 一种高精度实时微纳卫星集群导航算法 | |
CN108205151B (zh) | 一种低成本gps单天线姿态测量方法 | |
CN103383460B (zh) | 一种高精度高动态上行注入信号生成方法 | |
Farroni et al. | A comparison among different methods to estimate vehicle sideslip angle | |
CN112540389A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |