0

我在处理渐变和粗麻布的函数句柄时遇到问题。

我有以下代码:

syms x1 x2

x = [x1,x2];

% Define the function phi(x)
phi = @(x1,x2) 10*x1^4 - 20*x1^2*x2 + 10*x2^2 + x1^2 - 2*x1 + 5; 

% Define the gradient of the function phi(x)
gradphi = @(x1,x2) jacobian(phi,x).';

% Define the Hessian of phi(x)    
hessphi = @(x1,x2) jacobian(gradphi,x);

现在当我输入命令终端时:

phi(1,2)

我得到一些标量值。

但是当我输入

gradphi(1,2)
hessianphi(1,2)

我想要在这些点评估梯度的相应向量。

对于渐变,我只是得到

EDU>> gradphi(1,2)

ans =

2*x1 - 40*x1*x2 + 40*x1^3 - 2
            - 20*x1^2 + 20*x2

这只是 grad 向量函数。但我想要输入 x1 = 1 和 x2 = 2 的实际数值结果。

EDU>> hessphi(1,2)

返回错误。

我不确定为什么。

4

1 回答 1

2

对于数组的乘法和除法,您需要在 and 的定义中使用逐元素运算符.*,和。否则,Matlab 将尝试做矩阵乘法/除法/幂,这不会顺利。./.^hessianphigradphi

/aside:搜索错误消息的文本将显示最可能的错误原因。

于 2014-02-21T08:51:43.343 回答