Преобразование геодезической широты и долготы в косинусную матрицу направления
n = dcmecef2ned(lat, lon)
n = dcmecef2ned(lat, lon) вычисляет матрицу направления косинуса, n, для данного набора геодезической широты и долготы, lat, lon. lat является m массив геодезических широт. lon является m массив долгот. Значения широты и долготы могут быть любыми. Однако значения широты + 90 и -90 могут возвращать непредвиденные значения из-за сингулярности на полюсах .n возвращает значение 3-за-3-by-m матрица, содержащая m матрицы косинусов направления. n выполняет координатное преобразование вектора в осях с земным центром (ECEF) в вектор в осях с северо-восточным направлением вниз (NED). Геодезические широты и долготы вводятся в градусах.
Определите матрицу направления косинуса по геодезической широте и долготе:
lat = 45;
lon = -122;
dcm = dcmecef2ned(lat, lon)
dcm =
0.3747 0.5997 0.7071
0.8480 -0.5299 0
0.3747 0.5997 -0.7071Определите матрицу направления косинуса из нескольких геодезических широт и долгот:
lat = [45 37.5];
lon = [-122 -85];
dcm = dcmecef2ned(lat, lon)
dcm(:,:,1) =
0.3747 0.5997 0.7071
0.8480 -0.5299 0
0.3747 0.5997 -0.7071
dcm(:,:,2) =
-0.0531 0.6064 0.7934
0.9962 0.0872 0
-0.0691 0.7903 -0.6088