myjw.m

来自「此函数通过两点的大地经纬度（Ja,Wa）和（Jb,Wb）计算两点间距离 Dis」· M 代码 · 共 24 行

clear all
clc

R=6371000;
JJ1=109.4888;
WW1=17.4897;

dis=9573.4;
Ang=96.44;

jj1=JJ1*pi/180;
ww1=WW1*pi/180;
ang=Ang*pi/180;

x=dis*sin(ang);
y=dis*cos(ang);

djj=acos(1-((y/R)^2)/2);
dww=acos(1-((x/R)^2)/2);

jj2=jj1+djj;
ww2=ww1+dww;
JJ2=jj2*180/pi
WW2=ww2*180/pi