dalsi husty clanok z tej stranky
http://technet.idnes.cz/vojnicuv-rukopis-0z7-/veda.aspx?c=A130701_171454_veda_mla

#file candidate-target-language-words should contain tokens which possible occur in the transliterated version of the Voynich manuscript

#scripts exploits evolutionary algos to converge to a set of substitution rules which transliterate (project) the source Voynich sequences upon the candidate language

#Voynich sequence file specifies in the only command line argument of the script

#Daniel D. Hromada (GPL Licence)

open(S,"<candidate-target-language-words");

my @rexp;



my $origin;



my %slova;

while (<>) {

chomp;

$slova{$_}=true;

# $origin.=$_;

}



#should we reverse the text, Mr. Dee ??? ;-)

#$origin=reverse($origin);





while (<S>) {

chomp;

print $_;

push @rexp, qr/^$_/i;

}



my $N_popul=42;

$chars="abcdehiklmnoprstuvz";

my @population;

my @Fitness;



for(0 .. ($N_popul)) {

my $ind="";

for (0 .. 37) {

$ind.=substr($chars,int(rand(length($chars))),1);

}

print "$ind n";

push @population,$ind;

push @Fitness,0;

}



sub Roulette_Wheel {

my $sum=0.0;

my @upper_bound;

my @lower_bound;

my $k;

for my $i (0..$N_popul) {

$sum=$sum+$Fitness[$i];

}

for my $i (0..$N_popul) {

$upper_bound [$i]=$Fitness[$i]/$sum;

}

$lower_bound[0]=0;

for my $i (1..$N_popul) {

$lower_bound[$i]=$upper_bound[$i-1];

$upper_bound[$i]=$upper_bound[$i-1]+$upper_bound[$i];

}

my $aux=0.001+0.998*(rand());

my $i=1;

my $j=$N_popul;

do {

$k=int(($i+$j)/2);

if ($aux>$upper_bound[$k]) {

$i=$k+1;

} else {

$j=$k-1;

}

} until ($lower_bound[$k]<=$aux && $aux<=$upper_bound[$k]);



return $k;

}



my $gen=0;

while (1) {

$gen++;

print "Generation $gen n";

for (0 .. $N_popul) {

my $i=$_;

my $text=$origin;

#print $text;

my $ind=$population[$i];

# eval "$text=~tr/Q$oldlistE/Q$newlistE/";



# initialization

my $old = "abcedfghijklmnopqrstuvwxyzAEFHIKOPSTYZ";

my $new = $ind;

my %translit;



# map chars in $orig to chars in $replace

@translit{split //, $old} = split //, $new;



my %matched;

for my $re (@rexp) {

for my $slovo (keys %slova) {

$slovo =~ s/(.)/defined($translit{$1}) ? $translit{$1} : $1/eg;

next if $matched{$re};

if ($slovo=~/$re/) {

#print "n $slovo matched by $renn";

$Fitness[$i]++;

$matched{$re}=true;

next;

}

}

}

print "individuum $population[$i] Fitness $Fitness[$i] n";

my $min=777;

my $min_i;

my $max=0;

my $max_i;

my $sum=0;

for my $i (0 .. $N_popul) {

if ($Fitness[$i] < $min) {

$min=$Fitness[$i];

$min_i=$i;

} elsif ($Fitness[$i]>$max) {

$max=$Fitness[$i];

$max_i=$i;

}

$sum+=$Fitness[$i];

}



print "Total fitness $sum n";



#the best replaces the worst

$population[$min_i]=$population[$max_i];





#mutations

for (0 .. $N_popul) {

my $subj=$_;

for (0..37) {

my $recpos=$_;

if ($subj!=$max_i && rand() < 0.07) {

my $win=Roulette_Wheel();

substr($population[$subj],$recpos,1,substr($population[$win],$recpos,1));

} elsif (rand() <0.03 && $subj !=$max_i) {

substr($population[$subj],$recpos,1,substr($chars,int(rand length($chars)),1));

}

}

}

@Fitness=();

}


node backed-up from kyberia.cz