1、第3章 MATLAB 程序设计(1),主要内容：,关系运算和逻辑运算规则； 常用的逻辑函数； 函数文件的定义和调用； 数据文件的读写。,3.1关系运算,3.1.1 概述关系运算表达式的计算结果是一个由0和1组成的逻辑数组（logical array），在数组中用1表示真，0表示假。,3.1.2 关系运算符,6种关系运算符(大于)= (大于或等于)= (等于)= (不等于),3.1.3 关系运算规则,关系运算符的运算法则为： 1 当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0。 2 当参与比较的量是两个同型的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规
2、则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个与原矩阵同型的矩阵，它的元素由0或1组成。,3 当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个与原矩阵同型的矩阵，它的元素由0或1组成。 4 对于复数，、=，仅对参加比较变量的实部进行比较，= = 和 = ，则同时对实部和虚部进行比较。,例3-1 分析语句a=2+2=4的执行结果,例3-2 分析语句 a=fate;b=cake;result=a=b的执行结果。,例3-3 已知a=1 2 3 4;b1=3，b2=0.1 2.1 2.9 4.1，
3、分析a与b1 ，b2的六种关系运算的结果。,3.1.4 关系运算的应用,1 find函数 函数find函数的格式为：,indices = find(X)功能：返回数组X的非零元素的序号，如无非零元素，返回空数组。,i,j = find(X)功能：返回数组X的非零元素的行号和列号，此格式很适合稀疏数组。,i,j,v = find(X)功能：返回数组X的非零元素的行号、列号以及列向量。,例3-4 已知A=4,-65,-54,0,6;56,0,67, - 45,0，分析下列语句的功能。,find(A4) %找出大于4的元素的序号 i j=find(A,1) %找出一个不为零的元素的行号、列号 i j
4、=find(A,1,first) %找出第一个不为零的元素的行号、列号 i j=find(A,1,last) %找出最后一个不为零的元素的行号、列号,3.2 逻辑运算,3.2.1 逻辑运算符和运算法则,5种逻辑运算符：& (与)| (或) (非)& (先决与)| (先决或),逻辑运算的运算法则为：,1) 在逻辑运算中，确认非零元素为真，用1 表示，零元素为假，用0表示。 2) 设参与逻辑运算的是两个标量a和b，那么： a&b a,b全为非零时，运算结果为1，否则为0。 a|b a,b中只要有一个非零，运算结果为1。 a 当a是零时，运算结果为1；当a非零时，运算结果为0。,&(先决与) 当该逻
5、辑运算符的左边为1(真)时，才继续执行该符号右边的运算。 | (先决或) 当逻辑运算符的左边为1(真)时，就不需要继续执行该符号右边的运算，而立即得出该逻辑运算结果为1(真)；否则，就要继续执行该符号右边的运算。注意：&(先决与)和| (先决或)对标量计算无意义。,3) 若参与逻辑运算的是两个同型矩阵，那么运算将对矩阵相同位置上的元素按标量规则逐个进行。最终运算结果是一个与原矩阵同型的矩阵，其元素由1或0组成。 4) 若参与逻辑运算的一个是标量，一个是矩阵，那么运算将在标量与矩阵中的每个元素之间按标量规则逐个进行。最终运算结果是一个与矩阵同型的矩阵，其元素由1或0组成。 5) 逻辑非是单目运算
6、符，也服从矩阵运算规则。,水文应用：径流级别划分,3.2.2 逻辑函数,1 关系逻辑函数 1) anyany(x) 如果向量x中的有一个元素为非零值，那么返回1；否则，返回0。any(A) 对矩阵A进行列运算，根据相应列是否所有元素都为非零值，返回一个元素为1和0的行向量。 2）all all(x) 如果向量x中所有的元素都是非零值，那么返回1；否则，返回0。all(A) 对矩阵A进行列操作，根据相应列是否包含非零元素，返回元素为1和0的一个行向量。,isequal(A,B) 如果A和B是同型的，并且数组的元素也相同，则返回1；否则返回0。 isscalar(S) 如果S是1 x 1的矩阵（标
7、量），则返回1；否则返回0。 isvector(V) 如果V 是1 x n的矩阵（行向量）或n x 1的矩阵（列向量），则返回1；否则返回0。 isnumeric (A) 如果A是数值数组，则返回1；否则返回0。 isinteger(A) 如果A是整数型数值数组，则返回1；否则返回0。,isfloat(A) 如果A是浮点型数值数组，则返回1；否则返回0。 ischar(s) 判断s是否为字符串，返回1表示为字符串；否则，返回0。 isletter(s1) 在逐个字符的基础上比较字符的类别是否为字母，返回1表示为字母。 iscell(C) 判断C 是否为元胞数组，如果是，返回1；否则返回0。 i
8、sstruct(S) 如果S是构架数组，返回1；否则，返回0。,3.2.3 运算符优先级,MATLAB中各种运算符的优先级如下： (矩阵转置)、(矩阵幂)和.(数组转置)、.(数组幂)(逻辑非) *(乘)、/(左除)、(右除)和.*(点乘)、./(点左除)、.(点右除) +、-(加减) : (冒号) 、=、= &(逻辑与) |(逻辑或) &(先决与) |(先决或),3.3 M文件,3.3.1 M文件概述用MATLAB语言编写的程序，称为M文件。 M文件可以根据调用方式的不同分为两类：命令文件和函数文件。 两种工作方式：交互式的命令行工作方式和M文件的程序工作方式。,3.3.2 M文件的建立,1
9、 建立新的M文件为建立新的M文件，启动MATLAB文本编辑器有3种方法：1) 菜单操作。2) 命令操作。3) 命令按钮操作。,3.4 分支语句,3.4.1 if 结构,格式1： if 条件语句组end,格式2： if 条件语句组1else语句组2end,格式3： if 条件1语句组1elseif 条件2语句组2 elseif 条件m语句组melse语句组m+1end,例3-12 计算分段函数 的值。,3.4.2 switch-case结构,一般switch-case语句格式为：,switch 表达式 case 值1语句组1case 值2语句组2 case 值m语句组motherwise语句组m+1 end,例 3-13 求分段函数 的值。,
