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

佛自业障

! Q( ?' I3 |8 |9 w6 o（i）首先编写M文件mente.m定义目标函数f 和约束向量函数g，程序如下：
function [f,g]=mengte(x);
0 f' v" e$ V5 h) ~1 g! of=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)...
. t" {' R, p, x6 _   -x(4)-2*x(5);
g(1)=sum(x)-400;
  J- B1 \5 |/ [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;
g(4)=x(3)+x(4)+5*x(5)-200;
. R, x0 R2 ^# s* D0 [  S2 @1 g（ii）编写如下程序求问题的解：

rand('state',sum(clock));
3 A% u6 j" `5 z1 g! z) b0 Yp0=0;
tic
for i=1:10^5
   x=99*rand(5,1);
6 N5 y. C+ X, |* g/ V: A5 }6 hx1=floor(x);x2=ceil(x);
. _$ ]& R/ D5 _' j[f,g]=mengte(x1);
if sum(g<=0)==4
: c) O1 X' q: g) u   if p0<=f
      x0=x1;p0=f;
0 a# w, J1 x5 l$ `, A- [# {   end
end
[f,g]=mengte(x2);
if sum(g<=0)==4
   if p0<=f
      x0=x2;p0=f;
4 k- m! a6 A: k# g7 v   end
/ Z( E  ?" M, wend
end
/ A, r- S+ ?! v9 m7 k' Bx0,p0
0 l: ?: w- w; y3 V7 @; {toc

/ k$ e- R0 o% c9 o! y* C3 h