MATLAB.rar_OPT MATLAB_can matlab_ga sa tabu sworm opt_tabu-searc

clc; m=input('enter number of jobs'); n=input('enter number of operations for each job'); pop=input('enter no of chromosomes in population'); crp=input('enter crossover rate'); mrp=input('enter mutation rate'); o=m*n; a=zeros(pop,o); for i=1:pop x=zeros(1,m); j=1; while(j<=o) out=randi([1,m],1,1); x(1,out)=x(1,out)+1; if (x(1,out)<=n) a(i,j)=out; j=j+1; end; end; end; disp (a); p=zeros(m,n); for i=1:m for j=1:n p(i,j)=input('enter precedence matrix\bp(i,j)\b'); end; end; disp (p); t=zeros(m,n); for i=1:m for j=1:n t(i,j)=input('enter processing time matrix\bp(i,j)\b'); end; end; disp (t); b=a; count1=zeros(pop,m);% this gives number of operations completed on job count2=zeros(pop,n);% this gives number of operations completed on machine st=zeros(m,n);% In this matrix,each row gives starting times of a job ct=zeros(m,n);%In this matrix,each row gives completion times of a job ctm=zeros(n,m);%in this matrix, each row gives completion times of machine stm=zeros(n,m);%in this matrix, each row gives starting times of machine mak=zeros(1,pop); for j=1:pop for i=1:o if (count1(j,(b(j,i)))<=n) count1(j,b(j,i))=count1(j,b(j,i))+1; end; if ((count2(j,(p((b(j,i)),count1(j,b(j,i))))))<=m) count2(j,(p(b(j,i),count1(j,b(j,i)))))=count2(j,(p(b(j,i),count1(j,b(j,i)))))+1; end; if (count1(j,b(j,i))==1&&count2(j,(p((b(j,i)),count1(j,b(j,i)))))==1) st(b(j,i),count1(j,b(j,i)))=0; ct(b(j,i),count1(j,b(j,i)))=st(b(j,i),count1(j,b(j,i)))+t(b(j,i),count1(j,b(j,i))); else if(count1(j,b(j,i))==1) st(b(j,i),count1(j,b(j,i)))= ctm(p(b(j,i),count1(j,b(j,i))),(count2(j,(p(b(j,i),count1(j,b(j,i))))))-1); ct(b(j,i),count1(j,b(j,i)))=st(b(j,i),count1(j,b(j,i)))+t(b(j,i),count1(j,b(j,i))); else if (count2(j,(p(b(j,i),count1(j,b(j,i)))))==1) st(b(j,i),count1(j,b(j,i)))=max(ct(b(j,i),(count1(j,b(j,i))-1)),ctm(p(b(j,i),count1(j,b(j,i))),count2(j,(p(b(j,i),count1(j,b(j,i))))))); ct(b(j,i),count1(j,b(j,i)))=st(b(j,i),count1(j,b(j,i)))+t(b(j,i),count1(j,b(j,i))); else st(b(j,i),count1(j,b(j,i)))=max(ct(b(j,i),(count1(j,b(j,i))-1)),ctm(p(b(j,i),count1(j,b(j,i))),(count2(j,(p(b(j,i),count1(j,b(j,i))))))-1)); ct(b(j,i),count1(j,b(j,i)))=st(b(j,i),count1(j,b(j,i)))+t(b(j,i),count1(j,b(j,i))); end; end; end; if (count1(j,b(j,i))==1&&count2(j,(p((b(j,i)),count1(j,b(j,i)))))==1) stm((p(b(j,i),count1(j,b(j,i)))),(count2(j,(p(b(j,i),count1(j,b(j,i)))))))=0; ctm((p(b(j,i),count1(j,b(j,i)))),(count2(j,(p(b(j,i),count1(j,b(j,i)))))))=stm((p(b(j,i),count1(j,b(j,i)))),count2(j,(p(b(j,i),count1(j,b(j,i))))))+t(b(j,i),count1(j,b(j,i))); else if(count1(j,b(j,i))==1) stm((p(b(j,i),count1(j,b(j,i)))),(count2(j,(p(b(j,i),count1(j,b(j,i))))))) =ctm(p(b(j,i),count1(j,b(j,i))),(count2(j,(p(b(j,i),count1(j,b(j,i)))))-1)); ctm((p(b(j,i),count1(j,b(j,i)))),(count2(j,(p(b(j,i),count1(j,b(j,i)))))))=stm((p(b(j,i),count1(j,b(j,i)))),count2(j,(p(b(j,i),count1(j,b(j,i))))))+t(b(j,i),count1(j,b(j,i))); else if (count2(j,(p(b(j,i),count1(j,b(j,i)))))==1) stm((p(b(j,i),count1(j,b(j,i)))),(count2(j,(p(b(j,i),count1(j,b(j,i)))))))=ct(b(j,i),(count1(j,b(j,i))-1)); ctm((p(b(j,i),count1(1,b(j,i)))),(count2(j,(p(b(j,i),count1(j,b(j,i)))))))=stm((p(b(j,i),count1(j,b(j,i)))),count2(j,(p(b(j,i),count1(j,b(j,i))))))+t(b(j,i),count1(j,b(j,i))); else stm((p(b(j,i),count1(j,(b(j,i))))),(count2(j,(p(b(j,i),(count1(j,b(j,i))))))))=max(ct((b(j,i)),((count1(j,b(j,i)))-1)),ctm((p(b(j,i),count1(j,b(j,i)))),(count2(j,(p((b(j,i)),(count1(j,b(j,i))))))-1))); ctm((p(b(j,i),count1(j,b(j,i)))),count2(j,(p((b(j,i)),count1(j,b(j,i))))))=stm((p(b(j,i),count1(j,b(j,i)))),count2(j,(p(b(j,i),count1(j,b(j,i))))))+t(b(j,i),count1(j,b(j,i))); end; end; end; end; mak(1,j)=max(max(ct)); end; disp(mak); b(:,o+1)=mak; %disp(b); c=zeros(pop,o+1); %tournament selection for i=1:pop r1=randi([1,pop],1,1); r2=randi([1,pop],1,1); while(r1==r2) r1=randi([1,pop],1,1); r2=randi([1,pop],1,1); continue; end if(b(r1,o+1)<b(r2,o+1)) c(i,:)=b(r1,:); end c(i,:)=b(r2,:); end; end; display(c); %cross over operator d=zeros(1,2); e=zeros(2,o); j=1; accept=0; for i=1:pop ra=rand(0,1); if(ra<=crp) d(1,j)=i; j++; accept=accept+1 end if (accept==2) e(1,:)=c(d(1,1),:); e(2,:)=c(d(2,2),:); crossover(e(1,:),e(2,:)); end disp(d);