用matlab计算有障碍物时雷达的视距

%视距计算
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