[jisuan.m]
funciton s=jisuan(x)
J=1.25;              % 相對介電常數    %
c=3*10.^8;            % 光速    %
P=0.233;             % 開槽周期   %
L=1000;              % 電纜總長   %
f=450*10.^6;          % 工作頻率    %
k0=2*pi*f/c;         % 波數    %
b=k0*J.^0.5;          % 傳播常數    %
N=fix(1000/0.233);   % 開槽口總數  %
s=0;
for i=1:N
    s(i)=2/sqrt(2.^2+(x-i*P).^2)*exp(-j*k0*sqrt(2.^2+(x-i*P).^2)-j*b*(i-1)*P)/sqrt(2.^2+(x-i*P).^2);
    s=s+s(i)
end;