Matlab函数进阶：使用匿名函数和内嵌函数处理多变量传递问题.docxVIP

Matlab函数进阶：使用匿名函数(Anonymous Function)和内嵌函数(Nested Function)处理多变量传递问题(Matlab 7.0以上)问题：有一个多变量函数f(abcx),现需要分别在a=a0b=b0c=c0和a=a1b=b1c=c1的条件下对f(abcx)进行某一操作。此类问题常在数值积分时出现。解决方案：1．使用全局变量可在主调函数和被调函数中分别将a,b,c声明为全局变量(globala b c)，这时f通过全局变量传递abc,定义时可以只有一个参数x。2．使用anonymous function3．使用 nested function下面举例说明anonymous function和nested function的使用。例：对任意二次多项式 ax2+bx+c 进行数值积分（此处显然可以解析得到，此例使用anonymousfunction做演示）解：（1）使用匿名函数编写文件intpoly2.m如下functiony_int=intpoly2(a,b,c)y_int=quad(@(x)(poly2(a,b,c,x)), 0，1); %此处利用matlab内部函%数quad(fun, x0,xt)进行积分，%被积函数fun我们使用匿名函%数quot;@(x)(poly2(a,b,c,x))quot;以便%将自变量限制为x. function y=poly2(ai,bi,ci,x)%此处定义一个多变量的子函数 poly2( a,b,c,x)y=ai.*x.^2+bi.*x+ci； %此处使用矩阵元素运算(.* 和.^)等便于%被矩阵化的内部函数调用，且提高程序效率执行：保存该文件并将matlab切换至该文件目录下，命令行输入intpoly2(1,2,3)，便给出积分结果ans=4.33（2）使用nested function编写函数保存为intnest.m，内容如下functiony_int=intnest(a,b,c)y_int=quad(@poly2, 0，1); %此处利用matlab内部函%数quad(fun, x0,xt)进行积分，%被积函数fun我们使用内嵌函%数poly2(x)的句柄@poly2 function y=poly2(x)%此处定义一个内嵌函数 poly2( a,b,c,x) y=a.*x.^2+b.*x+c； %直接调用母函数中的变量a,b,cend% 结束内嵌函数poly2end% 结束母函数intpoly2保存后执行，同样效果。可见nested function只是将主调函数和被调函数封装到了一起以共享主调函数的变量。注意：一般函数不需要end，当出现nested function时，主调函数和被调函数都必须有end 来结束！附：M脚本文件和M函数文件、函数句柄m文件分两种：函数和脚本。1、matlab script file:脚本文件没有输入输出，对工作空间（workspace）中的变量进行操作。任何可执行的matlab命令都可以写入脚本文件。先来看一个简单的例子：例1: 产生一个20元素的一维随机数组并画出。解: 新建一个m文件命名为randplt.m(或者任何你喜欢的名字), 在其中加入如下两行代码:data=randn(1,20); plot(data);保存后选择debug/run菜单或按F5键运行,观察workspace中多出了一个变量data。脚本文件执行时, 就如同将文件中的每一条命令依次输入到matlab命令行中一样, 顺次执行。你可以尝试在command history窗口中按住ctrl键选择几行执行过的命令并单击鼠标右键在快捷菜单中选择creat M file来快速创建一个m文件。2、matlab function函数文件可以接受输入和给出输出,当然也可以没有,就像c语言的函数。matlab函数最大的书写特点(同C比起来)是它的输出变量定义在函数名前面。看一例:例2: matlab内部函数std给出数组的标准差std=sqrt(sum(X)/N),编写函数stderr求实验误差err=sqrt(sum(X)/(N*(N-1))), 代码如下：function err=stderr(arr)%实验误差（标准差估计）%caculate along each column of arrerr=std(arr)./sqrt(size(arr,1)-1);文件存为stderr.m, 运行时输入a=[1;2;3];建立一个列数组a, 然后输入stderr(a)则给出a的方差。(若输入行数组则发散,因N-1=0)此例中输入变量为arr, 函数名为stderr, 输出变量err, 此三者皆写在关键字fun