QQ登录

只需要一步,快速开始

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

[建模教程] 不看会后悔系列——国赛加分算法之遗传算法(下)

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

52

主题

12

听众

676

积分

  • TA的每日心情
    奋斗
    2021-6-27 15:42
  • 签到天数: 27 天

    [LV.4]偶尔看看III

    版主

    国际赛参赛者

  • TA的关系
  • 群组冬令营普通班

    群组Latex研学群

    群组2018美赛护航培训课程

    群组2018美赛冲刺培训

    群组2017科技论文写作

    跳转到指定楼层
    1#
    发表于 2018-8-9 00:04 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    本帖最后由 2336426014 于 2018-8-9 00:08 编辑
    $ J* {/ M* T0 p9 u0 a8 m" w/ B) x- a! U" ?' i0 b
            在每年的 国赛赛题中,几乎都会遇到多约束非线性规划问题 。 组建合理的多约束非线性规划数学模型不容易,如何精确求解这类同題也不容易 。 求解这类同题,目前的方法有很多种,比如 MATLAB优化工具箱甚至穷举法 。MATLAB优化工具箱求解此类难题精度不够;穷举法 效率太低。 目前的首选方法是 GA和1ingo和lindo软件。GA方法求解多约束非线性目标规划被广泛应用于各个领域。 例如,通用电气公司和 Rensselaer综合技术学院的组研究人员成功地将遗传算法用到一种商业客机使用的高函道比喷气发动机的涡轮的设计之中 。 这种涡轮是由多级静止和旋转的叶扇组藤安装在近似lll简形的願道里,是发动机开发计划的核心部分。 涡轮设计至少涉及1oo个变量 ,每个变量的取值范国各不相同,由此形成10^387以上个点的搜索空间。 涡轮设计方案的“适应度”取决于它满足一组限制的程度如何,这组限制有50个左右,如内壁和外壁的形状.函道内各点处燃气气流的压力、速度和扰动情况等。在一般情況下,一个工程师独立工作获得1个满意的设计结果大约要用几周的时间。 运行基于GA的发动机模拟软件和专家家系统有助于引导设计人员找出有意义的修改,工程师用这样的专家系统不到一天就可以完成一种设计。
    ( m8 o1 n2 C. R: S3 l) k' T
    2 Y% F- c. H  J0 u6 `* t2 E3 J        上次的例题及代码我们分析的是无约束目标函数最大值,而这次我们要做的是多约束非线性目标规划5 H$ N3 F4 i! G& s3 P

    " W  i* w; s( [        这次的操作比较简单,首先看一下目标函数及其三维图(我画的坐标系没有小于0部分):
    5 f: |! w& S& D3 U+ \7 i; p2 e 1.png
    8 Y0 c/ {* i+ Z  A, T 2.png
    8 Q+ w( z- ~6 ^+ J0 @/ M
    ; a. l( f& K0 W% E* [/ a        从这个角度看哈,首先可以看出它沿着某一个方向是递增的,考虑到而且要么同时正,同时负的时候函数可以取很大或者很小的值(st的限制看出来的,二条件同正下x1无穷大那函数无穷,同负下x1无穷小那就是负无穷,所以最小值应该就是图像所显示的,在以原点为中心,半径100的圆圈内)1 v* ?2 @, x, a$ o3 t4 L
    4 L0 ~" t) ]1 C, ]
           我们要找这个函数的极小值,也就是要找适应度最低的“染色体组”,所以接下来的事情,就是上个帖子的套路一样,编码,解码,适应度,迭代,新群体。
    # j: V8 _/ A) H, Y首先写出我们的新函数:
    0 @0 k1 z7 O- Z* v 3.png / B9 `2 `. b& f3 k0 K. ?
    其次是我们的主程序:* A' W: q1 Q7 ^2 C
    4.png : h8 s# T0 |: J6 ~& H
    ; N6 `- ~# ]" X# e) Y- R) E! G6 j
    程序结果呢就是最小值是0.0246,坐标位置为(-9.489,1.052),的确有那么点意思哈: f! o+ z: h# E4 d
    ) Q; n+ L5 \6 l# f% w6 v
           但是这远没有结束,我们的遗传算法貌似通过这种“爬山”的方式到达最低或者最高,但是我们只有说迭代次数越大,它是全局最优的概率才打,也就是说不管你迭代多少次,总是有几率证明它可能不是全局最优。比如因为迭代次数少陷入局部最优,或者爬山的能力差,所以呢,为了解决这些问题,我们就引入了遗传退火算法(当~当~当~当~~~~)。
    - A; R+ ]- I+ y2 C3 J8 O

    : B$ C, r9 u) e: C% A3 x操作流程如下:2 v' U0 @2 }8 l1 C% m
    5.png

    , q' ^- @9 r3 Q3 a8 X) y
    " \8 Z$ X9 \4 u考虑到这个需要懂退火算法,所以暂时不细说,只是给看一下流程图,与传统的遗传不同点在于:增加了对优秀个体的退火这一步骤,让优秀个体概率性的进行复制、变异。
    & A: L$ z) @7 ^# e1 H9 b' o7 a1 ~: \, t
           当然,我们论文里面就用传统的遗传算法其实也可以说足够了,为啥不怎么推荐退火呢,主要是 退货程序极其麻烦,占用内存,就怕一退就是一天,那岂不是凉凉,所以如果非专业研究,我们其实传统遗传退火就可以了,一个数无约束目标规划,另一个就是多约束目标规划。这样的话,国赛论文逼格高大上。拿奖分分钟!
    ' a0 G/ ]  g8 v( `9 [
    0 K2 x/ y7 @% l) U+ Y) f
    2 v3 m4 ~+ X. {% g; X9 f5 N5 W

    8 B5 ?$ K; M4 l' U) a& e. ~* d! W/ v- P9 {" W3 }

    多目标代码.rar

    698 Bytes, 下载次数: 19, 下载积分: 体力 -2 点

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏2 支持支持0 反对反对0 微信微信
    3963095 实名认证       

    1

    主题

    4

    听众

    88

    积分

    升级  87.37%

  • TA的每日心情
    开心
    2020-12-31 16:13
  • 签到天数: 20 天

    [LV.4]偶尔看看III

    邮箱绑定达人

    群组2016国赛备战群组

    66666666666666666666666666666666666666666666666666666666666666666666666666666
    2 e( T, U) A0 ?- J. R/ b; I! H4 U- z
    回复

    使用道具 举报

    0

    主题

    1

    听众

    2

    积分

    升级  40%

    该用户从未签到

    回复

    使用道具 举报

    ojbk        

    0

    主题

    3

    听众

    34

    积分

    升级  30.53%

  • TA的每日心情
    奋斗
    2018-11-25 01:10
  • 签到天数: 18 天

    [LV.4]偶尔看看III

    回复

    使用道具 举报

    fjdieb        

    0

    主题

    3

    听众

    4

    积分

    升级  80%

    该用户从未签到

    回复

    使用道具 举报

    1

    主题

    3

    听众

    38

    积分

    升级  34.74%

  • TA的每日心情
    奋斗
    2018-9-17 15:02
  • 签到天数: 12 天

    [LV.3]偶尔看看II

    群组2014国赛优秀论文解析

    回复

    使用道具 举报

    0

    主题

    3

    听众

    10

    积分

    升级  5.26%

  • TA的每日心情
    开心
    2018-8-9 16:02
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    哈哈哈哈哈哈哈

    邮箱绑定达人

    回复

    使用道具 举报

    0

    主题

    3

    听众

    18

    积分

    升级  13.68%

  • TA的每日心情

    2019-9-18 19:45
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    自我介绍
    一名热爱数学学生

    邮箱绑定达人

    回复

    使用道具 举报

    0

    主题

    3

    听众

    31

    积分

    升级  27.37%

  • TA的每日心情
    奋斗
    2019-1-24 13:00
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    群组全国大学生数学建模竞

    回复

    使用道具 举报

    0

    主题

    3

    听众

    31

    积分

    升级  27.37%

  • TA的每日心情
    奋斗
    2019-1-24 13:00
  • 签到天数: 8 天

    [LV.3]偶尔看看II

    群组全国大学生数学建模竞

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-6-11 09:41 , Processed in 0.451138 second(s), 106 queries .

    回顶部