太阳能小屋的设计

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