MATLAB实现蒙特卡罗方法.pdf

MATLAB实实现现蒙蒙特特卡卡罗罗⽅⽅法法

1.⾸先编写M⽂件mengte.m定义⽬标函数f和约束向量g

%%%蒙特卡洛法

%%定义⽬标函数和约束向量函数

function[f,g]=mengte(x);

f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)-*x(1)-2*x(2)-3*x(3)-x(4)-2*x(5);

g=[sum(x)-400

x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-00

2*x(1)+x(2)+6*x(3)-200

x(3)+x(4)+5*x(5)-200];

end

2.编写M⽂件mainint.m，求问题的

1%%%通过迭代求

2rand(state,sum(clock));%产⽣5个0-1之间的随机数

3p0=0;%p0代表⽬标值z，x0代表⽬标x

4tic

5fori=1:10^6%采集10^6个点

6x=99*rand(5,1);%0=x=99

7x1=floor(x);x2=ceil(x);%x1是x向下取整，x2是x向上取整

[f,g]=mengte(x1);%对x1求值

9ifsum(g=0)==4%若满⾜4个则符合条件

10ifp0=f

11x0=x1;

12p0=f;

13end

14end

15[f,g]=mengte(x2);%对x2求值

16ifsum(g=0)==4

17ifp0=f

1x0=x2;

19p0=f;

20end

21end

22end

23x0,p0

24toc