lcxeight.m

关于泄漏同轴电缆性质

clear;
a=9e-3;                        %内导体半径
b=22.8e-3;                       %外导体半径
epsir=1.25;                 %介质的介电常数 
sita=pi/4;             %开槽的倾角
l0=20e-3;                      %开槽长度的一半
w=5e-3;                       %开槽宽度
V0=0.01;                         %激励电压
f=24e8;                       %工作频率900MHz
womiga=2*pi*f;                %角频率
c=3e8;                        %自由空间电磁波传播速度
lamda=c/f;                    %波长
k0=2*pi*f/c;                  %自由空间波数   
mu0=4e-7*pi;
beta=k0*sqrt(epsir);           %电缆内电波传播常数       
beta1=beta*cos(sita);         %开缝方向传播常数分量
P=0.2;                          %开缝周期必须大于半波长



syms l;
Fi=double(int(exp(-j*k0*sqrt(l^2+(w/2)^2))/(2*pi*sqrt(l^2+(w/2)^2)),l,-20e-3,20e-3))


H0=V0*sqrt(epsir)/(b*log(b/a));       %外导体内侧的周向磁场

A=-0.5j*womiga*mu0*H0*sin(sita)/[(k0^2-beta1^2)*Fi*sin(2*k0*l0)];    %磁流的系数
   

 Nh=floor(k0*(1-sqrt(epsir))*P/(2*pi))
 Nl=ceil(-k0*(1+sqrt(epsir))*P/(2*pi))
 
   for n=Nl:Nh;
     syms z1;
    J(n-Nl+1)=double(1/P*int(2*A*cos(sita)*(sin((beta1-k0)*l0)*exp(j*k0*z1/cos(sita))-sin((beta1+k0)*l0)*exp(-j*k0*z1/cos(sita))...
        +sin(2*k0*l0)*exp(-j*beta1*z1/cos(sita)))*exp(-j*n*2*pi/P*z1),z1,-l0*cos(sita),l0*cos(sita)));
   end;                               %计算轴向磁流的傅里叶系数

   s=1;
   
for z=0:0.1:200;

   Z(s)=z; 
   Efi=0;
 
 for n=Nl:Nh;
  
     kesai(n-Nl+1)=sqrt(k0^2-(beta+2*pi*n/P)^2)
     besselh(1,2,kesai(n-Nl+1)*1)
     Efi=Efi+j/4*J(n-Nl+1)*kesai(n-Nl+1)*besselh(1,2,kesai(n-Nl+1)*2)*exp(-j*(beta+2*pi*n/P)*z);
 end;


Pr(s)=abs(Efi);
 
s=s+1;
end;
Z
Pr
plot(Z,Pr);hold on;