MATLAB程序文件

1、主讲教师主讲教师 穆志民穆志民 1、matlab基本操 作 2、高数中的 matlab 3、线代中的 matlab 4、概率方程中的 matlab 1.优化软件：优化软件： lingolingo 2.数据计算软件：数据计算软件：MatlabMatlab，MathematicaMathematica 3.文章的排版软件：文章的排版软件：wordword，latex(latex(数学数学) ) 4.数据统计处理软件：数据统计处理软件： excelexcel，spssspss，sassas 5.论文阅读软件：论文阅读软件：word(word(.doc.doc),Adobe ),Adobe Acrob

2、at (Acrobat (.pdf).pdf), ,超星阅览器超星阅览器rread(rread(电子书，电子书， ( (.pdg.pdg, ,.pdf.pdf),cajview(),cajview(.kdh.kdh, ,.pdf.pdf) ) 第第 一一 章章 u1节节.matlab简介及其基本内容简介及其基本内容 u2节节.高等数学中的高等数学中的matlab u3节节.线性代数中的线性代数中的matlab u4节节.概率统计中的概率统计中的matlab u5节节.方程方程(组组)中的中的matlab u6节节.数值处理与数值处理与matlab MATLAB是一个集数值计算、符号分析、图象显

3、示、是一个集数值计算、符号分析、图象显示、 文字处理于一体的大型集成化软件文字处理于一体的大型集成化软件. .它最初由美国的它最初由美国的 Cleve Moler博士所研制博士所研制. .其目的是为线性代数等课程中其目的是为线性代数等课程中 的矩阵运算提供一种方便可行的实验手段的矩阵运算提供一种方便可行的实验手段. .经过十几年的经过十几年的 市场竞争和发展，市场竞争和发展，MATLABMATLAB已发展成为在自动控制、生物已发展成为在自动控制、生物 医学工程、信号分析处理、语言处理、图像信号处理、医学工程、信号分析处理、语言处理、图像信号处理、 雷达工程、统计分析、计算机技术、金融界和数学界

4、等雷达工程、统计分析、计算机技术、金融界和数学界等 各各行各业中都有极其广泛应用的数学软件行各业中都有极其广泛应用的数学软件. .归纳归纳起起 来来,MATLAB,MATLAB具有以下几个特点：易学、适用范围广、功能具有以下几个特点：易学、适用范围广、功能 强、开放性强、网络资源丰富强、开放性强、网络资源丰富. .由于由于MATLABMATLAB的强大功能的强大功能, , 它能使使用者从繁重的计算工作中解脱出来它能使使用者从繁重的计算工作中解脱出来, ,把精力集中把精力集中 于研究、设计以及基本理论的理解上于研究、设计以及基本理论的理解上, ,所以所以,MATLAB,MATLAB已成已成 为在

5、校大学生、硕士生、博士生所热衷的基本数学软件为在校大学生、硕士生、博士生所热衷的基本数学软件. . 在此在此, ,我们把我们把MATLABMATLAB作为学习数学的工具介绍给读者作为学习数学的工具介绍给读者, ,希希 望能有利于读者今后的学习望能有利于读者今后的学习. . 由于由于MATLAB 拥有大量的常用数学函数、工拥有大量的常用数学函数、工 具箱，基本包括了现今数学、物理及工程应用领具箱，基本包括了现今数学、物理及工程应用领 域的常用函数，如能熟练掌握，无论在编程时间域的常用函数，如能熟练掌握，无论在编程时间 或程序量上，都远远低于或程序量上，都远远低于Fortran、C 等高级语等高级

6、语 言，尤其在处理矩阵和向量计算方面更具有无可言，尤其在处理矩阵和向量计算方面更具有无可 比拟的优势。比拟的优势。 我们在数学建模竞赛中，由于只有短短的三我们在数学建模竞赛中，由于只有短短的三 到四天，而论文的评判不仅注重计算的结果更注到四天，而论文的评判不仅注重计算的结果更注 重模型的创造性等很多方面，因此比赛中把大量重模型的创造性等很多方面，因此比赛中把大量 的时间花费在编写和调试程序上只会喧宾夺主，的时间花费在编写和调试程序上只会喧宾夺主， 是很不值得的，这时使用是很不值得的，这时使用MATLAB 可以很大程可以很大程 度上的方便计算、节省时间，使我们将精力更多度上的方便计算、节省时间，

7、使我们将精力更多 的放在模型的完善上，所以是较为理想的的放在模型的完善上，所以是较为理想的。 Matlab的安装的安装 Matlab有各种版本，早期有有各种版本，早期有Matlab 1.0 for 386的的DOS版本，版本， 后来逐步发展这里介绍的版本是后来逐步发展这里介绍的版本是 Matlab 6 6.x x for Windows因为它使用方便，界面美观，因为它使用方便，界面美观， 我们选择它作为主要讲解版本我们选择它作为主要讲解版本Matlab还有许多附加的部还有许多附加的部 分，最常见的部分称为分，最常见的部分称为Simulink，是一个用作系统仿真的，是一个用作系统仿真的 软件包，

8、它可以让您定义各种部件，定义各自对某种信号软件包，它可以让您定义各种部件，定义各自对某种信号 的反应方式及与其它部件的连接方式最后选择输入信号，的反应方式及与其它部件的连接方式最后选择输入信号， 系统会仿真运行整个模拟系统，并给出统计数系统会仿真运行整个模拟系统，并给出统计数 据据Simulink有时是作为有时是作为Matlab的一部分提供的，称为的一部分提供的，称为 Matlab with Simulink版本版本Matlab还有许多工具箱，它还有许多工具箱，它 们是根据各个特殊领域的需要，用们是根据各个特殊领域的需要，用Matlab自身的语言编写自身的语言编写 的程序集，使用起来非常方便您

9、可以视工作性质和需要的程序集，使用起来非常方便您可以视工作性质和需要 购买相应的工具箱常见的工具箱有：购买相应的工具箱常见的工具箱有： Signal Process信号处理信号处理System Identification 系统辨识系统辨识 Optimization优化优化Neural Network神经网络神经网络 Control System自动控制自动控制Spline样条样条 Symbolic Math符号代数符号代数Image Process图像处理图像处理 Nonlinear Control非线性非线性 控制控制 Statistics统计统计 常见的工具箱常见的工具箱 MATLABM

10、ATLAB的运行的运行 启动启动MATLABMATLAB 点击点击MATLAB图标图标, ,进入到进入到MATLAB 命令窗（命令窗（MATLAB Command WindowMATLAB Command Window）. .在命令窗内在命令窗内, ,可以输可以输 入命令、编程、进行计算入命令、编程、进行计算. . 学会使用学会使用helphelp命令命令 在命令窗内输入在命令窗内输入helphelp命令命令, ,再敲再敲 回车键回车键. .在屏幕上出现了在线帮助总览在屏幕上出现了在线帮助总览. .（注意：（注意：MATLABMATLAB 命令被输入后命令被输入后, ,必需敲回车键才能执行必需

11、敲回车键才能执行. .为行文方便为行文方便, ,以后以后 不再每次提醒不再每次提醒“敲回车键敲回车键”.”.学会使用学会使用helphelp命令命令, ,是学习是学习 MATLAB的有效方法的有效方法. .例如：要想知道例如：要想知道MATLAB中的基中的基 本数学函数有哪些本数学函数有哪些, ,可以在总览的第五行查到：可以在总览的第五行查到：MATLAB 中的中的“基本数学函数基本数学函数”用用elfunelfun(elementary function)(elementary function) 表示表示, ,于是于是, ,可进一步键入：可进一步键入：“help elfun”,“help

12、elfun”,屏幕上将出屏幕上将出 现现“基本数学函数基本数学函数”表表. .（注意：（注意：help elfunhelp elfun之间有空格之间有空格, , 以后不再每次提醒以后不再每次提醒. .）如果想了解）如果想了解sinsin函数怎样使用函数怎样使用, ,可进可进 一步键入一步键入help sin.help sin.在工具栏中点击在工具栏中点击helphelp按扭按扭, ,或点击？或点击？ 号按扭号按扭, ,与上面获取帮助信息的方法是等效的与上面获取帮助信息的方法是等效的. . 常用控制命令： clc %清屏 clear %清变量 save %保存变量 load %导入变量 一些演示

13、命令： logo surf(peaks(20)等 学会使用学会使用demodemo命令命令 在命令窗内输入在命令窗内输入demo命令命令,再敲回车再敲回车, 键屏幕上将出现演示窗口键屏幕上将出现演示窗口.（MATLAB Demo Window）一）一 共有三个窗口共有三个窗口,左边的窗口显示欲演示内容的大标题左边的窗口显示欲演示内容的大标题,选定选定 其中一项其中一项,右下方的小窗口显示欲演示的具体内容右下方的小窗口显示欲演示的具体内容,选中其选中其 中一栏中一栏,再点击再点击run按扭按扭,屏幕上将演示选定的演示程序屏幕上将演示选定的演示程序.右上右上 方的窗口显示关于大标题的一些说明方的窗

14、口显示关于大标题的一些说明.在命令窗内输入在命令窗内输入type (文件名文件名),将显示演示程序的将显示演示程序的M文件文件,仔细研究演示程序的仔细研究演示程序的M 文件文件,是学习是学习MATLAB的又一有效方法的又一有效方法. 进入演示窗还有另一方法：在工具栏中点击进入演示窗还有另一方法：在工具栏中点击Help栏栏,下拉式下拉式 菜单中点击菜单中点击examples and demos项项,即可进入演示窗口即可进入演示窗口. 退出退出 退出退出Matlab可键入可键入quit或或exit或选择相应的菜单中或选择相应的菜单中 止止Matlab运行会引起工作空间中变量的丢失运行会引起工作空间

15、中变量的丢失。 命命 令令含含 义义 help在线帮助在线帮助 helpwin在线帮助窗口在线帮助窗口 helpdesk在线帮助工作台在线帮助工作台 demo运行演示程序运行演示程序 ver版本信息版本信息 readme显示显示ReadmeReadme文件文件 who显示当前变量显示当前变量 whos显示当前变量的详细信息显示当前变量的详细信息 clear清空工作间的变量和函数清空工作间的变量和函数 pack整理工作间的内存整理工作间的内存 load把文件调入变量到工作间把文件调入变量到工作间 save把变量存入文件中把变量存入文件中 quit/exit退出退出MATLABMATLAB wha

16、t显示指定的显示指定的matlabmatlab文件文件 lookfor在在HELPHELP里搜索关键字里搜索关键字 which定位函数或文件定位函数或文件 path获取或设置搜索路径获取或设置搜索路径 echo命令回显命令回显 cd改变当前的工作目录改变当前的工作目录 pwd显示当前的工作目录显示当前的工作目录 dir显示目录内容显示目录内容 unix执行执行unixunix命令命令 dos执行执行dosdos命令命令 !执行操作系统命令执行操作系统命令 computer显示计算机类型显示计算机类型 MATLAB系统命令系统命令 1.1.2 1.1.2 变量、语句、矩阵与函数变量、语句、矩阵与

17、函数 变量变量，在在MATLABMATLAB中中, ,变量由字母、数和下划线变量由字母、数和下划线 组成组成. .第一个字符必须是字母第一个字符必须是字母. .一个变量最多由一个变量最多由3131 个字符组成个字符组成, ,并区分大小写并区分大小写. .下面是下面是MATLABMATLAB中表示中表示 特殊量的字符：特殊量的字符：pipi（圆周率）、（圆周率）、epseps（最小浮点（最小浮点 数）、数）、InfInf（正无穷大）、（正无穷大）、NaNNaN（表示（表示0/00/0或或inf-inf- infinf等不定值）、等不定值）、i,j,j（虚数单位）（虚数单位） 语句语句，Matla

18、bMatlab的表述语句、变量的类型说明由的表述语句、变量的类型说明由 MatlabMatlab系统解释和判断系统解释和判断MatlabMatlab语句通常形式为：语句通常形式为： 变量变量= =表达式表达式 或者使用其简单形式为：或者使用其简单形式为： 表达式表达式 表达式由操作符或其它特殊字符、函数和变量表达式由操作符或其它特殊字符、函数和变量 名组成表达式的结果为一个矩阵，显示在屏幕名组成表达式的结果为一个矩阵，显示在屏幕 上，同时保存在变量中以留用如果变量名和上，同时保存在变量中以留用如果变量名和 “=”=”省略，则具有省略，则具有ansans名名( (意思指回答意思指回答) ) 的变

19、量将自动建立例如：键入的变量将自动建立例如：键入1900/81 结果为：结果为：ans = 23.4568 需注意的问题有以下几点：需注意的问题有以下几点： * * 语句结束键入回车键，若语句的最后一个字符是分号，语句结束键入回车键，若语句的最后一个字符是分号， 即即“;”，则表明不输出当前命令的结果，则表明不输出当前命令的结果 * * 如果表达式很长，一行放不下，可以键入如果表达式很长，一行放不下，可以键入“ ” （三个点，但前面必须有个空格，目的是避免将形如（三个点，但前面必须有个空格，目的是避免将形如 “数数22 ”理解为理解为“数数2.”2.”与与“.”“.”的连接，从而导致的连接，从

20、而导致 错误），然后回车错误），然后回车 * * 变量和函数名由字母加数字组成，但最多不能超过变量和函数名由字母加数字组成，但最多不能超过6363 个字符，否则系统只承认前个字符，否则系统只承认前6363个字符个字符 * * Matlab变量字母区分大小写，如变量字母区分大小写，如A和和a不是同一个变不是同一个变 量，函数名一般使用小写字母，如量，函数名一般使用小写字母，如inv(A)不能写成不能写成 INV(A)，否则系统认为未定义函数，否则系统认为未定义函数 矩阵方式输入矩阵方式输入 A=2 ,3 5;3 5 7;4 5 8 %矩阵输入矩阵输入 （a为为2*3阶阶矩阵）矩阵） b=3;5;

21、7 %矩阵输入矩阵输入 det(A) %方阵行列式方阵行列式 inv(A) %方阵的逆方阵的逆 x=Ab %ax=b方程组的解方程组的解 y=inv(A)*b %与与x相同相同 disp(a,b,x) %显示矩阵显示矩阵 矩阵矩阵，MATLAB的变量及数组均是以向量或矩阵方式存储的变量及数组均是以向量或矩阵方式存储 MATLAB中矩阵的输入方法如下：中矩阵的输入方法如下： A=a11,a1n；am1,amn.逗号是同一行数逗号是同一行数 之间的分隔符（也可用空格代替）；分号是不同之间的分隔符（也可用空格代替）；分号是不同 行的分隔符行的分隔符. MATLAB中基本代数运算符中基本代数运算符 运

22、算运算 符号符号 举例举例 加法，加法，a+b + 5+3 减法，减法，a-b - 5-3 乘法，乘法，ab * 5*3 除法，除法，ab / or 48/4=448=12 乘幂，乘幂，ab 52=25 MATLAB中向量中向量(数组数组)、矩阵基本运算符、矩阵基本运算符 加法，A+B + 1 2+3 4 =1+3 2+4 减法，A-B - 1 2-3 4 = 1-3,2-4 乘法，A*B * 1,2*3;4=11 A.*B .* 1,2.*3,4=3,8 除法, A/B(AB) / 1,2;3 4/ 1,4;9,16=0.1,0.1; -0.6 0.4 A./B(A.B) . ./ 1,2;

23、3 4./1,2;3,2 =1,1;1,2 乘幂，方阵的幂 1,2;3 42=7 10;15,22 元素的幂 . 1,2;3 4.2=1,4;9,16 例：我们要计算例：我们要计算 可键入可键入：y=sin(pi/6)y=sin(pi/6)得得 y =0.5000y =0.5000 如果我们键入：如果我们键入： x=0,pi/6,pi/3,pi/2,2x=0,pi/6,pi/3,pi/2,2* *pi/3,5pi/3,5* *pi/6,pi;pi/6,pi; y=sin(x),y=sin(x),得得 y = 0 0.5000 0.8660 1.0000 0.8660 0.y = 0 0.500

24、0 0.8660 1.0000 0.8660 0. 5000 0.00005000 0.0000 这里这里, ,对于对于x x有更简洁的输入方法：有更简洁的输入方法：x=0 x=0：pi/6pi/6：pi,pi,此命令表此命令表 示示x x从从0 0开始开始, ,以以pi/6pi/6为步长变到为步长变到pipi为止为止. . 如果我们键入：如果我们键入： x=0,pi/6;pi/3,pi/2;x=0,pi/6;pi/3,pi/2; y=sin(x) y=sin(x) 得得 y = 0 0.5000y = 0 0.5000 0.8660 1.0000 0.8660 1.0000 其它的函数的用法

25、与此类似其它的函数的用法与此类似. . 6 sin MATLAB中 符号函数的 表示 函数函数，MATLAB提供了大量的函数提供了大量的函数.可以通过可以通过help查询查询.例如例如sqrt（开方）、（开方）、 log（常用对数）、（常用对数）、log10（以（以10为底的对数）为底的对数）、sin（正弦）等（正弦）等. MATLAB常用的基本数学函数常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度：纯量的绝对值或向量的长度 sqrt(x)：开平方：开平方 round(x)：四舍五入至最近整数：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数：无论正负，舍去小数至最近

26、整数 floor(x)：地板函数，即舍去正小数至最近整数：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数：天花板函数，即加入正小数至最近整数 rat(x)：将实数：将实数x化为分数表示化为分数表示 rats(x)：将实数：将实数x化为多项分数展开化为多项分数展开 sign(x)：符号函数：符号函数 (Signum function)。 当当x0，sign(x)=1。 rem(x,y)：求：求x除以除以y的余数的余数 gcd(x,y)：整数：整数x和和y的最大公因数的最大公因数 lcm(x,y)：整数：整数x和和y的最小公倍数的最小公倍数 exp(x)：自然

27、指数：自然指数 pow2(x)：2的指数的指数 log(x)：以：以e为底的对数，即自然对数或为底的对数，即自然对数或 log2(x)： 以以2为底的对数为底的对数 log10(x)：以：以10为底的对数为底的对数 MATLAB常用的三角函数常用的三角函数 sin(x)： 正弦函数正弦函数 cos(x)： 馀弦函数馀弦函数 tan(x)： 正切函数正切函数 asin(x)：反正弦函数：反正弦函数 acos(x)：反余弦函数反余弦函数 atan(x)：反正切函数：反正切函数 atan2(x,y)：四象限的反正切函数：四象限的反正切函数 sinh(x)： 超越正弦函数超越正弦函数 cosh(x)：

28、 超越馀弦函数超越馀弦函数 tanh(x)： 超越正切函数超越正切函数 asinh(x)： 反超越正弦函数反超越正弦函数 acosh(x)： 反超越馀弦函数反超越馀弦函数 atanh(x)： 反超越正切函数反超越正切函数 finverse(f) 反函数反函数 任何任何Matlab语句执行结果都可在屏幕上显示，同时赋给语句执行结果都可在屏幕上显示，同时赋给 指定的变量，没有指定变量时赋给指定的变量，没有指定变量时赋给ans数字显示格式可数字显示格式可 由由format命令来控制（命令来控制（Windows系统下的系统下的Matlab系统的数系统的数 字显示格式可以由字显示格式可以由Option菜

29、单中的菜单中的Numerical Format菜单菜单 改变）改变）format仅影响矩阵的显示，不影响矩阵的计算与仅影响矩阵的显示，不影响矩阵的计算与 存贮存贮(Matlab以双精度执行所有的运算以双精度执行所有的运算) 首先，如果矩阵元素是整数则矩阵显示就没有小数，如首先，如果矩阵元素是整数则矩阵显示就没有小数，如 x=-1 0 1，结果为：，结果为： x= -101 如果矩阵元素不是整数则输出形式有：（用命令：如果矩阵元素不是整数则输出形式有：（用命令： formatformat 格式进行切换）格式进行切换） 格式格式 中文解释中文解释 说明说明 format 短格式短格式 （缺省格式）

30、（缺省格式） Default. Same as SHORT format short 短格式短格式 （缺省格式）（缺省格式） Scaled fixed point format with 5 digits （只显示五位十进制数） format long 长格式长格式 Scaled fixed point format with 15 digits format short e 短格式短格式e方式方式 Floating point format with 5 digits format long e 长格式长格式e方式方式 Floating point format with 15 digits

31、 format short g 短格式短格式g g方式方式 Best of fixed or floating point format with 5 digi ts format long g 长格式长格式g g方式方式 Best of fixed or floating point format with 15 di gits format hex 16进制格式进制格式 Hexadecimal format format bank 银行格式银行格式 Fixed format for dollars and cents format rat 有理数格式有理数格式 Approximation

32、by ratio of small integers format compact 压缩格式压缩格式 Suppress extra line-feeds format loose 自由格式自由格式 Puts the extra line-feeds back in MATLAB中关系和逻辑运算中关系和逻辑运算 1 1关系运算符关系运算符 小于小于, , 小于等于小于等于, ,大于大于, ,大于等于大于等于, , 等于，等于， 不等于不等于 运算法则：运算法则：如果两个比较量a、b是标量,那么,当 a、b之间的关系成立时输出值为1；否则输出值为 0. 如果两个比较量a、b是相同维数的数组,那么就

33、 按标量的运算法则,对a、b的对应元素进行运算,最 后的输出结果为一个与a（或b）同维的01数组. 如果a是标量,b是数组,那么按标量的运算法则将 a与b的每个元素逐一比较,最后的输出结果为一个 与b同维的01数组. 在算术运算、关系运算中,算术运算优先. 2 2逻辑运算符逻辑运算符 a=y b=t 例：例：交换两个变量的值交换两个变量的值 myswap.m l 有两个输入参数和两个输出参数有两个输入参数和两个输出参数 clear clc a=3,b=4, X=a; a=b, b=X %交换交换a,b值值 函数文件函数文件 命令文件命令文件 命令文件的创建和运行命令文件的创建和运行： 在在MA

34、TLABMATLAB的命令窗口点击的命令窗口点击“新建新建”工具栏或在工具栏或在“file”“file” 下拉菜单中选下拉菜单中选“New”“New”中的中的 “M-file”“M-file”项项, ,进入编辑进入编辑/ /调试器调试器. . 例例 计算所有小于计算所有小于1000的的Fibonnaci数数. % 计算小于计算小于10001000的的Fibonnaci（斐波那契数列斐波那契数列）数数 clearclear clc clc f=1,f=1,2 2; i=1;1; whilewhile f( f(i)+f()+f(i+1)1000+1)1000 f(i+2)=f( f(i+2)=f

35、(i)+f(i+1);)+f(i+1); i= =i+1;+1; endend f, f,i 在在“file”“file”下拉菜单中选下拉菜单中选“Save”“Save”项项, ,依提示输入文件依提示输入文件 名名“fibno”“fibno”至此至此, ,完成了命令文件完成了命令文件fibno.mfibno.m的创建的创建. . function 输出形参列表输出形参列表=函数名函数名(输入形参列表输入形参列表) % 注释说明部分注释说明部分(可选可选) 函数体语句函数体语句(必须必须) l 第一行为第一行为引导行引导行，表示该，表示该 M文件是函数文件文件是函数文件 l 函数名的命名规则与变

36、量名相同函数名的命名规则与变量名相同 ( 必须以字母开头必须以字母开头 ) l 当当输出形参输出形参多于一个时，用多于一个时，用方括号方括号括起来括起来 l 以以百分号百分号开始的语句为开始的语句为注释语句注释语句 函数文件创建函数文件创建 q 函数文件由函数文件由 function 语句引导语句引导 l 函数函数文件名必须与函数名一致文件名必须与函数名一致 l 函数必须是一个函数必须是一个单独的单独的 M文件文件 functionf=ffibno(n) f=1,1;f=1,1; i=1;1; whilewhile f(f(i)+f()+f(i+1)n+1)n f(i+2)=f(f(i+2)=

37、f(i)+f(i+1);)+f(i+1); i= =i+1;+1; endend f 在在“file”“file”下拉菜单中选下拉菜单中选“Save”“Save”项项, ,依提示输入文依提示输入文 件名件名“ffibno”fibno”至此至此, ,完成了函数文件完成了函数文件ffibno.mfibno.m的创建的创建. . 执行执行ffibno.fibno. 在在MATLABMATLAB窗口中输入窗口中输入ffibnofibno(1000)并敲回车键即可并敲回车键即可. . 例如：一个自编的函数文件例如：一个自编的函数文件mean.m，用于求向量的（或矩阵按列的）平均值，用于求向量的（或矩阵按

38、列的）平均值 function y=mean(x) % MEAN Average or mean value，For Vectors, % MEAN (x) returns the mean value % For matrix MEAN (x) is a row vector % containing the mean value of each colu umn m,n=size(x); if m=1 m=n; end y=sum(x)/m; 磁盘文件中定义的新函数称为磁盘文件中定义的新函数称为mean函数，它与函数，它与Matlab函数一样使用，例如函数一样使用，例如z为从为从1到到99

39、 的实数向量：的实数向量： z=1:99; 计算均值：计算均值：mean(z) ans= 50 mean.m程序的说明：程序的说明： （1 1）第一行的内容：函数名，输入变量，输出变量，没有这行这个文件就是命令文件，）第一行的内容：函数名，输入变量，输出变量，没有这行这个文件就是命令文件， 而不是函数文件而不是函数文件 （2 2）：表明右边的行是说明性的内容注释前一小部分行来确定）：表明右边的行是说明性的内容注释前一小部分行来确定M文件的注释，文件的注释， 并在键入并在键入help mean后显示出来显示内容为连续的若干个右边的文字后显示出来显示内容为连续的若干个右边的文字 （3 3）变量）变

40、量m,n和和y是是mean的局部变量，在的局部变量，在mean运行结束后，它们将不在工作空间运行结束后，它们将不在工作空间z z中中 存在如果在调用函数之前有同名变量，先前存在的变量及其当前值将不会改变存在如果在调用函数之前有同名变量，先前存在的变量及其当前值将不会改变 MATLAB编程编程语句语句 1 1 控制语句控制语句 MATLABMATLAB也有控制流语句也有控制流语句, ,用于控制程序的流程用于控制程序的流程. .主要有主要有forfor 循环、循环、whilewhile循环、循环、ifif和和breakbreak三种控制语句三种控制语句. .虽然语句很少虽然语句很少, , 但功能很

41、强但功能很强. . for for循环语句循环语句 for循环语句的一般表达形式为：循环语句的一般表达形式为： i= =表达式表达式 可执行语句可执行语句1 1 可执行语句n 作用：重复执行命令集作用：重复执行命令集commands 格式： for i=n1:(step):n2 commands; end 例：求例：求S=1+2+3+50,S=1+2+3+50,可编程如下可编程如下： s=0;s=0; forfor k=1:50 k=1:50 s=s+k; s=s+k; endend whilewhile循环语句用来控制一个或一组语句在某逻辑条件下循环语句用来控制一个或一组语句在某逻辑条件下

42、重复预先确定或不确定的次数重复预先确定或不确定的次数. . while循环语句的一般表达形式为：循环语句的一般表达形式为： 表达式表达式 循环体语句循环体语句 例：例：求求S=1+2+3+50, ,可编程如下：可编程如下： S=0;k k=0; whilek51 S=S+k;k=k+1; % %当条件当条件k51kep k=k+1; x=exp(-x0); dx=abs(x-x0)/(1+abs(x); x0=x; end if if和和breakbreak语句语句 MATLABMATLAB中中ifif和和breakbreak语句的作用与使用方式同其它编程语句的作用与使用方式同其它编程 语言一

43、样语言一样, ,用来将控制流程用来将控制流程 进行分流与中断退出进行分流与中断退出. . 例：可以把上面的解方程的例子中的循环语句改写成例：可以把上面的解方程的例子中的循环语句改写成: : x0=0; while(1while(1) ) k=k+1; k=k+1; x=exp(-x0); x=exp(-x0); dx=abs(x-x0)/(1+abs(x); dx=abs(x-x0)/(1+abs(x); ifif dx=eps dx=eps breakbreak; ;endend x0=x; x0=x; end 程序中程序中 while(1)while(1)说明循环条件总是真说明循环条件总是

44、真, ,直到满足直到满足 dx=epsdx=eps条件跳出循环体条件跳出循环体. . if else end分支结构分支结构 分支结构有三种形式：分支结构有三种形式： if 表达式表达式 执行语句执行语句 end 如果表达式的值非0,则执行下面的语句.否则执行end 后面的语句. if 表达式 执行语句1 else 执行语句2 end if 表达式表达式1 1 执行语句执行语句1 1 elseif elseif 表达式表达式2 2 执行语句执行语句2 2 elseif elseif 表达式表达式3 3 执行语句执行语句3 3 else else （此句可以省略）（此句可以省略） 执行语句执行语句n n end 例：函数例：函数 10, x11-3x 10,x1 12 1, x )(x x xf ifif x1 x=1x-1; elseelse y=3y=3* *x-11;x-1