GDOP

利用matlab计算雷达网对目标定位的GDOP值

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