在遺傳算法中
執行交叉變異得到正確的子代結果
是十分重要的一步

為了更好的交叉變異
有人發明了Grefenstette編碼

Grefenstette的編碼定義如下：
既是：通過已有的順序（非常重要，編碼解碼都靠它），
產生一個訪問序列，
且每訪問完一個點，就需要把點從已有序列刪除，形成新的序列
作為新的訪問前置條件
例如已有順序
1，2，3，4，5
編碼集：4，2，1，2，1
通過編碼集得到的解碼集：
解碼集：4，2，1，5，3

注意編碼集最后一位永遠都是以1結束

以種群數量S 與 生存環境數量M為研究對象，

產生一個1——M的序列
加上s就是：
產生一個SxM的原始已有序列矩陣，每一行為一個種群對于的已有序列

再根據Grefenstette的編碼規則：
訪問生成序列下標并刪除對應點
得到編碼程序：

%Grefenstette的編碼與主函數 s=3;%種群數量M=5;%城市數量 pop=zeros(s,M); for i=1:s%產生1-t的隨機數列pop(i,1:M)=randperm(M); end fprintf('原始信息（一行為一個信息）'); pop % 編碼grePop=zeros(s,M); %原始種群Grefenstette編碼的結果%Grefenstette的編碼for i=1:s %一行代表一個生存環境temp=[1:1:M];for j=1:Mloc=find(pop(i,j)==temp);grePop(i,j)=loc(1);temp(loc(1))=[];endendfprintf('編碼結果：');grePop

又根據編碼結果來計算解碼集：
得到編碼程序：
grePop是已知的編碼集矩陣
s,M對應于編碼程序的s,M

%Grefenstette的解碼 function GrefenstetteDecode(grePop,s,M) %解碼for i=1:s %初始數據1--Mtemp=[1:1:M];for j=1:M%解碼集%temp(grePop(i,j))pop(i,j)=temp(grePop(i,j));%一行代表一個生存環境%去掉temp(grePop(i,j))temp(grePop(i,j))=[];endendpop

舉例：
2種群，4城市：

編碼截圖：

解碼截圖：

可以看到完全將隨機產生的原始信息還原回來了。

完。
加油，奧利給！

總結

以上是生活随笔為你收集整理的matlab-Grefenstette的编码与解码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。