Matlab求解规划问题之 fgoalattain函数

语法：[X,FVAL] = fgoalattain(fun,x0,goal,a,b,Aeq,Beq,LB,UB,nonlcon)；

    X 为最终解 ， FVAL为最终解对应的函数值
    *注意：求最大值时，结果FVAL需要取反*
    
    fun 是定义的决策函数，通常通过M文件或者匿名函数进行定义。
    *注意：当所求为最大值时，系数需要取反*
    
    x0 为初始值
    
    goal 为欲达到的目标，通常通过linprog函数先计算得到每个决策函数目标值
    
    a 为约束条件中不等式组的系数矩阵 ，a的列数等于f的列数
    *注意：当不等号为  > 或  ≥ 时，矩阵需要取反*
    
    b 为约束条件中不等式组右边的值
    *注意：当不等号为  > 或  ≥ 时，矩阵需要取反*
    
    Aeq 为约束条件中等式组的系数矩阵 ，Aeq的列数等于f的列数
    
    Beq 为约束条件中等式组右边的值
    
    LB、UB 是解的范围
    
    nonlcon 为定义的向量函数

例题：

matlab求解代码：

clear
close all
clc
a = [-1 -1 0 0;0 0 -1 -1;3 0 2 0;0 3 0 2];
b = [-30 -30 120 48]';
c1 = [-100 -90 -80 -70];
c2 = [0 3 0 2];
fun = @(x) [c1;c2]*x;
[x1,g1] = linprog(c1,a,b,[ ],[ ],zeros(4,1));
[x2,g2] = linprog(c2,a,b,[ ],[ ],zeros(4,1));

g3 = [g1,g2];
[x,fval] = fgoalattain(fun,rand(4,1),g3,abs(g3),a,b,[ ],[ ],zeros(1,4));

结果为:

x = [19.0652 10.9348 31.4023 0];
fval = [-5402.8 32.8045];
即，max Z1 = 5402.8, min Z2 = 32.8 。

linprog()函数用法