%視距計算
clear?all
R0=6370e3*4/3;%地球半徑，米
hr=[7];%雷達天線的高度，米
ht=75;%目標的高度，米
htf=ht/0.3;%目標的高度，英尺
ho=12;%障礙物的高度，米


Rsighta=（（R0+hr）.^2-R0^2）.^0.5+（（R0+ht）^2-R0^2）^0.5;%精確的視距計算
Rsightc=4.12*（hr.^0.5+ht^0.5）;%視距估算?
for?mm=1:length（ht）
dR=50:10:10000;%障礙物的距離
for?kk=1:length（dR）
theta=2*asin（dR（kk）/2/R0）;
k=（（R0+ho（mm））.*cos（theta）-（R0+hr））./（（R0+ho（mm））.*sin（theta））;
C=R0+hr;
%?dAC=（（（R0+ho（mm））.*sin（theta））.^2+（（R0+ho（mm））.*cos（theta）-（R0+hr））.^2）.^0.5;
x=100:100:6000000;
ff=x.^2+（k.*x+C）.^2-（R0+ht）^2;
%?plot（ff）grid?on
ind=find（ff>01）;
xd=x（ind）;
%?if?xd>=Rsightc
%?????xd=Rsightc;
%?end
yd=k*xd+C;
Rsight0=（xd^2+（yd