辐射方向图2.m

来自「关于泄漏同轴电缆性质」· M 代码 · 共 25 行

clear;
clf;
J=1.5;              % 相对介电常数    %
c=3*10.^8;            % 光速    %
P=0.35;             % 开槽周期   %
L=70;              % 电缆总长   %
f=590*10.^6;          % 工作频率    %
k0=2*pi*f/c;         % 波数    %
b=k0*J.^0.5;          % 传播常数    %
N=fix(L/P);   % 开槽口总数  %
s=1;
for x=0:0.001*pi:pi
    X(s)=x;
    d=0;    
for i=1:N
    d=d+exp(j*i.*(k0.*cos(x)-b).*P);
end;
y(s)=abs(d);
s=s+1;
end
polar(X,y);