MatLab中带约束的函数最大化

我想在MatLab中最大化这个函数- http://goo.gl/C6pYP

但是变量x, y和z只能是非负整数。有什么想法如何在MatLab中实现吗?

您想要整数的事实使得这个问题很难解决(即无法在合理的时间内解决)。

你将不得不使用一些通用的优化器，并通过蛮力尝试许多启动条件。你不能保证找到一个全局最大值。有关更多可能的优化器，请参阅Matlab的优化包。

您必须将问题表述为ILP(整数线性规划)。要解决ILP，您需要对matlab LP求解器的输入进行一些更改。您还可以从LP求解器中获得解决方案，然后将解决方案四舍五入为整数。解决方案可能不是最优的，但将是接近的。

还可以在文件交换站点使用混合整数线性规划求解器，而文件交换站点又使用LP求解器。对于二进制变量，您可以使用matlab二进制整数编程求解器。

好吧，幸运的是问题的大小很小，所以我们可以用蛮力来解决。

首先得到一些上限，下面是如何为x做的:

xmax= 0;
while 12546975*xmax+525*xmax^2<=4000000000
    xmax=xmax+1;
end

这给了我们所有三个变量的上限。现在我们可以看到，这些极限的乘积不是很多，所以我们可以尝试所有的解。

bestval = 0;
for x = 0:xmax
    for y = 0:ymax
        for z = 0:zmax
            val = 3*x+6*y+9*z;
            if val> bestval && 12546975*x+525*x^2+25314000*y+6000*y^2+47891250*z+33750*z^2<=4000000000
                bestval = val;
                best = [x y z];
            end
        end
    end
end
best, bestval

这可能不是最有效的方式，但它应该非常容易阅读。

y和z(152,79)的最大值不是很高，所以我们可以逐一检查以快速找到解(在我的笔记本电脑中只需0.040252秒)。

我的matlab代码:

function [MAX,x_star,y_star,z_star]=stackoverflow1
%maximize 3x+6y+9z
% s.t.    12546975x+525x^2+25314000y+6000y^2+47891250z+33750z^2<=4000000000
MAX=0;
y_max=solver(6000,25314000,-4000000000);
z_max=solver(33750,47891250,-4000000000);
for y=0:floor(y_max)
    for z=0:floor(z_max)
        x=solver(525,12546975,+25314000*y+6000*y^2+47891250*z+33750*z^2-4000000000);
        x=floor(x);
        if isnan(x) || x<0
            break;
        end
        if 3*x+6*y+9*z>MAX
            MAX=3*x+6*y+9*z;
            x_star=x;
            y_star=y;
            z_star=z;
        end
    end
end
end
function val=solver(a,b,c)
% this function solve equation a*x^2+b*x+c=0.
% this equation should have two answers,this function returns the bigger one only.
if b*b-4*a*c>=0
    val=(-b+sqrt(b*b-4*a*c))/(2*a);
else
    val=nan; % have no real number answer.
end
end

解决方案是:

MAX =

945年

x_star =

287年

y_star =

z_star =

相关内容

最新更新

热门标签：