【C】【已知经纬度求两点间距离】【三角函数】【函数调用】

题目描述

假设地球是一个标准的球体，其半径是R，忽略地形对距离的影响。球面上的点使用经纬度表示，规定东经为正、西经为负、北纬为正、南纬为负。设A点的经度是α1、纬度是β1；B点的经度是α2、纬度是β2，则A和B之间的距离可按下式近似计算：

同时约定：东经为正,西经为负；南纬为90°+地理纬度值、北纬为90°-地理纬度值
D=Rarccos(C)
其中：C=cos(β1)cos(β2)cos(α1-α2)+sin(β1)sin(β2)；
特别注意：arccosC的单位是弧度.
 

假设R=6371km，编写函数计算地球表面任意两点之间的距离。

输入

两个点的经纬度

输出

两点之间的距离

 实现代码如下：

#include<stdio.h>
#include<math.h>
#define pi 3.14159
#define R 6371
double D(double C){
     return R * acos(C);

}
int main(){
    double a1,a2,b1,b2;
    scanf("%lf %lf",&a1, &b1);
    scanf("%lf %lf",&a2,&b2);
    a1=a1*(pi/180);
    a2=a2*(pi/180);
    b1=b1*(pi/180);
    b2=b2*(pi/180);
    double C=cos(b1)*cos(b2)*cos(a1-a2)+sin(b1)*sin(b2);
    double d = D(C);
    printf("%lf",d);
    return 0;
}

注意：使用三角函数或反三角函数 要将度数转化为弧度

反三角函数：acos()、asin()、atan()