编写一个具体求黑塞矩阵的matlab代码，要具体的例子

黑塞矩阵（Hesse Matrix）也称为雅可比矩阵，它在数值计算中用于描述函数的一阶导数信息。对于给定的二元或多元函数f(x)，其黑塞矩阵是关于x的元素的Hessian矩阵，即：

\[ H(f) = \begin{bmatrix}
\frac{\partial^2 f}{\partial x_1^2} & \frac{\partial^2 f}{\partial x_1 \partial x_2} & \cdots \\
\frac{\partial^2 f}{\partial x_2 \partial x_1} & \frac{\partial^2 f}{\partial x_2^2} & \cdots \\
\vdots & \vdots & \ddots
\end{bmatrix}

在MATLAB中，你可以通过符号数学工具箱来计算黑塞矩阵。假设有一个简单的二次函数f(x,y) = x^2 + xy + y^2，我们可以编写如下代码来求解它的Hesse矩阵：

syms x y % 定义变量x和y为符号
f = x^2 + x*y + y^2; % 定义函数

% 计算一阶导数
dfdx = diff(f, x); % 对x的偏导数
dfdy = diff(f, y); % 对y的偏导数

% 计算二阶导数构成Hesse矩阵
hesse_matrix = [diff(dfdx, x), diff(dfdx, y); diff(dfdy, x), diff(dfdy, y)];

disp(hesse_matrix); % 输出黑塞矩阵

运行此代码后，你会看到函数f(x,y)的Hesse矩阵。