fx-5800P的矩阵计算
fx-5800P内值Mat A、Mat B、Mat C、Mat D、Mat E、Mat F六个矩阵及答案矩阵Mat Ans,矩阵的行列数可以由用户定义,大不应超过10行×10列。允许按矩阵运算规则进行矩阵的加、减、乘、行列式(det)、转置(Trn)、求逆(-1)、平方(2),求元素绝对值(Abs)、乘或除以标量等运算。
[例5-1] 已知某水准网的观测值为LT=(1.359 2.009 0.363 1.012 0. 567 0.238 -0.595)(m)
条件方程AV-W=0的系数矩阵A闭合差矩阵W分别为
(mm)
高差观测值的协因数阵为Q=diang(1.1 1.7 2.3 2.7 2.4 1.4 2.6)
试求高差观测改正数的小二乘解、观测值的平差值及单位权中误差。
[解] ① 定义与输入已知矩阵
按 键定义矩阵A为4×7阶,输入条件方程系数矩阵A。
按 键定义矩阵B为4×1阶,输入条件方程闭合差列矩阵W。
按 键定义矩阵C为7×7阶,输入协因数矩阵Q。
按 键定义矩阵F为7×1阶,输入观测值列矩阵L。
② 计算法方程联系数的解K=(AQAT)-1W。
执行矩阵表达式(Mat AMat CTrn(Mat A))-1Mat B→Mat E计算联系数向量K,结果见图5-1左上1-2图所示。
图5-1 在COMP模式下求解[例5-1]的操作过程
③ 计算观测值的改正数V=QATK。
执行矩阵表达式Mat CTrn(Mat A)Mat E→Mat B计算观测值改正数向量V,结果见图5-1右上1-2图所示。
④ 计算验后单位权中误差 ,式中r=4为多余观测数。
执行矩阵表达式Trn(Mat B)Mat C-1Mat B计算VTPV值,见图5-1左下图1-2所示。执行表达式√(Mat Ans[1,1]÷4)→S计算单位权中误差 的值,结果见图5-1下右2图所示。
⑤ 计算高差观测值的平差值
执行矩阵表达式Mat F+Mat B÷1000计算出 的值,结果见图5-1右下图所示。
|