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

佛自业障

% w+ d* D, I: w: d7 F（i）首先编写M文件mente.m定义目标函数f 和约束向量函数g，程序如下：
function [f,g]=mengte(x);
0 h& J: N, U) ?" A% l/ W2 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;
g(3)=2*x(1)+x(2)+6*x(3)-200;
g(4)=x(3)+x(4)+5*x(5)-200;
# I7 g) \' l1 _（ii）编写如下程序求问题的解：
+ W: N$ r. V1 T) L) ~4 z) ]' Y
rand('state',sum(clock));
9 P5 D, ?4 P' ]# R/ vp0=0;
! d. h" E+ y" K0 b% Rtic
/ {% t9 W# a0 r; [, }7 hfor i=1:10^5
7 y5 |8 n8 @+ x7 i- t   x=99*rand(5,1);
, G' f) K7 d5 I  W) {x1=floor(x);x2=ceil(x);
[f,g]=mengte(x1);
if sum(g<=0)==4
   if p0<=f
      x0=x1;p0=f;
/ Z4 q, ~, C. s* P* l7 g; X   end
9 t# O  s4 b! [% \2 wend
[f,g]=mengte(x2);
7 |3 L' {% a4 \/ [9 o# j9 {if sum(g<=0)==4
$ B8 D% g( M% f+ V   if p0<=f
8 x) x9 I& l- C2 P& M7 p; c      x0=x2;p0=f;
( T" T+ Z9 ^" \) _! T4 `  @   end
) `- z8 F! I8 h* Yend
end
x0,p0
toc
" t( \) F- P% v3 v7 R
- q% Q( [# q: G8 R: D! l
7 [' d; H7 @6 V) M- R6 k
! E; I7 J$ ^0 Q7 b
) L4 B' U: P3 g3 t8 i0 @