设为首页收藏本站

路桥隧 - 免费·分享·共赢

 找回密码
 注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

开启左侧

[路桥] 介绍一种利用Excel表格解算七参数的方法

    [复制链接]
发表于 2010-1-14 10:51:36 | 显示全部楼层 |阅读模式
本帖最后由 路桥吾爱资料员 于 2010-1-14 10:55 编辑



    目前计算机的操作系统下都安装有Office2000---2003,利用Office下的Excel表格和表格中提供的多种矩阵运算函数,可以解算在GPS RTK作业中经常用到的由地心坐标系转换到国家坐标系所使用的坐标转换参数或称作七参数。利用EXCEL表格解算七参数的方法非常方便和实用且便于掌握。下面把解算过程中用到的有关矩阵运算的函数列出如下:

1.
=Mmult(数组1,数组2) ――――――-两矩阵乘积函数


2.
=Transpose(数组)――――――――――矩阵转置函数

3.
=Minverse(数组)―――――――――― 矩阵求逆函数

4.
=数组1 ± 数组2――――――――――矩阵加减函数

    这里的数组其含义是:例如我们把一个矩阵存放到第3行~20行,第A列~G列,那么在进行矩阵运算时数组应输入A3:G20。

   下面给出解算七参数的数学模型如下:

详情请下载附件查看
如果t=3时,则B矩阵为9行×7列,L为9行×1列矩阵。

2. 求B矩阵的转置矩阵BT
  利用=Transpose(数组)转置矩阵函数进行转置运算,首先在表格中确定将要转置的这个矩阵的行数和列数并把其范围确定下来。然后在输入栏中输入转置矩阵函数,数组输入B矩阵的所在单元格的行数和列数例如,=TANSPOSE(A3:I20),
然后按CTRL+SHIFT+ENTER键,所得结果就显示在所划定的范围内。

3. 法方程组成(BTB)或(BTPB)
  利用矩阵乘积函数=Mmult(数组1,数组2)组成法方程,首先在表格中确定乘积矩阵的行数和列数并把其范围确定下来。然后在输入栏中输入矩阵乘积函数,例如: = Mmult(A24:R30,A3:G20),然后按CTRL+SHIFT+ENTER键,所得结果就显示在所划定的范围内。
4. 法方程常数项矩阵组成(BTL)或(BTPL)
  操作步骤与步骤3相同。
5. 法方程解算(矩阵求逆) (BTB)-1或(BTPB)-1
  利用矩阵求逆函数=Minverse(数组)解算法方程,具体操作同步骤3,
例如: = Minverse(A33:G39)。

6. 求转换参数Y=(BTB)-1(BTL)或(BTPB)-1(BTPL)
  操作步骤与步骤3相同。

精度评定:
7. 解算改正数矩阵V
=BY
- L
  利用矩阵加法或减法函数进行计算,例如:
=B55:B72-D55:D72

8. 求V的转置矩阵VT
  计算步骤与步骤2相同

9. 计算VTV或VTPV
  操作步骤与步骤3相同。利用公式(12)计算出单位权中误差σ20

10. 利用公式(13),计算转换参数的方差协方差矩阵,其中主对角线上的元素
   分别为所解算的七个参数的精度。

                             武汉大学索佳测绘技术服务中心  赵新维

                  
※※※ 完 ※※※

详情请下载附件查看

Excel表格解算七参数的方法.rar

14.71 KB, 下载次数: 1863

发表于 2010-3-15 20:37:09 | 显示全部楼层
呵呵好东西啊
发表于 2010-5-21 18:16:08 | 显示全部楼层
下下来研究下
发表于 2010-5-22 10:50:21 | 显示全部楼层
好东西 必须顶
发表于 2010-5-23 11:04:43 | 显示全部楼层
好程序 下来研究研究
发表于 2010-5-28 13:07:32 | 显示全部楼层
好东西,支持。。。。。
发表于 2012-5-4 15:36:14 | 显示全部楼层
下下来研究下
发表于 2013-5-12 07:42:07 | 显示全部楼层
呵呵,有点复杂啊啊,看不懂啊
发表于 2014-9-4 09:25:27 | 显示全部楼层
好贴,值得收藏,大家都顶












12V7A电源适配器
发表于 2014-9-20 00:16:48 | 显示全部楼层
谢谢楼主,真是太好了












19V4A电源适配器
您需要登录后才可以回帖 登录 | 注册

本版积分规则

QQ|【路桥隧】 ( 湘ICP备16018960号-1 )  QQ群:【路桥吾爱-lq52.com】

GMT+8, 2024-3-29 01:36

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表