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

佛自业障

（i）首先编写M文件mente.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)-8*x(1)-2*x(2)-3*x(3)...
. f7 A0 y! o% b2 R   -x(4)-2*x(5);
g(1)=sum(x)-400;
8 G0 b' J  }' c& eg(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;
$ N. A6 n# W( o% c6 S* wg(4)=x(3)+x(4)+5*x(5)-200;
（ii）编写如下程序求问题的解：
* G! r! L) O2 d6 {9 u: Y( U; ^+ H
rand('state',sum(clock));
p0=0;
4 e5 `, |, O3 o5 ~tic
0 ~* E8 x2 W* ^; Z) Efor i=1:10^5
   x=99*rand(5,1);
9 Z. t  r$ L. q. _; Q; R# \x1=floor(x);x2=ceil(x);
[f,g]=mengte(x1);
if sum(g<=0)==4
   if p0<=f
; V+ b6 l5 w) y; f3 c$ S# c0 N. B  L      x0=x1;p0=f;
   end
end
[f,g]=mengte(x2);
if sum(g<=0)==4
   if p0<=f
0 \9 Q* a% j, d" x, h6 y: Q5 a0 C      x0=x2;p0=f;
   end
end
end
x0,p0
# A/ P4 K4 Z: K1 u+ ^" r9 G& Ttoc


( `. }$ Z- c" c4 a% C- Y5 c) [' b* o