MATLAB程序设计及应用实例.ppt

5.1MATLAB程序设计
5.2应用实例5.1MATLAB程序设计5.1.1M文件1.建立新的M文件
(1)菜单操作：MATLAB命令窗口的File菜单→New菜单项→M-file命令。
(2)命令操作：在MATLAB命令窗口输入命令edit。
(3)命令按钮操作：单击MATLAB命令窗口工具栏上的“新建”命令按钮。
2.编辑已有的M文件
3.M文件的分类P1164.1.2内容
分类：脚本文件(ScriptFile)和函数文件(FunctionFile)。
主要区别。
M文件的基本结构1脚本文件
没有输入和输出
由一系列指令组成
可在命令窗口直接运行
产生的所有变量存储在workspace中。
例exam5_1.m。(1)函数文件的基本结构
函数文件由function语句引导
其基本结构为：
function输出形参表=函数名(输入形参表)
%注释说明部分
函数体语句函数文件的格式说明
(1)关于函数文件名:函数文件名与函数名也可以不相同（一般设为相同！）。当两者不同时，MATLAB将忽略函数名而确认函数文件名，因此调用时使用函数文件名。
(2)关于注释说明部分。注释说明包括三部分内容：①紧随函数文件引导行之后以%开头的第一注释行。②第一注释行及之后连续的注释行。③与在线帮助文本相隔一空行的注释行。函数调用的一般格式：
[输出实参表]=函数名(输入实参表)

注意：等号左边是方括号，右边输入参数用括号例5.1编写函数文件求半径为r的圆的面积和周长

function[s,p]=fcircle(r)
%FCIRCLEcalculatetheareaandperimeterofacircleofradiir
%r圆半径
%s圆面积
%p圆周长
s=pi*r*r;
p=2*pi*r;
以上函数文件以文件名fcircle.m存入c:\matlab7\work\mas下，然后在MATLAB命令窗口调用该函数：
[s,p]=fcircle(10)函数的递归调用
例：利用递归调用函数文件factor.m：
functionf=factor(n)
ifn<=1
f=1;
else
f=factor(n-1)*n;
end
在命令文件main5_1.m中调用函数文件factor.m求s=1!+2!+3!+4!+5!。函数所传递参数的可调性
在调用函数时，MATLAB用两个永久变量nargin和nargout分别记录调用该函数时的输入实参和输出实参的个数。只要在函数文件中包含这两个变量，就可以准确地知道该函数文件被调用时的输入输出参数个数，从而决定函数如何进行处理
例：nargin用法示例函数文件examp.m：
functionfout=charray(a,b,c)
ifnargin==1
fout=a;
elseifnargin==2
fout=a+b;
elseifnargin==3
fout=(a*b*c)/2;
end
命令文件main5_2.m：
x=[1:3];y=[1;2;3];
examp(x)
examp(x,y')
examp(x,y,3)全局变量用命令global定义。函数文件的内部变量是局部的，与其他函数文件及MATLAB工作空间相互隔离
例全局变量应用示例。
先建立函数文件wadd.m，该函数将输入的参数加权相加。
functionf=wadd(x,y)
globalALPHABETA
f=ALPHA*x+BETA*y;
在命令窗口中输入：
globalALPHABETA
ALPHA=1;
BETA=2;
s=wadd(1,2)5.1.3数据的输入输出指令input函数disp函数例求一元二次方程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)];
disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);

其调用格式为：
pause(延迟秒数)
如果省略延迟时间，直接使用pause，则将暂停程序，直到用户按任一键后程序继续执行。
若要强行中止程序的运行可使用Ctrl+C命令。5.1.4选择结构-if语句语句格式为：
if条件
语句组
end

当条件成立时，则执行语句组，执行完之后继续执行if语句的后继语句，若条件不成立，则直接执行if语句的后继语句
双分支if语句格式为：
if条件
语句组1
else
语句组2
end
当条件成立时，执行语句组1，否则执行语句组2，语句组1或语句组2执行后，再执行if语句的后继语句。例计算分段函数值
程序如下：
x=input('请输入x的值:'