Matlab基础及其应用（陈姿羽）6 Matlab程序设计.pptVIP

局部变量和全局变量 如果在命令窗口执行命令与函数，效率很低，程序调试不方便，输入内容不易保存。 MATLAB的函数文件中除了自变量、因变量以外，在计算过程中如果用到了中间变量或其它的变量，这些都将是局部变量。局部变量只在这一函数中有效，在其它函数或程序中无效。例如，将三角形面积S视为三条边长a,b,c的函数，则可以编缉一个函数文件(文件名为：heron.m)来实现海伦公式的计算： function s=heron(a, b, c) p=(a+b+c)/2; s=sqrt(p*(p-a)*(p-b)*(p-c)); 将这一函数文件保存在工作目录下，在MATLAB环境中，键入如下指令 heron(3, 4, 5) 计算机执行后，将显示出计算结果：ans= 6。在这一函数文件中，三角形的半周长p就是一个局部变量。调用这一函数时，计算机用实际的数据 3、4、5代替形式变量a、b、c计算出变量p的值，然后再计算出s的值，最后将数据结果输出到屏幕上。但是只要计算出函数值，则自变量的实际数据和中间变量都将被自动删除，计算机内存中只保留了函数值变量（缺省变量名为ans）。 如果想保留中间变量p的值，就必须将它定义为全局变量。定义p为全局变量用命令 global p 注意：这一命令必须出现在两个地方，一是函数文件中第一次使用变量p之前；二是在命令窗口第一次调用这一函数文件之前。 M文件的编辑与运行 例6-1 分别建立命令文件和函数文件，将华氏温度f转换为摄氏温度c。 例6-1 分别建立命令文件和函数文件，将华氏温度f转换为摄氏温度c。 M文件颜色属性 建立新的M文件 打开已有的M文件 函数文件 函数调用 函数调用的一般格式是： [输出实参表]=函数名(输入实参表) 函数参数的可调性 在调用函数时，MATLAB用两个永久变量nargin和nargout分别记录调用该函数时的输入实参和输出实参的个数。 6.2 程序控制结构 循环结构 程序如下： y=0; n=100; for i=1:n y=y+1/(2*i-1); end disp(y); for语句更一般的格式为： for 循环变量=矩阵表达式 循环体语句 end 写出下列程序的执行结果。 s=0; a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a s=s+k; end disp(s); 2、while语句 例3-11 从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 3．break语句和continue语句 与循环结构相关，它们一般与if语句配合使用。 比较 选择（条件）结构：if，switch 单分支if语句： 双分支if语句： 例: 计算分段函数的值。 程序如下： x=input(请输入x的值:); if x=0 y= (x+sqrt(pi))/exp(2); else y=log(x+sqrt(1+x*x))/2; end disp(y) 多分支if语句：语句用于实现多分支选择结构。 计算函数y=|x|+|x-2|的值。 switch语句 try语句：试探性执行语句 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。 程序如下： A=[1,2,3;4,5,6]; B=[7,8,9;10,11,12]; try C=A*B; catch C=A.*B; end disp(C); disp(lasterr); %显示出错原因 6.3 数据的输入输出 1、数据的输入函数：input 2、数据的输出函数：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(x); disp([x1=,num2str(x(1)),,x2=,num2str(x(2))]); 4、M数据文件的存储/加载(save/load) save: save save filename save filename X Y Z load: