clc;
clear?all;
close?all;
%?x1=-30;x2=-30;x3=60;xt=-60;
%?y1=51.96;y2=-51.96;y3=0;yt=0;
%?z1=0;z2=0;z3=0;zt=0.1;

xt=0;yt=0;zt=0;
x1=-12.99;y1=7.5;z1=-0.1;
x2=12.99;y2=7.5;z2=-0.2;
x3=0;y3=15;z3=-0.3;

z=15;
y=-300:1:300;x=-300:1:300;

for?i=1:length（x）
for?j=1:length（y）
m=x（i）;
n=y（j）;
r1=（（m-x1）.^2+（n-y1）.^2+（z-z1）.^2）.^（1/2）;
r2=（（m-x2）.^2+（n-y2）.^2+（z-z2）.^2）.^（1/2）;
r3=（（m-x3）.^2+（n-y3）.^2+（z-z3）.^2）.^（1/2）;
rt=（（m-xt）.^2+（n-yt）.^2+（z-zt）.^2）.^（1/2）;
c11=（m-x1）/r1;c21=（m-x2）/r2;c31=（m-x3）/r3;ct1=（m-xt）/rt;
c12=（n-y1）/r1;c22=（n-y2）/r2;c32=（n-y3）/r3;ct2=（n-yt）/rt;
c13=（z-z1）/r1;c23=（z-z2）/r2;c33=（z-z3）/r3;ct3=（z-zt）/rt;
c=[（ct1+c11）?（ct2+c12）?（ct3+c13）;（ct1+c21）?（ct2+c22）?（ct3+c23）;（ct1+c31）?（ct2+c32）?（ct3+c33）];

b=inv（c‘*c）*（c‘）;
b11=b（1）;
b12=b（4）;
b13

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1955??2012-11-22?11:28??GDOP\GDOP.asv

?????文件???????1967??2012-11-22?11:35??GDOP\GDOP.m

?????目錄??????????0??2013-06-05?14:19??GDOP

-----------?---------??----------?-----??----

?????????????????3922????????????????????3