私は原則的にsolve.QP(quadprog)で解決できるはずの非常に単純な二次問題を持っています:
max x^2+3*y^2 subject to x>=0, y>=0, x+y <=1
問題は、最大化を必要な最小化に変えると、行列が (明らかに) 負定値になることです。小さな固有値や丸めに関する問題ではありません。
solve.QP は一般的なソルバーだと思っていましたが、オンライン資料を読んで努力したにもかかわらず、solve.QP を使用してコンパクトな領域 (線形制約によって定義された) で正定値二次形式を最大化できないようです。
本当?
私はこれと同様の問題を他の関数で解決できることを知っていますが(正常constrOptim
に動作します)、ラグランジュ乗数をマキシマイザーに接続したいと思っていました。
(非常に効率的な)solve.QPを使用して上記の問題を解決し、正の定義性に関連する非対称の制限を克服する方法を提案できますか?