matlab编程代码求 2,附录2求解RDC的 MATLAB 编程源代码

《附录2求解RDC的 MATLAB 编程源代码》由会员分享，可在线阅读，更多相关《附录2求解RDC的 MATLAB 编程源代码(3页珍藏版)》请在人人文库网上搜索。

1、附录2求解RDC的 MATLAB 编程源代码附录：求解RDC的 MATLAB 编程源代码 程序用MATLAB 6.0 编制而成的，基于遗传算法的程序，实现了求解配送模型的过程。 1.适值函数 function sol，eval=fitness(sol，options) d=0 8 12 15 18 40 20 32 16； 8 0 13 8 20 10 15 22 20； 12 13 0 15 20 20 15 15 15； 15 8 15 0 20 10 18 18 30； 18 20 20 20 0 20 15 15 20； 40 10 20 10 20 0 14 18 15； 20 15。

2、 15 18 15 14 0 14 20； 32 22 15 18 15 18 14 0 20； 16 20 15 30 20 15 20 20 0； G=0 3 5 2 6 4 4 1 3； s=rand(size(G) T=0 2 2 1 2.5 1.5 2 0.5 1.5； A=0 6 6 7 8 7 7 7 9； B=0 9 8 9 12 11 10 10 12； sum2=0； sum3=0； sum4=0； for k=1:4 for i=1:9 forj=1:9 if (sum2=0)&(c2=0) fsum=0； for i=1:9 for j=1:9 fsum=fsum+d(。

3、i，j)*xijk； end end eval=fsum； else eval=-10000； end eval=-eval； 2.遗传算法求解的MATLAB 代码 %the matrix for distance d=0 8 12 15 18 40 20 32 16； 8 0 13 8 20 10 15 22 20； 12 13 0 15 20 20 15 15 15； 15 8 15 0 20 10 18 18 30； 18 20 20 20 0 20 15 15 20； 40 10 20 10 20 0 14 18 15； 20 15 15 18 15 14 0 14 20； 32 22。

4、 15 18 15 18 14 0 20； 16 20 15 30 20 15 20 20 0； G=0 3 5 2 6 4 4 1 3； T=0 2 2 1 2.5 1.5 2 0.5 1.5； A=0 6 6 7 8 7 7 7 9； B=0 9 8 9 12 11 10 10 12； bounds=2；2；2；2；2；2；2；2；2； flag=0； while flagelse ykj=0； end sum2=sum2+G(j)*ykj； sum4=sum4+ykj； end end end s(i)=rand(size(G) c1=s(i)-A(i)； c2=B(i)-s(i)； i。

5、f (sum2=0)&(c2=0) flag=flag+1； initPop(flag，:)=init； else Continue end end p endPop bpop traceJ=ga(bounds，fitness，initPop，1e-5 1 1，. maxGeomSelect，0.08，arithXover，20，nonUnifMutate，2 1 3)； 3.初始化函数的MATLAB 代码 function pop=initializega(num， bounds， evalFN，evalOps，options) defaults 1e-6 1 if nargin=maxGe。

6、n； 5.基于正态分布的序列选择函数MATLAB 代码 functionnewPop=normGeomSelect(oldPop，options) q=options(2)； e=size(oldPop，2)； n=size(oldPop，1)； newPop=zeros(n，e)； fitp2，bounds，Ops) c1=pl*a+p2*(1-a)； c2=pl*(1-a)+p2*a； 7.变异函数 function parent=nonUnifMutate(parent，bounds，Ops) cg=Ops(1)； mg=Ops(3)； b=Ops(4)； df=bounds(:，2)-。

7、bounds(:，1)； numVar=size(parent，2)-1； mPoint=round(rand*(num Var-1)+1； and=round(rand)； if md newValue=parent(mPoint)+delta(cg，mg，bounds(mPoint，2)-parent(mPoint)，b)； else newValue=parent(mPoint)-delta(cg，mg，parent(mPoint)-bounds(mPoint，1)，b)； end parent(mPoint)=newValue；%Make the child 8.主运行函数GA 的MA。

8、TLAB 代码 function x，endPop，bPop，traceInfo=ga(bounds，evalFN，evalOps，startPop，opts，. termFN，termOps，selectFN，selectOps，xOverFNs，xOverOps，mutFNs，mutOps) n=nargm； if n2|n=6|n=10|n=12 disp(lnsufficient arguements) end if n3 evalOps=； end if n5 opts=1e-6 1 0； end if isempty(opts) opts=1e-6 1 0； end if any(evalFN48) if opts(2)=1 e1str=x=c1； c1(xZomeLength)=， evalFN； e2str=x=c2； c2(xZomeLength)=， evalFN； else elstr=x=b2f(endPop(j，:)，bounds，bits)； endPop(j，xZomeLength)=，. eval FN ； end else。