ver84tover90.asv

本程序设计了常用的几个坐标系的转换

% WGS-84直角坐标转换成PE-90直角坐标
clear;
input('input vertical_84:');                                                  % 输入WGS-84坐标系下的直角坐标值
X_84=input('X_84=');                                                          % X轴分量（米）
Y_84=input('Y_84=');                                                          % Y轴分量（米）
Z_84=input('Z_84=');                                                          % Z轴分量（米）
X_90=1/(1+1.9e-6*1.9e-6)*[X_84+1.9e-6*(Y_84-2.5)];
Y_90=1/(1+1.9e-6*1.9e-6)*[-1.9e-6*X_84+(Y_84-2.5)];
Z_90=1/(1+1.9e-6*1.9e-6)*[(1+1.9e-6*1.9e-6)*Z_84];
X_90=fprintf('X_90=%.3f \n',X_90);                                           % 以指数形式输出实数，且输出10位小数，输出换行
Y_90=fprintf('Y_90=%.3f \n',Y_90);
Z_90=fprintf('Z_90=%.3f \n',Z_90);
%DeltaX=fprintf('DeltaX=%.3f \n',X_84-X_90);
%DeltaY=fprintf('DeltaY=%.3f \n',Y_84-Y_90);
%DeltaZ=fprintf('DeltaZ=%.3f \n',Z_84-Z_90);