已知地球上两点经纬度，如何求其相对方位？

已知地球上两点经纬度，如何求其相对方位？
我用近似公式求出其相对方位如下，例如：
如何知广州或深圳在吴川的什么方位？如已知经纬度A(a,b)和B(c,d)
，求其相对方位。设地球半径为R，则BC=2*Pi*R*(d-b)/(2*Pi)=(d-b)*R，CC'=R*cos(b),AC=2*pi*CC'*(c-a)/(2*pi)=(c-a)*cos(b)*R,
tg(X)=BC/AC=(d-b)/((c-a)*cos(b)),求出了夹角X，就可求出相对方位了。
吴川经纬度为(110.47,21.26)，广州经纬度为(113.14,23.08)，深圳经纬度为(114.07,22.33),中山(113.22,22.31)
tgx1=1.82/(2.67*0.9319)=0.7315,x=36.186,所以广州在吴川的寅方。
tgx2=1.07/(3.6*0.9319)=0.3189,x=17.69,所以深圳也在吴川寅方，但很接近卯方。
tgx3=1.05/(2.75*0.9319)=0.4097,x=22.28,所以中山也在吴川寅方，比深圳远离卯方些。
通过公式求出深圳和广州的夹角为41度，即：
tgx4=(-0.75)/(0.93*0.92)=(-0.75)/0.85556=--0.8766188,x=--41.24
所以深圳位于广州辰方而非巳方，但已接近巳方，考虑到计算是近似值的，所以深圳有可能位于广州的辰方和巳方的交界处。
这个问题似乎并非这样就完结的，如果用切线来表示相对方位，则又不同了，在电脑上画出图来就麻烦了//////////设球面上经纬度A(ja,wa),B(jb,wb),平面AOB与A的子午线平面夹角为AAA，则有公式：
tgAAA=(coswb*sin(jb-ja))/(sinwb*coswa-coswb*sinwa*cos(jb-ja))
这是我推导出的精确公式，花了我几天时间才推导出这精确公式的。//////////吴川经纬度为(110.47,21.26)，广州经纬度为(113.14,23.08)，深圳经纬度为(114.07,22.33),中山(113.22,22.31)
tgAAA=(0.92*0.0466)/(0.3920*0.9319-0.92*0.3626*0.9989)=0.0429/(0.3653-0.3332)=0.0429/0.0321=1.3364   
AAA=53     广州在吴川寅方
tgAAA=(coswb*sin(jb-ja))/(sinwb*coswa-coswb*sinwa*cos(jb-ja))
吴川经纬度为(110.47,21.26)，广州经纬度为(113.14,23.08)，深圳经纬度为(114.07,22.33),中山(113.22,22.31)
tgAAA=(0.925*0.0627)/(0.3799*0.9319-0.925*0.3626*0.998)=0.058/(0.354-0.3347)=0.058/0.0193=3.005   
AAA=71.6    深圳也在吴川寅方
tgAAA=(coswb*sin(jb-ja))/(sinwb*coswa-coswb*sinwa*cos(jb-ja))
吴川经纬度为(110.47,21.26)，广州经纬度为(113.14,23.08)，深圳经纬度为(114.07,22.33),中山(113.22,22.31)
tgAAA=(0.9251*0.048)/(0.3796*0.9319-0.9251*0.3626*0.9988)=0.0444/(0.3537-0.3350)=0.0444/0.01867=2.378   
AAA=67    中山也在吴川寅方 
tgAAA=(coswb*sin(jb-ja))/(sinwb*coswa-coswb*sinwa*cos(jb-ja))
吴川经纬度为(110.47,21.26)，广州经纬度为(113.14,23.08)，深圳经纬度为(114.07,22.33),中山(113.22,22.31)
tgAAA=(0.925*0.0162)/(0.3799*0.92-0.925*0.392*0.9999)=0.01499/(0.3495-0.3626)=0.01499/-0.0131=-1.1443   
AAA=-48.9    深圳在广州辰方//////////
高校导师都无法解答的题目：球面已知两点坐标与运动方向，求运动轨迹相交的两个点坐标
已经A、B两点的坐标，运动方向（即是与经线的夹角），求A，B两点运动轨迹交叉点的坐标。设球面上经纬度A(ja,wa),B(jb,wb),平面AOB与A的子午线平面夹角为AAA，则有公式：
tgAAA=(coswb*sin(jb-ja))/(sinwb*coswa-coswb*sinwa*cos(jb-ja))
这是我推导出的精确公式。以此来建立方程求解就很易了，
设交点经纬度为P(x,y),则有方程
tgAAA*sin(x-ja)*coswb+sinwa*cos(x-ja)*coswb=tgBBB*coswa*sin(x-jb)+sinwb*cos(x-jb)*coswa
用和角公式代入，得到一个关于tgx的一元一次方程，求解就可解出tgx
tgy=(tgAAA*sinx+sinwa*cos(x-ja))/coswa,则tgx,tgy都可求出矣