【正文】 ，即：，其所占的適應(yīng)度比例分別是，在保持這個(gè)比例的情況下對(duì)這個(gè)取值范圍放大1000倍. 按照順序在1～1000內(nèi)分別占不同的區(qū)間，當(dāng)隨機(jī)函數(shù)產(chǎn)生1～1000以內(nèi)的時(shí)，即使是適應(yīng)度比例最小的也有可能被選中，從而較好的保持了群體的多樣性. 由上所述則可選擇出適應(yīng)力強(qiáng)的，淘汰適應(yīng)力弱的個(gè)體從而得到總體適應(yīng)能力強(qiáng)的群體. 經(jīng)過選擇算子得出總體適應(yīng)能力強(qiáng)的A、B群體數(shù)據(jù)矩陣為，因數(shù)據(jù)量太大，文中就不給出具體數(shù)據(jù). （5）交叉重組依據(jù)選擇性集成思想 ,等概率使用兩點(diǎn)交叉法和區(qū)域交叉法這兩種差異性較大的交叉方法，擴(kuò)大遺傳算法的搜索范圍，避免過早收斂. 其中，兩點(diǎn)交叉法是先隨機(jī)設(shè)定兩個(gè)基因交叉位置，將父輩兩個(gè)個(gè)體在這兩個(gè)交叉點(diǎn)之間的基因鏈碼相互交換，從而形成新的個(gè)體；區(qū)域交叉法是隨機(jī)在染色體中選擇一個(gè)交叉區(qū)域，將第二條染色體的交叉區(qū)域加在第一條染色體的前面，第一條染色體的交叉區(qū)域加在第二條染色體的前面，在交叉區(qū)域后依次刪除與交叉區(qū)域相同的基因，得到最后的兩條子染色體. 將第（3）步得到的關(guān)于A，B種群的數(shù)據(jù)分別用兩種交叉算法來實(shí)現(xiàn)操作. 其中一半數(shù)據(jù)用兩點(diǎn)交叉法，另一半的數(shù)據(jù)用區(qū)域交叉法來進(jìn)行染色體的交叉重組. 其具體算法的MATLAB程序?qū)崿F(xiàn)見附件四. 經(jīng)過交叉重組得出的A、B群體數(shù)據(jù)矩陣為，因數(shù)據(jù)量太大，文中就不給出具體數(shù)據(jù). （6）染色體的變異為了避免過早收斂，依據(jù)選擇性集成思想選擇鄰居交換變異和兩點(diǎn)交換變異這兩種個(gè)性好且差異性較大的變異方法，等概率使用以擴(kuò)大搜索范圍. 其中，鄰居交換變異是產(chǎn)生一個(gè)隨機(jī)數(shù)，將該數(shù)對(duì)應(yīng)的基因和其后的基因交換；若該數(shù)對(duì)應(yīng)的基因是染色體中的最后一個(gè)基因，則將該基因與染色體的第一個(gè)基因交換；兩點(diǎn)交換變異是先產(chǎn)生兩個(gè)不同的隨機(jī)數(shù)，確定兩個(gè)交換點(diǎn)，然后對(duì)個(gè)體在此兩點(diǎn)的基因進(jìn)行交換. 經(jīng)過染色體變異得出的A、B群體數(shù)據(jù)矩陣為，因數(shù)據(jù)量太大，文中就不給出具體數(shù)據(jù). （7）種群交叉將兩個(gè)種群中的最優(yōu)解取出，再在每個(gè)種群中隨機(jī)選取n個(gè)染色體，將這n+1個(gè)染色體互換，進(jìn)入對(duì)方種群. 經(jīng)過種群交叉得出的A、B群體數(shù)據(jù)矩陣為，因數(shù)據(jù)量太大，文中就不給出具體數(shù)據(jù). （8）最佳個(gè)體保留法 要把群體中適應(yīng)度最高的個(gè)體不經(jīng)過配對(duì)交叉直接復(fù)制到下一代中，然后從新群體中淘汰一個(gè)適應(yīng)度最差的個(gè)體. 分別對(duì)A、B進(jìn)行獨(dú)立的操作. 經(jīng)過最佳保留法選擇后得出的A、B群體數(shù)據(jù)矩陣為，因數(shù)據(jù)量太大，文中就不給出具體數(shù)據(jù). （9）迭代的結(jié)束條件 在本文中，采用進(jìn)化的代數(shù)N作為循環(huán)迭代過程的結(jié)束，如果等于N，則算法結(jié)束，輸出適應(yīng)值最高的解；否則，繼續(xù)重復(fù)上述過程. 重復(fù)步驟（3），（4），（5），（6），（7），（8）依次進(jìn)行循環(huán)迭代，本題中設(shè)定迭代次數(shù)為N=1000. 最后得到30個(gè)近似的最優(yōu)解. 以上（2）—（9）流程的MATLAB程序?qū)崿F(xiàn)見附件二. （10）結(jié)果選出這30個(gè)近似最優(yōu)解中以時(shí)間作為權(quán)值最小的那一組解作為題設(shè)中要求的近似最優(yōu)解. 其中滿足要求的基因鏈碼（站點(diǎn)數(shù)）的順序即是顧客所需從起始點(diǎn)到終點(diǎn)站的路線 從起點(diǎn)站到終點(diǎn)站轉(zhuǎn)一次車從起點(diǎn)站到終點(diǎn)站轉(zhuǎn)一次車遺傳算法流程中基于站點(diǎn)序號(hào)的編碼，交叉重組，染色體的變異，種群交叉，迭代的結(jié)束條件和結(jié)果的原理與從起點(diǎn)站到終點(diǎn)站不轉(zhuǎn)車相同. 只有在適應(yīng)度函數(shù)與選擇流程中每一個(gè)個(gè)體所對(duì)應(yīng)的適應(yīng)度函數(shù)處有所改變，其函數(shù)為：其中5（分鐘）表示公汽換乘公汽一次耗時(shí)(其中步行時(shí)間2分鐘).除這之外，這一流程中的其他的原理沒變. 從起點(diǎn)站到終點(diǎn)站轉(zhuǎn)兩次車從起點(diǎn)站到終點(diǎn)站轉(zhuǎn)兩次車遺傳算法流程中基于站點(diǎn)序號(hào)的編碼，交叉重組，染色體的變異，種群交叉，迭代的結(jié)束條件和結(jié)果的原理與從起點(diǎn)站到終點(diǎn)站不轉(zhuǎn)車相同. 只有在適應(yīng)度函數(shù)與選擇流程中每一個(gè)個(gè)體所對(duì)應(yīng)的適應(yīng)度函數(shù)處有所改變，其函數(shù)為： 其中（分鐘）表示公汽換乘公汽兩次耗時(shí)(其中步行時(shí)間分鐘).除這之外，這一流程中的其他的原理沒變. 模型求解 由程序運(yùn)行最終得出：找不一條路線使從起點(diǎn)站到終點(diǎn)站不轉(zhuǎn)車. 從起點(diǎn)站到終點(diǎn)站轉(zhuǎn)一次車 。clc。 a2=data(find(data(:,1)==i1)+2,i2)。 (i4i2)==1 i3(a2,a1)=3。 a4=data(find(data(:,1)==i1)+3,i2)。 (i4i2)==1 i3(a3,a4)=3。 i3(j1,j2)~=3 i3(j1,j2)=inf。M=30。 %隨機(jī)產(chǎn)生一個(gè)由自然數(shù)1到num 組成的全排列endpopm2=zeros(M,num)。a3=[]。%群體2中，不同個(gè)體a1和a2的站點(diǎn)差enda5=3.*a3+3956*3。a7=(a5./sum(a5)).*1000。for k4=1:30 a11=a11+a7(k4)。for k5=1:30 a13=a13+a8(k5)。%種群2產(chǎn)生30個(gè)11000的均勻分布的隨機(jī)數(shù)a17=[]。elseif a12(2)(a15(k6)) amp。elseif a12(4)(a15(k6)) amp。elseif a12(6)(a15(k6)) amp。elseif a12(8)(a15(k6)) amp。elseif a12(10)(a15(k6)) amp。 elseif a12(12)(a15(k6)) amp。elseif a12(14)(a15(k6)) amp。elseif a12(16)(a15(k6)) amp。elseif a12(18)(a15(k6)) amp。elseif a12(20)(a15(k6)) amp。elseif a12(22)(a15(k6)) amp。elseif a12(24)(a15(k6)) amp。elseif a12(26)(a15(k6)) amp。elseif a12(28)(a15(k6)) amp。endenda18=[]。elseif a14(2)(a16(k7)) amp。elseif a14(4)(a16(k7)) amp。elseif a14(6)(a16(k7)) amp。elseif a14(8)(a16(k7)) amp。elseif a14(10)(a16(k7)) amp。 elseif a14(12)(a16(k7)) amp。elseif a14(14)(a16(k7)) amp。elseif a14(16)(a16(k7)) amp。elseif a14(18)(a16(k7)) amp。elseif a14(20)(a16(k7)) amp。elseif a14(22)(a16(k7)) amp。elseif a14(24)(a16(k7)) amp。elseif a14(26)(a16(k7)) amp。elseif a14(28)(a16(k7)) amp。endenda19=[,]。end a22=[,]。p1=min(fix(unifrnd(1,3957,1,2)))。 for j=p2p1+2:3957+p2p1+1 for i=1:p2p1+1 if(a1(j)==a1(i)) a1(j)=nan。 elseif r1~=1 a23(r12,:)=a1(~isnan(a1))。b22=[,]。p1=min(fix(unifrnd(1,3957,1,2)))。 for j=p2p1+2:3957+p2p1+1 for i=1:p2p1+1 if(a1(j)==a1(i)) a1(j)=nan。 elseif r2~=1 b23(r22,:)=a1(~isnan(a1))。w1=[,]。z3=a2(q,z1)。w1(q,:)=[a2(q,1:z11),z5,z3,a2(q,(z1+2):3957)]。for q=16:30z1=fix(unifrnd(1,3957,1,2))。z21=max(z2)。z4=b2(q,z22)。endw5=[w1。%a22=[,]。p1=min(fix(unifrnd(1,3957,1,2)))。 for j=p2p1+2:3957+p2p1+1 for i=1:p2p1+1 if(a1(j)==a1(i)) a1(j)=nan。 elseif r1~=1 a23(r12,:)=a1(~isnan(a1))。b22=[,]。p1=min(fix(unifrnd(1,3957,1,2)))。 for j=p2p1+2:3957+p2p1+1 for i=1:p2p1+1 if(a1(j)==a1(i)) a1(j)=nan。 elseif r2~=1 b23(r22,:)=a1(~isnan(a1))。%w1=[,]。z3=a2(q,z1)。w1(q,:)=[a2(q,1:z11),z5,z3,a2(q,(z1+2):3957)]。for q=16:30z1=fix(unifrnd(1,3957,1,2))。z21=max(z2)。z4=b2(q,z22)。endw5=[