QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3697|回复: 0
打印 上一主题 下一主题

复合形法求解约束最优化Fortran代码

[复制链接]
字体大小: 正常 放大

2620

主题

162

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2015-3-12 15:35
  • 签到天数: 207 天

    [LV.7]常住居民III

    社区QQ达人 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组第六届国赛赛前冲刺培

    群组国赛讨论

    群组2014美赛讨论

    群组2014研究生数学建模竞

    群组数学中国试看培训视频

    跳转到指定楼层
    1#
    发表于 2015-1-12 11:16 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    !* 说明: 1.本程序为复合形法                                                                          *
    ' Y5 L0 [( j3 ~!*                                                                                                        *  Z) e4 E* k, J& N9 P1 a
    !*       2.程序功能是求解约束最优化问题                                                              *8 A( y1 z4 {$ H, |8 G
    !*   max F(x1,x2,…,xn)                                                               */ c) U% l" p1 Z$ _- u
    !*        s.t.        Gi≤xi≤Hi                                                              *
    ; p: z, P4 D' E- T!*      其中 x1,x2,…,xn为独立自变量,                                                              *9 G' l8 c& |6 V1 Y: _% A
    !*       xn+1~xm为隐式变量,是x1,x2,…,xn的函数,Gi,Hi为下界和上界,它们可以是常数(显式约束),*# R4 v6 V. I* K6 W. D% b& J4 J- a+ Y
    !*  也可以是自变量的函数(隐式约束)                                                            *+ {  G3 b( S2 O7 \
    !* 主程序中提供自变量的初始值,输入已知参数及打印最后结果                                            *' Q1 V  W' J% t
    !*                                                                                                        *
    , _! K, [9 K( I9 H. z6 P!*       3.输入变量说明                                                                              *
    - O# a, v7 K) ?) y!*  N--显式自变量数                                                                           *5 S" ?4 X$ I( P0 ?7 k0 c
    !*  M--约束组数                                                                               *8 }) l% G$ S, W' h; c7 s
    !*  K--构成复合形的顶点数,常用N+1,可多取                                                      *& [, d" |! g& l+ K
    !*  ITMAX--允许最多迭代次数                                                                   *! U8 u" g  h# }. y: ^: Z. h7 f
    !*  IPRINT--打印控制参数,IPRINT=1,打印中间结果,IPRINT=0,不打印中间结果                        *
    3 d8 e  C# L$ e+ k% r: P* s!*  ALPHA--反射因子,常用1.3                                                                   *
    # _/ S: ]" e, ~3 L# P; _0 q!*  BETA--收敛参数,例如函数的幅值乘1E-4                                                       *: {3 o) o$ u. }' X+ k1 A( C0 v
    !*  GAMMA--收敛参数,整数,常用值为5                                                            *5 V# b  f) v, b" }* _, m
    !*  DELTA--显式约束违反校正,小正数,如X向量幅值乘1E-4                                          *
    6 n& ~$ _. D7 X( @% n/ S!*  X(1,J)--自变量初始可行点,=1,N                                                             *, _& E/ P: J8 w. Z3 s  O4 {+ N
    !*                                                                                                      *
    1 M" x- d+ U+ e7 F!*      4.输出变量说明                                                                               *$ v! V# F) J7 t4 t: |
    !*  F--目标函数最大值勤                                                                       *
    $ [, F3 g! a$ z4 x!*  X(I)--自变量最优值,I=1,N                                                                  *
    9 v% H) g; f0 B!*                                                                                                      *  u% E  q9 ^5 i$ m' Z8 Z3 Z
    !*      5.使用方法                                                                                   *
    , d( }. \/ K0 k- U/ ^6 s!*  1)用户按照待解问题修改主程序PARAMETER的N,M,K值                                            *  M0 I, H. \* k1 E8 V0 s' Z) ~1 ?
    !*  2)在子程序FUNC中给定目标函数                                                              *( c% x1 U# ^# f- S' O! f
    !*  3)在子程序CONST中给定Hi和Gi,显式约束必须放在隐式约束前面                                  *) i0 X0 }1 v) b; @" h9 x( F
    !*  4)编写一个数据文件COMPDAT,文件中的数据依次为下列变量的值,每个数之间用逗号分开             *   
    1 R: P0 E  B7 j!*      ITMAX,IPRINT,ALPHA,BETA,GAMMA,DELTA,X(1,1),X(1,2),…,X(1,N)                          *: `( i; p% V2 e; Q  S
    !*      6.程序内容                                                                                   *
    - v8 x0 I9 d+ |$ C8 t* L!*  1)本程序由一个主程序和六个子程序组成,主程序首先给定N,M,K值,然后定义数组维数             *
    ; O/ y- g- W! t1 J* m$ i!*   PARAMETER(N=,M=,K=)                                                               *5 t  K7 c& I4 e/ R; S2 Y
    !*   DIMENSION X(K,M),R(K,N),F(K),G(M),H(M),XC(N)                                      *
    0 v* N" f; T3 Q2 h!*  主程序中提供自变量的初始值,输入已知参数及打印最后结果面                                  *5 w" x2 A) e: i+ ^  f
    !*  2)各子程序的作用分别为了                                                                  *   # _" _, {5 |2 _% v
    !*      CONSX—这是一个主要子程序,调用其它于程序及输出中间结果                               *- Y/ L; {$ v1 z+ D: \
    !*      CHECK—检查所有的点是否满足约束条件,对违背约束的点进行校正                            *" c  [! r$ ]2 M- \9 b4 F
    !*      CENTR—计算中心点                                                                     *
    ; M0 V% J4 K# Q1 |- }! ^1 J# t4 o!*      FUNC —目标函数,由用户提供                                                           *, g- T; f/ p0 _  @3 G0 Y
    !*      CONST—规定显式和隐式约束                                                             *# F: _  o/ e- k/ G+ g2 ^5 C5 R
    !*      RANDU—产生随机数                                                                     *5 {8 k0 n8 R" m( a" f1 ~
    !*************************************************************************************************
    3 A) O  p1 V/ O1 ?  ^
    ! z- y$ M6 v& Q' }! H
    游客,如果您要查看本帖隐藏内容请回复

    4 K: a3 K1 i- I* i; z
    . |$ p) n1 K  Z
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-4 08:42 , Processed in 0.408237 second(s), 54 queries .

    回顶部