数学建模社区-数学中国

标题: 复合形法求解约束最优化Fortran代码 [打印本页]
作者: 建不了的模。    时间: 2015-1-12 11:16
标题: 复合形法求解约束最优化Fortran代码
!* 说明: 1.本程序为复合形法                                                                          *2 ~& U$ n; ?1 b+ r
!*                                                                                                        *
& K  H6 N7 G2 ]!*       2.程序功能是求解约束最优化问题                                                              *) c8 n' p3 W: I. h, {+ k
!*   max F(x1,x2,…，xn)                                                               *2 m. J7 T' c5 r8 m1 _2 t, N  t
!*        s.t.        Gi≤xi≤Hi                                                              *) h' o- h$ j( ~3 l- @
!*      其中 x1,x2,…，xn为独立自变量，                                                              *
- l, n" _) a/ R& P" ~7 X  {& i1 t!*       xn+1~xm为隐式变量，是x1,x2,…，xn的函数，Gi，Hi为下界和上界,它们可以是常数(显式约束),*4 f' E, M' }3 W3 w8 t. j: Z9 L
!*  也可以是自变量的函数(隐式约束)                                                            *
( ^5 }3 i# w! m- e( G9 b) X6 b" h!* 主程序中提供自变量的初始值，输入已知参数及打印最后结果                                            *
( e* |( G9 d& j$ Y( {. N" c!*                                                                                                        *9 ]6 N/ k: Y0 P, ~
!*       3.输入变量说明                                                                              *' l. g# H( w3 L0 {* c
!*  N--显式自变量数                                                                           *
& g7 P3 y9 U: B!*  M--约束组数                                                                               *
! O9 w. o, b- E( y  \!*  K--构成复合形的顶点数,常用N+1,可多取                                                      *
+ I- B; V$ z9 o2 S!*  ITMAX--允许最多迭代次数                                                                   *: A, X/ h' J2 c2 p; u% r, N
!*  IPRINT--打印控制参数,IPRINT=1,打印中间结果,IPRINT=0,不打印中间结果                        *
' c+ j  T' }4 Q. R/ g/ Z' k% ?!*  ALPHA--反射因子,常用1.3                                                                   *; {. ~# w) I3 h! W3 @% g. ~
!*  BETA--收敛参数,例如函数的幅值乘1E-4                                                       *# T, s# u" h% ?4 q* z- _
!*  GAMMA--收敛参数,整数,常用值为5                                                            *" }: b4 P% s3 F  _: L1 u; p
!*  DELTA--显式约束违反校正,小正数,如X向量幅值乘1E-4                                          *+ u+ V/ V. V0 ?5 c8 K% p
!*  X(1,J)--自变量初始可行点,=1,N                                                             *
* s( `& L7 q  `!*                                                                                                      *0 S5 M. s: j8 |. r8 |: \1 E& w  @& {6 b
!*      4.输出变量说明                                                                               *( w% G% c6 z  U# z( E9 `! }* D# z
!*  F--目标函数最大值勤                                                                       *
" ~' V- E/ K( q. o9 z* W!*  X(I)--自变量最优值,I=1,N                                                                  *
/ W% e* [% Q2 s3 O* h; G+ k% J!*                                                                                                      *, e4 e3 V2 s- j/ C
!*      5.使用方法                                                                                   *+ v. x# [' B/ v0 Q* M4 u
!*  1)用户按照待解问题修改主程序PARAMETER的N,M,K值                                            *. I. _8 z. A6 D% s
!*  2)在子程序FUNC中给定目标函数                                                              *8 o" X$ H' Y( {  Y) p
!*  3)在子程序CONST中给定Hi和Gi,显式约束必须放在隐式约束前面                                  *
+ n5 s5 L/ t- O0 e5 \! n!*  4)编写一个数据文件COMPDAT,文件中的数据依次为下列变量的值,每个数之间用逗号分开             *   , }8 b$ M- P7 D
!*      ITMAX,IPRINT,ALPHA,BETA,GAMMA,DELTA,X(1,1),X(1,2),…，X(1,N)                          *& G, X* a# J! _; r5 a
!*      6.程序内容                                                                                   *
  L: d1 d, w8 z9 ?8 B6 p. n!*  1)本程序由一个主程序和六个子程序组成，主程序首先给定N,M,K值，然后定义数组维数             ** O, w/ B. `/ x# I, s
!*   PARAMETER(N=,M=,K=)                                                               *8 A! ?; f, u8 i+ A( F
!*   DIMENSION X(K,M),R(K,N),F(K),G(M),H(M),XC(N)                                      *  Y% y' X' d; Y
!*  主程序中提供自变量的初始值，输入已知参数及打印最后结果面                                  */ I" i: r3 `4 ~9 w8 t9 @
!*  2)各子程序的作用分别为了                                                                  *   
$ d. @. ^$ {5 r0 z$ t' H, F/ \. n) H!*      CONSX—这是一个主要子程序，调用其它于程序及输出中间结果                               *
$ ]+ c+ u1 P, @) V! p!*      CHECK—检查所有的点是否满足约束条件,对违背约束的点进行校正                            *  \2 P% X2 w, P5 i; q
!*      CENTR—计算中心点                                                                     *! m0 o7 ]1 |' v( a: ]$ w" Q; W4 t
!*      FUNC —目标函数，由用户提供                                                           *; q7 ~' b; j/ Q( A$ H
!*      CONST—规定显式和隐式约束                                                             *
6 i' O, p( ~+ [* [; }* B( o$ e!*      RANDU—产生随机数                                                                     *$ x% r/ _9 @2 l5 S0 z
!*************************************************************************************************6 }  T' S1 p! W* s0 ?

" s: q9 z' D. Y7 L. Y
6 W: }+ f7 e& Z3 n& a! S2 Y. Y5 T6 X: y! J$ R. N- X




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