|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
- F5 D2 Z( ~4 n' b5 N3 c/ a' o I一、简介
; A) ?- Q) Q+ d& Z+ t7 F1 遗传算法基本理论9 ~5 @, p6 m5 m2 }4 c
遗传学认为,遗传是作为一种指令遗传码封装在每个细胞中,并以基因的形式包含在染色体中,每个基因有特殊的位置并控制某个特殊的性质。每个基因产生的个体对环境有一定的适应性。基因杂交和突变可能产生对环境适应性强的后代,通过优胜劣汰的自然选择,适应度值高的基因结构就保存下来。
* l3 L5 [7 L& E: N遗传算法借鉴“适者生存”的遗传遗传学理论,将优化问题的求解表示成“染色体”的“适者生存”过程,通过“染色体”群的一代代复制、交叉、变异的进化,最终得到的是最适应环境的个体,从而得到问题的最优解或者满意解。这是一种高度并行、随机和自适应的通用的优化算法。
2 _# X! ]" O+ `& W遗传算法的一系列优点使它近年来越来越受到重视,在解决众多领域的优化问题中得到了广泛的应用,其中也包括在交通领域的成功应用。# s M0 |* A% o6 o! v
2 遗传算法的特点
* X! C2 Y3 V9 W. t3 b遗传算法是模拟生物自然环境中的遗传和进化过程而形成的一种自适应全局优化概率搜索算法。是一类可用于复杂系统优化计算的鲁棒搜索算法,与其他一些优化算法相比,它具有很多特点。) [2 @+ Q+ a$ G# ?) g4 w# I
传统的优化算法主要有三种:枚举法、启发式算法和搜索算法:
6 y2 D9 u1 f2 T+ b4 i( F" Q5 y1.枚举法2 A* T. l7 d6 U/ N! y
枚举法在可行解集合内枚举所有可行解,以求出精确最优解。对于连续函数,该方法要求先对其进行离散化处理,这样就可能因离散处理而永远达不到最优解。此外,当枚举空间比较大时,该算法的求解效率非常低,极其耗时。
: ~- L: s' K C+ f2.启发式算法- Q: R; v9 D; |. W$ v
启发式算法寻求一种能产生可行解的启发式规则,以找到一个最优解或近似最优解。启发式算法的求解效率比较高,但对每一个需求解的问题必须找出其特有的启发式规则,这个启发式规则一般无通用性,不适合于其他问题。
1 @! T5 k" }7 d. e! k$ h3.搜索算法
) Y' K5 X6 M! w" l搜索算法在可行解集合的一个子集内进行搜索操作,以找到问题的最优解或者近似最优解。搜索算法虽然保证不了一定能够得到问题的最优解,但若适当的利用一些启发知识,就可在近似解的质量和效率上达到一种较好的平衡。* o' |6 J2 `- I z: C/ G# L
3 基本遗传算法的工作流程
k$ N( Z( M. H7 L2 n
: F. j; ?* X% h h! a1 L: _, R1 a& t
4 E; U) R) f! h0 H" |' m6 U4 适应度函数
9 D; `, m$ `; S* H, N; J- G6 D( u G- w) y: w+ E1 ]+ w
' S' S2 ^% Y& ]7 l3 |7 B
# @8 P$ |. E7 G: ^4 L+ I5 选择算子 ~/ |8 u# ?# Y! ]1 S# t
选择又称为复制,是在群体中选择生命力强的个体产生新的群体的过程。遗传算法使用选择算子来对群体中的个体进行优胜劣汰操作,根据每个个体的适应度大小选择,适应度较高的个体被遗传到下一代群体中的概率较大;反之亦然。这样就可以使得群体中个体的适应度值不断接近最优解。选择算子的确定的好坏,直接影响到遗传算法的计算结果。* {, ]5 \: o* a5 ? t- V/ L
下面介绍几种典型常用的选择算子:
4 ?/ ~8 C2 v6 o, D. h$ i0 ~& o: b1.轮盘赌选择
4 n) _0 v) h( r# F# ~6 Y S$ b2.随机竞争选择
' j* a& [: I, q$ `/ n, H: J3.随机遍历选择
4 d% g$ u/ I* _" ~/ D6 m4 B5 B4.排序选择) `$ U6 ]7 C: ^
5.联赛选择1 R6 v8 L* i. y
6 交叉算子5 [6 |$ P6 G' o( j
下面介绍几种适合于二进制编码个体或十进制编码个体的交叉算子。
* ~/ \* G* m( D( a5 H8 ^' }l.单点交叉
; J0 n9 L8 s1 n, p& F: T单点交叉(One-point Crossover)又称为简单交叉,是最常用和最基本的交叉操作算子。它以二值串中的随机选择点开始,对每一对相互配对的个体,依设定的交叉概率在其交叉点处相互交换两个个体的部分染色体,从而产生出两个新的个体。
$ ^" c' a8 B4 I8 }3 G( h2.两点交叉与多点交叉5 w+ j: K: w+ i8 U1 L
两点交叉(Two-point Crossover)是指在个体编码串中随即设置两个交叉点,然后再进行部分基因交换。两点交叉的具体过程是:% ? J& E# X) Y0 A9 g5 _# [
(1)在相互配对的两个个体编码串中随即设置两个交叉点。
% l5 U7 Y6 I6 l) _, Y( m(2)交换两个个体在所设定的两个交叉点之间的部分染色体。
' O- F! I* D; H7 t* P, P2 e3.均匀交叉
3 q6 S( Q* f6 ~0 e均匀交叉是指两个配对个体的每个基因座上基因都以相同的交叉概率进行交换,从而形成两个新的个体。其具体运算可通过设置一屏蔽字来确定新个体的各个基因如何由哪一个父代个体来提供。1 ^* k3 o& K& a2 t/ n: g7 `2 s
4.算术交叉4 ~5 I4 \+ p2 {$ f4 }% w
算术交叉是指由两个个体的线性组合而产生出两个新的个体。为了能够进行线性组合运算,算术交叉的操作对象一般是浮点数编码所表示的个体。
# k+ i/ v8 D- }- d8 O7 变异算子# E- q, ~2 u. e4 o8 E4 `( t" C
遗传算法中所谓的变异运算,是指将个体染色体编码串中的某些基因座上的基因值用该基因座的其他等位基因来替换,从而形成一个新的个体。变异是遗传算法生成新个体的主要方法之一,变异运算可以使算法在运行过程中维持种群的多样性,有效避免早熟,起到改善遗传算法局部搜索能力的作用。- u$ P; W; J! \) M4 o# `3 A# A* @8 ~
遗传算法中变异算子也应根据不同的要求进行选择和设计,下面是几种常用的变异算子。
/ a! ?8 o' }: n8 遗传算法的改进
" p. v, z4 e8 ^( L1 t" p0 B
7 ^* `0 E, Q, W: l2 k
5 z4 L* P2 h v; n7 \
' x2 S) K9 c# Y: r4 z, e二、源代码) Q- i( C) Z& j6 Z) A% K
%% GA
$ y, F; O% `2 c4 u& _! b/ h4 |%% 清空环境变量
* i# g* |3 K; T% Kclc,clear,close all
% D1 u' ~" F$ n; P% \warning off
5 F9 c# `1 I" B Y2 o' _4 U+ i% feature jit off8 d; Y/ p0 `# c$ f5 @
%% 遗传算法参数初始化
5 D( H- O' {( ]9 _7 W) t& Pmaxgen = 50; % 进化代数,即迭代次数
! N, F7 P' }/ n9 C1 k" `" \sizepop = 100; % 种群规模
8 h9 n8 o* e7 b6 {pcross = [0.7]; % 交叉概率选择,0和1之间
2 s+ P' g& Y6 H1 n7 W# }pmutation = [0.1]; % 变异概率选择,0和1之间
! \' D5 s: X. V* C9 V. P3 d%染色体设置
) H5 M, i% @# glenchrom=ones(1,3); % t1、t2、t3
5 P2 Z; Y5 q1 J/ a" @bound=[38,59;26,37;33,44;]; % 数据范围
5 D/ }1 F& N8 F4 [% [( I%---------------------------种群初始化------------------------------------
1 U, v6 D" j6 k! x# p* ]4 rindividuals=struct('fitness',zeros(1,sizepop), 'chrom',[]); %将种群信息定义为一个结构体
; X2 q4 z) Q; d. c0 o; n. y# Wavgfitness = []; %每一代种群的平均适应度
" D3 H: i9 u- r obestfitness = []; %每一代种群的最佳适应度
9 e+ E4 {* T0 {; ~) p% e, A* Rbestchrom = []; %适应度最好的染色体
" J8 x# F# c; S) K9 s" K7 H% p1 `) L! | ?! H
%% 初始化种群+ A" A, K- L A4 i5 X8 y1 ~ V
for i=1:sizepop
2 Q5 Y8 }+ N1 c7 C0 N % 随机产生一个种群
+ _8 e$ I! T( a" o% \( C8 w+ C individuals.chrom(i,:)=Code(lenchrom,bound); % 编码(binary和grey的编码结果为一个实数,float的编码结果为一个实数向量)& L, p5 t3 x9 e5 [2 W
x=individuals.chrom(i,:);
) E' w5 W: E+ X( [ % 计算适应度3 j' m: J* s4 u! ~9 r, M
individuals.fitness(i)=fun(x); % 染色体的适应度
[% X2 @! i( C/ C' u) O& tend) p2 R! J: q6 o. ]
& j+ x2 G J7 g* J( ]
%% 找最好的染色体
$ H7 D3 R6 V9 \8 M d2 {; @6 Z0 z[bestfitness bestindex] = min(individuals.fitness);, M: |4 Z3 D' p) b9 U( b* P
bestchrom = individuals.chrom(bestindex,:); % 最好的染色体
/ s; A# J r. R( T2 z- U) P% 记录每一代进化中最好的适应度和平均适应度& u; x8 a7 V% }; r/ f$ D9 R
trace = [bestfitness];
! [# F# t0 f/ F: |) ?: M1 I8 C+ R2 Q# H
%% 迭代求解最佳初始阀值和权值
% P5 W: }8 \" J- v. W% 进化开始* L+ p( u* `: L0 t$ Q! Q; z
for i=1:maxgen0 o" X( H, }+ k4 j+ F! h( C+ }
disp(['迭代次数: ',num2str(i)])8 u4 w1 ?- e! L9 o
% 选择7 y" f- P( V4 R: ]) G1 X+ p
individuals=Select(individuals,sizepop);
: i: q" A. R7 \1 X8 | % 交叉
' V o3 q, }0 m/ X; I. G individuals.chrom=Cross(pcross,lenchrom,individuals.chrom,sizepop,bound);
6 x: {6 }: M/ N7 P/ } % 变异, z( W. o5 U; s- s# U
individuals.chrom=Mutation(pmutation,lenchrom,individuals.chrom,sizepop,i,maxgen,bound);8 g* x: F i! w# J
0 z- \# z8 p# E+ x% y % 计算适应度9 J$ u E. ?& r$ O% E* q8 f5 }
for j=1:sizepop: @$ j% i( X! D
x=individuals.chrom(j,:); % 解码 m6 d4 N/ C3 X0 J$ z
individuals.fitness(j)=fun(x); % 染色体的适应度
/ a$ E) |" x f" E end
7 X7 `$ ?" C" E) j6 k%% 改进的GA1 M7 e1 I8 Z0 K3 @
%% 清空环境变量
. @6 \0 y& z+ ~% |( J% clc,5 o3 v. o: r# M/ v" ]8 A
clear % ,close all % 清除变量空间
- s$ ]3 u6 S l1 u! Kwarning off % 消除警告
+ p8 G$ }. _9 p4 g1 V% feature jit off % 加速代码执行
6 ]' R# S- m& q1 J3 ` H%% 遗传算法参数初始化
- L( r$ n4 a4 `7 V. Z1 e, pmaxgen = 50; % 进化代数,即迭代次数
8 I! F2 P4 \9 o, p. Esizepop = 100; % 种群规模
0 ~9 r# f) c" C7 rpcross = [0.7]; % 交叉概率选择,0和1之间
" N, ]3 K. b. T; b$ W( W5 opmutation = [0.01]; % 变异概率选择,0和1之间7 r ~" O+ f/ }! C F
delta = 0.1;
2 m: U b2 q. f- _5 ^/ t, Y# |% 城市交通信号系统参数& c: d R$ a" Y. G
C = 140;( z0 f' I1 i! h, ?) H
L = 10;
_. ?. S, [3 _/ E5 g% Gload('data.mat') % 包含交通流量q以及饱和流量xij
! Y$ T. D1 d$ }% P- v/ s* Nq = q./3600; % 转化为秒s
; j1 Z8 ?0 W! }5 @ Txij = xij./3600; % 转化为秒s8 l6 f) K8 g( v; T1 h$ M: N' \( O4 ^
%染色体设置
0 n0 P( H( ?% m: r7 B0 Z/ t9 v/ [lenchrom=ones(1,3); % t1、t2、t3- R9 i3 h" J m' @! l6 N
bound=[38,59;26,37;33,44;]; % 数据范围) ^" G- }4 Z: @( ^
%---------------------------种群初始化------------------------------------8 c& l+ d& W) }3 o* P
individuals=struct('fitness',zeros(1,sizepop), 'chrom',[]); %将种群信息定义为一个结构体! r3 i3 s2 G, p8 K' n5 G6 i: ?
avgfitness = []; %每一代种群的平均适应度
2 D9 |+ E, |7 R: Tbestfitness = []; %每一代种群的最佳适应度5 O" ~8 c' [3 M) p
bestchrom = []; %适应度最好的染色体9 s4 f* M7 f; ?1 e" [7 m3 N
8 ?; T9 k2 _( v! v: m
%% 初始化种群5 O& R2 _- y# w
for i=1:sizepop! Q* h2 k q5 u2 W I Y
% 随机产生一个种群
: _( L4 |/ q( j/ P: F individuals.chrom(i,:)=Code(lenchrom,bound); % 编码(binary和grey的编码结果为一个实数,float的编码结果为一个实数向量)- {; F' `$ g; y
x=individuals.chrom(i,:);
/ ~4 q+ p: t1 H % 计算适应度$ W. S8 L3 g% @, Z% Z k
individuals.fitness(i)=fun1(x); % 染色体的适应度
3 C# A! O: ^9 K; z3 P8 kend) U' s2 z: a7 A3 [0 T
. w4 p8 K2 Y. Z& d, r3 g
%% 找最好的染色体
9 z& z; _0 W- Q[bestfitness bestindex] = min(individuals.fitness);
/ w' F4 a* s ~( r4 i' \9 s# Y vbestchrom = individuals.chrom(bestindex,:); % 最好的染色体. ?3 }/ K+ ?# P' V4 ^, Y
% 记录每一代进化中最好的适应度和平均适应度
/ ~& ^5 b( p- s& otrace = [bestfitness];
, R' a/ S/ X5 P% i/ s4 l1 h: B/ \: A4 K& a" E
%% 迭代求解最佳初始阀值和权值7 d8 @) X" ~% _% s3 i+ |
% 进化开始$ ]) g% ~' {# B
for i=1:maxgen
& P; h! i0 O& l disp(['迭代次数: ',num2str(i)])
' Y9 @8 z' Q) I5 | % 选择
) O) J# s; W$ [4 _# \7 e4 x individuals=Select1(individuals,sizepop); # E# v# ?; Q& M8 A
% 交叉6 {8 \9 B- O" a* j" }
individuals.chrom=Cross1(pcross,lenchrom,individuals.chrom,sizepop,bound);
! p# n- q# e# P# ]! ]- c* t % 变异
* q0 W( w; M3 Y7 `: | individuals.chrom=Mutation1(pmutation,lenchrom,individuals.chrom,sizepop,i,maxgen,bound);& y2 N8 t- m! R0 l0 \4 D+ W
' \7 {; t4 P0 t5 K) G' b# m % 计算适应度( \& V- r3 k' P: N: Q
for j=1:sizepop
( a3 \7 P, r0 a% Q* n! k/ ]! M6 v) N x=individuals.chrom(j,:); % 解码
0 h* p0 a* u# M0 ]% w individuals.fitness(j)=fun1(x); % 染色体的适应度
" y |2 y Q& G/ v+ g( I+ n0 ^: n end
- }2 h0 t1 f" }; S fmax = max(individuals.fitness); % 适应度最大值
* p; u9 a/ y: f3 N fmin = min(individuals.fitness); % 适应度最小值; e( Q& [1 s0 v' s3 U S% u* L
favg = mean(individuals.fitness); % 适应度平均值- {" a7 S- T/ i+ t! }& J
individuals.fitness = (individuals.fitness + abs(fmin))./(fmax+fmin+delta); %适应度标定 F c1 o& K& K" d" j9 i
, p* D9 p7 p/ ~ i3 P0 z& ]/ ?. t! z
- {7 \2 U- J9 | [newbestfitness,newbestindex]=min(individuals.fitness);; @, h% n' k. e: s, |( N/ H
[worestfitness,worestindex]=max(individuals.fitness);
6 z' J* S; `/ U$ Z# a1 \2 m / Q9 K% q6 D1 ?) t& e
if bestfitness>newbestfitness
6 W/ O; ^9 L/ o( T3 U0 Z. [ bestfitness=newbestfitness;
9 V5 W# m( U, t: I* o7 B9 a bestchrom=individuals.chrom(newbestindex,:);! [; u- B5 e" t6 b+ S
end/ r: w" n# j8 \' ?% C6 i
# s9 X5 m+ ?' I4 O2 K8 [. |
2 E; e. r3 }- ]! U. k三、运行结果
) C0 e/ y8 E5 ?$ }2 T( g% K2 J
* X$ I! r$ {+ k; ^8 |) ] |
|