0

我试图用 Maple 中的线性等式约束最小化二次能量。目前我正在发布如下内容：

with(Optimization):
p := (t) -> c3*t^3 + c2*t^2;
m := Minimize(int(diff(p(t),t)^2,t=0..1),{eval(p(t),t=1)=1,eval(diff(p(t),t),t=1)=0});

但这似乎给了我一个带有浮点错误的数值优化解决方案：

m := [1.19999999999997, [c2 = 3.00000000000000, c3 = -2.00000000000000]]

（正确答案是 m:= [6/5,[c2=3,c3=-2]]）

有没有办法使用枫树象征性地计算解决方案？

我宁愿不必自己计算拉格朗日量。我希望有一个像 symbolic=true 这样的标志。

4

1 回答 1

1

逐步应用拉格朗日乘数。这很容易做到。

p:= t-> c3*t^3+c2*t^2:
Obj:= int(D(p)^2,  0..1):
Con:= {p(1)=1, D(p)(1)=0}:
L:= Obj - l1*lhs(Con[1]) - l2*lhs(Con[2]):
solve(Con union {diff(L,c3), diff(L,c2)});
              /                      12       -1\ 
             { c2 = 3, c3 = -2, l1 = --, l2 = -- }
              \                      5        5 / 
eval(Obj, %);
                               6
                               -
                               5
于 2013-08-21T17:47:36.620 回答