Question

我一直在尝试使用以下代码在输入为向量的特定点找到函数的渐变，函数返回标量。

以下是我尝试计算渐变的函数。

%fun.m    
function [result] = fun(x, y)
     result = x^2 + y^2;

这就是我称之为渐变的方式。

f = @(x, y)fun(x, y);
grad = gradient(f, [1 2])

但是我收到以下错误

octave:23> gradient(f, [1 2])
error: `y' undefined near line 22 column 22
error: evaluating argument list element number 2
error: called from:
error:    at line -1, column -1
error:   /usr/share/octave/3.6.2/m/general/gradient.m at line 213, column 11
error:   /usr/share/octave/3.6.2/m/general/gradient.m at line 77, column 38

如何解决此错误？

Answer 1

我的猜测是gradient无法在2D函数句柄上运行，因此我就这样做了。请考虑以下 lambda-flavoring 解决方案：

让fz成为您某些功能的函数句柄

fz = @(x,y)foo(x,y);

然后考虑这段代码

%% definition part:
only_x = @(f,yv) @(x) f(x,yv);  %lambda-like stuff, 
only_y = @(f,xv) @(y) f(xv,y);  %only_x(f,yv) and only_y(f,xv) are
                                %themselves function handles

%Here you are:
gradient2 =@(f,x,y) [gradient(only_x(f,y),x),gradient(only_y(f,x),y)];

用作

gradient2(fz,x,y);

最后一点测试：

fz = @(x,y) x.^2+y.^2
gradient2(f,1,2);

结果

octave:17> gradient2(fz,1,2)
ans =

    2   4

八度：在特定点计算多维函数的梯度

1 个答案: