GPS定位问题数学建模讲课讲稿
- 格式:doc
- 大小:575.50 KB
- 文档页数:12
数学建模GPS 定位问题摘要本次建模中要解决根据GPS 卫星位置来确定GPS 信号接收机位置的问题,在本次建立的模型中主要用到的是点定位的数学模型,用码伪距进行点定位。
再用Matlab 编程解得地点位置,最后转换成其经度和纬度。
对于问题一,我们采用GPS 定位中单点定位的方法(单点定位利用一点采集的观测数据和广播星历确定点的坐标)。
题目中假定了卫星所在的空间位置是准确值因此不考虑广播星历。
往往伪距方程解算的基本思路是将非线性观测方程进行Taylor 级数展开至一阶,忽略二阶及以上的高阶项,得到线性观测方程。
我们将上面的每两个非线性观测方程相减消去二阶及以上的高阶项可得到42C 个四元一次方程。
在此基础上派生出64C 个线性方程组并用2222R iz i y i x =++ 进行验证选择最符合的坐标,得到四个地点在地心空间直角坐标系的坐标是(-2179,4373,4081) ; (-2174,3,4381,4090);(-2169,4410.1,4123);(-2159,4382.4,4142.3);再转换成经度和纬度就是(40:08:38.58167N ,116:10:14.01669E);(40:05:39.12131N ,116:23:48.72859E);(40:10:46.58408N ,116:11:20.90291E); (40:29:04.29791N ,116:13:23.03773E)然后再在地图上标出各个点的位置对于问题二,由于添加了一个点,多出了一个数据,可以同样的继续采用上述方法,只是每两个非线性观测方程相减消去二阶及以上的高阶项可得到52C 个四元一次方程。
在此基础上派生出104C 个线性方程组并用2222R iz i y i x =++进行验证选择最符合的坐标(2129,4361,4125)转换成经纬度(40°33′05.71354″N ,116°01′10.64958″E)关键词:点定位 码伪距 钟差 单点定位 MATALAB 编程一、问题重述全球定位系统(GPS )是美国国防部研制的导航定位授时系统,由24颗等间隔分布在6个轨道面上20200公里高度的卫星组成。
Geometrical tolerancing - Tolerances of form ,orientation, location and run-out几何公差 形状、方向、位置和跳动公差标注Geometrical Product Specification (GPS )(IDT ISO 1101Geometrical tolerancing - Tolerances of form ,orientation, location and run-out GB/T 1182-2008:2004)王喜力中国航空综合技术研究所2010.5G P S 标准的发展和体系介绍 Geometrical Product Specification (GPS )Geometrical tolerancing - Tolerances of form ,orientation, location and run-outGeometrical tolerancing - Tolerances of form ,orientation, location and run-out2.2.3 GPS通用标准矩阵(General GPS standards metrix)标准链:链环1:产品图样表示;链环2:公差定义;链环3:实际要素特征或参数定义;链环4:工件误差评定(与公差极限比较);链环5:测量器具;链环6:测量器具标定和校准;确定规范之间关系Geometrical tolerancing - Tolerances of form ,orientation, location and run-out2.3.1 产品(工件)的表面模型产品(工件)的表面模型是产品功能描述、几何设计、加工制造及检验/验证过程中建立的几何表达模型,是规范产品几何定义和实际工件检验/验证的基础。
公称表面模型规范表面模型验证表面模型Geometrical tolerancing - Tolerances of form ,orientation, location and run-outGeometrical tolerancing - Tolerances of form ,orientation, location and run-out•分离(Partition )提取(Extraction )•滤波(Filtration )拟合(Association )•组合(Collection )构建(Construction ))Geometrical tolerancing - Tolerances of form ,orientation, location and run-outGeometrical tolerancing - Tolerances of form ,orientation, location and run-out 2009-03-1239组成要素(表面、轮廓)导出导出要素(中心点、中心线、中心面)几何要素定义之间相互关系的结构框图。
G P S定位问题数学建模数学建模GPS 定位问题摘要本次建模中要解决根据GPS 卫星位置来确定GPS 信号接收机位置的问题,在本次建立的模型中主要用到的是点定位的数学模型,用码伪距进行点定位。
再用Matlab 编程解得地点位置,最后转换成其经度和纬度。
对于问题一,我们采用GPS 定位中单点定位的方法(单点定位利用一点采集的观测数据和广播星历确定点的坐标)。
题目中假定了卫星所在的空间位置是准确值因此不考虑广播星历。
往往伪距方程解算的基本思路是将非线性观测方程进行Taylor 级数展开至一阶,忽略二阶及以上的高阶项,得到线性观测方程。
我们将上面的每两个非线性观测方程相减消去二阶及以上的高阶项可得到42C 个四元一次方程。
在此基础上派生出64C 个线性方程组并用2222R iz i y i x =++ 进行验证选择最符合的坐标,得到四个地点在地心空间直角坐标系的坐标是(-2179,4373,4081) ; (-2174,3,4381,4090);(-2169,4410.1,4123);(-2159,4382.4,4142.3);再转换成经度和纬度就是(40:08:38.58167N ,116:10:14.01669E); (40:05:39.12131N ,116:23:48.72859E); (40:10:46.58408N ,116:11:20.90291E); (40:29:04.29791N ,116:13:23.03773E)然后再在地图上标出各个点的位置对于问题二,由于添加了一个点,多出了一个数据,可以同样的继续采用上述方法,只是每两个非线性观测方程相减消去二阶及以上的高阶项可得到52C 个四元一次方程。
在此基础上派生出104C 个线性方程组并用2222R i z i y i x =++进行验证选择最符合的坐标(2129,4361,4125)转换成经纬度(40°33′05.71354″N ,116°01′10.64958″E)关键词:点定位 码伪距 钟差 单点定位 MATALAB 编程一、问题重述全球定位系统(GPS )是美国国防部研制的导航定位授时系统,由24颗等间隔分布在6个轨道面上20200公里高度的卫星组成。
GPS 用户从接收的GPS 信号可以得到足够的信息进行精密定位和定时。
卫星所在的空间位置由卫星的轨道参数确定,为简化问题,本题题目里假定它是准确值。
题目中为了简化问题,假定卫星所在的空间位置是准确值。
GPS 信号到达接收机的时间是由卫星上的时钟(铯原子钟)和地面接收机上的时钟(低成本钟)决定,钟差是未知的。
今给出了4颗卫星在地心空间直角坐标系上的坐标(地心空间直角坐标系就是将坐标系的原点O与地球质心重合,Z 轴指向地球北极,X 轴指向经度原点E ,Y 轴垂直于XOZ 平面构成右手坐标系),地球的半径,光速,以及4颗卫星的GPS 信号到达四个GPS 接收机地点处的时间。
根据4颗卫星在地心空间直角坐标系上的坐标见表A.1,以及4颗卫星的GPS 信号到达四个GPS 接收机地点处的时间见表A.2,求得四个GPS 接收机地点处在地心空间直角坐标系上的坐标。
然后将坐标转换成经度与纬度。
对于多于四颗卫星的问题,怎么建立一个更好的模型,才能精确的确定某个地点的位置,并将其转化为经纬度在图中标出。
对于多点定位问题,应该考虑周全,误差值,偏差值,并使得每颗卫星都能准确的将其地位。
二、问题分析在问题一的求解上,已知四颗卫星的地心直角坐标系位置,并且知道每颗卫星GPS 信号到达每个地点的时间,由于卫星在发送时有时延,并且GPS 接收仪在接收的过程中也会有时延,但这题题目中假定了卫星所在的空间位置是准确值,那我们就直接通过码伪距进行点定位,用派生出的线性方程组求解得到4个GPS 接收仪在地心空间直角坐标系上的坐标,并通过google 地球在地图上标明位置所在;对于第二问,在通常的情况下,地面的GPS 接收机能收到5—8颗卫星的信号,对于多于4颗卫星的情况,应该周全考虑派生出的104C 个线性方程组,因此在模型一的基础上,在码伪距的测量上对其时间误差进行考虑,从而得到)(*δ-=ij T c i k R ,并通过卫星坐标列出四个方程组成方程组,求得新加点的坐标,并在google 地球上表明。
特别要注意的是,在地球上,每个点的空间直角坐标x ,y ,z 都必须满足2222R i z i y i x =++ (R 为地球的半径=6371公里),我们将以这个式子来检验方程组解得的数据三、模型假设1假设每一颗卫星的发送时延都是一样的2 假设每一个GPS接收仪的接收时延都是一样的3 假设每颗卫星到达每个地点的时间值天气状况一样,所得的结果都是准确的4 忽略GPS信号在传播过程中所收到的干扰5 忽略大气层和电离层的残差对水平位置定位误差6 假定单点定位的精度不受广播星历误差和钟信息(包括选择可用性误差)的限制。
四、符号说明i D 接收机和卫星之间的实际距离R 地球半径五、模型分析、建立及求解 5.1用码伪距进行点定位:码伪距的观测方程可表示为:)(*δ-=ij T c i k R i R 为卫星i 与测站k 的码伪距观测值,c 为光速, δ为接收机钟差和卫星钟差之差,几何距离i D =2)(2)(2)(iz Zi i y Yi i x Xi -+-+- 并且满足:2222R i z i y i x =++根据以上可以列出地点位置关于卫星位置的关系等式:c iz i Z i y i Y i x i X =-+-+-2)(2)(2)(*(i j T -δ)可列出第一个地点的458.2997922)110100(2)115150(2)18747(=-+-+-z y x *(0.054354-δ) 458.2997922)15228(2)116898(2)19756(=-+-+--z y x *(0.0489226-δ) 458.2997922)117494(2)110100(2)1(=-+-+-z y x *(0.0491307-δ) 458.2997922)114284(2)17142(2)112370(=-+-+--z y x *(0.0489224-δ)并用2212121R z y x =++检验最后的结果。
往往伪距方程解算的基本思路是将非线性观测方程进行Taylor 级数展开至一阶,忽略二阶及以上的高阶项,得到线性观测方程。
我们将上面的每两个相减消去二阶及以上的高阶项在此基础上派生出来的方法有解线性方程组: a=[37006 -3496 9744 1.0808e+00917494 10100 -14788 1.0434e+00942234 16016 -8368 1.0805e+0095228 19512 -18112 -3.2355e+005];b=[-5.6125e+007-11415-5.4291e+007-7527-5.6109e+007-152111.5829e+004-3796];用MATLAB 求解可以得到第一个地点在地心空间直角坐标系的坐标(1,1,1z y x )大致为(-2179,4373,4081) ;同样的方法可以得到剩下的三个地点的坐标分别是(-2174.3 , 4381 , 4090);(-2169 , 4410.1 , 4123);(-2159 , 4382.4 , 4142.3);转换成经度和纬度,这四个地点就分别为40°08′16.90161″N 116°10′14.30207″E2. 40°05′39.1213″N , 116°23′48.72859″E3. 40°10′46.58408″N, 116°11′20.90291″E4. 40°29′04.29791″N, 116°13′23.03773″E5.2第二问对于多于4颗卫星的情况458.2997922)510100(2)515150(2)58747(=-+-+-z y x *(0.0547118-δ) 458.2997922)55228(2)516898(2)59756(=-+-+--z y x *(0.0489472-δ) 458.2997922)517494(2)510100(2)5(=-+-+-z y x *(0.0489068-δ) 458.2997922)514284(2)57142(2)512370(=-+-+--z y x *(0.0488635-δ) 458.2997922)510100(2)515723(2)57669(=--+-+--z y x *(0.0633407-δ) 用MATLAB 求解可以得到第五个地点在地心空间直角坐标系大致为 (-2129,4361,4125)转换成经度和纬度,这个点就为40°33′05.71354″N ,116°01′10.64958″E六、模型评价6.1缺点1.伪距定位法是利用全球卫星定位系统进行导航定位的最基本的方法,其基本原理是:在某一瞬间利用GPS接收机同时测定至少四颗卫星的伪距,根据已知的卫星位置和伪距观测值,采用距离交会法求出接收机的三维坐标和时钟改正数。
伪距定位法定一次位的精度并不高,存在着误差2.对于伪距的测量与计算中,某些时候忽略了时间的误差,造成数据的偏差3.对于点定位问题,往往考虑的更多的是点对点定位,却忽略了相互点位6.2优点对于点定位模型,定位速度快,很容易通过卫星的位置以及信号接收时间求得定位点的位置七、模型的改进与推广1.,在接收机多余于四个的求解过程中,如果考虑用到由最小二乘原理求未知数的未知量则最好不过;2. 由于轨道误差和电离层效应,基准站接收机直接测量的伪距不同于精确距离,两者之间的差异应该伪距该正数;3.对于此模型,可以推广到对于四点定位的求解上,在四点定位中能快速确定地点坐标。