matlab-整数规划(非线性规划之蒙特卡洛法（随机取样法）)

首先编写M 文件mengte.m 定义目标函数f 和约束向量函数g，程序如下：
function [f,g]=mengte(x);%定义目标函数f和约束向量函数g

f=x()^+x()^+*x()^+*x()^+*x()-*x()-*x()-*x()-x()-*x();

g=[sum(x)-

x()+*x()+*x()+x()+*x()-

*x()+x()+*x()-

x()+x()+*x()-];注:文件名一定要与定义的函数名一致！

编写M文件mainint.m如下求问题的解：
rand('state',sum(clock));

p0=;

tic  %tic表示计时的开始，toc表示计时的结束

for i=:^

    x=*rand(,)  %rand函数产生0-1之间的随机数，表明了0=<xi<=

x1=floor(x);  %floor朝负无穷方向取整

x2=ceil(x);   %ceil朝正无穷方向取整

[f,g]=mengte(x1);

if sum(g<=)==

    if p0<=f

        x0=x1;p0=f;

    end

end

[f,g]=mengte(x2);

if sum(g<=)==

    if p0<=f

        x0=x2;p0=f;

    end

end

end

x0,p0

toc

MATLAB中四个取整函数具体使用方法如下：

Matlab取整函数有: fix, floor, ceil, round.
fix
朝零方向取整，如fix(-1.3)=-1; fix(1.3)=1;
floor
朝负无穷方向取整，如floor(-1.3)=-2; floor(1.3)=1;
ceil
朝正无穷方向取整，如ceil(-1.3)=-1; ceil(1.3)=2;
round
四舍五入到最近的整数，如round(-1.3)=-1;round(-1.52)=-2;round(1.3)=1;round(1.52)=2。

巴特西

matlab-整数规划(非线性规划之蒙特卡洛法（随机取样法）)

最新文章

热门文章