- 在线时间
- 90 小时
- 最后登录
- 2018-12-27
- 注册时间
- 2016-4-22
- 听众数
- 17
- 收听数
- 0
- 能力
- 20 分
- 体力
- 23472 点
- 威望
- 2 点
- 阅读权限
- 200
- 积分
- 7535
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 126
- 主题
- 100
- 精华
- 2
- 分享
- 0
- 好友
- 6
升级   50.7% TA的每日心情 | 开心 2018-6-4 15:01 |
|---|
签到天数: 7 天 [LV.3]偶尔看看II
 群组: 2018年大象老师国赛优 群组: 高考备战 群组: 2018中小学数学建模冬 |
将要规划求解的问题运用各种原理写成要最大化或者最小化的数学公式后,运用MATLAB求解问题。针对非线性规划,MATLAB提供了如下命令:fimincon.
: [5 v7 R0 Z9 D2 C8 O2 ~" BMATLAB中非线性规划的数学模型可写成如下形式:
4 _5 u6 Q) l& o4 R$ K 6 L- c) v( Z5 x3 }$ L# q
f(x)为标量函数,A,B,Aeq,Beq是相应维数的矩阵和向量,C(x),Ceq(x)是非线性向量函数
0 V5 q M+ P, w$ ]$ Z![]()
) V7 w- v+ v g如,求解下列非线性规划: # a; S' S! \4 }$ `3 }" B. h9 k8 o
! I4 k: S) d' F9 @6 e- N/ |
先编写funl.m定义目标函数:! `! d, Q& ?, ?# d; u1 y: K
function f = fun1(x);f = sum(x.^2) + 8;
3 G$ [: F4 p5 `. F* v: P! I文件fun2.m定义了非线性约束条件:
' @+ f# b/ v& r/ Ffunction [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 J( L2 d9 ?: p/ G主程序文件example2.m如下:
" F( {6 ?$ ?7 D2 Q3 t, G4 H4 b3 [2 zoptions = optimset('largescale','off');[x,y] = fmincon('fun1',rand(3,1),[],[],[],[],zeros(3,1),[],... 'fun2',options)' I: y& W; S7 H8 H; `, \. X+ ]
结果为:
7 J5 t2 V1 y3 r9 K# R* V) w2 Q6 Mx = 0.5522 1.2033 0.9478y = 10.6511
7 S4 p+ h& {5 T- L1 A2 v& K
2 T4 W Y- Z9 g; \4 }' c5 `& N1 P, E+ n/ K& `
: c! `! Z) w' t" c, F# `
. [: ^! L, v/ `' K, W
. o$ n: W! E' O1 j; ? |
zan
|