测试 表面波、单模、多模.m

关于泄漏同轴电缆性质

clear;
q=pi/4;              % 开槽倾角  %
b=0.025 ;                 %电缆外径 %
c=3*10.^8;            %   光速%
f=600*10.^6;           %  工作频率  %       
k=2*pi*f/c;                  % 波数    %
L=0.08;                  % 开槽口长度    %
S=0.042;                  % 1/4个周期  %
sym t;
sita=pi/2;            %     %
H=-j/pi*int('exp(-j*cosh(cos(sita)))','t',-inf,inf)
s=1
for fai=0:0.01*pi:2*pi
   x(s)=fai;
   d=0;
   for  n=-6000:6000
       d=d+sin(q).^2.*(    
                         (exp(j*n*fai)*j.^n)./H.*
                                         (
                              
                                             exp(j*k*S*cos(sita)).* 
                                                      (
                                                        (cos(L/b*(n*sin(q)-k*b*cos(fai)*cos(q)))-cos(k*L))./ ((k*b).^2-(n*sin(q)-k*b*cos(sita)*cos(q)).^2)).*(cos(q)+n*cos(sita)*sin(q)./(k*b*sin(sita).^2))
                                                       )
                                               +exp(-j*k*S*cos(sita)).*
                                                      (
                                                         (cos(L/b*(n*sin(q)+k*b*cos(sita)*cos(q)))-cos(k*L))./((k*b).^2-(n*sin(q)+k*b*cos(sita)*cos(q)).^2)).*(cos(q)-n*cos(sita)*sin(q)./(k*b*sin(sita).^2))
                                         
                                                       )
                              
                                           )
                            )
                         );
               
               
           end
           f(s)=abs(d);
           s=s+1;
       end;
       polar(x,f);