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

佛自业障

3 a% s) Y9 F! z. b7 s2 h（i）首先编写M文件mente.m定义目标函数f 和约束向量函数g，程序如下：
/ v- ~; ~; D3 r# f: l9 \, b. T, vfunction [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)...
6 T: W0 ?, g, _: v+ K   -x(4)-2*x(5);
3 b  u7 F; r5 }. sg(1)=sum(x)-400;
% D7 d* ]. F; U6 b; b# D  A! ^g(2)=x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800;
7 i8 v7 ]' ]  m' f. f. f0 ~" d. O9 Mg(3)=2*x(1)+x(2)+6*x(3)-200;
g(4)=x(3)+x(4)+5*x(5)-200;
# y% m" @' J: |7 _" m; K' }0 @7 [（ii）编写如下程序求问题的解：
1 i2 m: L& G( m1 j8 J0 F/ \) B3 i1 \
rand('state',sum(clock));
$ |8 h  d( Y5 ~9 W1 Qp0=0;
- l# H8 Z" s( C/ A* l  F/ D+ Etic
for i=1:10^5
# j1 d# l9 R: e" D" n7 f   x=99*rand(5,1);
x1=floor(x);x2=ceil(x);
[f,g]=mengte(x1);
if sum(g<=0)==4
   if p0<=f
. h5 h- a( S7 p" }: \0 d& H      x0=x1;p0=f;
1 f( e% k5 }9 P" v7 ]( P( ]   end
; G! m4 h* J# A: ]/ Oend
: o  E0 U4 h7 q- h[f,g]=mengte(x2);
if sum(g<=0)==4
   if p0<=f
3 C4 i& f5 G  u0 b) M4 B* P& w      x0=x2;p0=f;
! S' b; x. N- L" j4 l& D9 q. a) g   end
8 s' k5 Z4 w8 S# kend
& W' P& M4 Y) V2 Tend
x0,p0
toc
- O8 B4 _9 f" ^  p! A: B' D& A+ J
) p& G8 p* T0 T+ k6 U

) ~/ X( H% X& \8 c5 v6 P
' M- _! d( V! P7 Q  a& W4 \