MATLAB程序设计95120.docVIP

MATLAB 程序设计 MATLAB提供了一个完善的程序设计语言环境，使用户能够方便地编制复杂的程序，完成各种计算。本节先介绍关系运算、逻辑运算，再介绍M-文件（即程序文件）的结构及MATLAB的程序控制流语句。 一、关系运算和逻辑运算 1．关系运算 （1）关系运算符： ； = ； ； = ；= = ；~ = （2）关系表达式： 用关系运算符将两个同类型的量（表达式）连接起来的式子。 【注】①关系运算本质上是标量运算，关系表达式的值是逻辑值（0-假1-真）； ②当作用于两个同样大小矩阵时，则分别对两个矩阵的对应元素运算，结果是一个0-1矩阵。 例1．对向量进行关系运算。 A=1:5,B=5:-1:1 % 输入向量 A = 1 2 3 4 5 B = 5 4 3 2 1 C=(A=4) % 对向量进行关系运算 C = 0 0 0 1 1 D=(A==B) % 对向量进行关系运算 D = 0 0 1 0 0 2．逻辑运算 （1）逻辑运算符： (and,与)、| (or,或)、~ (not,非) （2）逻辑表达式：用逻辑将两个逻辑量连接起来的式子。 【注】①逻辑运算本质上是标量运算，它将任何非零元素视为1（真）； ②当作用于两个同样大小矩阵时，则分别对两个矩阵的对应元素运算，结果是一个0-1矩阵。（真值表见P27） 例2．对向量进行逻辑运算。 a=1:9,b=9-a a = 1 2 3 4 5 6 7 8 9 b = 8 7 6 5 4 3 2 1 0 c=~(a4) % 非运算 c = 1 1 1 1 0 0 0 0 0 d=(a=3)(b6) % 与运算 d = 0 0 0 1 1 1 1 1 1 3．逻辑函数 any(x) 向量x中有非零元返回1，否则返回0。（向量函数） all(x) 向量x中所有元素非零返回1，否则返回0。（向量函数） isempty(x) 矩阵x为空矩阵返回1，否则返回0。 find(x) 返回向量x中非零元下标，若x是矩阵，则视其为一长列向量。 例3．对向量或矩阵进行逻辑运算，或调用逻辑函数。 a=[1 0 -5 0;-3 0 8 2] a = 1 0 -5 0 -3 0 8 2 b=all(a)，c=any(a) b = 1 0 1 0 c = 1 0 1 1 x=isempty(a) % 判定矩阵a是否为空，不能用关系式”a==[]”! x = 0 y=find(c) %返回向量c 中非零元下标 y = 1 3 4 z=find(a) % 1个输出参数 z = [1 2 5 6 8]’ [m,n]=find(a) % 2个输出参数 二、M-文件的结构 M-文件就是MATLAB程序文件,它是一个包含MATLAB语句组（命令序列）的普通ASCII文本文件，其扩展名为“.m”。M-文件有两类：脚本文件和函数文件。 1．脚本文件（亦称命令M-文件） 脚本文件的结构比较简单，它没有输入参数和输出参数，只是一些MATLAB命令行的组合。脚本文件中定义的变量都是全局变量。 例4．给定矩阵 , ,编写命令文件，输入矩阵A、B，并计算输出它们的和与差。 解：建立如下M-文件SY00504.m A=[1 2 3;4 5 6]; % 输入矩阵A B=[4 0 1;0 1 2]; % 输入矩阵B C=A+B % 计算输出A与B的和 D=A-B % 计算输出A与B的差 将上述文件以SY00504.m为名存盘，然后在命令窗口调用此M-文件，则有 SY00504 C = 5 2 4 D = -3 2 2 4 6 8 4 4 4 2．函数文件（亦称函数M-文件） 和命令文件相比,函数文件稍微复杂一些。通常函数文件包含以下几个部分： （1）函数