用户到卫星的距离可以先简单理解成通过光速×时间差得到

而用户测得的是包含各种误差影响在内的距离, 称之为伪距。

这是伪距

是用户机钟差,加上用户三维坐标,共四个未知量,所以需要至少同时看到四颗卫星才能实现定位。

设用户坐标是(X,Y,Z),三颗卫星位置根据星历可得是已知量,设为(Xi,Yi,Zi)

则用户对第i个卫星的伪距为:

后面的一串是各种误差,比如电离层影响等等,这些都根据模型可提前求得,算作已知量。

写完整:

提前写明,最终求解这个四元方程的方法是  迭代法

因此对于(X,Y,Z)会需要一个初值,可以理解为近似值,但实际你随便取无影响,只会对迭代系数有一点点的影响。

我们把这个初值定为(X0,Y0,Z0),令

然后将一开始伪距的表达式泰勒一阶展开,这里是三元的泰勒展开,anyway只有套用三元的一阶展开公式就行

三元函数一阶泰勒展开公式

得到(5)(左图ΔX=X-X0带入了)

引入这个ρ01,类似到ρ04,带入(5)并求偏导

可得    ,再带入(5),(5)就可用偏导表示 

完整版

,记

(6)写成矩阵表示:,因此

最终可解得ΔX,ΔY,ΔZ,CΔtu。通过更新用户坐标

如果超过四颗卫星的信息,前面都一样,相当于方程数多于未知量个数,此时只需要在最后加上最小二乘法的处理就可以得到最优解了。

考虑到涉及到的矩阵运算比较多,所以选用MATLAB来计算

北斗导航系统伪码定位原理,MATLAB迭代法求解相关推荐

  1. MATLAB秦九韶多项式求值算法的原理和迭代法求解的近似数值方法。

    1..熟悉常用的Matlab操作: 2.了解秦九韶多项式求值算法的原理和迭代法求解的近似数值方法. 秦九韶多项式求值算法: 迭代法求解的近似数值: x=2; for k=1:10x=(x+2/x)/2 ...

  2. 北斗信号服务器解算,北斗导航系统接收机定位解算设计与实现

    摘要: 随着北斗导航系统的建设不断推进,其应用范围越来越广,因此北斗接收机需求也越来越大.不同的应用场景的接收机结构和侧重点有所不同,但是其中的定位解算模块都是其关键部分.本文主要对北斗接收机的整体结 ...

  3. 北斗导航 | 卫星导航定位原理与应用——导航电文内容

    ===================================================== github:https://github.com/MichaelBeechan CSDN: ...

  4. 北斗系统基础知识1(北斗系统定位原理说明“图文详述”)

    需求说明:知识储备 内容       :北斗系统起源及工作原理 来自       :时间的诗 原文:http://www.bingdun.com/news/bingdun/3252.htm 从来没有那 ...

  5. 单片机+北斗模块实现定位

    本文原创,转载请注明出处.小编可能以后不会在简书简书地址上更新了,因为这里的markdown 编辑器太不好用了,没有CSDN 上面的好用,也希望粉丝们多多关注本人CSDN(一样的昵称呦). 导航是继移 ...

  6. 北斗导航系统部署完毕,为何没推出类似于谷歌地图高精度地图?

    许多人都认为地图是由北斗.GPS等卫星下发给我们的,但其实不是,定位数据只是电子地图中的很小的一部分,但它绝对是相当重要的部分,因为定位后,地图就可以实现导航的功能,所有的矢量数据.底图图层就活了起来 ...

  7. 非线性方程组牛顿迭代法matlab,matlab实现牛顿迭代法求解非线性方程组

    <matlab实现牛顿迭代法求解非线性方程组>由会员分享,可在线阅读,更多相关<matlab实现牛顿迭代法求解非线性方程组(5页珍藏版)>请在人人文库网上搜索. 1.matla ...

  8. 北斗系统基础知识2(北斗一代定位原理详述)

    需求说明:知识储备 内容       :北斗系统起源及工作原理 来自       :时间的诗 原文:http://bbs.tianya.cn/post-funinfo-3117078-1.shtml ...

  9. 用matlab求解jacobi,用jacobi迭代法求解线性方程,求助matlab大师

    问题描述: 用jacobi迭代法求解线性方程,求助matlab大师 用jacobi迭代法求解线性方程 9x1-x2+x3=10 -x1+10x2-2x3=7 -2x1+x2+10x3=6,设迭代初值为 ...

最新文章

  1. 内存泄漏的定位与排查:Heap Profiling 原理解析
  2. opencv图像处理常用完整示例代码总结
  3. Yii的Querybuild ActiveRecord
  4. python把汉字转换为二进制数_在Python中,如何将8位二进制数转换为ASCII字符?
  5. Java Persistence API:快速入门
  6. java spring mvc json ajax 优势_SpringMVC后台json数据前台ajax获取不到!!!急求解答!!!...
  7. dbscan聚类算法_一种视频人群流的轨迹聚类方法
  8. FairMOT: On the Fairness of Detection and Re-Identification in Multiple Object Tracking 效果展示
  9. 吴恩达机器学习笔记-非监督学习
  10. shell判断字符串为空
  11. css中标准盒模型和怪异盒模型的区别,如何将标准盒模型转换为怪异盒模型
  12. pureftpd 配置 mysql_Pure-ftp配置文件详解
  13. 大疆创新2019校招
  14. 使用EditPlus替换文件中的某一字符为换行符
  15. MPI并行编程: 矩阵乘法,积分计算,PSPR排序
  16. 直连附加存储服务器,直接连接存储
  17. 怎样设置二维码的尺寸
  18. C# CAD二次开发之基本图形
  19. 我有一个梦,袁隆平!
  20. sim900芯片—GPRS模块ppp拨号上网

热门文章

  1. 明白什么是“五险一金”
  2. 在路径为……的上下文中,servlet[....]的Servlet.service()引发异常
  3. 儿童背包亚马逊美国站CPC认证ASTMF963
  4. 在 D 天内送达包裹的能力
  5. Multisim_Analyses and Simulation
  6. 几何观止(经典几何对象)
  7. Python实现钱龙长线(LON)指标
  8. MPU6050的初始化
  9. Vue-教你实现一个订单列表倒数30分钟的功能
  10. 十一乐一乐 - 让程序员崩溃的瞬间(1)