raoshe.m

在matlab环境上 用几何绕射法计算缝隙的辐射场

clc
clear all
tic
% 初始化 定义常量
d1=0.61;d2=0.61;
f=1*10^9;
c=3*10^8;
lamda=c/f;
k=2*pi/lamda;
r=100*lamda;
E0=1;
EziP1=E0*exp(-j*k*d1)/d1;
EziP2=E0*exp(-j*k*d2)/d2;
% 程序主题部分
theta=0:pi/180:pi; % 角度的离散
f=inline('exp(-j*t.^2)','t');
for i=1:length(theta)  %主循环体
    X1(i)=sqrt(2*k*d1*(cos(theta(i)/2+pi/4)).^2); 
    X2(i)=sqrt(2*k*d2*(cos(theta(i)/2-pi/4)).^2); %费涅尔积分中的积分变量。
    F1(i)=quad(f,0,X1(i));
    F2(i)=quad(f,0,X2(i));
    F3=sqrt(pi)*exp(-j*pi/4)/2;
    K1(i)=(1/sqrt(pi))*exp(j*(X1(i).^2+pi/4)).*(F3-F1(i));
    K2(i)=(1/sqrt(pi))*exp(j*(X2(i).^2+pi/4)).*(F3-F2(i));%改型的费涅尔积分
    F11(i)=2*sqrt(j*pi)*X1(i).*K1(i);
    F22(i)=2*sqrt(j*pi)*X2(i).*K2(i);% F函数与改型的费涅尔积分的函数关系
    Dh1(i)=(-exp(-j*pi/4)/sqrt(2*pi*k))*(F11(i)/cos(theta(i)/2+pi/4));
    Dh2(i)=(-exp(-j*pi/4)/sqrt(2*pi*k))*(F22(i)/cos(theta(i)/2-pi/4)); %绕射系数
    r1(i)=r-d1*sin(theta(i));
    r2(i)=r+d2*sin(theta(i));%相位中远区场的计算
    if theta(i)<=pi/2
     EiP(i)=E0*exp(-j*k*r)*sin(theta(i))/r;
    else
        EiP(i)=0.0001;
    end                      %地面以下部分，直射场很小。
    E1dP(i)=0.5*EziP1*Dh1(i)*sqrt(d1)*exp(-j*k*r1(i))/r;
    E2dP(i)=0.5*EziP2*Dh2(i)*sqrt(d2)*exp(-j*k*r2(i))/r;
    EP(i)=EiP(i)+E1dP(i)+E2dP(i);  
end
  theta(181:361)=2*pi-theta(181:-1:1);
  EP(182:361)=EP(180:-1:1); %根据对称关系得出pi到2*pi的场值。
  EPm=max(abs(EP));
for i=1:361
    EPn(i)=20*log10(abs(EP(i))/EPm);%将场值转化为db表示。且归一化。
end
    polar(theta,EPn+40,'r')  
    hold on
  toc