地球的平均半径为6371千米，已知地球上两个城市A、B的经度和纬度，编程序求出这两个城市之间的地面距离

地球的平均半径为6371千米，已知地球上两个城市A、B的经度和纬度，编程序求出这两个城市之间的地面距离。~

#include
#include

using namespace std;

#ifndef M_PI
#define M_PI (3.141592653)
#endif

#define EARTH_R (6371.0)
#define MIN_DISTANCE (0.00001)//最小分辨度为1厘米

//求double类型的绝对值
double dabs(double value){
return value > 0.0 ? value : (0.0 - value);
}

//求double类型的较大数
double max(double a, double b){
return a > b ? a : b;
}
//求double类型的较小数
double min(double a, double b){
return a < b ? a : b;
}

double toRadian(double angle){
return angle / 360.0 * 2 * M_PI;
}

//已知半径、圆心角角度（弧度制），求圆心角对应的弦长
double getChordByRadiusAngle(double radius, double angle){
return 2.0 * radius * sin((angle / 2.0));//求得弦长
}

//已知半径、弦长，求弧对应的弧度（弧度制）
double getAngleByRadiusChord(double radius, double chord){
return 2.0 * asin(chord / 2.0 / radius);//求得弦对应的大圆圆心角，弧度制
}

//已知半径、弦长，求弦端点在球面上的距离
double getArcByRadiusChord(double radius, double chord){
double angle =getAngleByRadiusChord(radius, chord);//求得弦对应的大圆圆心角，弧度制

double distance = (M_PI * radius) * (angle / M_PI);
return distance;
}

//经度相同，求纬度不同的两点A、B间的弦长
//输入：A点的纬度（[90,-90]，北纬为正）、B点的纬度
//输出：点A、B间的弦长
double chord_Longitude(double latitudeA, double latitudeB){
double latitudeDiff = dabs(latitudeA - latitudeB);
double angle = toRadian(latitudeDiff);//求得大圆圆心角，弧度制

return getChordByRadiusAngle(EARTH_R, angle);//求得弦长
}

//纬度相同，求经度不同的两点A、B间的弦长
//输入：A点的经度（[180,-180]，东经为正）、B点的经度、纬度（[90,-90]，北纬为正）
//输出：点A、B间的弦长
double chord_Latitude(double longitudeA, double longitudeB, double latitude){
latitude = dabs(latitude);
double laR = EARTH_R * cos((latitude / 180.0) * M_PI);//求得纬线圈半径

double longitudeDiff = dabs(longitudeA - longitudeB);
longitudeDiff = longitudeDiff < 180.0 ? longitudeDiff : (360.0 - longitudeDiff);
double angle = toRadian(longitudeDiff);

return getChordByRadiusAngle(laR, angle);//求得弦长
}



//输入两点的经纬度，求出两点的距离
//输入：A点的经度（[180,-180]，东经为正）、A点的纬度（[90,-90]，北纬为正）、B点的经度、B点的纬度
//输出：点A、B的距离
//先求得两点的弦长，再求弦对应的大圆圆心角，最后得到大圆弧线长即为距离
double distanceOf2Points(double longitudeA, double latitudeA, double longitudeB, double latitudeB){
double chord = 0.0;
/*
if(longitudeA == longitudeB){//经度相同
chord = chord_Longitude(latitudeA, latitudeB);
}else if(latitudeA == latitudeB){//纬度相同
chord = chord_Latitude(longitudeA, longitudeB, latitudeA);
}else{*/
//求以AB为对角线，以过点A、B的两条纬线、经线的四个交点为顶点的等腰梯形的四条边长
double baselineA = chord_Latitude(longitudeA, longitudeB, latitudeA);//过点A的底线长
double baselineB = chord_Latitude(longitudeA, longitudeB, latitudeB);//过点B的底线长
double waistline = chord_Longitude(latitudeA, latitudeB);//腰线长
/*
if(baselineA == baselineB){//实为矩形
chord = sqrt(baselineA * baselineA + waistline * waistline);
std::cout<<"矩形: "<<chord<<endl;
}else{*/
double baselineS = min(baselineA, baselineB);//较短的底线
double baselineL = max(baselineA, baselineB);//较长的底线
double shortline = (baselineL - baselineS) / 2.0;//分割中间矩形后，边侧直角三角形的边

chord = sqrt(waistline * waistline - shortline * shortline + (baselineL - shortline) * (baselineL - shortline));
//}
//}

double distance = getArcByRadiusChord(EARTH_R, chord);
if(distance >= MIN_DISTANCE){
return distance;
}
return 0.0;
}

int main(){
double params[4];
int i;
double distance;

std::cout<<"输入格式：A点经度 A点纬度 B点经度 B点纬度"<<std::endl<<std::endl;

while(true){
i = 0;
distance = -1.0;
while(std::cin>>params[i]){//读入四个经纬值
i++;
if(4 == i){
break;
}
}
if(4 == i){
distance = distanceOf2Points(params[0], params[1], params[2], params[3]);
std::cout<<"A点("<<params[0]<<","<<params[1]<<")、B点("<<params[2]<<","<<params[3]<<") 的距离为："<<std::endl<<distance<<"(千米)"<<std::endl<<std::endl;
}else{
break;//输入非数字 字符时，结束程序
}
}
return 0;
}

我正在给你写，请不要采纳别人的答案

求出经纬度差，根据cosθ=cosθ1cosθ2求出两城到球心的角的弧度θ，再根据L = θ * R可求球面距离。
编程如下，仅是给个思路，考虑不周之处楼主请自行完善。

#include <stdio.h>
#include <math.h>
#include <string>
using namespace std;

const double R = 6371.0;
const double PI = 3.141592653;

struct Coordinate
{
string X;//经度
string Y;//纬度
};

double GetTheta(Coordinate& A,Coordinate& B)
{
int la,lb;
double a,b,dx,dy;
char cha,chb;
string ax,ay,bx,by;

ax = A.X;
ay = A.Y;
bx = B.X;
by = B.Y;

la = ax.length();
lb = bx.length();
cha = ax[la-1]; ax[la-1] = 0;
chb = bx[lb-1]; bx[lb-1] = 0;
sscanf(ax.data(),"%lf",&a);
sscanf(bx.data(),"%lf",&b);
if(cha == chb)
dx = abs(a-b);
else
{
dx = abs(a+b);
if(dx > 180) dx = 360 - dx;
}

la = ay.length();
lb = by.length();
cha = ay[la-1]; ay[la-1] = 0;
chb = by[lb-1]; by[lb-1] = 0;
sscanf(ay.data(),"%lf",&a);
sscanf(by.data(),"%lf",&b);
if(cha == chb)
dy = abs(a-b);
else
dy = abs(a+b);

dx = dx * PI / 180;//角度转弧度
dy = dy * PI / 180;
return acos(cos(dx)*cos(dy)) ;

}

int main()
{
double theta,dis;

Coordinate A,B;
A.X = "12.35E";
A.Y = "23.05N";

B.X = "30.25E";
B.Y = "20.28S";

theta = GetTheta(A,B);

dis = theta * R;

printf("%.2lfkm\n",dis);

return 0;
}

#include<stdio.h>
#include<math.h>
#define PI 3.141592653
double Yxdl_C(double a,double b,double c) /*余弦定理*/
{
return acos((a*a+b*b-c*c)/(2*a*b));
}
void main()
{
double R=6371.0,S1,S2,S3,S4,S5,DMJL;
double A_jd,A_wd,B_jd,B_wd,wdc,jdc,zjdc;
printf("规定东经为正，北纬为正；西经为负，南纬为负，请对应输入符号。\n");
lp:
printf("请输入A点经度：");
scanf("%lf",&A_jd);if(fabs(A_jd)>180.0){printf("输入有误，请重新输入!\n");goto lp;}
printf("请输入A点纬度：");
scanf("%lf",&A_wd);if(fabs(A_wd)>90.0){printf("输入有误，请重新输入!\n");goto lp;}
printf("请输入B点经度：");
scanf("%lf",&B_jd);if(fabs(B_jd)>180.0){printf("输入有误，请重新输入!\n");goto lp;}
printf("请输入B点纬度：");
scanf("%lf",&B_wd);if(fabs(B_wd)>90.0){printf("输入有误，请重新输入!\n");goto lp;}
printf("A点位置：");if(A_jd>0)printf("东");else printf("西");printf("经%.3f°，",fabs(A_jd));
if(A_wd>0)printf("北");else printf("南");printf("纬%.3f°\n",fabs(A_wd));
printf("B点位置：");if(B_jd>0)printf("东");else printf("西");printf("经%.3f°，",fabs(B_jd));
if(B_wd>0)printf("北");else printf("南");printf("纬%.3f°\n",fabs(B_wd));
if(fabs(A_jd-B_jd)>180)jdc=fabs(A_jd-B_jd)-180;else jdc=fabs(A_jd-B_jd); /*计算经度差*/
wdc=fabs(A_wd-B_wd); /*计算纬度差*/
S5=2*R*sin(PI*jdc/2/180)*cos(B_wd*PI/180);
zjdc=Yxdl_C(R,R,S5);
S1=R*sin(PI*wdc/180);
S2=R*sin(zjdc);
S3=fabs(R*(cos(PI*wdc/180)-cos(zjdc)));
S4=sqrt(S1*S1+S2*S2+S3*S3); /*两点间直线距离*/
DMJL=R*Yxdl_C(R,R,S4);
printf("A、B两点地面距离为：%lf千米。\n",DMJL);
goto lp;
}
/*终于弄好了，不懂的可以问我*/

这不是数学问题么，吧两城市的地面距离看成扇面。

---

白河县15090399313： 地球的平均半径为6371 - ？
骆备显克：[答案] 地球平均半径 6371.004千米 地球赤道半径 6378.140千米 地球极地半径 6356.755千米

白河县15090399313： 地球的平均半径约为6371千米,可用科学记数法表示为 米,此时有效数字为 个,精确到 位. - ？
骆备显克： 地球的平均半径约为6371千米,可用科学记数法表示为 6.371*10的6次方米,此时有效数字为 4个,精确到 千位.

白河县15090399313： 地球的平均半径是多少千米 - ？
骆备显克： 6371千米. 地球(Earth)是太阳系中由内及外的第三颗行星,是太阳系中直径、质量和密度最大的类地行星,地球自西向东自转,同时围绕太阳公转,距太阳约1.496亿千米.地球表面积5.1亿平方公里,其中71%为海洋,29%为陆地,在太空...

白河县15090399313： 已知地球的平均半径为6371千米,假设在地球的同一纬度上,有两个处于不同经度的城市 - ？
骆备显克：[答案] #include #define R 6371 main() { int i=0; int j=0; float x1=0; float x2=0; float distance=0; float angle=0; float Perimeter=0; //周长 printf(＂备注:经度度为东经,请输为正数(东经20度输入为20),西经请输为负数(西经20度输入为-20)\n＂); loop:...

白河县15090399313： 地球上同一纬度两地距离计算公式已知地球的平均半径为6371千米,假设在地球的同一纬度上,有两个处于不同经度的城市A、B,请给出两地距离计算公式- - - ？
骆备显克：[答案] 地球的半径是R,纬度的度数是a,两个城市的经度分别是b1、 b2. 城市所在纬度平面的半径(与赤道面平行)r = Rxcosa; 两城市的距离是纬度平面的弧长,D = (b1-b2)x 3.1416/180 x r.

白河县15090399313： 地球的平均半径为6371千米,已知地球上两个城市A、B的经度和纬度,编程序求出这两个城市之间的地面距离 - ？
骆备显克： 求出经纬度差,根据cosθ=cosθ1cosθ2求出两城到球心的角的弧度θ,再根据L = θ * R可求球面距离.编程如下,仅是给个思路,考虑不周之处楼主请自行完善.#include <stdio.h>#include <math.h>#include <string> using namespace std; const ...

白河县15090399313： 已知地球半径为6371千米求1)北纬40°和北纬80°分别所在的两个平面之间的距离2)北纬40°和南纬60°分别所在的两个平面之间的距离 - ？
骆备显克：[答案] 1)答:2179千米 6371*(sin80-sin40)=2179 2)答:9612.65千米 6371*(sin40+sin60)=9612.65

白河县15090399313： 已知地球的平均半径是6371km,求地球表面积( )km² - ？
骆备显克： S=4πR²=509805890.96km²

白河县15090399313： 已知地球半径约为6371千米.上海的位置约为东经121°、北纬31°,大连的位置约为东经121°、北纬39°,里斯本的位置约为西经10°、北纬39°.(1)若飞机以平... - ？
骆备显克：[答案] (1)∵上海与大连在同一经线上,∴它们在地球的同一个大圆上. 设地球的球心为O,上海、大连分别为点A、B. 由上海、大连的经、纬度知∠AOB=8°地球半径r≈6371千米 经计算得AB的弧长:6371*π*8°180°≈889.56(千米) 889.56÷720≈1.2(小时...

白河县15090399313： 地球平均半径是6371千米,赤道周长是多少千米? - ？
骆备显克：[答案] 赤道周长 =2*赤道半径*π =2*6371*3.14 =40 075.24千米 约4万千米