divergence.m

有关雷达的MATLAB仿真资料！！！！！！！！！！！！！！！

function [D] = divergence(r1, r2, ht, hr, psi)
% calculates divergence
%       inputs  %
%   r1  ground range between radar and specular point in KM
%   r2  ground range between specular point and target in Km
%   psi grazing angle in gegrees
%       parameters  %
%   re  4/3 earth radius 4/3 X 6375 Km
%   r = r1 + r2
psi = psi .* pi ./180; % psi in radians
re = (4/3) * 6375e3;
r = r1 + r2;
arg1 = re.* r .* sin(psi) .*cos(psi);
arg2 = ((2.*r1.*r2./cos(psi)) + re.*r.*sin(psi)) .* ...
    (1+hr./re) .* (1+hr./re);
D = sqrt(arg1 ./ arg2);
return