matlab飞行航程问题,飞行航程计算.ppt

*/10 飞行航程计算实验 问题背景与球面短程线 经纬度转换和实验任务 练习题和思考题 飞行航线演示程序 ? ? ? ? 航程计算的问题背景 当地时间9月27日15点30分，降落在纽约肯尼迪机场；“北京 ? 纽约” 新航线正式开通。 在此之前，“北京 ? 纽约” 航线为 北京时间2002年9月27日14点，国航CA981航班从首都国际机场起飞 ·········· 北京 ? 上海 ? 东京 ? 旧金山 ? 纽约 新航线飞越北极,途经俄罗斯、加拿大两国.加拿大航空管理部门官员称: 新航线至少节约了四小时飞行时间。 假设:飞机飞行高度为10km,飞行速度为每小时980公里;地球半径为6400km 新航线(飞越北极) 航程： 10914.77(km) 飞行时间：11.14(h) 航点 北京 纽约 纬度 N40o N41o 经度 E116o W76o 球面短程线长度计算公式 L = R×? 球面短程线计算公式 利用两城市的经纬度可算出球面上两点坐标 P1(x1，y1，z1)和P2(x2，y2，z2), (单位:弧度) 设地球半径 R =6400(km) ? x = R cosθ cosφ y = R cosθ sinφ z = R sinθ 经纬度转换为地心直角坐标公式 航点 北京 上海 东京 旧金山 纽约 纬度 N40o N31o N36o N37o N41o 经度 E116o E122 o E140 o W123 o W76o city=[40,116;41,-76]; %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); op=[x,y,z]; format bank Dmatrix=R*acos(op*op'/R^2) 两地飞行距离:10993.77 北京?纽约直飞航程计算(mlab21) 电子地图数据: 11008km 思 考 问 题 根据“北京时间2002年9月27日14点，航班从北京起飞，纽约时间9月27日15点30分，降落在纽约”这段消息，分析两地经差和时差，计算飞行时间。 2. 如果飞机时速为980公里/小时，则两条飞行航线的飞行时间差异是多大(设旧航线上中转站不停) 3.数学实验所用的计算方法是近似方法,其误差跟哪些参数有关?如何提高精度？ 4.在互联网上利用“谷歌地球”(google earth)，了解电子地图，寻找“北京、纽约”的经纬度等数据。 city=[40,116;31,122;36,140;37,-123;41,-76]; figure(1),sphere(24),colormap([1 1 1]) axis off,hold on p1=city(1,:);p2=city(5,:); skyway(p1,p2) figure(2),sphere(24),colormap([1,1,1]) axis off,hold on for k=1:4 p1=city(k,:);p2=city(k+1,:); skyway(p1,p2) end 5.下面程序(mlab22)模拟北京到纽约两城市之间飞行航线,试对程序中的语句作功能性注释 function skyway(p1,p2) city=[p1;p2]; theta=city(:,1)*pi/180;fai=city(:,2)*pi/180; x=cos(theta).*cos(fai); y=cos(theta).*sin(fai); z=sin(theta); t=linspace(0,1,20); xt=(1-t)*x(1)+t*x(2); yt=(1-t)*y(1)+t*y(2); zt=(1-t)*z(1)+t*z(2); r=sqrt(xt.*xt+yt.*yt+zt.*zt); xt=xt./r;yt=yt./r;zt=zt./r; plot3(x,y,z,'ro',xt,yt,zt,'b','linewidth',2) 《实验报告》用于表达数据和信息,基本要求： 实验结果与分析—— 获取了哪些数据，反映哪方面信息？ 实验结论—— 实验数据的应用价值？现有数据是否足以支持实验结论? 实验涉及的技术、原理(或数学模型)？ 数据和文字无涂改； 表达信息简洁清楚； 设计数据表格,反映实验结果 */10