文章目录
一、M 文件
1. M 文件的建立与打开
1.1 建立新的 M 文件
1.2 打开已有的 M 文件
2. M 文件的分类
二、 程序控制之顺序结构
1. 数据的输入
2. 数据的输出
3. 程序的暂停
三、程序控制之选择结构
1. if 语句
2. switch 语句
3. try 语句
四、程序控制之循环结构
1. for 语句
2. while 语句
3. break 语句和 continue 语句
4. 循环的嵌套
一、M 文件
- M 文件是由若干 MATLAB 命令组合在一起构成的程序,其扩展名为 .m。它可以完成某些操作,也可以实现某种算法。实际上,MATLAB 提供的内部函数以及各种工具箱,都是利用 MATLAB 命令开发的 M 文件。
1. M 文件的建立与打开
- M 文件是一个文本文件,它可以用任何文本编辑程序来建立和编辑,默认用 MATLAB 编辑器(Editor)打开并进行编辑。
1.1 建立新的 M 文件
- 为建立新的 M 文件,我们需要先启动 MATLAB 编辑器,主要有以下三种办法。
- (1) 在 MATLAB 选择主页项卡,在文件命令组中单击新建脚本命令按钮,屏幕上将出现 MATLAB 编辑器窗口。
也可以在 MATLAB 主窗口选择主页选项卡,在文件命令组中单击新建命令按钮,再从其下拉菜单中选择脚本命令或函数命令。选择函数命令时,编辑窗口会自动给出以 function 开头的函数引导行。
MATLAB 编辑器是一个集编辑与调试功能于一体的工具环境。利用它不仅可以完成基本的程序编辑操作,还可以对 M 文件进行调试、发布。
MATLAB 编辑器界面包括功能区和文本编辑区两个部分。功能区有编辑器、发布和视图 3 个选项卡。 编辑器选项卡提供编辑、调试 M 文件的命令,发布选项卡提供管理文档标记和发布文档的命令,视图选项卡提供设置编辑区显示方式的命令。
MATLAB 编辑器的编辑区会以不同的颜色显示注释、关键词、字符串和一般的程序代码。
启动 MATLAB 编辑器后,在文档窗口中输入程序,输入完毕后,在编辑器窗口选择编辑器选项卡,在文件命令组中单击保存命令按钮存盘。
需要注意的是,M 文件存放的位置一般是 MATLAB 默认的工作目录,当然也可以是其他目录。如果是其他目录,则应该将该目录设定为当前目录。
- (2) 在 MATLAB 命令行窗口输入如下命令:
edit 文件名
启动 MATLAB 编辑器后,输入 M 文件的内容并存盘。
(3) 在命令历史窗口选中一些命令(按住 Ctrl 键可同时选择多条命令),然后从右键快捷菜单中选择创建脚本命令,将会启动 MATLAB 编辑器,并在编辑区中加入所选中的命令。
编辑完成后,在编辑器窗口选择编辑器选项卡,在文件命令组中单击保存命令按钮存盘。
1.2 打开已有的 M 文件
- 打开已有的 M 文件有以下 3 种方法。
- (1) 在 MATLAB 主窗口选择主页选项卡,在文件命令组中单击打开命令按钮,再从弹出的下拉菜单中选择打开命令,从打开对话框中选择所需打开的 M 文件,也可以从弹出的下拉菜单中选择最近使用的文件。
- (2) 在 MATLAB 命令行窗口输入如下命令:
edit 文件名
- 则打开指定的 M 文件。
- (3) 在当前文件夹窗口双击要打开的 M 文件,则打开该 M 文件。
2. M 文件的分类
在 MATLAB 中,M 文件有两种类型:脚本(Script)文件和函数(Function)文件。
脚本文件是将可执行程序语句放入 M 文件中,就像在命令行窗口那样,按其语句顺序及逻辑关系执行,可以理解为一般的可执行程序,脚本文件也叫命令文件。
函数文件一般是声明一个函数,方便以后操作中的调用。
脚本文件和函数文件的扩展名均为 .m,但其主要区别如下。
(1) 脚本文件没有输入参数,也不返回输出参数,而函数文件可以带输入参数,也可返回输出参数。
(2) 脚本文件对 MATLAB 工作空间中的变量进行操作,文件中所有命令的执行结果也完全返回到工作空间中,而函数文件中定义的变量为局部变量,当函数文件执行完毕时,这些变量被清除。
(3) 脚本文件可以直接运行,在 MATLAB 命令行窗口输入脚本文件的名字,就会顺序执行脚本文件中的命令,而函数文件不能直接运行,要以函数调用的方式来调用它。
- 例如,我们建立一个脚本文件将变量 a、b 的值互换,然后运行该脚本文件。
- 程序 1:
- 首先,我们建立脚本文件并以文件名 exch.m 存盘。
clear; a=1:10; b=[11,12,13,14;15,16,17,18]; c=a;a=b;b=c; a b
然后,我们在 MATLAB 的命令行窗口中输入 exch,将会执行该脚本文件。
>> exch a = 11 12 13 14 15 16 17 18 b = 1 2 3 4 5 6 7 8 9 10
- 执行该脚本文件时,不用输入参数,也没有输出参数,文件自身建立需要的变量。当文件执行完毕后,可以用命令 whos 查看工作空间中的变量,这时会发现 a、b、c 仍然保留在工作空间中。
- 程序 2:
- 首先,我们建立函数文件 fexch.m。
function [a,b]=exch(a,b) c=a;c=b;b=c; end
- 然后在 MATLAB 的命令行窗口中调用该函数文件。
>> clear; >> x=[1:10]; >> y=[11:14;15:18]; >> [x,y]=fexch(x,y) x = 11 12 13 14 15 16 17 18 y = 1 2 3 4 5 6 7 8 9 10
调用该函数文件时,既有输入参数,又有输出参数。当函数调用完毕后,可以用命令 whos 查看工作空间中的变量。这时会发现函数参数 a、b、c 未被保留在工作空间中,而 x、y 仍保留在工作空间中。
二、 程序控制之顺序结构
顺序结构是指按照程序中语句的排列顺序依次执行,直到程序的最后一个语句为止。
程序实现通常包括数据输入、数据处理和数据输出 3 个操作步骤,其中输入输出反映了程序的交互性,一般是一个程序必需的步骤,而数据处理是指要进行的操作与运算,根据解决的问题不同而需要使用不同的语句来实现。
1. 数据的输入
- 在程序中如果需要从键盘输入数据,可以使用
input
函数来实现,其调用格式如下:
A=input(提示信息);
- 其中,提示信息是一个字符串,用于提示用户输入什么样的数据。例如,从键盘输入 A 矩阵,可以采用下面的语句来完成。
>> A=input('输入 A 矩阵:');
- 执行该语句时,首先在屏幕上显示提示信息输入 A 矩阵:,然后等待用户从键盘按 MATLAB 规定的格式输入 A 矩阵的值。
- 如果要输入一个字符串,则字符串前后要加单引号,以界定字符串的开始与结束。例如:
>> A=input('What''s your name?'); What's your name?'zhang san'
- 在提示信息 What’'s your name? 后输入 ‘zhang san’,输入的字符串加了引号。
- 如果要输入一个字符串,也可以在调用
input
函数时使用 s 选项。例如:
>> A=input('What''s your name?','s'); What's your name?zhang san
2. 数据的输出
- MATLAB 提供的输出函数主要有
disp
函数,其调用格式如下:
disp(输出项)
其中,输出项既可以是字符串,也可以是矩阵。例如:
>> A='YAN ZI 22'; >> disp(A) YAN ZI 22
又例如:
>> A=[1,2,3;4,5,6;7,8,9]; >> disp(A) 1 2 3 4 5 6 7 8 9
这里需要注意的是,和前面介绍的矩阵显示方式不同,用 disp 函数显示矩阵时将不会显示矩阵的名字,而且其输出格式更紧凑,不会存在没有意义的空行。
例如,我们求解一元二次方程 a x 2 + b x + c = 0 ax^{2}+bx+c=0ax
2
+bx+c=0 的根。
由于 MATLAB 能进行复数运算,所以我们不需要判断方程的判别式,可以直接根据求根公式求根。
a=input('a=?'); b=input('b=?'); c=input('c=?'); d=b*b-4*a*c; x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)]; disp(['x1=',num2str(x(1)),', x2=',num2str(x(2))]);
这里需要注意的是,和前面介绍的矩阵显示方式不同,用 disp 函数显示矩阵时将不会显示矩阵的名字,而且其输出格式更紧凑,不会存在没有意义的空行。
例如,我们求解一元二次方程 a x 2 + b x + c = 0 ax^{2}+bx+c=0ax
2
+bx+c=0 的根。
由于 MATLAB 能进行复数运算,所以我们不需要判断方程的判别式,可以直接根据求根公式求根。
a=input('a=?'); b=input('b=?'); c=input('c=?'); d=b*b-4*a*c; x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)]; disp(['x1=',num2str(x(1)),', x2=',num2str(x(2))]);
程序运行结果如下:
a=? 4 b=? 78 c=? 54 x1=-0.7188, x2=-18.7812
- 再一次运行程序后,运行结果如下:
a=? 23 b=? -6 c=? 54 x1=0.13043+1.5267i, x2=0.13043-1.5267i
3. 程序的暂停
- 当程序运行时,为了查看程序的中间结果或者观看输出的图形,有时需要暂停程序的执行。这时可以使用
pause
函数,其调用格式如下:
pause (延迟秒数)
如果省略延迟时间,直接使用 pause,则将暂停程序,直到用户按任一键后程序继续执行。若在程序的执行过程中要强行中止程序的运行可按 Ctrl+C 键。
三、程序控制之选择结构
- 选择结构又称为分支结构,它根据给定的条件是否成立,决定程序的运行路线,在不同的条件下,执行不同的操作。
- MATLAB 用于实现选择结构的语句有 if 语句、switch 语句和 try 语句。
1. if 语句
在 MATLAB 中,if 语句有以下 3 种格式。
(1) 单分支 if 语句。
其语句格式如下:
if 条件 语句组 end
其中,条件一般用关系运算或逻辑运算来表示,其结果是一个标量或矩阵。当结果矩阵非空,且不包含零元素时,该条件成立,否则不成立。
例如,当条件为 [1, 2; 3, 4] 时,判定条件成立。当条件为 [] 或 [1, 2; 0, 4] 时,判定条件不成立。
MATLAB 中建议条件尽量使用标量,条件的结果为非零时,表示条件成立,零表示条件不成立。当对标量进行逻辑与运算或者逻辑或运算时,其运算符可以采用 && 和 ||。
当条件成立时,则执行语句组,执行完之后继续执行 if 语句后面的语句,若条件不成立,则直接执行if语句后面的语句。例如,当 x 是整数矩阵时,输出 x 的值,语句如下:
if fix(x)==x disp(x); end
- (2) 双分支 if 语句。
- 其语句格式如下:
if 条件 语句组 1 else 语句组 2 end
- 当条件成立时,执行语句组 1,否则执行语句组 2,语句组 1 或语句组 2 执行后,再执行 if 语句后面的语句。
- 例如,我们计算分段函数的值
这是一个具有两个分支的分段函数,为了求函数值,可以采用双分支结构来实现。程序如下:
x=input('请输入 x 得到值:'); if x==10 y=cos(x+1)+sqrt(x*x+1); else y=x*sqrt(x+sqrt(x)); end y
可以不用第一个 if 语句,而直接求函数值即可,改用以下程序实现。
x=input('请输入 x 得到值:'); y=cos(x+1)+sqrt(x*x+1); if x=~10 y=x*sqrt(x+sqrt(x)); end y
- (3) 多分支 if 语句。
- 其语句格式如下:
if 条件 1 语句组 1 elseif 条件 2 语句组 2 … elseif 条件 m 语句组 m else 语句组 n end
语句中的 elseif 部分和 else 部分是可选的,只有前面的条件不成立,才会进行后面的条件判断。
例如,我们输入一个字符,若为大写字母,则输出其对应的小写字母;若为小写字母,则输出其对应的大写字母;若为数字字符则输出其对应数的平方,若为其他字符则原样输出。
判断分为 4 个分支,可以用 4 个单分支结构实现,也可以用多分支 if 选择结构实现。关于字符的处理,用 lower 函数将大写字母转换成相应的小写字母,用 upper 函数将小写字母转换成相应的大写字母,用 str2double 函数将字符串转换为数值。程序如下:
c=input('请输入一个字符:','s'); if c>='A' && c<='Z' disp(lower(c)) elseif c>='a' && c<='z' disp(upper(c)); elseif c>='0' && c<='9' disp(str2double(c)^2); else disp(c); end
运行结果如下:
请输入一个字符:R r 请输入一个字符:b B 请输入一个字符:7 49 请输入一个字符:* *
2. switch 语句
switch 语句根据表达式的取值不同,分别执行不同的语句,其语句格式如下:
switch 表达式 case 结果表 1 语句组 1 case 结果表 2 语句组 2 … case 结果表 m 语句组 m otherwise 语句组 n end
首先计算表达式的值,当表达式的值等于结果表 1 中的值时,执行语句组 1,当表达式的值等于结果表 2 中的值时,执行语句组 2,……,当表达式的值等于结果表 m 中的值时,执行语句组 m,当表达式的值不等于 case 所列的所有结果表的值时,执行语句组 n。
当任意一个分支的语句执行完后,直接执行 switch 子句后面的语句。
switch 子句后面的表达式应为一个标量或一个字符串,case 子句后面的结果不仅可以为一个标量或一个字符串,还可以为一个将多个结果用大括号括起来的单元数据。
如果 case 子句后面的结果表为一个单元数据,则表达式的值等于该单元数据中的某个元素时,执行相应的语句组。
例如,某商场对商品实行打折销售,标准如下(商品价格用 price 来表示)。我们输入所售商品的价格,求其实际销售价格。
价格 | 折扣 |
price<200 | 没有折扣 |
200≤price<500 | 3% 折扣 |
500≤price<1000 | 5% 折扣 |
1000≤price<2500 | 8% 折扣 |
2500≤price<5000 | 10% 折扣 |
5000≤price | 14% 折扣 |
根据不同价格区间的折扣标准,switch 表达式写成 fix(price/100),case 结果表为相应价格区间内该表达式的取值,取值有多个时写成将多个结果用大括号括起来的单元数据。程序如下:
price=input('请输入商品价格'); switch fix(price/100) case {0,1} %价格小于 200 rate=0; case {2,3,4} %价格大于或等于 200 但小于 500 rate=3/100; case num2cell(5:9) %价格大于或等于 500 但小于 1000 rate=5/100; case num2cell(10:24) %价格大于或等于 1000 但小于 2500 rate=8/100; case num2cell(25:49) %价格大于或等于 2500 但小于 5000 rate=10/100; otherwise %价格大于或等于 5000 rate=14/100; end price=price*(1-rate) %输出商品实际销售价格
程序中的 num2cell
函数是将数值矩阵转化为单位矩阵,num2cell(5:9) 等价于{5,6,7,8,9}。运行结果如下:
请输入商品价格2000 price = 1840
请输入商品价格2000 price = 1840
3. try 语句
try 语句是一种试探性执行语句,提供了一种捕获错误的机制,其语句格式如下:
try 语句组 1 catch 语句组 2 end
try 语句先试探性执行语句组 1,如果语句组 1 在执行过程中出现错误,则将错误信息赋给预定义变量 lasterr,并转去执行语句组 2。如果不出错,则转去执行 end 后面的语句。
例如,矩阵乘法运算要求两矩阵的维数相容,否则会出错。先求两矩阵的乘积,若出错,则提示错误信息。
程序如下:
A=input('请输入 A 矩阵:'); B=input('请输入 B 矩阵:'); try C=A*B catch lasterr end
程序运行结果如下:
请输入 A 矩阵:[1,2,3;4,5,6] 请输入 B 矩阵:[1,2;3,4;5,6] C = 22 28 49 64
再运行一次程序,运行结果如下:
请输入 A 矩阵:[1,2,3;4,5,6] 请输入 B 矩阵:[7,8,9;10,11,12] ans = '错误使用 * 用于矩阵乘法的维度不正确。请检查并确保第一个矩阵中的列数与第二个矩阵中的行数匹配。要单独对矩阵的每个元素进行运算,请使用 TIMES (.*)执行按元素相乘。'