第5章 MALAB程序设计.docVIP

第5章 MATLAB程序设计 1 5.1脚本文件和函数文件 1 5.1.1 M文本编辑器 2 5.1.2 M文件的基本格式 3 5.1.3 M脚本文件（Script File） 6 5.1.4 M函数文件(Function File) 9 5.2程序流程控制 11 5.2.1 for ... end循环结构 11 5.2.2 while ... end循环结构 14 5.2.3 If…else…end条件转移结构 15 5.2.4 switch…case开关结构 17 5.2.5 try... catch... end试探结构 19 5.2.6流程控制语句 21 5.3函数调用和参数传递 26 5.3.1子函数和私有函数 26 5.3.2局部变量和全局变量 28 5.3.3函数的参数 30 5.3.4程序举例 37  第5章 MATLAB程序设计 MATLAB除了可以在命令窗口中编写行命令外，作为一种高级应用软件还可以生成自己的程序文件。 5.1脚本文件和函数文件 MATLAB程序代码所编写的文件通常以“.m”为扩展名，因此这些文件称为M文件。 M文件是一个ASCII码文件，可以用任何文本处理软件来编写，M文件在第一次运行时由于逐句解释运行程序，故速度比编译型的慢，单M文件一经运行就将编译代码放在内存中，再次运行的速度就大大加快。 M文件有两种形式：M脚本文件和M函数文件。M函数文件是MATLAB程序设计的主流，一系列工具箱的内部函数都是M函数文件。 5.1.1 M文本编辑器 MATLAB的M文件是通过M文件编辑／调试器窗口(Editor／Debugger)来创建的。 单击MATLAB桌面上的图标，或者单击菜单“File”——“New”——“M-file”，可打开空白的M文件编辑器，也可以通过打开已有的M文件来打开M文件编辑器。如图5.1所示为打开已创建的M文件。 5.1.2 M文件的基本格式 M脚本文件和M函数文件的格式不同。 下面介绍绘制二阶系统时域曲线的M文件，欠阻尼系统的时域输出y与x的关系为，【例5.1】为M脚本文件，【例5.2】为M函数文件。 【例5.1】用M脚本文件绘制二阶系统时域曲线。 %EX0501 二阶系统时域曲线 %画阻尼系数为0.3的曲线 x=0:0.1:20; y1=1-1/sqrt(1-0.3^2)*exp(-0.3*x).*sin(sqrt(1-0.3^2)*x+a cos(0.3)) plot(x,y1,r) 【例5.2】创建一个画二阶系统时域曲线的函数，阻尼系数zeta为函数的输入参数。 function y=Ex0502(zeta) % EX0502 Step response of quadratic system. % 二阶系统时域响应曲线 % zeta 阻尼系数 % y 时域响应 % % copyright 2003-08-01 x=0:0.1:20; y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta)) plot(x,y) M函数文件的基本格式： 函数声明行 H1行(用%开头的注释行) 在线帮助文本(用%开头) 编写和修改记录(用%开头) 函数体? 说明： 函数声明行是M函数所必须的，M脚本文件没有函数声明；函数名与文件名保持一致，当不一致时，MATLAB以文件名为准，上例中，Ex0502函数保存为Ex0502.m文件。 H1行通常包含大写的函数文件名，可以提供给help和lookfor关键词使用。 在线帮助文本通常包含函数输入、输出变量的含义、格式说明。 编写和修改记录一般在空1行后，记录作者、日期和版本记录，用于软件档案管理。 例如，在命令窗口输入help和lookfor命令查看帮助信息： help Ex0502 EX0502 Step response of quadratic system. 二阶系统时域响应曲线 zeta 阻尼系数 y 时域响应 lookfor 二阶系统时域响应 Ex0502.m: %二阶系统时域响应 Help命令显示M文件的第一个连续注释块，lookfor命令显示第一行注释。 函数体由实现M函数文件功能的MATLAB命令组成。 5.1.3 M脚本文件（Script File） 脚本文件（命令文件）：实际上是一串指令的集合，与在命令窗口逐行执行文件中的所有指令，其结果是一样的。没有输入输出参数。 脚本文件的特点： (1) 脚本文件中的命令格式和前后位置，与在命令窗口中输入的没有任何区别。