第二階段題特等1013程序附件

1GM(1,1)pp>0.95C<0.351pp>0.80C<0.502pp>0.70C<0.653pp<=0.70C>=0.654%%%%function%x=[3.7113.7233.7163.721%x=[2.673.133.253.363.56%x=[2.973.233.293.463.59%x=[43.4547.0552.7557.1462.64%%7.328.53%%%IIi=[5678910% %formatlongforfori=1:n-fori=2:n%預(yù)測第I%XI1=xua*exp((- %51fori=1:nfor%53計算絕對誤差及相對誤差序列fprintf(%fprintf('%f%%%forfprintf(%fprintf('%f%%fprintf('%f%%fprintf('%f%for%檢驗預(yù)測精度fprintf('pp>0.95或C<0.35為1\npp>0.80或C<0.50為2級(合格)\npp>0.70或C<0.65為3級(勉強合格)\npp<=0.70或C>=0.65為4if pp>0.95ifpp<=0.95&pp>0.80 pp<=0.95&pp>0.80ifpp<=0.80&pp>0.70 pp<=0.80&pp>0.70if pp<=0.70ifC<0.35 if if if ifpdd==1pddstr='1（好 fprintf('1（好）\n\n');ifpddstr='2（合格）' fprintf('2（合格）\n\n');ifpddstr='3（勉強合格 fprintf('3（勉強合格）\n\n');ifpddstr='4（不合格 fprintf('4（不合格 fori=1:n1 fprintf('預(yù)測結(jié)果：第%d個值(原數(shù)據(jù)的第%d個序列號)的預(yù)測結(jié)果為： for%fprintf('%f%fprintf(%fprintf('%f%fork=0:n-%%for%%endfprintf('');fork=1:n1 fprintf('預(yù)測結(jié)果：第%d個值(原數(shù)據(jù)的第%d個序列號)的預(yù)測值為：%四說明:fga.m為遺傳算法的主程序;采用二進制Gray編碼,采用基于賭法的非線性選擇,均勻交叉,變異操作,而且還引入了倒位操作!%%Findsaumofafunctionofseveral%fmaxgasolvesproblemsofthe maxF(X)subjectto:LB<=X<=% -最優(yōu)的群體即為最優(yōu)的% -最佳所對應(yīng)的目標(biāo)函數(shù)% % % % % % % % ifnargin<3,error('FMAXGArequiresatleastthreeinputarguments');endifnargin==3,ns=[01e-4];endifnargin==4,ifnargin==5,pCross=0.8;pMutation=0.1;pInversion=0.15;options=[0ifnargin==6,pMutation=0.1;pInversion=0.15;options=[01e-4];endifnargin==7,pInversion=0.15;options=[01e-4];end find((LB- globalmnNewPopchildren1children2 bounds(:,1))'./precision));%由設(shè)定精度劃分區(qū)間forj=1:mpMutation=pm0+(i^4)*(pCross/3-holdon;plot(I,MaxFval,'*');text(I+5,MaxFval,['FMAX='str1=sprintf('%d%s%f\n應(yīng)是：%s',I,num2str(X),MaxFval,num2str(BestPop(I,:)));ififelapsed_time(5)<0)-end%像這種程序當(dāng)然不考慮運行上小時啦str2=sprintf('%d%d%.4ffunction[initpop]=InitPopGray(popsize,bits)initpop=zeros(popsize,len);%Thewholezeroencodingindividualfori=2:popsize-1pop=mod(([0pop]+[pop0]),2);initpop(popsize,:)=ones(1,len);%Thewholeoneencoding%function[fval]= bounds scale=(bounds(:,2)-bounds(:,1))'./(2.^bits-1);%Therangeofthecs=[0cumsum(bits)];fori=1:numV[ywhilenewIn<=mifrNums(newIn)<newfit(fitIn)

functionglobalmnNewPopiflen>2&mod(len,2)==1%如果用來進行交叉的的條數(shù)為奇數(shù)，將其調(diào)整ififopts==0functionglobalnchildren1children2globalnChildren1Childre