背包问题基于量子免疫算法求解背包问题matlab源码
Posted Matlab走起
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了背包问题基于量子免疫算法求解背包问题matlab源码相关的知识，希望对你有一定的参考价值。
一、简介
1 量子免疫克隆
-
clear; C=[253 245 243 239 239 239 238 238 237 232 231 231 230 229 228 227 224 217 213 207 203 201 195 194 191 187 187 177 175 171 169 168 165 164 161 160 158 150 149 147 141 140 139 136 135 132 128 126 122 120 119 116 116 114 111 110 105 105 104 103 93 92 90 79 78 77 76 76 75 73 62 62 61 60 60 59 57 56 53 53 51 50 44 44 42 42 38 36 34 28 27 24 22 18 12 10 7 4 4 1]; W=[253 245 243 239 239 239 238 238 237 232 231 231 230 229 228 227 224 217 213 207 203 201 195 194 191 187 187 177 175 171 169 168 165 164 161 160 158 150 149 147 141 140 139 136 135 132 128 126 122 120 119 116 116 114 111 110 105 105 104 103 93 92 90 79 78 77 76 76 75 73 62 62 61 60 60 59 57 56 53 53 51 50 44 44 42 42 38 36 34 28 27 24 22 18 12 10 7 4 4 1]; V=6666; a=2.40; q=(sqrt(1/2)).*ones(20,100); [N,L] = size(q); ger =500; pm=0.0088; pc=0.88; fat=0.1; disp(sprintf('Number of generations: %d',ger)); disp(sprintf('Population size: %d',N/2)); % General parameters & Initial operations vmfit = []; it = 1; vx = []; T=[]; vw=[]; f='dot(v,C1,2)'; updatef=0; % Generations t0 = clock; c1=0; while it <= ger itemp=[1 3 5 7 9 11 13 15 17 19]; v=observe(q(itemp,:)); rate=C./W; [a,ind] = sort(rate); for j=1:N/2 i1=1; %repair knoverfilled=0; if dot(v(j,:),W)>6666 knoverfilled=1; end while knoverfilled==1 v(j,ind(i1))=0; if dot(v(j,:),W)<=6666 knoverfilled=0; end i1=i1+1; end v(j,:)=v(j,:); end %and/or Crossver for i=1:N/2 cind(i)=i; end for i=1:N/2 point=unidrnd(N/2-i+1); temp=cind(i); cind(i)=cind(i+point-1); cind(i+point-1)=temp; end for i=1:2:N/2 p=rand(1); if(p<pc) c1=and(v(cind(i),:),v(cind(i+1),:)); c2=or(v(cind(i),:),v(cind(i+1),:)); v(cind(i),:)=c1; v(cind(i+1),:)=c2; end end for j=1:N/2 i1=1; %repair knoverfilled=0; if dot(v(j,:),W)>6666 knoverfilled=1; end while knoverfilled==1 v(j,ind(i1))=0; if dot(v(j,:),W)<=6666 knoverfilled=0; end i1=i1+1; end v(j,:)=v(j,:); end C1=ones(N/2,1)*C; fit=eval(f); [sol,indb] = max(fit); maxf=sol; max1=v(indb,:); minf=min(fit); if updatef>=sol maxf=updatef; max1=updatec; else maxf=sol; max1=v(indb,:); end %量子变异 for i=1:2:size(q,1) for j=1:size(q,2) if v((i+1)/2,j)==max1(1,j) if max1(1,j)==0 s=0; angle=0; else if q(i,j)*q(i+1,j)>0 s=1; angle=0.015; end if q(i,j)*q(i+1,j)<0 s=-1; angle=0.015; end if q(i,j)==0 s=0; angle=0.015; end if q(i+1,j)==0 s=1; angle=0.015; end end elseif v((i+1)/2,j)==0 if maxf>fit((i+1)/2) s=0; angle=0; else if q(i,j)*q(i+1,j)>0 s=-1; angle=0.015; end if q(i,j)*q(i+1,j)<0 s=1; angle=0.015; end if q(i,j)==0 s=1; angle=0.015; end if q(i+1,j)==0 s=0; angle=0.015; end end else if maxf>fit((i+1)/2) if q(i,j)*q(i+1,j)>0 s=-1; angle=0.015; end if q(i,j)*q(i+1,j)<0 s=1; angle=0.015; end if q(i,j)==0 s=1; angle=0.015; end if q(i+1,j)==0 s=0; angle=0.015; end else if q(i,j)*q(i+1,j)>0 s=1; angle=0.015; end if q(i,j)*q(i+1,j)<0 s=-1; angle=0.015; end if q(i,j)==0 s=0; angle=0.015; end if q(i+1,j)==0 s=1; angle=0.015; end end end o=s*angle*pi; u=[cos(o) -sin(o);sin(o) cos(o)]; temp1=[q(i,j);q(i+1,j)]; temp2=u*temp1; q(i,j)=temp2(1,1); q(i+1,j)=temp2(2,1); end三、运行结果
四、备注
完整代码或者仿真咨询添加QQ1575304183
以上是关于背包问题基于量子免疫算法求解背包问题matlab源码的主要内容，如果未能解决你的问题，请参考以下文章
背包问题基于matlab遗传算法结合贪婪算法求解背包问题含Matlab源码 791期
背包问题基于matlab带权重的贪心萤火虫算法求解0-1背包问题含Matlab源码 045期
背包问题基于matlab带权重的贪心萤火虫算法求解0-1背包问题含Matlab源码 045期
背包问题基于matlab萤火虫算法求解背包问题含Matlab源码 1440期