clcclear??all
tic
sigma_r=0.1;?%距離的均方差
sigma_alpha=0.3*pi/180;?%測角的均方差
%?sigma_beta=1*pi/180;
%?sigma_x=0.007;?sigma_y=sigma_x;%站址誤差
%?sigma_z=0.001;
%站址坐標?
x1=50;?y1=-86;z1=0;%千米
x2=100;?y2=0;z2=0;%千米
x3=50;?y3=86;z3=0;%千米
%?x4=-50;?y4=86;z4=0;
%?x5=-100;?y5=0;z5=0;
%?x6=-50;?y6=-86;z6=0;
%目標范圍
x=-151:2:151;%千米
y=-151:2:151;%千米
z=1:2:21;????%千米
for?i=1:length（x）
????for?j=1:length（y）
????????for?l=1:length（z）
d1=（x（i）-x1）^2+（y（j）-y1）^2;?R1=（x（i）-x1）^2+（y（j）-y1）^2+（z（l）-z1）^2;
d2=（x（i）-x2）^2+（y（j）-y2）^2;?R2=（x（i）-x2）^2+（y（j）-y2）^2+（z（l）-z2）^2;
d3=（x（i）-x3）^2+（y（j）-y3）^2;?R3=（x（i）-x3）^2+（y（j）-y3）^2+（z（l）-z3）^2;
%?d4=（x（i）-x4）^2+（y（j）-y4）^2;?R4=（x（i）-x4）^2+（y（j）-y4）^2+（z（l）-z4）^2;
%?d5=（x（i）-x5）^2+（y（j）-y5）^2;?R5=（x（i）-x5）^2+（y（j）-y5）^2+（z（l）-z5）^2;
%?d6=（x（i）-x6）^

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????5213??2014-04-14?22:20??xiaoneng_GDOP.m