基于syms实现含60个动态变量函数的Hessian矩阵自动求解

哈哈，手动定义60个变量确实是个噩梦！我之前做优化问题的时候也遇到过类似需求，用Matlab的符号工具箱完全可以搞定动态变量生成和自动求Hessian，给你一套灵活的方案，以后改变量数或者函数形式都超方便~

1. 动态创建符号变量 #

不用逐个敲syms x1 x2 ... x60，直接用符号变量数组就能一键生成任意数量的变量：

n = 60; % 变量总数，想改多少直接改这个数
syms x(1:n) % 自动生成x₁到x₆₀的符号变量数组

不管你以后要改成30个还是100个变量，只需要调整n的值就行，省超多时间。

2. 灵活定义目标函数f #

你可以根据自己的需求随便写函数形式，用数组索引就能访问对应的变量。举个示例（你直接替换成自己的函数就行）：

% 示例函数：平方和 + 前30个变量与后30个变量的乘积和
f = sum(x.^2) + sum(x(1:30).*x(31:60));

要是以后想修改函数逻辑，比如调整乘积项的变量范围，直接改索引或者表达式就行，完全不用碰变量定义的部分。

3. 自动计算Hessian矩阵 #

用Matlab内置的hessian函数一步到位，直接传入函数f和变量数组x就行：

H = hessian(f, x);

这时候H就是一个60×60的符号矩阵，所有二阶偏导数都自动计算好了。

小技巧：优化矩阵显示 #

因为矩阵太大，直接打印可能看不全，你可以用这些方法优化显示：

• 用simplify(H)简化矩阵中的表达式
• 用pretty(H)以更友好的格式输出
• 用disp(H)快速打印矩阵

完整示例代码 #

把所有步骤整合起来，直接运行就能得到结果：

% 动态生成符号变量并求解Hessian矩阵
clear; clc;

% 1. 定义变量数量并生成符号变量
n = 60;
syms x(1:n)

% 2. 定义目标函数（可按需修改）
f = sum(x.^2) + sum(x(1:floor(n/2)).*x(floor(n/2)+1:n));

% 3. 计算Hessian矩阵
H = hessian(f, x);

% 可选：简化并显示矩阵
H_simplified = simplify(H);
disp('简化后的Hessian矩阵：');
pretty(H_simplified);

这个方案的核心就是利用符号变量数组实现动态生成，不管变量数量怎么变，函数形式怎么改，都能快速适配，完全不用手动重复劳动~

内容的提问来源于stack exchange，提问作者Naomi Li