三种工作模式的比较.m

关于泄漏同轴电缆性质

clear;
x=500;              % 电缆500m处的场强    %            
J=1.25;              % 相对介电常数    %
c=3*10.^8;            % 光速    %
P=0.233;             % 开槽周期   %
L=1000;              % 电缆总长   %
f1=400*10.^6;     % 表面波工作频率    %
f2=800*10.^6;     % 单模工作频率    %
f3=1200*10.^6;     % 多模工作频率    %
k1=2*pi*f1/c;         % 波数 1   %
k2=2*pi*f2/c;         % 波数 2   %
k3=2*pi*f3/c;         % 波数 3  %
b1=k1*J.^0.5;          % 传播常数1    %
b2=k2*J.^0.5;          % 传播常数2    %
b3=k3*J.^0.5;          % 传播常数3    %
N=fix(L/P);   % 开槽口总数  %
n=1;
for r=0:0.1:12;
    R(n)=r;
    d1=0;
    d2=0;
    d3=0;
for i=1:N
d1=d1+r./sqrt(r.^2+(x-i.*P).^2).*exp(-j.*k1.*sqrt(r.^2+(x-i.*P).^2)-j.*b1.*(i-1).*P)./sqrt(r.^2+(x-i.*P).^2);
d2=d2+r./sqrt(r.^2+(x-i.*P).^2).*exp(-j.*k2.*sqrt(r.^2+(x-i.*P).^2)-j.*b2.*(i-1).*P)./sqrt(r.^2+(x-i.*P).^2);
d3=d3+r./sqrt(r.^2+(x-i.*P).^2).*exp(-j.*k3.*sqrt(r.^2+(x-i.*P).^2)-j.*b3.*(i-1).*P)./sqrt(r.^2+(x-i.*P).^2);
end;
y1(n)=abs(d1);
y2(n)=abs(d2);
y3(n)=abs(d3);
n=n+1;
end
plot(R,y1,'b--');
hold on;
plot(R,y2,'r:');
hold on;
plot(R,y3,'k-');
axis([0,12,0,10]);