calculateobjectvalue.m

来自「对于天线的设计和应用进行了仿真,并进行了优化处理」· M 代码 · 共 45 行

function[Fitness]=CalculateObjectValue(X,M,N)
Max=25;
Index=zeros(1,M);
Fmin=ones(1,M).*Max;
Value=zeros(M,101);
s=0;
for j=1:M
    s=0;
    for u=0:0.01:1
        s=s+1;
        Value(j,s)=1+exp(i*100*pi*u);
        for n=1:N
            Value(j,s)=Value(j,s)+exp(i*(pi*X(j,n)*u));
        end
        Value(j,s)=abs(Value(j,s));
        if Fmin(j)<Value(j,s)
            Index(j)=u;
            break;
        end
        Fmin(j)=Value(j,s);
    end
end

Rmax=zeros(1,M);
for j=1:M
    s=0;
    for u=Index(j):0.01:1
        s=s+1;
        Value(j,s)=1+exp(i*100*pi*u);
        for n=1:N
            Value(j,s)=Value(j,s)+exp(i*(pi*X(j,n)*u));
        end
        Value(j,s)=abs(Value(j,s));
        if Rmax(j)<Value(j,s)
            Rmax(j)=Value(j,s);
        end
    end
end
Fitness=20*log10(Rmax./Max);