上机练习五 函数文件编写与程序调试.docVIP

上机练习五 函数文件编写与程序调试一、上机目的1、掌握函数文件的创建、编辑、运行和调用方法；2、掌握全局变量的定义和使用方法；3、通过练习进一步理解函数文件的编程方法。二、上机仪器设备1、计算机；2、MATLAB软件环境。三、上机内容和步骤1、函数文件的结构、创建与运行打开Matlab M编辑器，输入以下代码：function y=average(x)% 函数average(x)用以计算向量元素的平均值。% 输入参数x 为输入向量，输出参数y为计算的平均值。% 非向量输入将导致错误。% 代码行。m,n=size(x); % 判断输入量的大小if(m=1)|(n=1)| (m=1)& (n=1) % 判断输入是否为向量error(必须输入向量。)endy=sum(x)/length(x); %计算向量x 所有元素的平均值将以上代码以文件名average.m保存，尝试以下操作：1）在命令窗口输入： z=1:99; average(z)观察运行结果，体会函数文件的调用方法。2）在命令窗口输入： lookfor average观察运行结果，体会函数文件H1行的作用。3) 在命令窗口输入： help average观察运行结果，体会函数文件注释行的作用。4) 将文件名改成average_1.m，分别在命令窗口输入以下指令： average(1 2 3) average_1(1 2 3)观察运行结果，总结函数文件名与函数名不同时的函数调用方法。2、子函数的使用打开Matlab M编辑器，输入以下代码：function avg,med = newstats(u) %主函数% NEWSTATS 计算均值和中间值n = length(u);avg = mean(u,n); %调用子函数med = median(u,n); %调用子函数 function a = mean(v,n) %子函数 %计算平均值a = sum(v)/n; function m = median(v,n) % 子函数 %计算中间值w = sort(v);if rem(n,2) = = 1 m = w(n+1)/2);else m = (w(n/2)+w(n/2+1)/2;end 将以上代码以文件名newstats.m保存，尝试以下操作。1）在命令窗口输入： x = 1:11; mean,mid = newstats(x)观察运行结果，体会含有子函数的函数文件的运行过程。2）在命令窗口分别输入以下指令： help newstats help newstatsmedian观察运行结果，总结含有子函数的函数文件H1行的作用。3、函数参数的可调性建立以下函数文件，并按照以下方式运行之，观察运行结果，体会nargin和nargout变量的使用方法。function apple = test(a,b,c)if nargin = =1apple n = a;elseif nargin = =2apple n = a+b;elseif nargin = =3apple n = a+b+c;end test(4) test(4,5) test(4,5,6)4、全局变量与局部变量1）建立以下函数文件（local.m）：function local%LOCAL 查看局部变量的例子x = rand(2,2);y = zeros(2,2);z = 函数中的变量;xu = x,y,z;disp(z)whos在命令窗口运行该文件 local在当前工作空间观察函数文件变量是否存在，并解释运行结果。 whos2）将上述函数文件改成（global_1.m）：function global_1%GLOBAL 查看全局变量的例子Global XX= rand(2,2);y = zeros(2,2);z = 函数中的变量;Xu = X,y,z;disp(z)whos在命令窗口运行该文件 global_15、函数的递归调用先用函数的递归调用定义一个函数文件求，然后调用该函数文件求四、上机要求1、完成以上上机任务；2、总结函数文件