clear;
a=9e-3;                        %內(nèi)導(dǎo)體半徑
b=22.8e-3;                       %外導(dǎo)體半徑
epsir=1.25;                 %介質(zhì)的介電常數(shù) 
sita=pi/4;             %開(kāi)槽的傾角
l0=20e-3;                      %開(kāi)槽長(zhǎng)度的一半
w=5e-3;                       %開(kāi)槽寬度
V0=0.01;                         %激勵(lì)電壓
f=24e8;                       %工作頻率900MHz
womiga=2*pi*f;                %角頻率
c=3e8;                        %自由空間電磁波傳播速度
lamda=c/f;                    %波長(zhǎng)
k0=2*pi*f/c;                  %自由空間波數(shù)   
mu0=4e-7*pi;
beta=k0*sqrt(epsir);           %電纜內(nèi)電波傳播常數(shù)       
beta1=beta*cos(sita);         %開(kāi)縫方向傳播常數(shù)分量
P=0.2;                          %開(kāi)縫周期必須大于半波長(zhǎng)



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));       %外導(dǎo)體內(nèi)側(cè)的周向磁場(chǎng)

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

 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;                               %計(jì)算軸向磁流的傅里葉系數(shù)

   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;