测量数据处理 MATLAB等工具软件及其应用 Matlab工具软件及其应用.docxVIP

PAGE PAGE 1 Matlab工具软件及其应用 案例导入： 如图2-6所示的水准网中，和是已知高程的水准点，并设这些点的已知高程值无误差，图中、、是待定点。已知高程、观测高差及相应水准路线的长度见表2-1，试用MATLAB按间接平差法求： （1）各待定点的平差高程； （2）至点间高差平差值的中误差。 图 图2-6 知识准备： 一、基本知识介绍 1．启动7.1程序后的界面如图2-7所示。 图2-7 2．命令窗口（Command Window） 命令窗口是的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。命令窗口不仅可以内嵌在的工作界面，如图2-7所示，还可以以独立窗口的形式浮动在界面如图2-8所示。 图2-8 命令窗口中的“”为命令提示符，表示正在处于准备状态。在提示符后输入命令并按下回车键后，会解释执行所输入的命令，并在后面给出计算结果。如在命令提示符下输入“version”，则会显示当前使用的版本号。 3． 命令介绍 命令是查询函数语法的最基本方法，查询信息直接显示在命令窗口。在命令窗口中输入加函数名，用来显示该函数的帮助说明。例如，为了显示求矩阵的逆阵的函数“inv”的使用方法与说明，可在后输入“inv”，屏幕将显示帮助信息，如图2-9所示。 二、数据及其运算 矩阵是最基本、最重要的数据对象，的大部分运算或命令都是在矩阵运算的意义下执行的。 向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个元素都可以作为矩阵的特例来处理。 1．中的变量与赋值 在7.1中，变量名是以字母开头，后接字母、数字或下划线的字符序列，最多63个字符；变量名区分字母的大小写，即“myexample”、“Myexample”表示不同的变量。此外，提供的标准函数名以及命令名必须用小写字母，如求矩阵A的逆用函数“inv（A）”，不能写成“Inv（A）”或“INV（A）”，否则会出错。 的赋值语句有两种格式： （1）变量名=表达式 （2）表达式 其中表达式是用运算符将有关运算量连接起来的式子，其结果是一个矩阵。在第一种语句形式下，将右边表达式的值赋给左边的变量，而在第二种语句形式下，将表达式的值赋给的预定义变量“ans”。 一般地，运算结果在命令窗口中显示出来。如果在语句的最后加分号“；”，那么仅仅执行赋值操作，不显示运算的结果。所以，当不需要显示运算结果时，则可以在赋值语句的后面加上“；”。 在语句后面或前面可以加上注释，用于解释或说明语句的含义，对语句处理结果不产生任何影响。注释以“%”开头，后面是注释的内容。 2．中矩阵的建立与运算 （1）直接输入法建立矩阵 将矩阵的元素用“[ ]”括起来，按行的顺序输入各元素，同一行的各元素之间用空格或用“，”分隔，不同行的元素之间用“；”分隔。 如输入A矩阵中各元素后按回车键，则显示A矩阵如图2-10所示。 图2-10 （2）特殊矩阵的建立 还可以用通过特殊矩阵函数产生特殊矩阵，如0矩阵，单位矩阵等。以产生阶单位矩阵为例，用函数“eye（4）”实现，如图2-11所示。 图2-11 （3）构造对角矩阵 在平差过程中，当观测值之间不相关时，观测值的权阵是对角矩阵。 设为具有个元素的向量，将产生一个对角矩阵，其主对角元素即为向量V的元素。 设，构造对角矩阵如图2-12所示。 图2-12 （4）矩阵的转置 矩阵的转置是单撇号“’”。 例：设，则，在中的执行情况如图2-13所示。 图2-13 （5）方阵的逆阵 直接调用函数inv(A)，即可。 例：设，求A的逆阵B，如图2-14所示。 图2-14 可以验证，，如图2-15所示。 图2-15 案例解答： 运用MATLAB前，先根据平差问题，计算多余观测数；根据参数的近似值，列出误差方程；平差值函数的表达式： 表2-1 路线编号 观测高差/m 路线长度/km 已知高程/m 1 1.359 1.1 2 2.009 1.7 3 0.363 2.3 4 1.012 2.7 5 0.657 2.4 6 0.238 1.4 7 -0.595 2.6 1）计算多余观测数，选择参数并计算其近似值 由图易知，本例中多余观测数，设、、点高程平差值为、、，相应的近似值取为： 2）列误差方程 将观测数据等代入观测方程，得误差方程的矩阵形式： 则可以得到误差方程的系数矩阵与常数项矩阵： ， 3）列至点间高差平差值的权函数关系式： 式中，为与求精度无关的常数，故省略具体值的计算。 即权函数式的系数： 第一种方法：在命令窗口中分步骤使用MATLAB 为应用方便起见，将单位权中误差用字母“”表示，将所求函数的中误差用“”表示，待定点高程平差值用“”表示，待定点高程近似值用 “”表示，及其改正数平差值用“”表示。则