1、 1楊金煥.J.太陽能學報.1992,13(1):86-92.2楊剛,CAD計算方法， HYPERLINK /view/4d6935e9aeaad1f346933fc4.html /view/4d6935e9aeaad1f346933fc4.html,2012年9月8日.3陳仕軍，曹炬. HYPERLINK /view/a99c32f5f61fb7360b4c655e.htmL2012 /view/a99c32f5f61fb7360b4c655e.htmL2012年9月8日.4沈輝，曾祖勤.M.北京：,2005.5中國建筑標準研究院.S.北京,2007.附錄二%下面舉例說明遺傳算法%2.1初始

2、化(編碼)functionpop=initpop(popsize,chromlength)pop=round(rand(popsize,chromlength);%2.2計算目標函數值functionpop2=decodebinary(pop)px,py=size(pop);%求pop行和列數fori=1:pypop1(:,i)=2.人(py-i).*pop(:,i);endpop2=sum(pop1,2);%求popl的每行之和%2.2.2將二進制編碼轉化為十進制數(2)%遺傳算法子程序%Name:decodechrom.m%將二進制編碼轉換成十進制functionpop2=decodech

3、rom(pop,spoint,length)pop1=pop(:,spoint:spoint+length-1);pop2=decodebinary(pop1);%2.2.3計算目標函數值functionobjvalue=calobjvalue(pop)temp1=decodechrom(pop,1,10);%將pop每行轉化成十進制數x=temp1*10/1023;%將二值口中的數轉化為變量域的數objvalue=10*sin(5*x)+7*cos(4*x);%計算目標函數值%2.3計算個體的適應值functionfitvalue=calfitvalue(objvalue)globalCmi

4、n;Cmin=0;px,py=size(objvalue);fori=1:pxifobjvalue(i)+Cmin0temp=Cmin+objvalue(i);elsetemp=0.0;endfitvalue(i)=temp;endfitvalue=fitvalue;%2.4選擇復制functionnewpop=selection(pop,fitvalue)totalfit=sum(fitvalue);%求適應值之和fitvalue=fitvalue/totalfit;%單個個體被選擇的概率fitvalue=cumsum(fitvalue);%如fitvalue=1234,貝Ucumsum(f

5、itvalue)=13610px,py=size(pop);ms=sort(rand(px,1);%從小至U大排歹Ufitin=1;newin=1;whilenewin=pxif(ms(newin)fitvalue(fitin)newpop(newin)=pop(fitin);newin=newin+1;elsefitin=fitin+1;endend%2.5交叉functionnewpop=crossover(pop,pc)px,py=size(pop);newpop=ones(size(pop);fori=1:2:px-1if(randpc)cpoint=round(rand*py);ne

6、wpop(i,:)=pop(i,1:cpoint),pop(i+1,cpoint+1:py);newpop(i+1,:)=pop(i+1,1:cpoint),pop(i,cpoint+1:py);elsenewpop(i,:)=pop(i);newpop(i+1,:)=pop(i+1);endend%2.6變異functionnewpop=mutation(pop,pm)px,py=size(pop);newpop=ones(size(pop);fori=1:pxif(randpm)mpoint=round(rand*py);ifmpointbestfitbestindividual=pop(

7、i,:);bestfit=fitvalue(i);endend%2.8主程序clearclfpopsize=20;%群體大小chromlength=10;%字符串長度(個體長度)pc=0.6;%交叉概率pm=0.001;%變異概率pop=initpop(popsize,chromlength);%隨機產生初始群體fori=1:20%20為迭代次數objvalue=calobjvalue(pop);%計算目標函數fitvalue=calfitvalue(objvalue);%計算群體中每個個體的適應度newpop=selection(pop,fitvalue);%復制newpop=crossov

8、er(pop,pc);%交叉newpop=mutation(pop,pc);%變異bestindividual,bestfit=best(pop,fitvalue);%求出群體中適應值最大的個體及其適應值y(i)=max(bestfit);n(i)=i；pop5=bestindividual;x(i)=decodechrom(pop5,1,chromlength)*10/1023;pop=newpop;endfplot(10*sin(5*x)+7*cos(4*x),010)holdonplot(x,y,r*)holdoffzindex=max(y);%計算最大值及其位置x5=x(index)%

9、計算最大值對應的x值y=z附錄三model:sets:guangfuA/1.6/:xi,yi,ni,mi;guangfuB/1.7/:xj,yj,nj,mj;guangfuC/1.11/:xk,yk,nk,mk;endsetsdata:xi,yi,ni=1580808-237.3361956991-392.291580808313.7811651992-294.5581650991-270.9871956991-354.474;xj,yj,nj=1650991344.50271956991383.34971482992616.31561640992322.01741956992778.0116

10、1956992505.697616681000324.7081;xk,yk,nk10671321711521.143814141114939.163514001100808.84714001100898.95531035535.91661518036.24261535571.8476920355107.3233818355107.86891645712618.0148;S=26980000;enddatamax=sum(guangfuA(i):mi*ni)+sum(guangfuB(j):mj*nj)+sum(guangfuC(k):mk*nk)7S=SUM(guangfuA(i):xi*yi*ai*mi)+sum(guangfuB(j):xj*yj