数学建模社区-数学中国
标题:
MATLAB数学建模(3)-非线性规划
[打印本页]
作者:
佛自业障
时间:
2018-10-31 09:12
标题:
MATLAB数学建模(3)-非线性规划
将要规划求解的问题运用各种原理写成要最大化或者最小化的数学公式后,运用MATLAB求解问题。针对非线性规划,MATLAB提供了如下命令:fimincon.
2 T# }# p6 s: z- n" E& u4 p
MATLAB中非线性规划的数学模型可写成如下形式:
C. f5 x, S# q/ b
, L5 s% M+ q# b- v+ H
f(x)为标量函数,A,B,Aeq,Beq是相应维数的矩阵和向量,C(x),Ceq(x)是非线性向量函数
7 a6 F9 |# P( @1 v( Y
. |- a+ s3 @' v4 I$ v% |
如,求解下列非线性规划:
( k( J# y9 I M% y( Q
9 g% i& X9 I% }* \) G4 k
先编写funl.m定义目标函数:
6 q) j# b( q- J: q
function f = fun1(x);
f =
sum
(
x.
^
2
) +
8
;
: ^* v; T: R6 l* B
文件fun2.m定义了非线性约束条件:
# ^) K9 b( N0 O% _. b
function [g,h] = fun2(x);g = [(-x(
1
)
^2
+ x(
2
) - x(
3
)
^2
);(x(
1
) +x(
2
)
^2
+ x(
3
)
^2
-
3
)];h = [(-x(
1
) - x(
2
)
^2
+
2
);(x(
2
) +
2
*x(
3
)
^2
-
3
)];
1 ]4 M5 N- W x3 [( y, O) c
主程序文件example2.m如下:
4 N) C8 a& S4 ~8 y
options = optimset(
'largescale'
,
'off'
);[x,y] = fmincon(
'fun1'
,rand(
3
,
1
),[],[],[],[],zeros(
3
,
1
),[],
...
'fun2'
,options)
$ Y9 Y( O6 ?( B# F8 o ^
结果为:
0 U/ Y* c" g+ @: J- Q7 j3 d! z
x
= 0.5522 1.2033 0.9478
y
= 10.6511
, U+ b$ @: L, l
9 T' d; e m6 s \
6 C% p$ k4 @1 \" e$ X$ y4 n, o# j
. v3 [2 r" c/ c' j" {
% _1 E+ X# W6 u
v* ?4 G* v3 a# r8 v4 H
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5