测量数据处理与程序设计.ppt

测量数据处理与程序设计 束蝉方 shuchanfang@ 南京工业大学测绘学院 课程概况 课程主要内容和学时安排 第一章 matlab语言介绍 一、MATLAB概述 MATLAB － MATrix LABoratory 20世纪七十年代后期：时任美国新墨西哥大学计算机科学系主任的Cleve Moler教授用FORTRAN设计了一组程序接口，此即萌芽状态的MATLAB。 1984年成立了MathWorks公司，并把MATLAB正式推向市场。从这时起，MATLAB的内核采用C语言编写，而且除原有的数值计算能力外，还新增了数据图视功能。 进入20世纪九十年代后，MATLAB已经成为国际上最常用的科学计算软件之一。 3M：Matlab、Mathematica、Maple 二、MATLAB语言特点 二、MATLAB语言特点 二、MATLAB语言特点 二、MATLAB语言特点 二、MATLAB语言特点 三、MATLAB数值运算 1.1 变量的命名 在MATLAB中，变量名是以字母开头，后接字母、数字或下划线的字符序列，最多19个字符。 在MATLAB中，变量名区分字母的大小写。MATLAB提供的标准函数名以及命令名必须用小写字母。 1.2 赋值语句 赋值就是把数赋予代表常量或变量的标识符。 MATLAB赋值语句有两种格式： (1) 变量=表达式（或数） (2) 表达式 矩阵的值放在“[ ]” 中，同一行元素之间以“,” 或空格分开，行与行之间用“;”号分开； 句尾用“, ”或回车，此时显示结果；若不希望显示结果以“;”号结尾再回车。 变量的元素用圆括号“( )”中的数字（下标）表示。若要显示整行或整列，则可以用(:)冒号来表示。冒号(:)代表矩阵中行(ROWS)或列(COLUMNS)的全部。 可把matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。 1、save —— 将工作空间中的变量存到*.mat文件中。 save data.mat a b； 如果赋值元素的下标超出原来矩阵的大小，矩阵的行列会自动扩展 把a的第2行和第4行及第1列和第3列交点上的元素提出，构成一个新矩阵b 要抽去a中的第2行，第4行，第5行，可用空矩阵[ ] 的概念。 1、矩阵加、减（＋,－）运算 规则： ? 相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。 ? 允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。 2、矩阵乘（?）运算 规则： A矩阵的列数必须等于B矩阵的行数 标量可与任何矩阵相乘。 a=[1 2 3;4 5 6;7 8 0]; b=[1;2;3]; c=a*b c =14 32 23 3、矩阵除法 ? 恰定方程组的解 方程ax=b(a为非奇异) x=a-1 b 矩阵逆 两种解: x=inv(a)?b — 采用求逆运算解方程 x=a\b — 采用左除运算解方程 ?超定方程组的解 方程 ax=b ,mn时此时不存在唯一解。 方程解 (a a)x=a b x=(a a)-1 a b —— 求逆法 x=a\b —— matlab用最小二乘法找一 个准确地基本解。 ?欠定方程组的解 当方程数少于未知量个数时,即不定情况,有无穷多个解存在。 matlab可求出两个解： 用除法求的解x是具有最多零元素的解 是具有最小长度或范数的解，这个解是基于伪逆pinv求得的。 4、矩阵的乘方 4、矩阵的乘方 M文件的语法类似于c语言，但又有其自身特点，它只是一个简单的ASCII码文本文件，执行程序时逐行解释运行程序，Matlab是解释性的编程语言。 M文件有两类： ?独立的m文件 — 称脚本文件(Script File) ?可调用m文件 — 称函数文件(Function File ) 函数m文件第一行必须以单词function作为引导词，必须遵循如下形式： function [因变量