Scattering by Conducting Cylinders.rar_多目标_多目标优化_遗传算法

function [curr]=scat(a,N) % a is the size of the square % % N is the number of the segment, it should be times of 8 % % the origin is located in the center of the square % segment=8*a/N; freq=8e+6; omiga=2*pi*freq; c=3e+8; mu=4*pi*(1e-7); namda=c/freq; k=2*pi/namda; k*a; x=ones(N,1); y=ones(N,1); for index=1:N if index<=(N/8) x(index)=a; y(index)=(index-1/2)*segment; elseif index>(N/8) & index<=(3*N/8) x(index)=a-(index-N/8-1/2)*segment; y(index)=a; elseif index>(3*N/8) & index<=(5*N/8) x(index)=-a; y(index)=a-(index-3*N/8-1/2)*segment; elseif index>(5*N/8) & index<=(7*N/8) x(index)=-a+(index-5*N/8-1/2)*segment; y(index)=-a; else x(index)=a; y(index)=-a+(index-7*N/8-1/2)*segment; end end l=zeros(N,N); for m=1:N for n=1:N dist=sqrt((x(m)-x(n)).^2+(y(m)-y(n)).^2); if m~=n l(m,n)=omiga*mu*besselh(0,2,k*dist)*segment/4; else l(m,n)=omiga*mu*segment*(1-j*2*log(1.781*k*segment/(4*2.718))/pi)/4; end end end g=zeros(N,1); for m=1:N g(m)=exp(-j*k*x(m)); end cur=inv(l)*g; plot([1:N],abs(cur));