电话:028-83223983 65819088
                87711546
手机:13618037680
传真:028-83223983
Email:cdxin.rui@163.com
   
[技术]CASIO新型编程计算器fx-5800P系列讲座(二)
作者:成都新瑞测绘  来源:覃辉   浏览:1573次  更新:2009/7/11 14:04:21

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右下图所示。

 
返回列表 | 打印本页 | 返回顶部 
地址:成都市营门口路18号附44号 电话:028-83223983 87711546 传真:028-83223983 Email:cdxin.rui@163.com
版权所有 © 成都新瑞测绘仪器有限公司 技术支持:仕航软件 备案号:蜀ICP备09012819号-1