Matlab范例_飞行航程计算

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