这个例子解决了典型的线性规划问题
加载sc50b.mat
文件,其中包含矩阵和向量一个
,Aeq
,b
,说真的
,f
,下界磅
.
负载sc50b
这个问题有48个变量,30个不等式,20个等式。
disp(大小(A))
30 48
disp(大小(Aeq))
20 48
属性的选项对偶单纯形
算法和迭代显示。
选择= optimoptions (@linprog,“算法”,对偶单纯形的,“显示”,“通路”);
这个问题没有上界,所以有个集合乌兰巴托
来[]
.
乌兰巴托= [];
通过打电话来解决问题linprog
.
[x, fval exitflag、输出]=...Aeq linprog (f, A, b,说真的,磅,乌兰巴托,选项);
LP预处理去除2个不等式,16个等式,16个变量,26个非零元素。Iter Time Fval Primal feas Dual feas 0 0.023 0.00000000e +00 0.00000000e +00 1.305013e+00 8 0.043 -1.587073e+02 3.760622e+02 0.00000000e +00 33 0.045 -7.000000e+01 0.00000000e +00 0.00000000e +00找到最佳解决方案。
检查退出标志、解决方案的目标函数值和所使用的迭代次数linprog
为了解决这个问题。
exitflag、fval output.iterations
exitflag = 1
fval = -70
ans = 33
您还可以在迭代显示中找到目标函数值和迭代次数。