garadareval2.m

运用遗传算法的工具箱解决宽带雷达吸波结构的优化问题

function [sol,val]=gaRadarEval(sol,options)
for l=1:3
 n(l)=round(sol(2*l-1));  
 b(l)=sol(2*l);
end
b=0.001*b;
d=sum(b);
t=zeros(1,3);
if d<0.002
  t=b;  
end 
t;
%反射系数随频率的变化关系
for p=1:20
    f=0.1+1*p;
     for i=1:3       
 switch n(i)
    case 1
        u(1)=1+j*0;  e(1)=10+j*0;
        kk(i)=u(1);  tt(i)=e(1);
    case 2
        u(2)=1+j*0;  e(2)=50+j*0;
        kk(i)=u(2);  tt(i)=e(2);
    case 3
        u(3)=5/f.^0.974-j*10/f.^0.961;  e(3)=15+j*0;
        kk(i)=u(3);  tt(i)=e(3);
    case 4
        u(4)=3/f.^1.000-j*15/f.^0.975; e(4)=15+j*0;
        kk(i)=u(4);  tt(i)=e(4);
    case 5
        u(5)=7/f.^1.000-j*12/f.^1.000;  e(5)=15+j*0;
        kk(i)=u(5);  tt(i)=e(5);
    case 6
        u(6)=1+j*0;  e(6)=5/f.^0.861-j*8/f.^0.569;
        kk(i)=u(6);  tt(i)=e(6);
    case 7
        u(7)=1+j*0;  e(7)=8/f.^0.778-j*10/f.^0.682;
        kk(i)=u(7);  tt(i)=e(7);
    case 8
        u(8)=1+j*0; e(8)=10/f.^0.778-j*6/f.^0.861;
        kk(i)=u(8);  tt(i)=e(8);
    case 9
        u(9)=35*0.64/(f.^2+0.64)-j*35*0.8*f/(f.^2+0.64);  e(9)=15+j*0;
        kk(i)=u(9);  tt(i)=e(9);
    case 10
        u(10)=35*0.25/(f.^2+0.25)-j*35*0.5*f/(f.^2+0.25);  e(10)=15+j*0;
        kk(i)=u(10);  tt(i)=e(10);
    case 11
        u(11)=30*1.00/(f.^2+1.00)-j*30*1.00*f/(f.^2+1.00);  e(11)=15+j*0;
        kk(i)=u(11);  tt(i)=e(11);
    case 12
        u(12)=18*0.25/(f.^2+0.25)-j*18*0.5*f/(f.^2+0.25);   e(12)=15+j*0;
        kk(i)=u(12);  tt(i)=e(12);
    case 13
        u(13)=20*2.25/(f.^2+2.25)-j*20*1.5*f/(f.^2+2.25);   e(13)=15+j*0;
        kk(i)=u(13);  tt(i)=e(13);
    case 14
        u(14)=30*6.25/(f.^2+6.25)-j*30*2.5*f/(f.^2+6.25);   e(14)=15+j*0;
        kk(i)=u(14);  tt(i)=e(14);
    case 15
        u(15)=30*4.00/(f.^2+4.00)-j*30*2.0*f/(f.^2+4.00);   e(15)=15+j*0;
        kk(i)=u(15);  tt(i)=e(15);
    case 16
        u(16)=25*12.25/(f.^2+12.25)-j*25*3.5*f/(f.^2+12.25);  e(16)=15+j*0;
        kk(i)=u(16);  tt(i)=e(16);
    otherwise
        input error
     end 
    end
    u=kk;
    e=tt;
    u;
    e;
N=3;
    for k=1:N    
    r(k)=j*20*pi*f*sqrt(u(k)*e(k))/3;
     zo=120*pi;
    zcm(k)=zo*sqrt(u(k)/e(k)); 
end
r;
zcm;
for i=1:N
   if i<2
   Zin(i)=zcm(i)*tanh(r(i)*t(i));
   else
   Zin(i)=zcm(i)*(Zin(i-1)+zcm(i)*tanh(r(i)*t(i)))/(zcm(i)+Zin(i-1)*tanh(r(i)*t(i)));
end
end
Zin;
    R(p)=(Zin(N)-zo)/(Zin(N)+zo);
    y(p)=1-abs(R(p));
end
val=max(y);