QQ登录

只需要一步,快速开始

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

组合优化算法-现代优化算法(四):改进的遗传算法

[复制链接]
字体大小: 正常 放大
浅夏110 实名认证       

542

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2020-11-14 17:15
  • 签到天数: 74 天

    [LV.6]常住居民II

    邮箱绑定达人

    群组2019美赛冲刺课程

    群组站长地区赛培训

    群组2019考研数学 桃子老师

    群组2018教师培训(呼伦贝

    群组2019考研数学 站长系列

    跳转到指定楼层
    1#
    发表于 2020-5-22 15:31 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    无人机航路规划问题实际上是一个组合优化问题,是优化理论中的 NP—hard 问题。 因为其解空间不连续,解邻域表达困难,所以难以用通常的算法求解。遗传算法作为现代优化算法之一[1],其主要特点是对非线性极值问题能以概率 1 跳出局部最优解,找 到全局最优解。而遗传算法这种跳出局部最优寻找全局最优特性都基于算法中的交叉和 变异。在传统遗传算法的结构中,变异操作在交叉操作基础上进行,强调的是交叉作用, 认为变异只是一个生物学背景机制[2]。在具体交叉操作中,人们通常采用断点交叉法 (段交叉)多点交叉与均匀交叉,其中断点交叉是指随机地在基因序列中选择一个断点, 然后交换双亲上断点右端的所有染色体。在变异操作中,变异算子一般是用 Guassian 分 布的随机变异来实现[3,4]。近年来,也有学者尝试用 Cauchy 分布的随机序列来实现变异[5], 希望通过 Cauchy 分布宽大的两翼特性实现更大范围的变异,以利于找到全局最 优解。[6]从理论上分析了采用 Cauchy 分布随机变异进化算法的局部收敛性。[7]进一 步把二者结合起来, 采用两种分布的线性叠加,但仿真结果显示,算法改进效果并不十 分明显。文献[8]将生物进化看成是随机性加上反馈,并指出其中的随机性主要是由系 统的内在因素所引起,而不是由外部环境的随机扰动所造成。而混沌系统在其混沌域中 表现为随机性,它是确定系统内部随机性的反映,不同于外在的随机特性。
    # I4 D% Z5 @3 d7 V3 l- X& A
    ! C2 R1 k% M; h5 j本文根据以 上特点对基于求解航路规划的遗传算法进行改进,首先将变异操作从交叉操作中分离出 来,使其成为独立的并列于交叉的寻优操作,在具体遗传操作中,混沌与遗传操作联系 在一起,在交叉操作中,以“门当户对”原则进行个体的配对,利用混沌序列确定交叉 点,实行强度最弱的单点交叉,以确保算法收敛精度,削弱和避免寻优抖振问题;在变异操作中,利用混沌序列对染色体中多个基因进行变异,以避免算法早熟。
    6 v9 n0 Y" }' x$ R
    # ?# a( h" ]- z6 B6 P下面我们研究 1.2 中同样的问题。, b: }1 b5 p. \% w3 K; G- z% Y- I- Z
    $ x' H3 {8 N+ e, G) @

    7 T/ i6 s  |3 \) \) m3 T$ j' N" N
    ' ]2 A& K5 |! K: [+ t" e- O# p' |
    $ ?! o8 o" L( u% o* N" O" Z3 q& z+ M: L" N% b! s) u
    我方有一个基地,经度和纬度为(70,40)。假设我方飞机的速度为 1000 公里/小时。 我方派一架飞机从基地出发,侦察完敌方所有目标,再返回原来的基地。在敌方每一目 标点的侦察时间不计,求该架飞机所花费的时间(假设我方飞机巡航时间可以充分长)。
    # A: K1 e6 Y/ v. {# o
    3 k6 P1 Z! s2 }* d( M4 Q& [2 q4 w4 X/ u
    ; M7 z/ e8 B8 D& |- G/ z& T5 }1 W8 ?+ D- \3 p: R' H; E! h; r
    2)我方有三个基地,经度、纬度分别为(70,40),(72,45),(68,48)。假设我方 所有无人侦察机的速度都为 1000 公里/小时。三个基地各派出一架飞机侦察敌方目标, 怎样划分任务,才能使时间最短,且任务比较均衡。; [, K# b+ y" H. r4 {1 }

    " R! q) w$ `5 f* A* F2 模型及算法
    $ u8 D9 K5 G  A1 i% S与标准的遗传算法相比,我们做了如下的两点改进。3 k  J  ^4 S! b  L+ A

    # c0 b7 [0 \" ~- `) \( M$ E(1)交叉操作
    6 G7 j2 r' s  {2 K
    4 g+ G7 }* M4 ^( L3 A- z& y2 _% l$ D9 {
    ( s( b. i( f% h
    / p, T- @1 \+ l2 A9 V. m  @* ^

    & J3 q' E3 R7 X6 [( s5 f很明显这种单点交叉对原来的解改动很小,这可以削弱避免遗传算法在组合优化应用中 产生的寻优抖振问题,可以提高算法收敛精度。6 G- j; z6 e; `1 R

    1 @' r% v- W6 |+ j7 T(2)变异操作0 M' A6 z+ j5 H' o5 J3 R  g6 \+ M$ M
    变异也是实现群体多样性的一种手段,是跳出局部最优,全局寻优的重要保证。在 本文具体变异算子设计如下,首先根据给定的变异率(本文选为 0.02),随机地取两个在 2 到 101 之间的整数,对这两个数对应位置的基因进行变异,具体变异以当前的基因值 为初值利用混沌序列 x(n +1) = 4x(n)(1− x(n)) 进行适当次数的迭代,得到变异后新的基因 值,从而得到新的染色体。! n5 L* q0 X1 s0 ?) F( U. v

    - R2 ]- `- C* F3 仿真结果对比及算法性能分析+ y7 ?1 i& {9 V6 k. W9 N; d# |
    计算的 MATLAB 程序如下:  \/ E- V6 e) U2 y  N& z* L. m

    % @9 l2 o' X% g6 x! Stic; [# l; M+ T4 @' q* N
    clc,clear- S& N' G3 p7 H/ {, }1 G
    load sj.txt %加载敌方 100 个目标的数据
    - B3 X+ m8 W7 U4 g- ~# _! G1 ^* @x=sj(:,1:2:8);x=x(;$ X; v1 y  o  c  D: W& H
    y=sj(:,2:2:8);y=y(;$ X* w: {0 e5 D9 o$ U8 y9 O% j
    sj=[x y];
    1 P% B# k1 r7 S5 L, V) gd1=[70,40];; r! m, i: W8 h* x; R
    sj=[d1;sj;d1];3 E% r. U6 `9 k; I- I5 w
    %距离矩阵 d! z( ~, h, u" T$ Q: y- g1 z* b
    sj=sj*pi/180;' z9 q2 s2 e% {
    d=zeros(102);: c* t& \% [0 ~" S  M
    for i=1:101/ S% g7 ?) P5 `* C/ M4 f
        for j=i+1:102' }3 M# @/ r( k1 C8 s% D
            temp=cos(sj(i,1)-sj(j,1))*cos(sj(i,2))*cos(sj(j,2))+sin(sj(i,2))*sin(sj(j,2));) S! h6 a5 b0 V* E& R" r# h' X
            d(i,j)=6370*acos(temp);8 ^$ D% a0 m; N
        end
    ! o+ b' s5 v3 L9 H/ r6 Zend! T& K: G5 n5 f% e, q
    d=d+d';L=102;w=50;dai=100;
    $ m  |& f) y5 C%通过改良圈算法选取优良父代 A
    ) }+ d! _( f% {9 D2 \for k=1:w
    1 ^% a( Y7 Z+ D$ q. l% i5 h    c=randperm(100);9 a& B& ]1 g% V. M# g3 B$ l
        c1=[1,c+1,102];
    1 a: w/ O% D4 k8 K& d% |    flag=1;
    - P, j. O% u5 q) S! e  z  a, V3 N    while flag>0+ k% x& y* \( Q* [6 w! `: O2 q
            flag=0;
    8 |$ A" E2 H- m& Q8 d) X, E        for m=1-3
    & g' W& g. `, l' ^3 L5 X$ ?            for n=m+2-11 S9 I, W0 }  u# k0 R& l
                    if d(c1(m),c1(n))+d(c1(m+1),c1(n+1))<d(c1(m),c1(m+1))+d(c1(n),c1(n+1))
    0 f* N( D+ {; D& u/ p1 v; g                    flag=1;
    , ]0 M! {/ [0 p$ L) s                    c1(m+1:n)=c1(n:-1:m+1);
    ; K& Q, ?* h4 D& Z" g                end
    ) P6 v4 j, D5 N8 r" h; e6 u            end
    3 P  f" s3 ^  _/ O& Z# i8 Z" l, q: Z        end
    & V2 A' K" L' Q. i$ e( h     end- f* r2 L1 H# ~) Z  w
        J(k,c1)=1:102;
    1 P$ ^. y  J7 Eend) R( |+ @* h" t
    J=J/102;2 n; F) a5 |  ]* x
    J(:,1)=0;J(:,102)=1;
    3 z( P9 s7 Y# L: s4 j* ^- |rand('state',sum(clock));" o; I3 y' w0 R2 o' i. D1 r3 @- }
    %遗传算法实现过程. v1 o$ Y* M, P3 ^- x5 X, i% d
    A=J;
    0 N7 u) E/ H/ S9 m! K) b) ?1 n3 nfor k=1:dai %产生 0~1 间随机数列进行编码
    $ \6 q$ T9 a- q2 f  v    B=A;
    # s  r  d5 E* A: U, F% D %交配产生子代 B
    7 a3 M; q& R- X+ Y- i* }$ X    for i=1:2:w7 U' _# @+ p  p& U8 y8 R; C/ Q- l. ?: M* h
            ch0=rand;ch(1)=4*ch0*(1-ch0);; a+ h7 I0 v9 P' v, q8 ?0 G' i3 n
            for j=2:505 }, R% L. U: h  \
                ch(j)=4*ch(j-1)*(1-ch(j-1));6 H7 U, A7 F9 A( W5 ?
            end
    * @6 j. S/ j0 l        ch=2+floor(100*ch);
    - }' S3 p" {% Z! a8 J1 T        temp=B(i,ch);' E! g" A9 f4 Z" u- U% Y
            B(i,ch)=B(i+1,ch);( R; K5 i4 n/ S9 L, P9 j
            B(i+1,ch)=temp;' I% l! t- h& S
        end7 ~: B- \- ~$ b2 R, F
    %变异产生子代 C" M1 q8 Y2 ?, d- M7 F
        by=find(rand(1,w)<0.1);! h3 a  a9 Z" k% T  \
        if length(by)==0
    5 {! K+ \2 t& P5 N( e        by=floor(w*rand(1))+1;    6 h, N' `" F1 P9 t0 v2 d# ]
        end
    7 L$ U& q5 d" N! L+ f* k# \    C=A(by,;
    * x$ l) |: b5 S. j    L3=length(by);8 x* Q) [; d3 F- l5 d
        for j=13
    ; |! |( q* T+ v     bw=2+floor(100*rand(1,3));; A& j6 D. L- a3 d7 J1 f
         bw=sort(bw);
    1 ~  i* F" z! }     C(j,=C(j,[1:bw(1)-1,bw(2)+1:bw(3),bw(1):bw(2),bw(3)+1:102]);
    0 d7 z8 ~% p- \2 x. ^7 Zend
    8 F* a+ a5 j* W# P4 Z6 U    G=[A;B;C];. c9 z! y. W. c/ P$ B
        TL=size(G,1);
    % Y8 m" }4 I7 J+ X4 d) n %在父代和子代中选择优良品种作为新的父代0 B; z/ ^/ h! m' [& K# y7 r
        [dd,IX]=sort(G,2);temp(1:TL)=0;0 ]2 d8 Z8 ?( O( R5 f; `$ a+ b
        for j=1:TL
    5 O- q  G, J+ u- H' A        for i=1:1019 R. |: G: _5 V+ A- ^
                temp(j)=temp(j)+d(IX(j,i),IX(j,i+1));
    6 {. S" I" s' [+ v2 U* g        end
    ; N9 q/ s. p9 |# {: k    end- Z/ E# S: N5 ?4 V* W! r
        [DZ,IZ]=sort(temp);
    ) h- y7 Z# G: K4 [2 P" f    A=G(IZ(1:w),;
    7 T/ R) n7 L/ ?6 Xend
    ' `+ ?, C6 |, c- Gpath=IX(IZ(1),9 ]8 C" t9 M* c/ D$ V3 M$ N
    long=DZ(1): Z' k- t! ^' y$ s( U# I# O- Q
    toc 8 l- e+ q8 f8 |6 |, r
    0 K; ?# H0 R  Q! u4 y$ r0 @5 k
    : C# x' e3 e$ v, p" G: H& C
    在仿真试验中,我们对文中航路规划问题分别利用断点交叉和换位变异结合的遗传 算法,多点交叉和移位变异结合的遗传算法[10]和文中提出的改进算法进行求解比较。 表 1 是各种算法种群规模( M = 50)和迭代次数(G = 100 )都相同时连续 20 次求 解的平均值(公里),算法平均运算时间(秒)。
    * h: D! W" |# L: y$ k- L
    ' Y7 p: a& }2 i& q9 I; I$ Y  B2 u% i/ F: C1 N/ y) K0 s
    ( X8 W0 t3 h; N0 h" ]; B
    & J- j( H+ g- g' }6 @
    本文从算法结构到具体的遗传操作都进行了改进,其中变异操作从交叉操作中分离 出来,使得遗传算法也可以通过并行计算实现,提高算法实现效率。其次改进后的算法, 分别采用变化强度不同的交叉操作和变异操作,其中交叉操作采用强度最弱的单点交 叉,保证了算法收敛精度,削弱和避免算法因交叉强度大而产生的寻优抖振问题。当然 单一的单点交叉很容易使算法早熟,文中采用较大强度的多个基因变异正好解决早熟问 题。从仿真结果可以看到改进后的算法效果较为明显。, d4 b  C% L/ y7 N- B3 o' E/ T

    2 T$ f9 m% t4 m, O& d1 \3 ?! ^2 T8 c& M2 l5 K; I6 J
    ; q% G  V& f% |0 S. [
    ————————————————- |  ?* E0 @/ Z: D, a' i" v* C# s& Q
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    * v" ]& {4 B: |, B原文链接:https://blog.csdn.net/qq_29831163/article/details/89672986
    0 S" ^7 [$ f* y; v/ b1 _& ^% r" T, c6 h' H: i( U9 N" p

    4 x- f: h/ U+ q, q  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-11 06:46 , Processed in 0.446518 second(s), 51 queries .

    回顶部