利用Excel表格解算七参数的方法.docx

介绍一种利用 Excel 表格解算七参数的方法目前计算机的操作系统下都安装有 Office2000---2003，利用 Office 下的Excel 表格和表格中提供的多种矩阵运算函数，可以解算在 GPS RTK 作业中经常用到的由地心坐标系转换到国家坐标系所使用的坐标转换参数或称作七参数利用 EXCEL 表格解算七参数的方法非常方便和实用且便于掌握下面把解算过程中用到的有关矩阵运算的函数列出如下：1. ＝Mmult(数组 1，数组 2) ――――――－两矩阵乘积函数2. ＝Transpose(数组)――――――――――矩阵转置函数3. ＝Minverse(数组)―――――――――― 矩阵求逆函数4. ＝数组 1 ± 数组 2――――――――――矩阵加减函数这里的数组其含义是：例如我们把一个矩阵存放到第 3 行～20 行，第 A列～G 列，那么在进行矩阵运算时数组应输入 A3：G20下面给出解算七参数的数学模型如下：Xi Xi Xi 0 -Z Y Εx/ρ ⊿XYi ＝ Yi ＋ Yi K×10- 6 ＋ Z 0 -X ΕY/ρ ＋ ⊿Y Zi Zi Zi -Y X 0 ΕZ/ρ ⊿ZⅡ Ⅰ Ⅰ ⅠXi XiLi ＝ Y i － Y i （2）Zi ZiⅡ Ⅰ 角标Ⅱ表示地心空间大地直角坐标角标Ⅰ表示参心空间大地直角坐标Li ＝B iY （3）Y ＝[ ⊿X ⊿Y ⊿Z K Εx ΕY Ε Z ] （4）Vi ＝B iY － Li ( i = 1，… t ) t 为公共点个数 （5）V ＝BY － L 所有误差方程 （6）Y ＝( B TB)- 1 BTL （7）上述求解过程认为两种坐标系统的坐标是等权的，或不知道它们的精度信息的情况下所采用的解算公式。

若已知了公共点上两套坐标的精度信息，就可以确定其相应的权矩阵，来确定转换参数若设σ 2xi 对称 σ 2xi 对称ΣⅡi ＝ σ YX σ 2Yi ΣⅠi ＝ σ YX σ 2Yi (8)σ ZX σ ZY σ 2Zi σ ZX σ ZY σ 2Zi ΣLi ＝ΣⅡi ＋ ΣⅠi (9)PLi ＝(ΣLi) - 1 (10)转换参数的最小二乘解为：Y ＝ ( B TPB)- 1 BTPL (11)精度评定：σ 20 =( V TPV) / (3t－7) (12)σ 2Y ＝σ 20 ( BTPB)- 1 (13)下面介绍利用 Excel 表格计算的过程如下：1. 首先确定误差方程的系数矩阵 B，B 矩阵由下式确定：1 0 0 Xi×10-6 0 －Zi/ρ Yi/ρBi ＝ 0 1 0 Yi×10-6 Zi/ρ 0 －Xi/ ρ (14)0 0 1 Zi×10-6 －Yi/ρ Xi/ρ 0i ＝1…t, ρ＝206265组成常数项矩阵 L XⅡi －XⅠiLi ＝ Y Ⅱi －YⅠi i＝1…tZⅡi －ZⅠi如果 t＝3 时，则 B 矩阵为 9 行×7 列，L 为 9 行× 1 列矩阵。

2. 求 B 矩阵的转置矩阵 BT利用=Transpose( 数组)转置矩阵函数进行转置运算，首先在表格中确定将要转置的这个矩阵的行数和列数并把其范围确定下来然后在输入栏中输入转置矩阵函数，数组输入 B 矩阵的所在单元格的行数和列数例如，=TANSPOSE(A3：I20) ，然后按 CTRL+SHIFT+ENTER 键，所得结果就显示在所划定的范围内3. 法方程组成（B TB）或（ BTPB）利用矩阵乘积函数=Mmult（数组 1，数组 2）组成法方程，首先在表格中确定乘积矩阵的行数和列数并把其范围确定下来然后在输入栏中输入矩阵乘积函数，例如： = Mmult（A24：R30，A3 ：G20） ，然后按CTRL+SHIFT+ENTER 键，所得结果就显示在所划定的范围内 4. 法方程常数项矩阵组成（B TL）或（B TPL）操作步骤与步骤 3 相同5. 法方程解算（矩阵求逆） （B TB） -1 或（B TPB） -1利用矩阵求逆函数=Minverse（数组）解算法方程，具体操作同步骤 3，例如： = Minverse（A33：G39） 6. 求转换参数 Y=（B TB） -1（B TL）或（B TPB） -1（B TPL）操作步骤与步骤 3 相同。

精度评定：7. 解算改正数矩阵 V ＝BY － L 利用矩阵加法或减法函数进行计算，例如： =B55：B72－D55 ：D728. 求 V 的转置矩阵 VT计算步骤与步骤 2 相同9. 计算 VTV 或 VTPV操作步骤与步骤 3 相同利用公式（12）计算出单位权中误差 σ 2010. 利用公式（13） ，计算转换参数的方差协方差矩阵，其中主对角线上的元素分别为所解算的七个参数的精度 武汉大学索佳测绘技术服务中心 赵新维※※※ 完 ※※※。