matlab:quadprog提示功能不对称时

Question

当我使用给定的功能 F matlab 输出运行 quadprog 时:

Warning: Your Hessian is not symmetric.
Resetting H=(H+H')/2. 

但是，检查泛函和它的转置之间的区别:

>> max(max(abs(F-F')))

ans =

   (1,1)     7.1054e-015

表明它们实际上是一样的。即使函数因机器错误而接近对称，quadprog 是否会输出此警告？

Answer 1

用于在 quadprog 源中存在反对称性的测试是

norm(H-H',inf) > eps

即它检查无穷范数是否大于机器 epsilon。由于

7.1054e-15 > eps = 2.22e-16

您的 Hessian 矩阵显示为非对称矩阵。