ftpe.m

来自「计算电波传播的一个MATLAB程序,用抛物线法求解的..还是初始阶段,共享一下.」· M 代码 · 共 34 行

%The program is used to illustrate the loss of wavefront propagating forward 
%hs-height of the source antenna in meter
%f-frequency of source in MHz
function FTPE(hs,f)

R=6400000;
k0=pi*f/150;
deltax = 0.18;
deltah = 0.5;
range = eps:deltax:200;
height = eps:deltah:1000;

for k=1:length(height)
    for m=1:length(range)
        costheta(k,m) = range(m)/sqrt((hs-height(k))^2 + range(m)^2);
    end
end

psip = 2*real( exp( -((height').^2.*log(2)) ./(2*k0^2*(sin(pi/12))^2) ) .*exp(-i.*height'.*hs) );
psi = ifft(psip);
result = psi;
clear psip;

for m=2:length(range)
    psi = exp(i*k0*deltax*(height')./R) .*ifft( exp(i*deltax*k0*(costheta(:,m) - 1)) .*fft(psi) );
    result = [result psi];
end

clear psi;
clear costheta;
clear k0 k m R hs f;
keyboard;
figure;
meshc(range,height,-20*log10(abs(result)))