Drake中调用GetSolution时1D数组与取反数组返回形状为何不一致？

问题复现代码 #

a = prog.NewContinuousVariables(1, "a")
neg_a = -a

调用result.GetSolution(a)返回形状为(1,)，而调用result.GetSolution(neg_a)返回形状为(1, 1)，二者形状存在差异，容易引发后续广播错误。

原因解析 #

这类符号式优化框架的运算逻辑里，原始变量a是一维向量类型（通过(1,)参数创建），但当你对它执行取负这类一元运算时，框架会默认按照矩阵运算规则处理：将一维向量转换为二维矩阵（单元素矩阵(1,1)）。

这是因为框架底层的符号计算系统以矩阵为核心单元，统一用矩阵维度来处理所有运算，哪怕是单元素的简单操作，目的是避免不同维度运算时的逻辑歧义，但这种处理会导致经过运算后的变量和原始变量的维度类型不一致。

解决办法 #

• 手动扁平化运算后的结果：比如用np.ravel(result.GetSolution(neg_a))将(1,1)形状转为(1,)
• 创建变量时明确使用二维形状：把prog.NewContinuousVariables(1, "a")改成prog.NewContinuousVariables(1, 1, "a")，这样原始变量和取负后的变量形状都会是(1,1)，保持一致

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