矩量法求解无限长理想导体圆柱的电流分布和双站RCS

题目：矩量法求解TE极化平面波0度入射半径a=1个波长的无限长理想导体圆柱的电流分布和双站RCS程序：%MOM___TE.m%矩量法求解TE极化平面波0度入射半径a=1个波长%的无限长理想导体圆柱的电流分布和双站RCS%矩形基函数和点匹配法clear;clc;t1=cputime;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%1.常数定义%Incident_Angle=0;%取入射波与x轴夹角为0度N=128;%将圆柱圆周分成N段c0=3.0e+8;%光速f=3.0e+8;%入射波频率300Mradius=c0/f;%二维无限长圆柱半径为1个波长step=2*pi*radius/N;%每段长度%eta=120*pi;%真空中的波阻抗k=(2*pi/c0)*f;%波数Z_coeff=j*k*step/4.0;%矩阵元素前系数,除对角元素%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2.计算所分各段中点对应的坐标(圆心为原点)fori=1:Nx(i)=radius*cos((i-0.5)*2*pi/N);y(i)=radius*sin((i-0.5)*2*pi/N);end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%3.计算阻抗矩阵(根据公式可以推得矩阵是对称的，并且主对角元素值相等）for(m=1:N)for(n=1:m)ifm==nZ(m,n)=1/2;elser=sqrt((x(m)-x(n))^2+(y(m)-y(n))^2);Z(m,n)=Z_coeff*(x(n)*(x(m)-x(n))+y(n)*...(y(m)-y(n)))/(radius*r)*besselh(1,2,k*r);Z(n,m)=Z(m,n);endendend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%4.计算激励向量和未知电流系数向量form=1:NH(m)=-exp(-j*k*x(m));end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5.计算未知电流系数向量a=Z\H';%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%6.计算双站RCSScattering_Angle=linspace(0,360,N);%将散射角从0-360°取N个点radian=Scattering_Angle*pi/180+pi;M=length(radian);fort=1:M%求不同散射角对应的RCSsum=0;forp=1:Nff=x(p)*cos(radian(t))+y(p)*sin(radian(t));sum=sum+a(p)*ff*exp(j*k*ff)*step/radius;endsum(t)=(abs(sum))^2;rcs(t)=k/4*sum(t);rcs(t)=10*log10(rcs(t));end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%7.作图：表面电流；双站RCSsubplot(211);plot(Scattering_Angle,abs(a));title('无限长二维理想导体圆柱表面电流与角度φ关系');legend('MOM');xlabel('角度φ(°)');ylabel('电流密度');subplot(212);plot(Scattering_Angle,rcs);title('无限长二维理想导体圆柱RCS与散射角φ关系');legend('MOM');xlabel('散射角φ(°)');ylabel('RCS/dbm');t2=cputime;t=t2-t1matlab结果图：