现代算法 基因算法

Posted zero27315

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了现代算法 基因算法相关的知识,希望对你有一定的参考价值。

  2019-03-10 16:37:59

  转眼间就到了第三周,刚考完Python二级,江苏省的题目有点太那个了吧,什么都有,第一道程序题给的最后结果都是错的?

花了我好长时间在那里,一直在找设置精度的函数,然鹅,回到宿舍用matlab计算也是python上跑的结果.呵呵,今天终于开始看算法了

第一个算法就是基因算法.话不多说开始.

  基因算法原理就不重述了,直接上代码

技术图片

 

  程序设计中选择的二进制编码,种群规模二 20,染色体杂交时的概率为 0.25,变异时概 率为 0.01,适应度函数为函数本身,终止条件为达到预定最大迭代次数。 

Objfunction & initial group

%n - 繁衍n代 不加上第一代,初始化父代
function gene(n)
f = @(x1,x2)21.5+x1*sin(4*pi*x1)+x2*sin(20*pi*x2);

% initial group
%种群数为20个
V =cell(1,20);
for i = 1:20
    temp1 = round(rand(1,33));
    %33 = 18+15=5+13+2+13
    temp2 = vec2str(temp1);
    V{i} = temp2;
end
clear temp1 temp2

record = adapt(V,f);
maxrec = maxrecord(record,V);
%自然选择,保存最适合
tic
for i=1:n
    V = chfather(V,record);
    V = opcrossover(V);
    V = variation(V);
    record = adapt(V,f);
    temp = maxrecord(record,V);
    if (temp(1,1)>maxrec(1,1))
        maxrec = temp;
        mark = i;
    end
end
time = toc;
format long
fprintf(       x1: %f
,maxrec(2));
fprintf(       x2: %f
,maxrec(3));
fprintf(f(x1,x2) : %f
,maxrec(1));
fprintf( from: %1.0f(th) generation 
,mark);
fprintf(    time: %f s
,time);

    

 

  

vec2str
%Char2str
%将向量转化成字符串
%bin2dec()中的参数为字符串!
function str = vec2str(vec)
len = length(vec);
str = ‘‘;
for i=1:len
    temp = num2str(vec(i));
    str = strcat(str,temp);
end
end

 

 
adapt
function record = adapt(V,f)
record = zeros(20,1,double);
for i=1:20
    temp1 =V{i};
    temp2 = temp1(1:18);
    temp3 = temp1(19:end);
    record(i,1) = f(bin_x(temp2,1),bin_x(temp3,2));
end
end

 

maxrecord
function maxrec = maxrecord(record,V)
[m,i] = max(record);
maxrec(1,1) = m;
temp1 = V{i};
temp2 = temp1(1:18);
temp3 = temp1(19:end);
maxrec(2,1) = bin_x(temp2,1);
maxrec(3,1) = bin_x(temp3,2);
end

 

 chfather
%选择父代,似乎染色体进行了异位
function V = chfather(V,record)
F = sum(record);
%pk-平均,qk-累计
pk = record/F;
qk = zeros(20,1,double);
for i=1:20
    qk(i) = sum(pk(1:i));
end
%模拟转盘转动20次
r = rand(20,1);
for i = 1:20
    k = 1;
    while(r(i)>qk(k))
        k = k+1;
    end
    r(i) = k;
end
temp = V;
for i =1:20
    V{i} = temp{r(i)};
end
end

 

 opcorssover
function V = opcorssover(V)
%交换概率?pc?
%one point corssover
pc = 0.25;
l = 1;
while(l==1)
    r = rand(20,1);
    mk = find(r<pc);
    l = numel(mk);
end
if mod(numel(mk),2) == 1
    mk = mk(1:end-1);
end
ri = randi([1,32],1,numel(mk)/2);
for i = 1:numel(mk)/2
    [V{mk(2*(i-1)+1)},V{mk(2*(i-1)+2)}] = onecross(V{mk(2*(i-1)+1)},...
        V{mk(2*(i-1)+2)},ri(i));
end
end

 

variation
function V = variation(V)
pm = 0.01;
for i = 1:20
    r = rand(33,1);
    k = find(r<pm);
    for j = 1:numel(k)
        V{i} = vari(V{i},k(j));
    end
end
end

 

Vari
function gnew = vari(g_old,pos)
gnew = g_old;
if gnew(pos) ==1
    gnew(pos) =0;
else
    gnew(pos) =1;
end
end

 

 onecorss
function [g1,g2]=onecross(gen1,gen2,pos)
len = numel(gen1);
g1 = gen1(1:pos);
g2 = gen2(1:pos);
g1(pos+1:len) = gen2(pos+1:end);
g2(pos+1:len) = gen1(pos+1:end);
end

 

bin_x
function num = bin_x(bin,opt)
switch opt
    case 1
        num = -3.0+bin2dec(bin)*((12.1-(-3.0))/(2^18-1));
        %12+3.0 = 16 ,0,1,2,3,4
        %6-4 = 2 ,0,1,2,
        %精度为0.5^13
    case 2
        num = 4.1+bin2dec(bin)*((5.8-4.1)/(2^15-1));
end
end

 


 

 

以上是关于现代算法 基因算法的主要内容,如果未能解决你的问题,请参考以下文章

MTSP遗传算法解决

如何捕获相似基因(两个相似哈希算法分析)

编程内功修炼 - 算法

c语言实现*/遗传算法改进BP神经网络原理和算法实现怎么弄

(基因分析黑科技)基于基数排序算法实现的samtools sort

有人可以解释啥是 SVN 平分算法吗?理论上和通过代码片段[重复]