数学建模常用模型20 ：非线性规划之蒙特卡洛法（优化模型）

佛自业障

& _9 H1 G7 U+ I6 o" J  @（i）首先编写M文件mente.m定义目标函数f 和约束向量函数g，程序如下：
( e3 w" k. F; k  Q# |% Zfunction [f,g]=mengte(x);
f=x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)-8*x(1)-2*x(2)-3*x(3)...
   -x(4)-2*x(5);
g(1)=sum(x)-400;
g(2)=x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800;
g(3)=2*x(1)+x(2)+6*x(3)-200;
) Y7 T5 z5 I! l- U6 E& u4 y  B! Yg(4)=x(3)+x(4)+5*x(5)-200;
（ii）编写如下程序求问题的解：
4 ]0 J3 y* ?# B$ @/ i
rand('state',sum(clock));
p0=0;
tic
for i=1:10^5
9 S: H( ^. a& j" b   x=99*rand(5,1);
x1=floor(x);x2=ceil(x);
[f,g]=mengte(x1);
if sum(g<=0)==4
2 \! i0 {/ x9 z   if p0<=f
      x0=x1;p0=f;
6 b3 a. h1 y/ `  T2 M0 i7 c   end
end
# G  C# R9 \5 e7 Q! v+ t5 J[f,g]=mengte(x2);
if sum(g<=0)==4
1 L: a. e( o$ g   if p0<=f
) c. |& t) t+ B" g/ a4 H      x0=x2;p0=f;
6 P2 u+ z  s( `7 `   end
+ h1 j3 o/ T8 e0 Iend
end
, v* @5 K8 ?5 V- F; o* Kx0,p0
" p! M+ X0 r5 C1 N' N# B: Btoc


; H" D3 p: Q5 o, g: {
2 k- J6 I3 j9 z$ F0 D