a1.asv

对于天线的设计和应用进行了仿真,并进行了优化处理

  end
end
Index
Lmin=MaxValue;
LIndex=zeros(M,1);
Value=zeros(M,1);
Value1=ones(M,1);
% find the first minimal value on the left
for j=1:M
    for f=Index(j,1):-0.01:-1
           for n=1:N
                Value1(j,1)=Value1(j,1)+exp(i*(k*f*n+Phi(j,n)));
            end
            
            Value(j,1)=abs(Value1(j,1));
            if Value(j,1)<=Lmin(j,1)
               Lmin(j,1)=Value(j,1);
           else
               break;
           end
            Value1(j,1)=1;
            LIndex(j,1)=f;          
    end
end


Value=ones(M,1);
Lmax=zeros(M,1);
for j=1:M
    for f=-1:0.01:LIndex(j,1)
        for n=1:N
            Value(j,1)=Value1(j,1)+exp(i*(k*f*n+Phi(j,n)));
        end
        Value(j,1)=abs(Value(j,1));
        if Lmax(j,1)<Value(j,1)
            Lmax(j,1)=Value(j,1);
        end
        Value(j,1)=1;
    end
end

Value=ones(M,1);
Rmax=zeros(M,1);
for j=1:M
    for f=-1:0.01:RIndex(j,1)
        for n=1:N
            Value(j,1)=Value1(j,1)+exp(i*(k*f*n+Phi(j,n)));
        end
        Value(j,1)=abs(Value(j,1));
        if Rmax(j,1)<Value(j,1)
            Rmax(j,1)=Value(j,1);
        end
        Value(j,1)=1;
    end
end

% the Fitness value
LV=Lmax./MaxValue;
RV=Rmax./MaxValue;
dBLV=10*log10(LV);
dBRV=10*log10(RV);
V=ones(M,1);
Fitness=V.*28-((max(0,(dBLV+17))+max(0,(dBRV+11))));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%