Matlab实训6-程序设计.ppt

《Matlab实训6-程序设计.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、本章详细讲解利用M语言进行编程的方法讲述的主要内容：脚本文件流程控制函数文件子函数M文件的调试实训六MATLAB编程基础（二）MATLAB编程基础（续）6.1概述6.2脚本文件6.3流程控制6.4函数文件6.5M文件调试6.6M文件性能分析6.4函数文件M语言函数文件能够接受用户的输入参数，进行计算，并将计算结果作为函数的返回值返回给调用者在MATLAB中具有不同类型的函数文件内建函数文件系统M函数文件系统MEX函数文件用户自定义MEX函数文件用户自定义的M函数文件函数文件（续）MATLAB的内建函数文件MATLAB自定义的函数文件称内建函数文件。调用内建函数的方法：使用函数名并给出相

2、应的输入、输出参数即可。例如：sin.m函数调用格式：y=sin(x)例：>>x=0:2*pi/180:2*pi;>>y=sin(2*x);>>plot(x,y)用户自定义的M函数文件是由用户自己利用M语言编写的文件系统M函数文件是由MATLAB提供的M语言函数文件函数文件（续）本节介绍6.4.1基本结构6.4.2输入输出参数6.4.3子函数6.4.4局部变量和全局变量6.4.5函数执行规则函数文件（续）M语言函数文件具有下面的不同部分函数定义行在线帮助注释行M语言代码6.4.1基本结构函数m文件的格式：function返回变量=函数名（输入变量）注释说明语句段程序语句段特定规则：1）

3、函数m文件第一行必须以单词function作为引导词，必须遵循如下形式：function<因变量>=<函数名>(<自变量>)2）程序中的变量均为局部变量，不保存在工作空间中，其变量只在函数运行期间有效。基本结构（续）001functiony=average(x)002%AVERAGE求向量元素的均值003%语法：004%Y=average(X)005%其中，X是向量，Y为计算得到向量元素的均值006%若输入参数为非向量则出错007008%代码行009[m,n]=size(x);010%判断输入参数是否为向量011if(~((m==1)

4、(n==1))

5、(m==1&n==1))012%若

6、输入参数不是向量，则出错013error('Inputmustbeavector')014end015%计算向量元素的均值016y=sum(x)/length(x);例函数文件示例函数定义行001functiony=average(x)包括关键字function函数输出参数y函数名称average函数输入参数x函数名称定义要求必须以字符开头，后面可以用字符、数字和下划线的组合构成函数名称MATLAB对函数名称的长度有限定函数的M文件名称最好和函数名称保持一致，若不一致，则调用函数时需要使用文件名称而非函数名称。基本结构（续）在线帮助M函数文件的在线帮助为紧随在函数定义行的注释行002%

7、AVERAGE求向量元素的均值003%语法：004%Y=average(X)005%其中，X是向量，Y为计算得到向量元素的均值006%若输入参数为非向量则出错基本结构（续）若在MATLAB命令行窗口中键入指令helpaverage则>>helpaverageAVERAGE求向量元素的均值语法：Y=average(X)其中，X是向量，Y为计算得到向量元素的均值若输入参数为非向量则出错基本结构（续）注释行008%代码行010%判断输入参数是否为向量012%若输入参数不是向量，则出错015%计算向量元素的均值注释行不会显示在在线帮助中，主要原因是这些注释行没有紧随在H1帮助行的后面基本结构（

8、续）M语言代码008%代码行009[m,n]=size(x)010%判断输入参数是否为向量011if(~((m==1)

9、(n==1))

10、(m==1&n==1))012%若输入参数不是向量，则出错013error('Inputmustbeavector')014end015%计算向量元素的均值016y=sum(x)/length(x)代码行需要完成具体的算法，实现用户的具体功能基本结构（续）MATLAB在定义输入输出参数时不需要指出变量的类型，而是将参数默认为双精度型MATLAB在定义参数时，没有确定输入参数的维数或者尺寸M语言的函数文件不仅可以有一个输入参数和一个返回值，还可以为M语言

11、函数文件定义多个输入参数和多个输出参数6.4.2输入输出参数例多个输入输出参数的M函数function[avg,stdev,r]=ourstats(x,tol)%OURSTATS多输入输出参数示例%该函数计算处理矩阵，得到相应的均值、标准差和矩阵的秩[m,n]=size(x);ifm==1m=n;end%Averageavg=sum(x)/m;%Sandaddeviationstdev=sqrt(sum(x.^2)/m-avg.^2);%Ran