close all
clear all
clc
% set parameter
M=80;
M2=10;
M3=M-M2;
T=500;
Pc=0.6;
Pm=0.004;
N=7;
Len=10;
Lenth=N*Len;
Phi=zeros(M,Lenth);
BestValue=zeros(1,T);
% Generater the first generation of the population
Ab1=round(rand(M,Lenth));
Phi=DecodeChro(Ab1,M,Lenth,N,Len);
Phi=pi*Phi/2046;
figure(1)
mesh(Phi)
MinFitness=zeros(1,T);
AveFitness=zeros(1,T);
for Gen=1:T
    
    
% Calculate Object fitness
Phi=DecodeChro(Ab1,M,Lenth,N,Len);
Phi=pi*Phi/2046;
Fitness=CalculateObjectValue(Phi,M,N);
Sum=sum(Fitness);
AveFitness(Gen)=Sum/M;
% select operation
[Cost,ind]=sort(Fitness);
MinFitness(Gen)=Fitness(ind(1));
Ab2=Ab1(ind(1:M2),:);
Ab4=Ab1(ind(1:M3),:);
% clone operation
Ab3=[Ab2;Ab2;Ab2;Ab2;Ab2;Ab2;Ab2;Ab2;Ab2;Ab2];
Ms=10*M2;
% mutation operation
Pm=-(0.7/T)*Gen+0.8;
for j=1:Ms
        for k=1:Lenth
        Randp=rand(1);
        if Randp<Pm
            if Ab3(j,k)==1
                Ab3(j,k)=0;
            else
                Ab3(j,k)=1;
            end
        end
    end
end
Phi3=DecodeChro(Ab3,Ms,Lenth,N,Len);
Phi3=pi*Phi3/2046;
Fitness3=CalculateObjectValue(Phi3,Ms,N);
Sum3=sum(Fitness3);
% select from ab3
[Cost3,ind3]=sort(Fitness3);
Ab5=Ab3(ind3(1:M2),:);
% generate new colony
Ab1=[Ab4;Ab5];
end
Phi=DecodeChro(Ab1,M,Lenth,N,Len);
Phi=pi*Phi/2046;
figure(2)
mesh(Phi)
figure(3)
hold on
plot( MinFitness,'r'),
plot(AveFitness,'b');
hold off

BestPhi=zeros(1,N);
BestPhi=Phi(ind(1),:);
dr=0.8;
k=2*pi*dr;
max=0;
s=0;
Value=ones(1,201);
    for f=-1:0.01:1
        s=s+1;
        for n=1:N
            Value(1,s)=Value(1,s)+exp(i*(k*f*n+BestPhi(1,n)));
        end
        Value(1,s)=abs(Value(1,s));
        if max<Value(1,s)
            max=Value(1,s);
        end
    end
    figure(4)
    f=-1:0.01:1;
    Value=Value./max;
    Value=20*log10(Value);
    plot(f,Value),grid;