matlab实验二飞行航程计算.ppt

*/10 实验二: 飞行航程计算实验 问题背景与球面短程线 经纬度转换和实验任务 练习题和思考题 飞行航线动态演示程序 ? ? ? ? 航程计算的问题背景 当地时间9月27日15点30分，降落在纽约肯尼迪机场；“北京 ? 纽约” 新航线正式开通。 在此之前，“北京 ? 纽约” 航线为 北京时间2002年9月27日14点，国航CA981航班从首都国际机场起飞 ·········· 北京 ? 上海 ? 东京 ? 旧金山 ? 纽约 新航线飞越北极,途经俄罗斯、加拿大两国.加拿大航空管理部门官员称: 新航线至少节约了四小时飞行时间。 假设:飞机飞行高度为10km,飞行速度为每小时980公里;地球半径为6400km 新航线(飞越北极) 航程： 10914.77(km) 飞行时间：11.14(h) W76o E116o 经度 N41o N40o 纬度 纽约 北京 航点 球面短程线长度计算公式 L = R×? 其中,? 是OP1与OP2之间夹角(单位:弧度) 球面短程线 半径为 R 的球面上 P1(x1，y1，z1)和P2(x2，y2，z2), x = R cosθ cosφ y = R cosθ sinφ z = R sinθ θ是P点处球面法线和赤道面的夹角(– 90o ~ +90o).向北 取正为北纬,向南取负为南纬. φ是P点与地球自转轴所 在平面与起始子午面的夹角(– 180o ~ +180o).由起始子 午线起算,向东取正为东经,向西取负为西经。 经纬度转换为地心直角坐标公式 city=[40,116;31,122;36,140;37,-123;41,-76]; R=6400+10; theta=city(:,1)*pi/180; fai=city(:,2)*pi/180; x=R*cos(theta).*cos(fai); y=R*cos(theta).*sin(fai); z=R*sin(theta); op1=[x(1),y(1),z(1)]; op5=[x(5),y(5),z(5)]; Dnew=R*acos(op1*op5 /R^2) mlab31 Dnew = 10994.38 北京?纽约直飞航程计算(mlab61) 电子地图数据11008km W76o W123 o E140 o E122 o E116o 经度 N41o N37o N36o N31o N40o 纬度 纽约 旧金山 东京 上海 北京 航点 航程差 直飞航线航程 太平洋航线航程