MATLAB程序设计
Matlab程序设计打印版
实验四 Matlab程序设计[实验目的]1．掌握字符串数组的创建和构造方法及常用字符串函数的使用。
2．熟练掌握MATLAB 控制流的使用方法。
3．熟悉M 脚本文件、函数文件的编写方法和技巧。
[实验原理]与数值数组相比，串数组在MATLAB 中的重要性较小，但不可缺少。
如果没有串数组及相应的操作，那么数据可视化、图形用户界面的制作将会遇到困难。
字符串与数值数组是两种不同的数据类，它们的创建方式也不同。
字符串的创建方式是：将待建的字符放在“单引号对”中。
注意，“单引号对”必须是在英文状态下输入，其作用是MATLAB 识别送来内容“身份”所必需的，如A=’This is an example!’；就创建了一个字符串A。
注意创建带单引号的字符串时，每个单引号符用“连续2 个单引号符”标识。
字符串的标识同数值数组同，而且也可以使用size 指令观察串数组的大小。
串数组的ASCII 码可以通过指令abs 和double 来获取，而用char 指令可以把ASCII 码变为串数组，另外，MATLAB 可以很好的支持中文字符串数组。
对于复杂串数组的创建，一是可以直接创建，但是要保证同一串数组的各行字符数相等，即保证各行等长，不推荐，太繁琐。
二是可以利用串操作函数创建多行数组，比如char, str2mat, strvcat 等，具体操作自己通过帮助体会。
另外还可以通过转化函数产生数码字符长，比如A_str=int2str(A) 就是把整数数组A 转换成串数组，如果是非整数将被四舍五入后再转换，类似的函数还有num2str（把非整数数组转换为串数组，常用于图形中数据点的标识）、mat2str（把数值数组转换成输入形态的串数组，常与eval 指令配用）。
假如想灵活运用MATLAB 去解决实际问题，想充分调动MATLAB——科学技术资源，想理解MATLAB 版本升级所依仗的基础，那么掌握M 脚本文件合函数的编写规则将十分有用。
MATLAB程序设计
MATLAB程序设计MATLAB提供了一个完善的程序设计语言环境，使用户能够方便地编制复杂的程序，完成各种计算。
本节先介绍关系运算、逻辑运算，再介绍M-文件（即程序文件）的结构及MATLAB的程序控制流语句。
一、关系运算和逻辑运算1．关系运算（1）关系运算符：< ；< = ；> ；> = ；= = ；~ =（2）关系表达式：用关系运算符将两个同类型的量（表达式）连接起来的式子。
【注】①关系运算本质上是标量运算，关系表达式的值是逻辑值（0-假1-真）；②当作用于两个同样大小矩阵时，则分别对两个矩阵的对应元素运算，结果是一个0-1矩阵。
例1．对向量进行关系运算。
>> A=1:5,B=5:-1:1 % 输入向量A = 1 2 3 4 5B = 5 4 3 2 1>> C=(A>=4) % 对向量进行关系运算C = 0 0 0 1 1>> D=(A==B) % 对向量进行关系运算D = 0 0 1 0 02．逻辑运算（1）逻辑运算符：& (and,与)、| (or,或)、~ (not,非)（2）逻辑表达式：用逻辑将两个逻辑量连接起来的式子。
【注】①逻辑运算本质上是标量运算，它将任何非零元素视为1（真）；②当作用于两个同样大小矩阵时，则分别对两个矩阵的对应元素运算，结果是一个0-1矩阵。
（真值表见P27）例2．对向量进行逻辑运算。
>> a=1:9,b=9-aa = 1 2 3 4 5 6 7 8 9b = 8 7 6 5 4 3 2 1 0>> c=~(a>4) % 非运算c = 1 1 1 1 0 0 0 0 0>> d=(a>=3)&(b<6) % 与运算d = 0 0 0 1 1 1 1 1 13．逻辑函数any(x) 向量x 中有非零元返回1，否则返回0。
（向量函数） all(x) 向量x 中所有元素非零返回1，否则返回0。
第二讲MATLAB的程序设计和M文件
第二讲MATLAB的程序设计和M文件在MATLAB中，程序设计主要是通过编写和运行M文件来实现。
M文件是MATLAB的主要代码文件，用于实现不同功能的程序，可以包含变量定义、算法实现、函数调用和图形绘制等。
下面将介绍MATLAB的程序设计和M文件的基本知识。
MATLAB的程序设计主要包括如下几个方面：1.变量和数据类型：在M文件中，可以使用不同的变量类型存储数据，如整型、浮点型、字符型等。
不同的变量类型在MATLAB中有不同的表示方法和功能。
2.运算和表达式：在M文件中，可以使用常见的数学运算符（如加减乘除、指数等）进行计算。
同时，也可以使用逻辑运算符（如与、或、非等）进行逻辑运算。
3. 条件语句：在M文件中，可以使用条件语句实现根据不同的条件执行不同的操作。
MATLAB中的条件语句主要有if语句和switch语句，通过判断条件的真假来决定执行路径。
4. 循环语句：在M文件中，可以使用循环语句实现对一段代码的重复执行。
MATLAB中的循环语句主要有for循环、while循环和do-while循环。
5.函数的定义和调用：在M文件中，可以通过定义函数来实现特定功能的封装。
函数可以包含输入参数和输出参数，通过参数的传递来实现数据的交互和函数的调用。
1.M文件的命名和保存：M文件的命名应该具有描述性，以体现文件中代码的功能。
M文件的保存格式是以.m为文件扩展名。
2.M文件的结构：在M文件中，一般会包含变量定义、函数定义和主程序等部分。
变量定义部分用于声明和初始化变量，函数定义部分用于定义自定义函数，而主程序部分用于调用函数和执行主要功能。
3.代码注释：为了提高代码的可读性和可维护性，需要在M文件中添加注释。
注释可以用于解释代码的目的和思路，以及描述变量、函数和算法等的作用和实现方法。
4.代码风格：为了代码的一致性和可读性，需要遵循一定的代码风格规范。
例如，可以在操作符周围留有空格，使用一致的缩进和命名规则，避免使用不必要的缩写等。
MATLAB程序设计及应用实例
注意：等号左边是方括号，右边输入参数用括号
第14页，共67页。
例1 编写函数文件求半径为r的圆的面积和周长
function [s,p]=fcircle(r)
%FCIRCLE calculate the area and perimeter of a circle of radii r
disp函数
调用格式为
disp(输出项) 其中输出项既可以为字符串，也可以为矩阵。 注意：用disp函数显示矩阵时将不显示矩阵的 名字，而且其格式更紧密，且不留任何没有意 义的空行。
第22页，共67页。
例 求一元二次方程ax2+bx+c=0的根。
程序如下： 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)];
c=input('请输入一个字符','s'); if c>='A' & c<='Z'
disp(setstr(abs(c)+1)); elseif c>='a'& c<='z'
disp(setstr(abs(c)-1)); elseif c>='0'& c<='9'
disp(abs(c)-abs('0')); else
end 命令文件main5_2.m：
matlab 程序设计
【例5.8】用try... catch... end结构来进行矩阵相乘运算. 例
% EX0508 try结构 n=4; a=magic(n); m=3; b=eye(3); try c=a*b catch c=a(1:m,1:m)*b end lasterr
5.2.6 流程控制语句
break, continue, return, pause, keyboard, input 1. break命令 命令 break命令可以使包含break的最内层的for或while 语句强制终止,立即跳出该结构,执行end后面的命令, break命令一般和if结构结合使用.
5.2.4 switch…case开关结构
语法: switch 开关表达式 case 表达式 表达式1 语句段1 语句段 case表达式2 case表达式2 表达式 语句段2 语句段 ... otherwise 语句段n 语句段 end
说明: (1) 将开关表达式依次与case后面的表达式进行比较,如 果表达式1不满足,则与下一个表达式2比较,如果都不 满足则执行otherwise后面的语句段n;一旦开关表达式 与某个表达式相等,则执行其后面的语句段. (2) 开关表达式只能是标量或字符串. (3) case后面的表达式可以是标量,字符串或元胞数组, 如果是元胞数组则将开关表达式与元胞数组的所有元素 进行比较,只要某个元素与开关表达式相等,就执行其 后的语句段.
(2) 将函数文件保存为"Ex0502.m". (3) 在MATLAB命令窗口输入以下命令,则会出现f的计算值 注意: 注意:M脚本文件和M函数文件的文件名及函数名的命名规 和绘制的曲线:f=Ex0502(0.3) 则与MATLAB变量的命名规则相同.
MATLAB语言程序设计基础
3.2 matlab语言基本运算及输入输出
3.2.5 输入与输出语句
input A=input(提示字符串）要求输
入矩阵
A=input(提示字符串，‘s’） 要求字符串eg:
n=input('how much')
n=input('ho第w23页m/共4u1页ch','s')
3.2 matlab语言基本运算及输入输出
关系运算和逻辑运
算
表3-6 关系运算和逻辑运算函数
函数 any all find
exist isnan
意义 逻辑条件任何一个
逻辑条件全部 寻找逻辑值的向量元素下 标
检查某变量是否存在 检查非数值量
函数 finite isempty isstr
strcmp
随机数元素矩阵 设三维绘图基底坐
第5页/共41页
单位矩阵
3.1.3构造多维数组
cat( ) a=cat(n,a1,a2,….) n:多维函数的维数 n=1:
cat(a1,a2,a3..)=[a1;a2;a3…] n=2:
cat(a1,a2,a3..)=[a1,a2,a3…] n=3: 图3－1示
重新定义维数
end 表示某一维末尾元素下标
2、复数矩阵： b=[1 2;3 4]+i*[5 6;7 8] b=[1+5i 2+6i;3+7i 4+8i]
第3页/共41页
3、空矩阵［］
0×0阶
与clear不同之处：clear删除变量
［］删除矩阵中的元素
A(:,[2,3])=[] 第2,3列元素删除
函数 abs angle sqrt real imag conj round fix
Matlab程序设计
19
例：计算2个数的和，根据输入的参数个数不同使用不同 的运算表达式。 function [sum]=Ex0514(x,y) if nargin==1 sum=x+0; else if nargin==0 sum=0; else sum=x+y end end [sum]=Ex0514(1,2)
可以在某一目录中建立一个private目录来存放相关函数，具有性质：在 私有目录下的私有函数，只能被其父目录的M函数文件调用，对其他目 录的文件私有函数是不可见的，可以和其他目录下的函数重名；私有函 数父目录的M脚本文件也不可调用私有函数；在函数调用搜索时，私有 函数优先于其他MATALB路径上的函数。
局部变量和全局变量
（1）局部变量是在函数体内部使用的变量，其影响范围只 能在本函数内，仅存在于函数的工作空间内，只在函数执行 期间存在，函数执行完变量消失。
（2）全局变量是可以在不同函数工作空间和MATLAB工作 空间中共享使用的变量。在使用前必须用global定义，而且 每个要共享全局变量的函数和工作空间，都必须逐个用 global对变量加以定义。
m文件模式
m文件编辑器缺省设置是(通过菜单Preferences) c:\windows\notepad.exe 还可设成：c:\windows\write.exe c:\windows\command\ 注意：所创建的m文件必须纳入matlab搜索路 径后，才能在matlab命令窗口运行。
实战matlab并行程序设计
实战matlab并行程序设计MATLAB并行程序设计是一种利用MATLAB的并行计算工具箱来提高计算效率的技术。
在处理大规模数据集或复杂算法时，传统的串行计算方式可能效率低下，而并行程序设计可以显著提升处理速度。
以下是MATLAB并行程序设计的一些关键概念和步骤。
1. 并行计算基础并行计算指的是同时使用多个处理器或核心来执行计算任务。
MATLAB并行程序设计主要依赖于MATLAB的Parallel Computing Toolbox。
2. 环境配置在开始并行程序设计之前，需要确保MATLAB安装了Parallel Computing Toolbox，并且计算机上安装了MATLAB的并行计算服务器。
3. 并行池的创建和管理并行池是执行并行任务的工作单元集合。
MATLAB提供了`parpool`和`delete(gcp('nocreate'))`函数来创建和管理并行池。
```matlab% 创建并行池parpool;% 删除并行池delete(gcp('nocreate'));```4. 并行计算函数MATLAB提供了多种并行计算函数，如`parfor`、`parfeval`、`parfor`等。
`parfor`是最常见的并行循环，用于并行化for循环。
```matlab% 并行for循环parfor i = 1:N% 执行并行任务end```5. 数据分配在并行程序设计中，数据需要在不同的工作单元之间分配。
MATLAB提供了`distribute`函数来帮助实现数据的自动分配。
```matlab% 数据分配data = distribute(largeArray);```6. 并行变量在并行环境中，MATLAB支持两种类型的变量：共享变量和私有变量。
共享变量在所有工作单元中共享，而私有变量则在每个工作单元中独立。
```matlab% 声明共享变量sharedVar = 0;% 声明私有变量parfor i = 1:NlocalVar = rand;end```7. 并行任务的同步在并行程序设计中，有时需要同步不同工作单元的状态。
《MATLAB程序设计》实验指导书
三、
实验仪器和设备
1、 计算机一台。 2、 MATLAB7.0 以上集成环境。
四、
预习要求
2
《MATLAB 程序设计》实验指导书
1、 复习 MATLAB 的启动与Байду номын сангаас出，熟悉 MATLAB 运行环境。 2、 复习 MATLAB 中矩阵的生成以及矩阵运算的基本原理。
五、
实验内容及步骤
实验内容：
1、 求下列表达式的值 1） z1
5
《MATLAB 程序设计》实验指导书
2） 函数调用 函数文件编制好后，就可调用函数进行计算了。函数调用的一般格式为 [输出实参表]=函数名（输入实参表） 注意：函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个 数一致，否则会出错。函数调用时，实参先传递给形参，然后再执行函数功能。 5、 选择结构 1） if 语句 a） 单分支 if 语句 if 条件 语句块 end b） 双分支 if 语句 if 条件 语句块 1 else 语句块 2 end c） 多分支 if 语句 if 条件 1 语句块 1 elseif 条件 2 语句块 2 …… elseif 条件 n 语句块 n else 语句块 n+1 end 2） switch 语句 switch 表达式 case 结果表 1 语句块 1 case 结果表 2 语句块 2 ……
二、
实验原理
1、 M 文件 用 MATLAB 语言编写的程序，称为 M 文件，它们的扩展名均为.m。M 文件根据 调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file） 。 2、 建立新的 M 文件 启动 MATLAB 文本编辑器有 3 种方法： 1） 单击工具栏上的“New M-File”命令按钮。 2） 从 MATLAB 主窗口的“File”菜单中选择“New”菜单项，再选择“M-file”命令。 3） 在 MATLAB 命令窗口输入命令“edit”。 3、 打开已有的 M 文件 1） 在当前目录窗口选中要打开的 M 文件，双击鼠标左键。 2） 单击 MATLAB 主窗口工具栏上的“Open File”命令按钮，再从弹出的对话框中 选择所需打开的 M 文件。 3） 从 MATLAB 主窗口的“File”菜单中选择“Open”命令，在“Open”对话框中选中 所需打开的文件。 4） 在 MATLAB 命令窗口输入命令“edit 文件名” 。 4、 函数文件 1） 函数文件的基本结构 函数文件由 function 语句引导，其基本结构为： function 输出形参表=函数名（输入形参表） 注释说明部分 函数体语句 注意：函数名的命名规则与变量名相同。当输出形参多于一个时，应用方括号 括起来。
matlab程序课程设计
matlab程序 课程设计一、课程目标知识目标：1. 理解Matlab程序的基本结构和编程语法；2. 掌握Matlab中数据类型、运算符和流程控制的使用方法；3. 学会使用Matlab进行数据处理、图像绘制和仿真计算；4. 了解Matlab在实际工程和科研中的应用。
技能目标：1. 能够运用Matlab编写简单的程序，解决实际问题；2. 能够运用Matlab进行数据分析和图像绘制，展示实验结果；3. 能够运用Matlab进行简单的科学计算和模型仿真；4. 能够在团队项目中，运用Matlab编程技能，完成所承担的任务。
情感态度价值观目标：1. 培养学生对Matlab编程的兴趣和热情，激发学生学习编程的积极性；2. 培养学生具备良好的编程习惯，注重代码的规范性和可读性；3. 培养学生面对编程问题时的耐心和毅力，克服困难，解决问题；4. 培养学生具备团队协作精神，能够在项目中发挥各自优势，共同完成任务。
课程性质：本课程为实践性较强的学科，注重培养学生实际操作能力和编程思维。
学生特点：学生具备一定的数学基础和计算机操作能力，但对Matlab编程尚处于入门阶段。
教学要求：结合学生特点，通过讲解、演示、实践等方式，引导学生掌握Matlab编程技能，并在实际应用中提高学生的编程水平。
教学过程中，注重培养学生的学习兴趣和团队协作能力，为后续相关课程和实际工作打下基础。
将课程目标分解为具体的学习成果，以便进行教学设计和评估。
二、教学内容1. Matlab基础知识- Matlab简介与安装- Matlab界面与基本操作- Matlab编程规范与技巧2. 数据类型和运算符- 基本数据类型：数值、字符串、逻辑等- 运算符：算术、关系、逻辑运算符3. 程序控制结构- 顺序结构- 选择结构：if-else-end- 循环结构：for、while循环4. 函数与脚本- 函数的定义与调用- 输入输出参数的使用- 脚本文件的编写与执行5. 数据处理与分析- 数组与矩阵操作- 数据可视化：绘图函数与图形美化- 数据分析：统计、插值、拟合等6. 模型仿真与优化- 模型建立与仿真- 参数优化与模型求解- 应用案例：物理、生物、经济等领域7. 综合应用与项目实践- 综合案例：结合实际问题的Matlab程序设计- 项目实践：团队协作完成指定项目，提高编程能力教学内容安排与进度：共安排16个课时，按照以下进度进行教学：1-4课时：Matlab基础知识与基本操作5-8课时：数据类型、运算符和程序控制结构9-12课时：函数与脚本、数据处理与分析13-16课时：模型仿真与优化、综合应用与项目实践教学内容与教材关联性：以上教学内容与教材章节相对应，涵盖了Matlab编程的主要知识点，旨在帮助学生系统地掌握Matlab编程技能。
第4章MATLAB程序设计
例 求一个值n，使n！最大但小于1050。 MATLAB程序如下： %初始化 r=1;k=1; 执行后结果如下： %做循环 The 41！Is while r<1e50 3.345253e+049 r=r*k; k=k+1; end k=k-1;r=r/k;k=k-1; % 输出结果 disp(['The',num2str(k),'!is',num2str(r)])
4.2.3 循环结构 for 循环：可完成指定次重复的循环，且在循环开 始之前，就知道代码重复的次数。
while 循环：while语句可完成不定次重复的循 环，与for语句不同，每次循环前要先判断条件 是否满足，再决定循环是否进行。
例 利用rand函数产生10个随机数，然后利用嵌套 for循环进行从小到大排序。
2 1.5 1 0.5 0 150 100 50
然后通过调用函数fun计算出z， 并利用mesh绘制出网格线。 编写的脚本文件如下： % 坐标系准备 x=[0:.02:2]; y=[-2:.02:0]; % 调用函数计算函数值 z=fun(x,y); % 绘图 figure mesh(z)
150 100 50
4.3.2 函数的调用 函数调用的一般格式如下： [输出实参表] = 函数名(输入实参表) 要注意的是，函数调用时各实参出现的顺 序、个数，应与函数定义时形参的顺序、 个数一致，否则会出错。函数调用时，先 将实参传递给相应的形参，从而实现参数 传递，然后再执行函数的功能。
例 对于函数
z ( x 1)2 ( y 1)2
fprintf函数中，format字符中常常会使用到一些特 殊字符来灵活的实现一些显示功能 format string 结果 %d %e %f %g 把值作为整数来处理 用科学记数法来显示数据 用于格式化浮点数，并显示这个 数 用科学记数格式，或浮点数格式， 根据那个表示的数位短，显示那 个 转到新的一行
第七章 MATLAB程序设计
1、输入程序： %对两个矩阵相 加 clc clear a=[1,2,3]; b=[4,5,6]; c=a+b
3、在命令窗口输入 myadd，回车 或打开文件，点击运 行按钮 给出结果： c=
5
7
9
18
打开已有的M文件，也有3种方法： (1) 菜单操作。从MATLAB主窗口的File菜单 中选择Open命令，则屏幕出现Open对话框， 在Open对话框中选中所需打开的M文件。在 文档窗口可以对打开的M文件进行编辑修改， 编辑完成后，将M文件存盘。 (2) 命令操作。在MATLAB命令窗口输入命令： edit 文件名，则打开指定的M文件。 (3) 命令按钮操作。单击MATLAB主窗口工具 栏上的Open File命令按钮，再从弹出的对话 框中选择所需打开的M文件。
直接访问work space 通过函数形参传递数据 中所有变量 直接选取MATLAB中 精心设计完成指定功能 执行的语句 重复操作 MATLAB功能扩展
7
1.3 M文件的创建
m文件的类型是普通的文本文件，我们可以使用 系统认可的文本文件编辑器来建立m文件。如 MATLAB下的m文件编辑器、windows的记事 本和word等。
s=pi*r*r; p=2*pi*r;
24
特定规则：
1.函数m文件第一行必须以单词function作为引导词，表示该M文件是一 个函数文件，必须遵循如下形式： function <因变量>=<函数名>(<自变量>)
2.函数名的命名规则与变量名相同，m文件的文件名必须是 <函数名> .m。 3.输入形参为函数的输入参数，输出形参为函数的输出参数。当输出形 参多于一个时，则应该用方括号括起来。 4.程序中的变量均为局部变量，不保存在工作空间中。其变量只在函数 运行期间有效。
matlab程序课程设计
matlab程序课程设计一、教学目标本课程的教学目标是使学生掌握MATLAB编程的基本知识和技能，能够运用MATLAB解决简单的数学和工程问题。
具体目标如下：1.理解MATLAB的基本概念，如变量、矩阵、数组等。
2.掌握MATLAB的基本语法，如运算符、函数、循环和条件语句等。
3.了解MATLAB的绘图功能，能够绘制基本的图形。
4.能够使用MATLAB进行简单的数学计算和数据分析。
5.能够编写MATLAB脚本程序，解决实际问题。
6.能够使用MATLAB的绘图功能，进行数据可视化。
情感态度价值观目标：1.培养学生的计算机编程思维，提高解决问题的能力。
2.培养学生团队合作的精神，提高沟通与协作能力。
3.培养学生对科学研究的兴趣，提高创新意识。
二、教学内容根据教学目标，本课程的教学内容主要包括以下几个方面：1.MATLAB概述：介绍MATLAB的发展历程、功能特点和应用领域。
2.MATLAB基本语法：变量、矩阵、数组、运算符、函数、循环和条件语句等。
3.MATLAB绘图功能：基本图形绘制、图形编辑和格式设置等。
4.MATLAB编程实践：数学计算、数据分析、实际问题解决等。
5.第1周：MATLAB概述和基本语法。
6.第2周：MATLAB绘图功能。
7.第3周：MATLAB编程实践（数学计算和数据分析）。
8.第4周：MATLAB编程实践（实际问题解决）。
三、教学方法为了达到教学目标，本课程将采用以下教学方法：1.讲授法：讲解MATLAB的基本概念、语法和绘图功能。
2.案例分析法：分析实际问题，引导学生运用MATLAB编程解决。
3.实验法：上机操作，让学生亲手编写MATLAB程序，巩固所学知识。
4.小组讨论法：分组完成项目任务，培养团队合作和沟通能力。
四、教学资源为了支持教学内容和教学方法的实施，本课程将采用以下教学资源：1.教材：《MATLAB程序设计》。
2.参考书：提供一些MATLAB编程的参考书籍，供学生自主学习。
MATLAB实验四：MATLAB程序设计 参考答案
function y=f(x) y=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14)
河南财经政法大学数学与信息科学学院
6
实验报告
f(1)*f(2)+f(3)
function f=f(a,b,c) f=g(a)*g(b)+g(c)^2; function g=g(x) g=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14);
a=rand(1,100); [b,c]=sort(a); index=c(1); i=1; while i<=7 if i==8 disp('You Lost!'); i=i+1; else d=input('please input the number you guess:'); if d>index disp('High'); elseif d<index disp('Low'); elseif d==index disp('You won!'); i=9;
syms i j k l n i=1;k=0; n=input('Input n:'); for j=1:n i=i*j; k=k+j; end fprintf('%d！=%d \n Sum(1,...,%d)=%d\n',n,i,n,k);
4、用 while-end 循环语句求不超过 1000 的偶数之和，并求显示出最大值。 s=0,n=0; while s<=1000 n=n+2; s=s+n; if s>1000 break end
实验五-MATLAB的程序设计
实验五MATLAB的程序设计一、实验目的1.熟悉MATLAB 的编程环境；2.熟悉和掌握脚本文件和函数文件建立与使用方法；3.掌握MATLAB程序设计的流程控制语句；4.掌握函数调用和参数传递；5.熟悉MATLAB 程序调试方法。
二、实验仪器1.计算机2.MATLAB 软件环境三、实验指导1. M 文本文件介绍一个M 文件通常包含五个部分，函数定义语句、H1 帮助行、帮助文本、函数体或者脚本文件语句和注释语句。
一个完整的M 文件的结构为：function f = fact(n) 函数定义语句% Compute a factorial value. H1 行% FACT(N) returns the factorial of N, 帮助文本% usually denoted by N!% Put simply, FACT(N) is PROD(1:N). 注释语句f = prod(1:n); 函数体函数定义语句只在函数文件中存在，定义函数名称、输入输出参数的数量和顺序。
函数定义语句的格式为：function y = fact(x)2. 函数流程控制顺序结构顺序结构是最简单的程序结构，系统在编译程序时，按照程序的物理位置顺序执行。
这种程序容易编制，但是结构单一，能够实现的功能有限。
判断语句（if…else…end）if…endif…else…endif…elseif…else…end1）if…end此时的程序结构如下：if 表达式执行代码块end这是最简单的判断语句，只有一个判断语句，其中的表达式为逻辑表达式，当表达式为真时，执行相应的语句，否则，直接跳到下一段语句。
2）if…else…end当程序有两个选择时，可以选择if…else…end 结构，此时程序结构为：if 表达式执行代码块1else执行代码块2end当判断表达式为真时，执行代码块1，否则执行代码块2。
3）if…elseif…else…end上面的两种形式中，分别包含一个选择和两个选择，当判断包含多个选择时，可以采用elseif 语句，结构为：if 表达式1执行代码块1elseif 表达式2执行代码块2elseif ..................else执行代码块end其中可以包含任意多个elseif 语句。
第四章 MATLAB程序设计
第二节 选择结构
2.1 if 语句
在MATLAB中，if 语句有三种格式
1、单分支if 语句
语句格式为： if 条件 语句组 end 当逻辑表达式的值为真时，执行该结 构中的执行语句，执行完之后继续向 下进行；若为假，则跳过结构中的内 容，向下执行。
实例分析
例4.4 当x是整数矩阵时，输出x的值。
程序3：
x=input('输入x='); y=cos(x+1)+sqrt(x*x+1);
u=cos(x+1);
v= sqrt(x*x+1); y=u+v end if x~=10
if x~=10
y=x*sqrt(x+sqrt(x));
end
y
y=x*sqrt(x+sqrt(x))
end
第二节 选择结构
if 条件1 语句组1 elseif 条件2 语句组2 …… elseif 条件m 语句组m else if-elseif的执行方式为：如果逻辑表达式1 的值为真，则执行语句1；如果为假，则判 语句组m+1 断逻辑表达式2，如果为真，则执行语句2， end 否则向下执行。
3、多分支if语句
实例分析
y
实例分析
例4.11 求
ye
0.5 x
sin( x
6
)
定积分。
a=0;b=3*pi;n=1000; h=(b-a)/n; x=a:h:b; f=exp(-0.5*x).*sin(x+pi/6);
for i=1:n
s(i)= (f(i)+f(i+1))*h/2; end s=sum(s)
第6章 MATLAB程序设计
如果用户调用toc时不指定输出参数t， 例如： >> tic >> toc elapsed_time = 4.0160 函数在命令窗口显示函数工作空间变量 elapsed_time 的值，但在MATLAB命令工作空间里不给输出参数t赋任 何值，也不创建变量t。
如果如果用户调用toc时指定输出参数t，例如： >> tic >> out=toc out = 2.8140 则以变量out的形式返回到命令窗口，并在 MATLAB命令工作空间里创建变量out。
function y = linspace(d1, d2, n) %LINSPACE Linearly spaced vector. if nargin == 2 n = 100; end n = double(n); y = [d1+(0:n-2)*(d2d1)/(floor(n)-1) d2];
function t=toc % TOC Read the stopwatch timer.
Global TICTOC If nargout<1 elapsed_time=etime(clock,TICTO C) else t=etime(clock,TICTOC); end
（8）全局变量可以为编程带来某些方便，但却破坏了函 数对变量的封装，所以在实际编程中，无论什么时候都应 尽量避免使用全局变量。如果一定要用全局变量，建议全 局变量名要长、采用大写字母，并有选择地以首次出现的 M文件的名字开头，使全局变量之间不必要的互作用减至 最小。 （9）MATLAB以搜寻脚本文件的同样方式搜寻函数文件。 例如，输入>>cow，MATLAB首先认为cow是一个变 量；如果它不是，那么MATLAB认为它是一个内置函数； 如果还不是，MATLAB检查当前cow.m的目录或文件夹； 如果仍然不是，MATLAB就检查cow.m在MATLAB搜 寻路径上的所有目录或文件夹。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验四第3章 MATLAB 程序设计Matlab 作为一种广泛应用于科学计算的工具软件，不仅具有强大的数值计算能力和丰富的绘图功能；可以人机交互式的命令行的方式工作；作为一种高级语言，同时也可以与 C 、FORTRAN 等高级语言一样进行程序设计.利用 Matlab 的程序控制功能，将相关 Matlab 命令编成程序存储在一个文件中（M 文件），然后在命令窗口中运行该文件，Matlab 就会自动依次执行文件中的命令，直到全部命令执行完毕.■ 例1 用 mesh 绘制半径为 3 的球命令行方式： 编程方式：新建一个M 文件 qiu.m 如下：保存后，在命令窗口输入 qiu ，即可执行该 M 文件.>> u=[0:pi/60:2*pi]; >> v=[0:pi/60:pi]; >> [U,V]=meshgrid(u,v); >> R=3;>> X=R*sin(V).*cos(U); >> Y=R*sin(V).*sin(U); >> Z=R*cos(V); >> mesh(X,Y,Z); >> axis equal;u=[0:pi/60:2*pi]; v=[0:pi/60:pi]; [U,V]=meshgrid(u,v); R=3;X=R*sin(V).*cos(U); Y=R*sin(V).*sin(U); Z=R*cos(V); mesh(X,Y,Z); axis equal;第一节M 文件一、M文件介绍●用Matlab 语言编写的程序称为M 文件●M 文件以.m 为扩展名●M 文件是由若干Matlab 命令组合在一起构成的，它可以完成某些操作，也可以实现某种算法●文件的命名规则与变量相同！文件名应尽量与程序要表达的意义相符合，以方便今后调用.（如例1）二、M文件的建立、打开与保存M文件是文本文件，可以用任何文本编辑器来建立和编辑，通常使用Matlab自带的M文件编辑器.①新建一个M文件●菜单操作( Fil e→New→M-File )●命令操作( edit M 文件名)●命令按钮( 快捷键)②打开已有的M 文件●菜单操作( File →Open )●命令操作( edit M 文件名)●命令按钮( 快捷键)●双击M文件(在当前目录窗口)③保存M 文件●菜单操作( File →Save )●命令按钮( 快捷键)三、M 文件分类（根据调用方式的不同）●Script：脚本文件/命令文件（可以直接运行的M文件）命令文件就是命令行的简单叠加，matlab会自动按顺序执行文件中的命令.这样就解决了用户在命令窗口中运行许多命令的麻烦，还可以避免用户做许多重复性的工作.（如例1）●Function：函数文件函数文件主要用以解决参数传递和函数调用的问题.(1) 第一行必须指定函数名、输入变量(参数)和输出变量(参数).输入参数是从MATLAB的工作空间复制到函数工作空间的变量.第一行举例如下：function [输出形参表] = name(输入形参表)(2) 一个函数可以有0个、一个或几个输入参数和返回值.当输出形参表的参数个数大于2时，[ ]不可缺省！(3) 建议函数名和文件名一样.调用时所用的变量并不需要与函数文件中定义的变量有相同的名字.■例2 比较下列两个程序，注意命令文件（以dd1命名）与函数文件（以dd2命名）的区别与联系.x=input(‘输入初值x=’);n=input(‘输入迭代次数=’); y(1)=x ;for k=1:nx=1/(x+1); y(k)=x;end function y=dd2(x,n) y(1)=x ;for k=1:nx=1/(x+1); y(k)=x;四、 两类文件的区别● 函数文件的第一行必须包含字 function ，命令文件没有这种要求.因此，没有这样第一行的M 文件是命令文件. ● 命令文件没有输入参数，也不返回输出参数，而函数文件可以带输入参数，也可以返回输出参数；● 命令文件对matlab 工作空间中的变量进行操作，文件中所有命令的执行结果也完全返回到工作空间中，而函数文件中定义的变量为局部变量，当函数文件执行完毕时，这些变量被清除.● 命令文件可以直接运行.在MATLAB 命令窗口输入命令文件的名字，就会顺序执行命令文件中的命令，而函数文件不能直接运行，而要以函数调用的方式运行. ● 函数文件可以建立函数，而命令文件不能。
■例3 建立函数23569(3)24()x xx x f x -+-++=解：程序如下：function y=f (x)56)/((3)(3)2)9/(4)(x x x x x x y *--*-++***+=；以文件名f.m 保存。
五、M文件的调用对已存在的M文件●命令文件在命令窗口直接输入该文件名即可；如例1>>qiu●函数文件调用的一般格式[输出实参数表]=函数名(输入实参数表)①函数调用时各实参数出现的顺序、个数应与函数定义时形式参数的顺序、个数一致，否则出错；②函数调用时，先将实参数传递给形式参数，从而实现参数的传递，然后再执行函数的功能.如例2，例3>>y=dd2(1,10) ， >>y=f(1)③当输出实参数表的参数个数大于2时，[ ]不可缺省！练习：将例1、例2和例3的程序在matlab中建立、保存、运行.第二节程序控制结构程序控制结构有三种：顺序结构、选择结构、循环结构.任何复杂的程序都可以由这三种基本结构构成.Matlab提供了实现控制结构的语句，利用这些语句可以编写解决实际问题的程序. 一、顺序结构顺序结构是指按照程序中语句的排列顺序从上到下依次执行，直到程序的最后一个语句如例1的qiu文件和例2的dd1文件.这是最简单的一种程序结构.一般涉及数据的输入、计算或处理、数据的输出等内容.1. 数据的输入通过input命令来接收从终端输入的内容，它也可以显示文本和提示.命令集14 输入命令input(out ,in)在屏幕上显示出字符串out的文本并等待终端的输入.如果变量in是's'，则输入的内容以字符串的形式进行保存，通常MATL AB在保存前要尽可能地求出表达式的值.如果使用格式控制符号如'\n'，字符串out可以是若干行.■例4 注意in的输入方式>> xm=input('What''s your name?','s')What's your name?syxxm =syx>> xm=input('What''s your name?')What's your name?syx??? Error using ==> inputUndefined function or variable 'syx'.What's your name?'syx'xm =syx2. 数据的输出可以通过简单地输入变量的名字来显示数字矩阵或者字符串向量的内容，结果将显示出变量的名字和内容.另一种显示变量的值就是使用命令disp.使用它只显示出变量的内容，这是有用的，特别是在字符串的应用中.命令集15显示命令disp(A)显示矩阵A的内容，如果A是字符串，则显示出它的文本.■例5 比较不同的输出方法及输出结果的表达方式>> A=magic(3);disp(A); >> disp(ones(3))8 1 6 1 1 13 5 7 1 1 14 9 2 1 1 1>> A >> A='Good afternoon';disp(A)A = Good afternoon8 1 6 >> disp('Hi,Tom')3 5 7 Hi,Tom4 9 2练习题1：编写M文件，对任意的a,b,c求方程20ax++=的根。
bx c二、循环结构循环结构是按照给定的条件，重复执行指定的语句。
Matlab 用于实现循环结构的语句有for — end 语句和while — end 语句。
1．for — end 语句 for 语句的格式为：for 循环变量=表达式1 : 表达式2 : 表达式3 循环体语句 end其中表达式1的值为循环变量的初值，表达式2的值为步长，表达式3的值为循环变量的终值。
步长为1时，表达式2可以省略。
■例6 1. 求1+2+…+100 ; 2. 求1+2+…+n . 1. 解： 2. 解■例7 1. 求初始值为1，任意等差数列的和？2. 求初始值任意，任意等差数列的和？1. 解：2. 解s=0;for i=1:100s=s+i; end sfunction s=mysum(n) s=0; for i=1:n s=s+i; end■ 例8 已知222111123y n =++++，求y 的值。
解：程序如下：运行并比较两程序。
循环语句应用时应注意： ● 循环语句可以嵌套使用；● 不能在 for — end 循环体内改变循环变量的值；● 为了提高代码的运行效率，应尽可能提高代码的向量化程度，避免 for — end 循环的使用，如例8；● 如果预先就知道循环的次数，则可以采用 for — end 循环；否则，如果预先无法确定循环的次数，则可以使用 while — end 循环。
n=input(‘n=’);y=0 ; for i=1:n y=y+1/i/i; end yn=input(‘n=’); i=1:n; y=sum(1./i.^2); y% d 为公差，n 为求和项数 function s=dc(n,d) s=0; for i=1:d:n s=s+i; end % c 为首项，d 为公差， % n 为求和项数function s=mysumdc(c,n,d) s=0; for i=c:d:n s=s+i; end练习题2：阅读下列程序，写出结果；并在计算机上运行验证：(1) bljh.mx=input(‘Input x please:’);y= input(‘Input y please:’);z=x; x=y; y=z;disp(x); disp(y)；（其中x=[ 12，3，3，4；43，6，3，-54 ]，y=1 : 8 ）(2) xlsum1.mn=input('n=');k=1:2:2*n-1;y=sum(1./k)(3) xlsum2.mn=input('n=');y=0;for k=1:2:2*n-1y=y+1/k;endy(4) Vander1.mt=[-1;2;1;3;5];n=length(t);A(: , n)=ones(n,1);for j= n -1: -1:1A(: , j)=t.*A(: , j+1);endA(5) Vander2.mt=[-1;2;1;3;5];n=max(size(t));A(: , n)=ones(n,1);for j= 1: nfor i=1: nA(i , j)=t(i)^(n-j);endA(6) ftx.mfunction s=ftx(a,b,n)h=(b-a)/n;x=a:h:b;f=exp(-0.5*x).*sin(x+pi/6);s=(2*sum(f)-f(1)-f(n+1))*h/2;练习题：3. 利用秦九韶算法编写M 文件，实现任意的多项式0()nk n k k P x a x ==∑的求值。