xyz坐标转换ybc_GNSS仰角和方位角的计算及代码，XYZ转BLH坐标的代码及原理

function [E,A]= Get_EA(sx,sy,sz,x,y,z)

%GET_EA Summary of this function goes here

%sx,sy,sz:站点的XYZ坐标，x,y,z:卫星的XYZ坐标

% Detailed explanation goes here

[sb,sl]=XYZtoBLH(sx,sy,sz);

T=[-sin(sb)*cos(sl) -sin(sb)*sin(sl) cos(sb);

-sin(sl) cos(sl) 0;

cos(sb)*cos(sl) cos(sb)*sin(sl) sin(sb)];%transition matrix(XYZ to NEU)

deta_xyz=[x,y,z]-[sx,sy,sz];

NEU=T*(deta_xyz)';

E=atan(NEU(3)/sqrt(NEU(1)*NEU(1)+NEU(2)*NEU(2)));

A=atan(abs(NEU(2)/NEU(1)));

if NEU(1)>0

if NEU(2)>0

else

A=2*pi-A;

end

else

if NEU(2)>0

A=pi-A;

else

A=pi+A;

end

end

end

计算仰角\(E\)和方位角\(A\)的公式：

\[E=arctan\left(\frac{cos(\phi_2-\phi_1)\times cos\beta-0.15}{\sqrt{1-\left[cos(\phi_2-\phi_1)\