- 在线时间
- 1084 小时
- 最后登录
- 2015-9-10
- 注册时间
- 2014-4-18
- 听众数
- 162
- 收听数
- 1
- 能力
- 10 分
- 体力
- 43953 点
- 威望
- 6 点
- 阅读权限
- 255
- 积分
- 15243
- 相册
- 0
- 日志
- 0
- 记录
- 1
- 帖子
- 3471
- 主题
- 2620
- 精华
- 1
- 分享
- 0
- 好友
- 513
升级 0% TA的每日心情 | 开心 2015-3-12 15:35 |
---|
签到天数: 207 天 [LV.7]常住居民III
群组: 第六届国赛赛前冲刺培 群组: 国赛讨论 群组: 2014美赛讨论 群组: 2014研究生数学建模竞 群组: 数学中国试看培训视频 |
!* 说明: 1.本程序为复合形法 *
+ g9 I4 C$ Y8 W4 ~1 q" u!* *# h+ ~* J% `( W5 D
!* 2.程序功能是求解约束最优化问题 *, U' k. u9 _$ K ]+ A& @8 P
!* max F(x1,x2,…,xn) *
# W! s/ N; }- R i) u# j7 k- v% ~!* s.t. Gi≤xi≤Hi *3 l7 o' Y- a" C5 |: K+ I+ n/ o
!* 其中 x1,x2,…,xn为独立自变量, *
1 U f* _+ U n!* xn+1~xm为隐式变量,是x1,x2,…,xn的函数,Gi,Hi为下界和上界,它们可以是常数(显式约束),*0 m' _4 m1 Z" y, b
!* 也可以是自变量的函数(隐式约束) *+ Y0 M/ D9 t, c/ N! Y
!* 主程序中提供自变量的初始值,输入已知参数及打印最后结果 *
; |. S( a5 W; `5 v+ q!* *+ d& i) N7 v) H/ _& K
!* 3.输入变量说明 *; o; e2 m, H. @/ y
!* N--显式自变量数 *2 D. j6 Y5 P9 o; u/ E
!* M--约束组数 *# H- i1 l: o3 `" H( A3 \* Q% U6 n
!* K--构成复合形的顶点数,常用N+1,可多取 *" U# j* s- y$ ^1 }
!* ITMAX--允许最多迭代次数 *
# n2 x4 ?. ]/ Y# z' I2 j$ i!* IPRINT--打印控制参数,IPRINT=1,打印中间结果,IPRINT=0,不打印中间结果 *& i" a7 v! m3 n# W. j; ~* T# e
!* ALPHA--反射因子,常用1.3 *
# E' i% Y2 I; @6 r* a# ]0 P!* BETA--收敛参数,例如函数的幅值乘1E-4 *2 U0 s: [9 i, \; i
!* GAMMA--收敛参数,整数,常用值为5 *) z a* V8 ]5 q. j5 |# r
!* DELTA--显式约束违反校正,小正数,如X向量幅值乘1E-4 *) ^! U2 J; G& y+ k
!* X(1,J)--自变量初始可行点,=1,N *
( C% v4 W% a: ^) s: J!* *
; Y# Y7 p8 O; Z1 h/ k" ~/ @$ Z!* 4.输出变量说明 *
$ _% ?) ?% `& j) R" V) w!* F--目标函数最大值勤 *( Z0 I- d1 K4 G/ T2 K0 V n
!* X(I)--自变量最优值,I=1,N *# ~ U, |9 W) t& B0 ~
!* *2 c# q, h; i4 a) Q% X5 u
!* 5.使用方法 *
: K% J" V# w* ~1 J!* 1)用户按照待解问题修改主程序PARAMETER的N,M,K值 ** A8 u+ X4 \5 @. }( m
!* 2)在子程序FUNC中给定目标函数 *
, D3 P# j" k8 Z. Y!* 3)在子程序CONST中给定Hi和Gi,显式约束必须放在隐式约束前面 *
1 ]. T' Y8 Y' e6 `7 j2 ?. r* H!* 4)编写一个数据文件COMPDAT,文件中的数据依次为下列变量的值,每个数之间用逗号分开 *
2 ]+ b% R- A. A+ x k!* ITMAX,IPRINT,ALPHA,BETA,GAMMA,DELTA,X(1,1),X(1,2),…,X(1,N) */ Z- i7 X: L2 z$ I5 Y
!* 6.程序内容 *( U2 ~3 E0 K- L# Q: ^
!* 1)本程序由一个主程序和六个子程序组成,主程序首先给定N,M,K值,然后定义数组维数 *
5 E) G1 U+ f* w!* PARAMETER(N=,M=,K=) *7 J- ]$ K; E9 w a) i3 Z
!* DIMENSION X(K,M),R(K,N),F(K),G(M),H(M),XC(N) *
- T) p ~8 m: m!* 主程序中提供自变量的初始值,输入已知参数及打印最后结果面 *) S) e( _. c# N# I
!* 2)各子程序的作用分别为了 *
1 ^3 u$ T% R6 n+ B" v/ J7 f! I!* CONSX—这是一个主要子程序,调用其它于程序及输出中间结果 *; ~' L& I$ J* s' C _0 g1 P
!* CHECK—检查所有的点是否满足约束条件,对违背约束的点进行校正 *
0 t) @/ I F9 b; l: K0 F( h!* CENTR—计算中心点 * V9 X) a3 i; |7 E6 D
!* FUNC —目标函数,由用户提供 *
?7 D | J7 v' ^!* CONST—规定显式和隐式约束 *
- \* G& g" ~. C' g0 k6 _' Y!* RANDU—产生随机数 *, u1 {/ M: E. ?
!*************************************************************************************************
+ }' Z; t8 O* v; M* t
- B1 \- d9 [# \ k. |& u* `4 [8 B) v) K. U0 m/ D
: t% E( h f) H |
zan
|