Desktop_矩量法求偶极子表面电流分布_

clc; clear; %基本参数设置 c=3*10^8; epslong=(-338.48-25.2j)*10^(-9)/(36*pi);% miu=1; lambda=0.693; omiga=2*pi*c/lambda; k=2*pi*lambda; L=0.5*lambda;%0.1450 a=0.006; xishu=j/(4*pi*omiga*epslong); %半波阵子天线的网格划分 Nz=5; dz=L/Nz; Ez_in=zeros(Nz,1); b=zeros(Nz,1); z=-L/2+dz/2:dz:L/2-dz/2; %电压激励项，且只在缝隙处才有电场 V=1; Ez_in((Nz+1)/2)=V/dz; b(1:Nz)=Ez_in(1:Nz); %自适应数值积分，阻抗矩阵h1的计算? for m=1:Nz for n=1:Nz x=z(m); y1=@(y)f(x,y,a,k).*cos((2.*n-1) *pi*y./L); y2=integral(y1,-L/2,L/2);%quad h1(m,n)=xishu*y2; end end %h1的LU分解 for p=1:Nz h(p)=det(h1(1:p,1:p)); end h2=h(1:Nz); for i=1:Nz if h(1,i)==0 disp h2; return end end if h2(1,i)~=0 for j=1:Nz U(1,j)=h1(1,j); end for k=2:Nz for i=2:Nz for j=2:Nz L1(1,1)=1; L1(i,i)=1; if i>j L1(1,1)=1; L1(2,1)=h1(2,1)/U(1,1); L1(i,1)=h1(i,1)/U(1,1); L1(i,k)=(h1(i,k)-L1(i,1:k-1)*U(1:k-1,k))/U(k,k); else U(k,j)=h1(k,j)-L1(k,1:k-1)*U(1:k-1,j); end end end end h2;U;L1; end %MOM方法系数的计算 U1=pinv(U); L2=pinv(L1); I=U1*L2*b; %求解半波振子的电流分布 y3=0; M=101; s=L/M; zd=-L/2+s/2:s:L/2-s/2; for m=1:M for i=1:Nz y4=I(i)*cos((2*i-1)*pi*zd(m)/L); y3=y3+y4; end II(m)=y3; y3=0; end current=abs(II); current_max=max(current); current=current/current_max; figure(1) % plot(zd,current,'r.-') %半波振子电路分布解析解 Dianliu=zeros(M,1); beta=2*pi/lambda; for m=1:M dianliu(m)=current_max*sin(beta*(lambda/4-abs(zd(m)))); end hold on plot(zd*1000,dianliu,'k-') title('偶极子天线Jz分布') % % legend('计算解','解析解'); xlabel('z(m)') ylabel('归一化电流') set(gca,'XLim',[0 180],'XTick',0:15:180); grid on