[求助]矩量法求解对称天线上的电流积分方程

V/)3d  
clear;clc;tic R(2MI}T  
x=sym('x'); V3_qqz}`r  
h=7.022e-03; ^.M_1$-  
digits(8); w_YY~Af  
}'\M}YM  
for n=1:3 E8o9ufj3  
    r1=(((n-1)/(8*h)-x)^2+1)^0.5; Y3xEFqMU  
    r2=(((n-1)/(8*h)-x)^2+1)^0.5; !9^GkFR6n  
  f1(n)=cos(2*pi*h*x)*exp(-j*2*pi*h*r1)/r1; 3ep L'My$  
  f2(n)=sin(4*pi*h*abs(x))*exp(-j*2*pi*h*r1)/r1;   Si|8xq$E;  
end 0qBXL;sE  
for n=1:3 k,y#|bf,Y  
    A(n)=int(f1(n),'x',-1/(4*h),1/(4*h)); eXdH)|l,\  
    B(n)=int(f2(n),'x',-1/(4*h),1/(4*h)); th=45y"C  
    A(n)=vpa(A(n)); %Ip=3($Ku[  
    B(n)=vpa(B(n)); ?hW(5]p|  
end R~bLEo  
mHc2v==X\-  
H1=[A(1),B(1),1;A(2),B(2),1/sqrt(2);A(3),B(3),0]; 7VJf~\%1j  
H2=[0;-j/(60*sqrt(2));-j/60]; ^1ks`1  
H3=inv(H1)*H2;  %对应系数值 )' 2vUt`_7  
H3=double(H3) V.[b${  
N]|U-fN\  
z=[0:0.01:1/4];  %源点坐标变化 
#|lVQ@=  
I1=real(H3(1)*sin(2*pi*(1/4-abs(z))))+real(H3(2)*sin(4*pi*(1/4-abs(z))));  %实部 =OR&,xt  
plot(I1,z) +@PZ3 [s  
hold on [{J1b  
I2=imag(H3(1)*sin(2*pi*(1/4-abs(z))))+imag(H3(2)*sin(4*pi*(1/4-abs(z))));  %虚部 o#IWH;ck.  
plot(I2,z,'r') l`];CALA4  
hold off !p)cP"fa  
grid on jw]IpGTt  
+7K]5p;!~  
legend('实部电流','虚部电流')< .. 'eoI~*}3WQ  
Q?%v b  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册