功能: 计算两个经纬度坐标间的距离

格式: 
CoorDistForLLA(pa, pb, ra, re)
CoorDistForLLA(pa, pb, ra)
CoorDistForLLA(pa, pb)

pa : 一个1 * 3的矩阵变量,里面存储第1个坐标经纬高
pb : 一个1 * 3的矩阵变量,里面存储第2个坐标经纬高
ra : 地球椭圆长半径(米),默认为6378137.0
re : 地球椭球第一偏心率,默认为0.0818191908426215

说明:
1. a,b参数里的经纬度单位为角度,经度参数范围[-180,180],纬度参数范围[-90,90],高的单位为m
2. 这里假设地球为椭球, 设地球赤道半径为6378137米,地球基准椭球体的极扁率为1/298.257223563, 则计算两个坐标点在椭球体上最短弧线距离.
3. 注意,这里求的是弧线距离,并不是两点的直线距离。

例子:
a = [110, 33,  389];
b = [111, 33.5,400];
s = CoorDistForLLA(a, b)//回车得到如下结果
s =
[ 108446.989848337 ]