このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。
mpsread
LP および MILP 最適化データの MPS ファイルの読み取り
説明
problem = mpsread(mpsfile)intlinprog または linprog ソルバーが受け入れる構造体にデータが返されます。
problem = mpsread(mpsfile,'ReturnNames',true)problem 構造体に、mpsfile 内に変数と制約の名前が格納された variableNames フィールドと constraintNames フィールドを追加します。
例
mps ファイルを読み込み、記述されている問題を解きます。
公開リポジトリから eil33-2.mps ファイルを読み込みます。問題のタイプを表示します。
gunzip("https://miplib.zib.de/WebData/instances/eil33-2.mps.gz") problem = mpsread("eil33-2.mps")
Running HiGHS 1.7.0 (git hash: n/a): Copyright (c) 2024 HiGHS under MIT licence terms
problem = struct with fields:
f: [4516×1 double]
Aineq: [0×4516 double]
bineq: [0×1 double]
Aeq: [32×4516 double]
beq: [32×1 double]
lb: [4516×1 double]
ub: [4516×1 double]
intcon: [4516×1 double]
solver: 'intlinprog'
options: [1×1 optim.options.Intlinprog]
problem.intcon が空ではなく、problem.solver が 'intlinprog' であることがわかります。これは、整数線形計画問題です。
反復表示を抑制し、ソルバーの進行状況に応じてプロットを生成するように、オプションを変更します。
options = optimoptions("intlinprog",... Display="final",PlotFcn=@optimplotmilp); problem.options = options;
intlinprog を呼び出して問題を解きます。
[x,fval,exitflag,output] = intlinprog(problem);
Optimal solution found. Intlinprog stopped because the objective value is within a gap tolerance of the optimal value, options.AbsoluteGapTolerance = 1e-06. The intcon variables are integer within tolerance, options.ConstraintTolerance = 1e-06.
mps ファイルを読み込んで、その変数名と制約名を取得します。
公開リポジトリから eil33-2.mps ファイルを読み込みます。返された問題構造体を表示します。
gunzip("https://miplib.zib.de/WebData/instances/eil33-2.mps.gz") problem = mpsread("eil33-2.mps",ReturnNames=true)
Running HiGHS 1.7.0 (git hash: n/a): Copyright (c) 2024 HiGHS under MIT licence terms
problem = struct with fields:
f: [4516×1 double]
Aineq: [0×4516 double]
bineq: [0×1 double]
Aeq: [32×4516 double]
beq: [32×1 double]
lb: [4516×1 double]
ub: [4516×1 double]
intcon: [4516×1 double]
solver: 'intlinprog'
options: [1×1 optim.options.Intlinprog]
variableNames: [4516×1 string]
constraintNames: [1×1 struct]
各タイプの最初の数個の名前を表示します。
problem.variableNames(1:4)
ans = 4×1 string
"x1"
"x2"
"x3"
"x4"
problem.constraintNames.eqlin(1:4)
ans = 4×1 string
"c1"
"c2"
"c3"
"c4"
問題に不等式制約は含まれていません。
problem.constraintNames.ineqlin
ans = 0×1 empty string array
