MATLAB 之 M 文件和程序控制结构(顺序、选择、循环)-1

简介: M 文件是由若干 MATLAB 命令组合在一起构成的程序,其扩展名为 .m。它可以完成某些操作,也可以实现某种算法。实际上,MATLAB 提供的内部函数以及各种工具箱,都是利用 MATLAB 命令开发的 M 文件。

文章目录


一、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 编辑器窗口。
  • 994bb128eb7a4f33951e12aa791b4e48.png



也可以在 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 语句后面的语句。
  • 例如,我们计算分段函数的值1.png

这是一个具有两个分支的分段函数,为了求函数值,可以采用双分支结构来实现。程序如下:

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 部分是可选的,只有前面的条件不成立,才会进行后面的条件判断。

aaa6f17e00d74bc983c9846ccf456597.png



例如,我们输入一个字符,若为大写字母,则输出其对应的小写字母;若为小写字母,则输出其对应的大写字母;若为数字字符则输出其对应数的平方,若为其他字符则原样输出。

判断分为 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 子句后面的语句。

2ec5ef8201844b239847bdf75706b900.png



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 (.*)执行按元素相乘。'


相关文章
|
3月前
|
JSON 数据格式 Python
【2023最新】Matlab 保存JSON数据集文件,并用Python读取
本文介绍了如何使用MATLAB生成包含数据和标签的JSON格式数据集文件,并展示了用Python读取该JSON文件作为训练集的方法。
123 1
|
27天前
|
监控 算法 数据安全/隐私保护
基于三帧差算法的运动目标检测系统FPGA实现,包含testbench和MATLAB辅助验证程序
本项目展示了基于FPGA与MATLAB实现的三帧差算法运动目标检测。使用Vivado 2019.2和MATLAB 2022a开发环境,通过对比连续三帧图像的像素值变化,有效识别运动区域。项目包括完整无水印的运行效果预览、详细中文注释的代码及操作步骤视频,适合学习和研究。
|
2月前
|
监控 算法 安全
基于颜色模型和边缘检测的火焰识别FPGA实现,包含testbench和matlab验证程序
本项目展示了基于FPGA的火焰识别算法,可在多种应用场景中实时检测火焰。通过颜色模型与边缘检测技术,结合HSV和YCbCr颜色空间,高效提取火焰特征。使用Vivado 2019.2和Matlab 2022a实现算法,并提供仿真结果与测试样本。FPGA平台充分发挥并行处理优势,实现低延迟高吞吐量的火焰检测。项目包含完整代码及操作视频说明。
|
3月前
Matlab批量修改指定文件下文件名
Matlab批量修改指定文件下文件名
174 1
|
3月前
|
存储 算法 Serverless
【matlab】matlab基于DTW和HMM方法数字语音识别系统(源码+音频文件+GUI界面)【独一无二】
【matlab】matlab基于DTW和HMM方法数字语音识别系统(源码+音频文件+GUI界面)【独一无二】
|
3月前
|
存储 Serverless
【matlab】matlab实现倒谱法基音频率检测和共振峰检测(源码+音频文件)【独一无二】
【matlab】matlab实现倒谱法基音频率检测和共振峰检测(源码+音频文件)【独一无二】
|
5月前
|
算法 计算机视觉 异构计算
基于FPGA的图像一维FFT变换IFFT逆变换verilog实现,包含tb测试文件和MATLAB辅助验证
```markdown ## FPGA 仿真与 MATLAB 显示 - 图像处理的 FFT/IFFT FPGA 实现在 Vivado 2019.2 中仿真,结果通过 MATLAB 2022a 展示 - 核心代码片段:`Ddddddddddddddd` - 理论:FPGA 实现的一维 FFT/IFFT,加速数字信号处理,适用于高计算需求的图像应用,如压缩、滤波和识别 ```
|
5月前
|
存储 算法 计算机视觉
m基于FPGA的FIR低通滤波器实现和FPGA频谱分析,包含testbench和滤波器系数MATLAB计算程序
在Vivado 2019.2平台上开发的系统,展示了数字低通滤波器和频谱分析的FPGA实现。仿真结果显示滤波效果良好,与MATLAB仿真结果一致。设计基于FPGA的FIR滤波器,利用并行处理和流水线技术提高效率。频谱分析通过离散傅里叶变换实现。提供了Verilog核心程序以示例模块工作原理。
48 4
|
5月前
|
机器学习/深度学习 数据可视化 算法
探索MATLAB世界:掌握基础知识与实用技能(1. MATLAB环境与基本操作 2. 数据类型与变量 3. 条件与循环,1. 数据分析与统计 2. 图像处理与计算机视觉 3. 信号处理与控制系统)
探索MATLAB世界:掌握基础知识与实用技能(1. MATLAB环境与基本操作 2. 数据类型与变量 3. 条件与循环,1. 数据分析与统计 2. 图像处理与计算机视觉 3. 信号处理与控制系统)
53 0
|
5月前
|
机器学习/深度学习 算法
m基于PSO-GRU粒子群优化长门控循环单元网络的电力负荷数据预测算法matlab仿真
摘要: 在MATLAB 2022a中,对比了电力负荷预测算法优化前后的效果。优化前为&quot;Ttttttt111222&quot;,优化后为&quot;Tttttttt333444&quot;,明显改进体现为&quot;Tttttttttt5555&quot;。该算法结合了粒子群优化(PSO)和长门控循环单元(GRU)网络,利用PSO优化GRU的超参数,提升预测准确性和稳定性。PSO模仿鸟群行为寻找最优解,而GRU通过更新门和重置门处理长期依赖问题。核心MATLAB程序展示了训练和预测过程,包括使用&#39;adam&#39;优化器和超参数调整,最终评估并保存预测结果。
54 0