【转载】MATLAB函数句柄和匿名函数

转自：http://bbs.seu.edu.cn/pc/pccon.php?id=950&nid=14857

下面代码创建一个内联函数a_humps：
>> a_humps = inline('1./((x-.3).^2+0.01)+1./((x-.9).^2+0.04)-6','x')
a_humps =
     Inline function:
     a_humps(x) = 1./((x-.3).^2+0.01)+1./((x-.9).^2+0.04)-6
上例中，函数inline从一个字符串创建一个函数，并以x为输入变量。要在一个函数中调用内联函数，只要将该内联函数的名字作为输入参数传递给函数即 可。例如，要将quad(Fun,low,high)中的Fun换为上面的内联韩式a_humps，只要按下面的方式调用即 可：quad(a_humps,low,high)。
要验证一个由字符串表示的函数或一个内联函数，可以使用feval函数。下面的代码验证了正弦函数和前面创建的a_humps函数
>> y = feval_r('sin',pi*(0:4)/4)
y =
         0    0.7071    1.0000    0.7071    0.0000
>> z = feval_r(a_humps,[1- 0 1])
z=
        -5.1378      5.1765     16
除了字符串函数和内联函数外，还由一种函数类型：匿名函数，并用函数句柄表示它。在应用中并不鼓励用户使用前两种方法，而是要尽量使用匿名函数句柄来引用函数。下面代码给出了一个匿名函数的例子：
b_humps = @ (x) 1./((x-.3).^2+0.01)+1./((x-.9).^2+0.04)-6;
其中，@符号意味这等号左边是一个函数句柄。@后面的（x）定义了函数的输入参数，最后一部分是函数表达式。我们同样可以利用feval函数来验证匿名函数，例如，可以使用下面代码验证b_humps：
>> z = feval_r(b_humps,[-1 0 1])
z =
   -5.1378    5.1765   16.0000
其实，用户根本没有必要利用feval函数来验证匿名函数，因为匿名函数可以使用自己的函数句柄直接进行验证，例如，上面的例子可以简写为：
>> z = b_humps([-1 0 1])
z =
   -5.1378    5.1765   16.0000
匿名函数在定义过程中可以调用任何Matlab函数（包括用户自定义的函数），也可以使用当时Matlab工作区中存在的任何变量。例如，下例中的匿名函数b_humpsab在定义时就是用了Matlab工作区中的变量a和b：
>> a = -.3;b=-.9;
>> b_humpsab = @ (x) 1./((x+a).^2+0.01)+1./((x+b).^2+0.04)-6;
>> b_humpsab([-1 0 1])
ans =
   -5.1378    5.1765   16.0000
我们看到：b_humpsab在定义时引用了前面定义的变量a和b。应当注意，当a或b的值变化时，匿名函数并不改变，这是因为函数句柄值捕捉它创建时刻的变量的值，并不随变量的变化而变化。
>> a = 0;   %changing the value of a does not chang the function.
>> b_humpsab([-1 0 1])          %evaluate again , get the same results
ans =
   -5.1378    5.1765   16.0000
我们也可以针对一个内置函数或一个M文件函数创建匿名函数句柄。下面函数分别给予M文件函数humps和Matlab内置函数cos创建了两个文件句柄：
>> c_Mfile = @humps      %function handle for M-file function
c_Mfile =
    @humps
>> c_Mfile(1)      %evaluate humps(1)
ans =
    16
>> c_builtin = @cos
c_builtin =
    @cos
>> c_builtin(pi)
ans =
    -1
由上面可知，要创建一个内置函数或一个M文件函数的句柄也很容易，只要在等号右边使用@符号，并在该符号后紧跟内置函数名或M文件函数名即可。
我们还可以利用单元数组同时创建多个内置函数和M文件函数的句柄，验证这些函数时，只要引用该函数所在单元即可。例如，下面的代码将上面两个独立创建的句柄利用一个单元数组c_dan创建在一起：
>> c_dan = {@humps @cos}
c_dan =
    @humps    @cos
>>c_dan{1}(1)              %evaluate humps(1)
ans =
    16
>> c_dan{2}(pi)             %evaluate cos(pi)
ans =
    -1
Matlab还专门提供了一些函数来处理和应用句柄。例如，函数functions将返回一个句柄的详细信息，如下例所示：
>> functions(c_Mfile)
ans =
    function: 'humps'
        type: 'simple'
        file: 'C:\Program Files\MATLAB\R2009a\toolbox\matlab\demos\humps.m'
>> functions(c_builtin)
ans =
    function: 'cos'
        type: 'simple'
        file: ''
>> functions(b_humps)
ans =
     function: '@(x)1./((x+a).^2+0.01)+1./((x+b).^2+0.04)-6'
         type: 'anonymous'
         file: ''
    workspace: {[1x1 struct]}
需要注意的是：fuctions函数通常只在调试程序时使用，因为它的返回值很容易发生变化。