数学建模社区-数学中国

标题: 复合形法求解约束最优化Fortran代码 [打印本页]
作者: 建不了的模。    时间: 2015-1-12 11:16
标题: 复合形法求解约束最优化Fortran代码
!* 说明: 1.本程序为复合形法                                                                          *3 K, B. o! q2 i& A5 A/ T5 \( A
!*                                                                                                        *
- g9 J; m& G& `7 K5 t9 \) F!*       2.程序功能是求解约束最优化问题                                                              *
- B1 V: L. R+ _. Z8 i!*   max F(x1,x2,…，xn)                                                               *
* S% F! ^0 e) Y6 X) a!*        s.t.        Gi≤xi≤Hi                                                              *0 D$ C+ P( t! g$ p9 r) i- N# z* \
!*      其中 x1,x2,…，xn为独立自变量，                                                              *
' I, H6 m8 t3 u2 c% Y8 t6 z!*       xn+1~xm为隐式变量，是x1,x2,…，xn的函数，Gi，Hi为下界和上界,它们可以是常数(显式约束),** c# V( W1 r' O, h
!*  也可以是自变量的函数(隐式约束)                                                            *1 B# j1 b5 N% V  |2 [3 J6 v7 Z- g
!* 主程序中提供自变量的初始值，输入已知参数及打印最后结果                                            *+ m" {* M6 A" P' r
!*                                                                                                        *; e% b- p) T2 `1 V! ^8 W# u- ?. T
!*       3.输入变量说明                                                                              *
: z7 s8 S4 J' u( ^. m9 V( s9 V!*  N--显式自变量数                                                                           *
" A( @, |' ~5 s$ b: w# `!*  M--约束组数                                                                               *0 b$ r1 A$ g  L8 r) x" N
!*  K--构成复合形的顶点数,常用N+1,可多取                                                      *, j. P% W, P' \% s- {9 s1 O- ~4 V& y9 p
!*  ITMAX--允许最多迭代次数                                                                   *
) C- Z9 ?6 c6 X& [* w!*  IPRINT--打印控制参数,IPRINT=1,打印中间结果,IPRINT=0,不打印中间结果                        *
4 W* e- K7 Z/ _' E7 n!*  ALPHA--反射因子,常用1.3                                                                   *
1 Q: {* v% B8 ~/ i! h!*  BETA--收敛参数,例如函数的幅值乘1E-4                                                       *
( v5 X; c4 v( b* ~& T!*  GAMMA--收敛参数,整数,常用值为5                                                            ** C1 k! b/ B/ v2 Z; {) h
!*  DELTA--显式约束违反校正,小正数,如X向量幅值乘1E-4                                          *
4 U  `0 |- I; T+ x!*  X(1,J)--自变量初始可行点,=1,N                                                             *. u9 }2 K* Q" C0 j- t
!*                                                                                                      *
# ?* w2 t* I, f/ ]6 P! [. ~7 ^1 M# r!*      4.输出变量说明                                                                               *9 c; O9 e4 s, L* C5 a' h
!*  F--目标函数最大值勤                                                                       *; T; ~4 b6 q8 D* O& `2 U+ x
!*  X(I)--自变量最优值,I=1,N                                                                  *
# K+ _3 b2 I# T/ S!*                                                                                                      *
' ]/ |2 o2 c- b4 B1 {# e; T!*      5.使用方法                                                                                   *6 E. I2 m1 Q. @1 h0 e) o
!*  1)用户按照待解问题修改主程序PARAMETER的N,M,K值                                            *5 y+ D! p5 j, L6 D" A
!*  2)在子程序FUNC中给定目标函数                                                              *
3 ]- S. A" B. I0 |3 F) ]!*  3)在子程序CONST中给定Hi和Gi,显式约束必须放在隐式约束前面                                  *3 h+ M: b. c! p+ k/ s/ v
!*  4)编写一个数据文件COMPDAT,文件中的数据依次为下列变量的值,每个数之间用逗号分开             *   
+ b% b; @- ^5 v% f9 _' E!*      ITMAX,IPRINT,ALPHA,BETA,GAMMA,DELTA,X(1,1),X(1,2),…，X(1,N)                          *& N7 r+ h% `) U, T
!*      6.程序内容                                                                                   ** w2 B( ?7 G( F8 I5 f
!*  1)本程序由一个主程序和六个子程序组成，主程序首先给定N,M,K值，然后定义数组维数             *. k/ |$ M( v' F; @! ~: D- ]
!*   PARAMETER(N=,M=,K=)                                                               *
' G/ U/ I3 w! ^) A) O; U!*   DIMENSION X(K,M),R(K,N),F(K),G(M),H(M),XC(N)                                      *( U3 O) ]1 W% K( K% R4 {
!*  主程序中提供自变量的初始值，输入已知参数及打印最后结果面                                  *# ^  X9 I9 J2 @" ~' D6 O( @
!*  2)各子程序的作用分别为了                                                                  *   ' D( D& E. q% z$ b+ M3 J
!*      CONSX—这是一个主要子程序，调用其它于程序及输出中间结果                               *, J$ P$ c3 ^: l. p
!*      CHECK—检查所有的点是否满足约束条件,对违背约束的点进行校正                            *
  K$ H7 [# x) n& h6 W0 O$ U!*      CENTR—计算中心点                                                                     *, ?( ^- }8 h1 m7 d; P
!*      FUNC —目标函数，由用户提供                                                           *
9 U4 p2 I- O* `! j!*      CONST—规定显式和隐式约束                                                             *
5 @. k9 `5 S5 P# F& T  U!*      RANDU—产生随机数                                                                     *5 ~: @1 v+ t! p7 o
!*************************************************************************************************
) |6 y" J  K- X+ q/ t, b8 d8 U& n) e& C
* U8 {; H  C" B; c% N3 z% i) U0 W

8 P5 J4 d9 h/ H" m



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5