Làm điều đó bằng ngôn ngữ PHP:
// giả sử độ cao = 0
$ earthR = 6371; // trong km (= 3959 trong dặm)
$ LatA = 37.418436;
$ LonA = -121.963477;
$ DistA = 0,265710701754;
$ LatB = 37.417243;
$ LonB = -121,961889;
$ DistB = 0,234592423446;
$ LatC = 37.418692;
$ LonC = -121,960194;
$ DistC = 0,0548954278262;
/ *
#use hình cầu tự động
#if sử dụng ellipsoid, bước này hơi khác
# Chuyển đổi trắc địa Lat / Dài sang ECEF xyz
# 1. Chuyển đổi Lat / Long sang radian
# 2. Chuyển đổi Lat / Long (radian) sang ECEF
* /
$ xA = $ earthR * (cos (deg2rad ($ LatA)) * cos (deg2rad ($ LonA)));
$ yA = $ earthR * (cos (deg2rad ($ LatA)) * sin (deg2rad ($ LonA)));
$ zA = $ earthR * (sin (deg2rad ($ LatA)));
$ xB = $ earthR * (cos (deg2rad ($ LatB)) * cos (deg2rad ($ LonB)));
$ yB = $ earthR * (cos (deg2rad ($ LatB)) * sin (deg2rad ($ LonB)));
$ zB = $ earthR * (sin (deg2rad ($ LatB)));
$ xC = $ earthR * (cos (deg2rad ($ LatC)) * cos (deg2rad ($ LonC)));
$ yC = $ earthR * (cos (deg2rad ($ LatC)) * sin (deg2rad ($ LonC)));
$ zC = $ earthR * (sin (deg2rad ($ LatC)));
/ *
CÀI ĐẶT, DỰNG LÊN:
sudo lê cài đặt Math_Vector-0.7.0
lê cài đặt Math_Matrix-0.8.7
* /
// Bao gồm PEAR :: Math_Matrix
// /usr/share/php/Math/Matrix.php
// include_path = ".: / usr / local / php / pear /"
request_once 'Math / Matrix.php';
allow_once 'Math / Vector.php';
allow_once 'Math / Vector3.php';
$ P1vector = new Math_Vector3 (mảng ($ xA, $ yA, $ zA));
$ P2vector = new Math_Vector3 (mảng ($ xB, $ yB, $ zB));
$ P3vector = new Math_Vector3 (mảng ($ xC, $ yC, $ zC));
#from wikipedia: http://en.wikipedia.org/wiki/Trilateration
#transform để có được vòng tròn 1 tại điểm gốc
#transform để có được vòng tròn 2 trên trục x
// CALC EX
$ P2minusP1 = Math_VectorOp :: sub trích ($ P2vector, $ P1vector);
$ l = Math_Vector mới ($ P2minusP1);
$ P2minusP1_length = $ l-> length ();
$ Norm = new Math_Vector3 (mảng ($ P2minusP1_length, $ P2minusP1_length, $ P2minusP1_length));
$ d = $ định mức; // lưu calc D
$ ex = Math_VectorOp :: split ($ P2minusP1, $ Norm);
// echo "ex:". $ ex-> toString (). "\ n";
$ ex_x = floatval ($ ex -> _ tuple-> getData () [0]);
$ ex_y = floatval ($ ex -> _ tuple-> getData () [1]);
$ ex_z = floatval ($ ex -> _ tuple-> getData () [2]);
$ ex = new Math_Vector3 (mảng ($ ex_x, $ ex_y, $ ex_z));
// CALC tôi
$ P3minusP1 = Math_VectorOp :: sub trích ($ P3vector, $ P1vector);
$ P3minusP1_x = floatval ($ P3minusP1 -> _ tuple-> getData () [0]);
$ P3minusP1_y = floatval ($ P3minusP1 -> _ tuple-> getData () [1]);
$ P3minusP1_z = floatval ($ P3minusP1 -> _ tuple-> getData () [2]);
$ P3minusP1 = new Math_Vector3 (mảng ($ P3minusP1_x, $ P3minusP1_y, $ P3minusP1_z));
$ i = Math_VectorOp :: dotSản phẩm ($ ex, $ P3minusP1);
// tiếng vang "i = $ i \ n";
// CALC EY
$ iex = Math_VectorOp :: scale ($ i, $ ex);
// tiếng vang "iex =". $ iex-> toString (). "\ n";
$ P3P1iex = Math_VectorOp :: sub khu vực ($ P3minusP1, $ iex);
// tiếng vang "P3P1iex =". $ P3P1iex-> toString (). "\ n";
$ l = new Math_Vector ($ P3P1iex);
$ P3P1iex_length = $ l-> length ();
$ Norm = new Math_Vector3 (mảng ($ P3P1iex_length, $ P3P1iex_length, $ P3P1iex_length));
// echo "Norm:". $ Norm-> toString (). "\ n";
$ ey = Math_VectorOp :: split ($ P3P1iex, $ Norm);
// echo "ey =". $ ey-> toString (). "\ n";
$ ey_x = floatval ($ ey -> _ tuple-> getData () [0]);
$ ey_y = floatval ($ ey -> _ tuple-> getData () [1]);
$ ey_z = floatval ($ ey -> _ tuple-> getData () [2]);
$ ey = new Math_Vector3 (mảng ($ ey_x, $ ey_y, $ ey_z));
// KÍNH
$ ez = Math_VectorOp :: crossSản phẩm ($ ex, $ ey);
// tiếng vang "ez =". $ ez-> toString (). "\ n";
// CALC D
// làm điều đó trước
$ d = floatval ($ d -> _ tuple-> getData () [0]);
// tiếng vang "d = $ d \ n";
// CALC J
$ j = Math_VectorOp :: dotSản phẩm ($ ey, $ P3minusP1);
// tiếng vang "j = $ j \ n";
# từ wikipedia
#plug và chug sử dụng các giá trị trên
$ x = (pow ($ DistA, 2) - pow ($ DistB, 2) + pow ($ d, 2)) / (2 * $ d);
$ y = ((pow ($ DistA, 2) - pow ($ DistC, 2) + pow ($ i, 2) + pow ($ j, 2)) / (2 * $ j)) - (($ i / $ j) * $ x);
# chỉ một trường hợp hiển thị ở đây
$ z = sqrt (pow ($ DistA, 2) - pow ($ x, 2) - pow ($ y, 2));
// tiếng vang "x = $ x - y = $ y - z = $ z \ n";
#triPt là một mảng với ECEF x, y, z của điểm trilateration
$ xex = Math_VectorOp :: scale ($ x, $ ex);
$ yey = Math_VectorOp :: scale ($ y, $ eye);
$ zez = Math_VectorOp :: scale ($ z, $ ez);
// CALC $ triPt = $ P1vector + $ xex + $ yey + $ zez;
$ triPt = Math_VectorOp :: add ($ P1vector, $ xex);
$ triPt = Math_VectorOp :: add ($ triPt, $ yey);
$ triPt = Math_VectorOp :: add ($ triPt, $ zez);
// echo "triPt =". $ triPt-> toString (). "\ n";
$ triPt_x = floatval ($ triPt -> _ tuple-> getData () [0]);
$ triPt_y = floatval ($ triPt -> _ tuple-> getData () [1]);
$ triPt_z = floatval ($ triPt -> _ tuple-> getData () [2]);
#convert trở lại lat / long từ ECEF
# chuyển sang độ
$ lat = rad2deg (asin ($ triPt_z / $ earthR));
$ lon = rad2deg (atan2 ($ triPt_y, $ triPt_x));
tiếng vang $ lat. ','. $ lon;