2011第4讲_Matlab程序设计解析.pptVIP

第四讲 Matlab程序设计 主讲人：祁荣宾 第四讲 主要内容 M文件 程序控制结构 函数文件 程序调试 1、M文件 Matlab命令执行方式 1、M文件 1.1 M文件的分类 1、M文件 （1）命令文件与函数文件的区别 1、M文件 例：实现变量a、b的值互换 1、M文件 1.2 M文件的建立 1、M文件 1.3 打开已有的M文件 2、程序控制结构 2.1 顺序结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 2、程序控制结构 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 3、函数文件 4、程序调试 4、程序调试 4、程序调试 4、程序调试 4、程序调试 4、程序调试 4、程序调试 如果一个循环结构的循环体又包括一个循环结构，就称为循环的嵌套，或称为多重循环结构。 （4） 循环的嵌套 例 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m); end end 结果为6，28，496 function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 （1）函数文件的基本结构 以function开头的一行为引导行，表示该M文件是一个函数文件。 函数名的命名规则与变量名相同。 输入形参为函数的输入参数，输出形参为函数的输出参数。 当输出形参多于一个时，则应该用方括号括起来。 函数文件名通常由函数名加上扩展名.m组成，函数名要与函数文件名一致； function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 （1）函数文件的基本结构 紧随函数文件引导行之后以%开头的第一注释行，一般包括大写的函数文件名和函数功能简要说明，供lookfor关键词查询和help在线帮助时使用； 第一注释行及之后连续的注释行，通常包括函数输入、输出参数的含义及调用格式说明等信息，构成全部在线帮助文本； 与在线帮助文本相隔一空行的注释行，包括函数文件编写和修改的信息，如作者、修改日期、版本等内容，用于软件档案管理。 function 输出形参表=函数名(输入形参表) 注释说明部分 函数体语句 （1）函数文件的基本结构 如果在函数中插入return语句，则执行到该语句就结束函数的运行，程序流程转到调用该函数的位置； 在函数文件中，也可以不使用return语句，这时被调用的函数执行完成后自动返回。 例 编写函数文件，求半径为r的圆的面积和周长。 function [s,p]=fcircle(r) %FCIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 %2004年7月30日编 s=pi*r*r; p=2*pi*r; 执行下面语句，观察结果 [s,p]=fcircle(10) help fcircle lookfor fcircle （2）函数调用 [输出实参表]=函数名(输入实参表) 要注意的是，函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个数一致，否则会出错。 函数调用时，先将实参传递给相应的形参，从而实现参数传递，然后再执行函数的功能。 例利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之间的转换。 函数文件tran.m： function [rho,theta]=tran(x,y) rho=sqrt(x*x+y*y); theta=atan(y/x); 调用tran.m的命令文件main1.m： x=input(Please input x=:); y=input(Please input y=:); [rho,the]=tran(x,y); rho the 例 利用函数的递归调用，求