matlab 求解目标规划的序贯算法_matlab 求解目标规划的序贯式算法

求解目标规划的序贯算法

clc,clear;
x=optimvar('x',2,'LowerBound',0);
dp=optimvar('dp',4,'LowerBound',0);
dm=optimvar('dm',4,'LowerBound',0);
p=optimproblem('ObjectiveSense','min');
p.Constraints.con1=2*sum(x)<=12;
con2=[200*x(1)+300*x(2)+dm(1)-dp(1)==1500
    2*x(1)-x(2)+dm(2)-dp(2)==0
    4*x(1)+dm(3)-dp(3)==16
    5*x(2)+dm(4)-dp(4)==15];
p.Constraints.con2=con2;
goal=10000*ones(3,1);
mobj=[dm(1);dp(2)+dm(2);3*dp(3)+3*dm(3)+dp(4)];
for i=1:3
    p.Constraints.cons3=[mobj<=goal];
    p.Objective=mobj(i);
    [sx,fval]=solve(p);
    fprintf('第%d级目标计算结果如下：\n',i)
    fval,xx=sx.x,sdm=sx.dm,sdp=sx.dp
    goal(i)=fval;
end
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

result

Solving problem using linprog.

Optimal solution found.

第1级目标计算结果如下：

fval =

 0
1

xx =

    1.875
     3.75
1
2

sdm =

        0
        0
      8.5
        0
1
2
3
4

sdp =

        0
        0
        0
     3.75
1
2
3
4

Solving problem using linprog.

Optimal solution found.

第2级目标计算结果如下：

fval =

 0
1

xx =

    1.875
     3.75
1
2

sdm =

        0
        0
      8.5
        0
1
2
3
4

sdp =

        0
        0
        0
     3.75
1
2
3
4

Solving problem using linprog.

Optimal solution found.

第3级目标计算结果如下：

fval =

       29
1

xx =

        2
        4
1
2

sdm =

 0
 0
 8
 0
1
2
3
4

sdp =

      100
        0
        0
        5
1
2
3
4

optimvar使用版本为2017b以后