QQ登录

只需要一步,快速开始

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

华为杯研赛数学建模之遗传算法

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

1178

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2023-7-31 10:17
  • 签到天数: 198 天

    [LV.7]常住居民III

    自我介绍
    数学中国浅夏
    跳转到指定楼层
    1#
    发表于 2021-10-13 18:18 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    遗传算法简介8 x5 g) p9 D7 f. b  G( T
    维基百科上的概念如下$ M5 @8 A1 n1 P
    遗传算法(英语:genetic algorithm (GA) )是计算数学中用于解决最优化的搜索算法,是进化算法的一种。进化算法最初是借鉴了进化生物学中的一些现象而发展起来的,这些现象包括遗传、突变、自然选择以及杂交等。
    * K# f; c2 }, W( S5 t; \. @7 B& m) \
    遗传算法通常实现方式为一种计算机模拟。对于一个最优化问题,一定数量的候选解(称为个体)可抽象表示为染色体,使种群向更好的解进化。传统上,解用二进制表示(即0和1的串),但也可以用其他表示方法。进化从完全随机个体的种群开始,之后一代一代发生。在每一代中评价整个种群的适应度,从当前种群中随机地选择多个个体(基于它们的适应度),通过自然选择和突变产生新的生命种群,该种群在算法的下一次迭代中成为当前种群。
    7 _5 y% v1 E' s( A! I2 k: C, y5 j" X
    3 |8 x2 _2 u, K6 M7 E, ^* `这看起来是不是很抽象?没关系,我们举一个例子  k2 ~& w8 ^& \* Z+ f; A* |
    比如你家住的那一块地方有一群老鼠,有一些跑得快而且也聪明,有一些傻得大白天在你面前跳舞,然后理所当然就被你捕获然后进行土葬。故那群跑得快而且也聪明自然活下来的概率比傻的大的多,我们称老鼠是否活下来为适应度,适应度高的更容易在你手上活下来,而低的更容易被你土葬。活下来的老鼠们,它们中大多都是适应度高的,它们进行交配繁殖,老鼠子类大多都是继承了老鼠父母的特点,当然也存在一定几率变异,变异数除以总体老鼠为变异率,变异率一般都不是很大的常数。但总体来说,由于淘汰机制,故所产生的鼠类比起它们的父类,平均来说更为聪明。通过一代代的优化,使得老鼠越来越聪明,而你越来越抓不到老鼠。这就是遗传算法的核心。/ w# t. `8 q0 L; d; c2 c
    9 l- v2 z8 k7 A
    遗传算法的步骤6 ]1 t; L/ M  W* ?/ L
    1.编码
    + |! m) R/ x3 U; t) l编码有多种,如二进制编码,浮点数编码,格雷码等。最常见的是二进制编码和浮点数编码。7 O+ S6 W9 c6 V
    , q/ @6 r, U) R# u/ {% ^2 p
    二进制编码' _' M% W" E3 Y3 j3 m+ C# J) x
    二进制编码方法是使用二值符号集{0,1},它所构成的个体基因型是一个二进制编码符号串。二进制编码符号串的长度与问题所要求的求解精度有关。
    1 Q; N% c+ ^0 |; W' S优点:, r5 s7 \, Z  I
    1 {3 x- u8 c" u" k; y
    简单。无论是编码还是解码操作都非常方便和快捷。4 S" y  H4 f# N+ J7 B) L3 L: \
    方便交叉和变异。
    8 M+ p9 {) k0 z9 `* S* `- f! n符合最小字符集编码原则。
    1 O# O. W9 r4 e4 d% {- d缺点:
    # i) G) `+ Y( m- G
    - ?( P. a, P8 @9 V不适合连续函数的优化问题,局部搜索能力差。
    4 z" Y  e1 f& `# f连续的数值之间有时候存在距离大的问题。例如63和64对应的二进制分别是0111111和1000000。(连续数值对应的二进制数7位全都不同)
    0 J# k6 G- A3 m$ ~8 S- A& V对于高精度的问题,变异后可能会出现远离最优解的情况,表现型不稳定。/ p* s* c: M: B( H3 n
    案例:( `; G7 _# g; }
    假设有f(x),x∈[0,1023],采用定长二进制编码,串0010101111就代表了175,编码精度为1.3 Q: C! U4 h, X' D: K( K9 O% l; n# `

    ) d3 S: s+ b  t5 c$ U浮点数编码" b8 K; @2 W1 E3 z7 z7 ]2 n
    定义:
    : s4 |- C' S9 I0 J个体基因值用某范围内的一个实数来表示。编码长度等于决策变量的个数。
    5 L( I9 L( |( K优点:& P' _1 s) `7 Z; U+ H
    9 u/ k1 L4 S7 `# x
    精度高,适用于连续变量问题。避免了海明悬崖问题。& d1 [8 i8 O. i, n5 x4 f% I- C
    适用于表示范围比较大的数值,适合空间较大的一串算搜索
    + t7 E' m2 e- J! }, u降低了计算复杂性,提升效率
    7 ?  d* u7 a8 |3 b/ y6 f9 X便于遗传算法与经典优化方法的混合使用: c2 D2 J4 M5 f/ l! o; ^
    便于设计针对问题的专门知识的知识型遗传算子# x2 t# ~& k5 b3 Z+ K
    便于处理复杂的决策变量约束条件,适合于组合优化问题
    5 l6 `+ @! p9 Q) K0 X. d8 N案例:
    7 _6 P* L9 X2 C' C: x. Y假设某优化问题有五个变量,每个变量的变化范围都不同。其中X={5.30,5.20,4.70,3.40,4.80}就是一个基因型,对应的表现型是x={5.30,5.20,4.70,3.40,4.80}。+ D7 i# H8 V4 e& t  F: k1 p! Z" {) K7 a

    $ d3 f0 b+ u/ A" B# C2.解码8 G. w2 X, y' _7 G
    以二进制为例。
    3 w- v. ?' m& A" U8 k. ^解码的目的就是将不直观的二进制数据串还原成十进制。
    1 E5 u  J0 Z& `0 T  w6 ]3 W! O5 i. d& G2 D
    3.交配(交叉)
      U6 \# n+ o6 E6 e1 O1 {以二进制为例。
    5 B# j4 b3 p, a“交配运算”是使用单点或多点进行交叉的算子。首先用随机数产生一个或多个交配点位置,然后两个个体在交配点位置互换部分基因码,形成两个子个体。例如,两条染色体S1=01001011,S2=10010101。交换其后4位基因,得S1’=01000101,S2’=10011011可以被看做原染色体S1和S2的子代染色体。
    . G& s$ n$ A" M8 M* S
    7 p# z+ q' j, t# E+ n# ~4 _4.变异! u% i6 a9 Y" D7 Z! G$ W
    突变- V- N. N2 u& Q  Z: O! `- L
    突变是指基因突变。例如对于S1=010110011,第三位0突变成1,那么我们得到S1‘=011110011。9 @: s1 N7 Z* O3 n5 P
    倒位
    , @. k. y' `$ a/ H$ U5 U倒位是指一个染色体某区段正常排列顺序发生180°的颠倒,造成染色体内的DNA序列重新排列。例如对于S1=1010100010110110101001进行倒位时得到S1‘=010101101000110101001。
    3 Z; C  y1 A! b其他
    % c' c6 `0 S+ I' A* q' T0 d不一一举例,大家自行查阅理解。  m3 q1 ~% e5 D* I0 S2 j( A

    3 Z4 u$ b" c; s; b% w3 T! \- y: |" W6.适应度评估
    # i- _* C1 ~$ q( L6 H- t遗传算法依照与个体适应度成正比的几率决定当前种群中各个个体遗传到下一代群体的机会。个体适应度大的个体更容易被遗传到下一代。通常情况下,求目标函数最大值的问题可以直接把目标函数作为检测个体适应度大小的函数
    ' i7 S) g7 T( f
    8 M. Y" [! T9 {/ h7.选择' }$ t9 b6 R* j8 X
    选择运算时根据个体适应度大小决定其下代遗传的可能性。设种群中个体总数为N,个体i的适应度为fi,则个体i被选取的几率为:
    $ q: A! v) a) _! L* |$ @/ ^+ j/ i, d7 G6 s* F0 E" I

    5 E+ o' i) d9 M遗传算法伪代码3 S. E& o( @* I% W
    BEGIN+ L$ d) h/ E, K
        t = 0;                       %遗传代数
    ; P6 l+ W6 ]. E$ q+ c& C0 [    初始化P(t);                  %初始化种群或者染色体* N  Z1 T3 y3 W9 l
        计算P(t)的适应值;
    . [6 P# q: w" f& W$ r& u" [: O9 E    while(不满足停止准则) do% F/ e( Y7 }+ H8 w( R) w
            begin( A' z6 w" P) s, r& q
            t = t+1;
    1 l3 h+ k8 ]0 c: g        从P(t-1)中选择P(t);       %选择2 G1 k& Y% H& v: e
            重组P(t);                 %交叉或变异6 ~% y% a: K) S, |+ d$ Q' `2 P
            计算P(t)的适应值;
    5 [( @3 d" a! U" U         end
    # Q! q& e1 ~4 d- t7 l9 G    end4 V; c2 z9 M8 m4 e5 Z
    END2 m3 T- a3 H2 y2 m. j
    遗传算法工具箱$ p, [6 N; t2 Q6 g/ P+ K
    如果前面没看懂,没关系!!!为了省略艰深难懂的遗传算法,MATLAB软件做成了专门的遗传算法工具箱GA Toolbox,方便用户调用。(但有时候遗传工具并不是万能的,很多的情况下更需要具体问题具体分析)
    - ~, P# j0 U+ p+ O我们在MATLAB中帮助文档查看ga函数的定义
    # h3 V) D, T2 x' j$ G, H. |; yFind minimum of function using genetic algorithm
    " \3 x  A( X8 A$ C; m& v意思就是通过使用遗传算法找到函数的最小值。
    % i& i; A% c- a& H它的用法,如下:" \# j5 k( y8 H2 B
    x = ga(fitnessfcn,nvars)
    ; G1 {! L' v$ ]1 F8 h5 h1 x& fx = ga(fitnessfcn,nvars,A,b)
    % M3 N, C. I* {1 m& ?* M' Gx = ga(fitnessfcn,nvars,A,b,Aeq,beq)6 r  g# y3 F7 n! y
    x = ga(fitnessfcn,nvars,A,b,Aeq,beq,LB,UB)
    # r# p- b5 Y' G. wx = ga(fitnessfcn,nvars,A,b,Aeq,beq,LB,UB,nonlcon)
    ( c2 }" r2 j- Xx = ga(fitnessfcn,nvars,A,b,Aeq,beq,LB,UB,nonlcon,options)1 F' H5 \2 V0 ]$ a' c
    x = ga(fitnessfcn,nvars,A,b,[],[],LB,UB,nonlcon,IntCon)  r% i( j% |+ h! l/ Q- Z2 y
    x = ga(fitnessfcn,nvars,A,b,[],[],LB,UB,nonlcon,IntCon,options)
    . D! e1 I7 A, f( T* N" P, w: ]x = ga(problem)8 H1 g* P5 M$ f* t
    [x,fval] = ga(fitnessfcn,nvars,…)9 z% p$ f$ m5 E+ |2 w
    [x,fval,exitflag] = ga(fitnessfcn,nvars,…)
    ( T9 d7 T9 [5 ^[x,fval,exitflag,output] = ga(fitnessfcn,nvars,…)! f# P6 L; @; N6 V& D
    [x,fval,exitflag,output,population] = ga(fitnessfcn,nvars,…)& |( Y& i1 H3 i% h; R0 }0 Y
    [x,fval,exitflag,output,population,scores] = ga(fitnessfcn,nvars,…)
    4 C- Y' [% _5 ]( F1 H8 tx为经过遗传进化以后的自变量最佳染色体返回值;fval为最佳染色体的适应度,exitflag为算法停止的原因;output为包含每一代输出的结构以及有关算法性能的其他信息;population返回矩阵,群体其行是最终总体;scores返回最终群体的得分。
    4 e$ P0 T" x# hfitnessfcn为适应度函数,nvars为目标函数自变量的个数,A,b,Aeq,beq为约束条件(与线性规划类似)LB,UB为自变量的上界和下界;函数nonlcon接受x并返回矢量C和Ceq,分别表示非线性不等式和等式;IntCon要求中列出的变量采用整数值;options为默认优化参数由options中的值替换。( N5 a8 \: V0 d
    我们用的更多的是这种1 |/ r: D  A& A+ {) Y1 q
    [x,fval] = ga(fitnessfcn,nvars,options)
      y9 W" G& A+ y$ v令options=gaoptimset()函数,gaoptimset函数用法如下:
    ( l& Y3 b- ~0 h0 b/ S0 z4 P
    9 e- S' R6 U6 q3 T属性如下:* t1 V8 S1 p/ U4 x. ]

    , y8 @+ u0 U. M% R
    ! K% W7 `# {$ ?; M
      B$ c" e# v5 J$ @) I
    / U' g- g/ v7 V7 l/ n
    $ ]. M+ b, q/ C3 w
    0 _% y  }) G+ Z+ L* I, R) J0 o3 x' B
    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-14 14:58 , Processed in 0.420776 second(s), 51 queries .

    回顶部