根据数学公式求解非线性规划 - 数学建模社区-数学中国

%MC搜索5 Q& h: D7 s# j$ {2 h+ b+ {) L
%复杂度低随机性强6 O5 M6 K' L$ S; G# e
r1=unifrnd(0,10,100000,1); %产生x1的n*1随机矩阵
( s' z; I0 O f5 s& F- |$ y
r2=unifrnd(0,10,100000,1); %产生x2的n*1随机矩阵
( n! ^# f4 i7 a7 ~7 a
sol=[r1(1) r2(1)];
2 @& y8 y+ R/ A# P
z0=-inf; %z0初始化
i$ m3 _! T1 A h
f=inline('-2*x(1)^2-x(2)^2+x(1)*x(2)+8*x(1)+3*x(2)','x'); %目标函数" t5 e- Q3 z; ^! I8 p# x
for i=1:100000
4 u7 M/ T6 z0 ?# t
x1=r1(i);6 e- w! ]9 y) K' A3 Q
x2=r2(i);/ D8 w1 Z8 z8 P) F
y=yueshu([x1 x2]);
4 ` j' W1 m+ v+ Y" K/ g: q6 M
if y==1 %当满足约束条件时# o8 B% R, I. c3 U# s
z=f([x1 x2]);
; M% [2 K4 }& S0 _0 ~- O+ i
if z>=z0 %求最大值; w/ N8 J* s( u& K% p9 f
z0=z;1 K2 L I, k1 F8 u
sol=[x1 x2]; %最值解
! [8 f& u" |% @6 T3 }4 m! _
end
& u6 D! p( J- _& x# i1 V7 F) O
end; n8 Q1 J* T' _& |3 k# e
end
" O9 U6 Q! g# ]: Z E
sol
/ @4 v* P3 C1 J# v
z0

复制代码