N3=200;
phifxt=0;
dth=pi/（2*N3）;?????????????????????
dph=2*pi/（2*N3）;???????????????????????
Eth=zeros（2*N3+1）;
Eph=zeros（2*N3+1）;
f=10e9;??????????????????????????????
lamda=3*10^8/f;????
k=2*pi/lamda;???????????????????????????
M=21;??????????????????????????????????
N=15;??????????????????????????????????
r0=（（0.7*M-0.35）^2+（0.7*N）^2）^0.5*lamda;
r=0.535;
M1=（M-1）/2;
N1=（N-1）/2;
dx=0.7*lamda;
dy=0.7*lamda;????????%先設(shè)置所需參數(shù)???????????????????
theta0=0;??????????????????????????????????????????????????????????????????????????????????????
m0=M;???????????????????????????????????
m1=（m0-1）/2;
%求解近場數(shù)據(jù)
for?p=0:2*N3
????theta=p*dth+0.0001;
???for?q=0:2*N3
???????phi=q*dph;
???????[xyz]=r2x（rthetaphi）;
Ex1=zeros（MN）;
Ey1=zeros（MN）;
Ez1=zeros（MN）;
for?N2=1:21