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

佛自业障

（i）首先编写M文件mente.m定义目标函数f 和约束向量函数g，程序如下：
# @  @! p9 p- S. T* F# ~function [f,g]=mengte(x);
6 o/ S$ V  j& E) V6 sf=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;
0 M3 a8 w, ^" g9 m- [! m" Tg(3)=2*x(1)+x(2)+6*x(3)-200;
g(4)=x(3)+x(4)+5*x(5)-200;
（ii）编写如下程序求问题的解：

* ?6 f" g! F* l; e/ r' ]% t* vrand('state',sum(clock));
6 I1 M+ u# k, F  C# h  s( T- bp0=0;
. }. T4 t: o% B, J0 X9 Wtic
6 K/ N- l& X' w0 w' Sfor i=1:10^5
' P) U4 r0 w2 C0 O- _   x=99*rand(5,1);
x1=floor(x);x2=ceil(x);
[f,g]=mengte(x1);
if sum(g<=0)==4
   if p0<=f
$ `6 P- \2 w! H( v+ i8 k: Q      x0=x1;p0=f;
. a0 J0 s) v# V! G3 K' A5 A( L' D   end
# U5 W5 L6 M1 `1 O" _end
' E; z$ v- d: N, A[f,g]=mengte(x2);
3 {3 h' d4 B* v9 }8 q/ O  bif sum(g<=0)==4
   if p0<=f
      x0=x2;p0=f;
8 _" i) ^; p) w2 b# m   end
end
7 R- G7 z. a; h3 zend
5 i) t% q+ |9 ?- v! L3 Wx0,p0
toc