線形計画問題に対する内点法 in Python - matsulibの日記内点法にもいろいろある。ここでやるのはもちろん一番簡単なやつ。 ・主アフィンスケーリング法 - 経営工学専攻 - 東京工業大学 http://www.me.titech.ac.jp/~mizu_lab/text/PDF-IP/IP2A-affine.pdf を見てPythonでナイーブに実装した。 ・行列演算はライブラリ任せ（これは速度面でも好ましい） ・本体はコメントを除けば20行足らず ・例外処理なし ・バグ放置 こんな感じで簡単に線形計画問題の最適解が求まります。 例題 以前の記事で商用ソルバを使って解いた問題をまた解く。 GurobiとPythonで数理最適化 - matsulibの日記 とりあえず前処理として標準形に変形する。 加えて、手続きの出発地点となる初期内点x0を与える必要がある。 下のコードでは x0=[20, 10, 2200, 1400, 670] としている。 ソ
