- 在线时间
- 1084 小时
- 最后登录
- 2015-9-10
- 注册时间
- 2014-4-18
- 听众数
- 162
- 收听数
- 1
- 能力
- 10 分
- 体力
- 43980 点
- 威望
- 6 点
- 阅读权限
- 255
- 积分
- 15251
- 相册
- 0
- 日志
- 0
- 记录
- 1
- 帖子
- 3471
- 主题
- 2620
- 精华
- 1
- 分享
- 0
- 好友
- 513
升级   0% TA的每日心情 | 开心 2015-3-12 15:35 |
|---|
签到天数: 207 天 [LV.7]常住居民III
 群组: 第六届国赛赛前冲刺培 群组: 国赛讨论 群组: 2014美赛讨论 群组: 2014研究生数学建模竞 群组: 数学中国试看培训视频 |
!* 说明: 1.本程序为复合形法 *; m! t9 }; }7 f# B
!* *3 C8 O% M% Y, ?: n: N: k
!* 2.程序功能是求解约束最优化问题 *2 O5 H* A. i& |) P6 L' N- {; r
!* max F(x1,x2,…,xn) *
( n5 i6 x9 v+ D; S6 G. i!* s.t. Gi≤xi≤Hi *$ m. y+ w2 o) q
!* 其中 x1,x2,…,xn为独立自变量, *
/ p1 H" h" o3 E5 K! ]!* xn+1~xm为隐式变量,是x1,x2,…,xn的函数,Gi,Hi为下界和上界,它们可以是常数(显式约束),*
3 a/ n$ c8 ~9 _: X% o!* 也可以是自变量的函数(隐式约束) *
7 a& r H+ b3 y" [) S!* 主程序中提供自变量的初始值,输入已知参数及打印最后结果 *- y5 w4 t. D: M+ x" E6 s
!* *7 b2 \5 _; }. R, W
!* 3.输入变量说明 *
" ?9 z) v* K1 w1 D/ Q0 T! q!* N--显式自变量数 *% g0 g# g- r# u# V) G( L. @1 L
!* M--约束组数 *
, _ \9 e0 B6 ?; X t" k& J- Q- G!* K--构成复合形的顶点数,常用N+1,可多取 *( W Y! q. K, y
!* ITMAX--允许最多迭代次数 *
& J. F- O6 Y, R!* IPRINT--打印控制参数,IPRINT=1,打印中间结果,IPRINT=0,不打印中间结果 *& W; o+ }+ l; J% U
!* ALPHA--反射因子,常用1.3 *
. n* S8 X) v" U- G!* BETA--收敛参数,例如函数的幅值乘1E-4 *
* D" a) i: S9 P6 N4 ^!* GAMMA--收敛参数,整数,常用值为5 *: k; R9 ?1 K& \: Z( |0 y2 |5 J
!* DELTA--显式约束违反校正,小正数,如X向量幅值乘1E-4 *
% T8 o5 [# z! P0 d!* X(1,J)--自变量初始可行点,=1,N *, w' B! r/ `9 S4 z( q
!* *
5 R/ v i7 l( w!* 4.输出变量说明 *9 J4 `' e( m& q" C0 A, F3 c: t% |
!* F--目标函数最大值勤 *4 A$ k% ~% ?# @5 c" \
!* X(I)--自变量最优值,I=1,N *% o) p+ {. d6 K. [& p }6 F6 @
!* *
* N( M' J" H+ W7 k!* 5.使用方法 *! U: K! R" Z/ n$ }; X% F9 d1 d$ X
!* 1)用户按照待解问题修改主程序PARAMETER的N,M,K值 *
. y* ~4 O. Y) Y, k$ H) L!* 2)在子程序FUNC中给定目标函数 *
! O1 u- |! Z+ W3 \5 I+ n& }!* 3)在子程序CONST中给定Hi和Gi,显式约束必须放在隐式约束前面 *& F; t5 L* e8 I, L4 Y5 h# ~% ]
!* 4)编写一个数据文件COMPDAT,文件中的数据依次为下列变量的值,每个数之间用逗号分开 *
7 q9 ^# I8 u/ \" k5 u$ {- R5 b0 q. P!* ITMAX,IPRINT,ALPHA,BETA,GAMMA,DELTA,X(1,1),X(1,2),…,X(1,N) *
* a. S# M v: [3 f. x!* 6.程序内容 *
5 z# G8 N( T. {!* 1)本程序由一个主程序和六个子程序组成,主程序首先给定N,M,K值,然后定义数组维数 *
T) c C% i4 ]( a$ }!* PARAMETER(N=,M=,K=) *
; p B* N+ c# N1 P!* DIMENSION X(K,M),R(K,N),F(K),G(M),H(M),XC(N) *
0 t4 X) O# t1 |. Z!* 主程序中提供自变量的初始值,输入已知参数及打印最后结果面 *8 X# ]& X2 e2 g2 ?- `) P& B8 x
!* 2)各子程序的作用分别为了 *
* ]& i5 E; m9 }0 l) J2 k!* CONSX—这是一个主要子程序,调用其它于程序及输出中间结果 *
5 e+ A2 ^/ [1 Z& W* \' k!* CHECK—检查所有的点是否满足约束条件,对违背约束的点进行校正 */ V6 ~, x* o9 ^+ l
!* CENTR—计算中心点 *# M0 z( C& ], m f# K' F7 Z
!* FUNC —目标函数,由用户提供 *
! [% @9 m q3 Z+ J9 w!* CONST—规定显式和隐式约束 *: @4 ]+ x- \0 n( K' j c" h
!* RANDU—产生随机数 *) c- c+ A# k1 {
!************************************************************************************************** H. x* }- a6 P# W' s# y
8 S) H" Y$ q( Y) @
: ?* h/ H& P. w4 U1 b6 S Z2 S
5 i3 d" Q# ~/ [& H
|
zan
|