资源描述：
《数模论文范文.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、海南大学第四届数学建模竞赛承诺书我们仔细阅读了海南大学数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从
2、A/B/中选择一项填写）：C.1我们的参赛报名号为（按报名顺序填写的号码）：1参赛队员：1.队长姓名：王磊所在班级：06食质联系方式：136989419022.队员姓名：王许所在班级：06经管联系方式：136989423653.队员姓名：所在班级：联系方式：日期：2009年4月14日一、问题重述在约10,000米高空的某边长160公里的正方形区域内，经常有若干架飞机作水平飞行。区域内每架飞机的位置和速度向量均由计算机记录其数据，以便进行飞行管理。当一架欲进入该区域的飞机到达区域边缘时，记录其数据后，要立即
3、计算并判断是否会与区域内的其它飞机发生相撞。如果发生相撞，则应计算如何调整各架（包括新进入的）飞机的飞行方向角，以避免碰撞。现假设条件如下：(1)不相撞的标准为任意两架飞机的距离大于8公里；(2)飞机飞行方向角调整的幅度不应超过30度；(3)所有飞机的飞行速度均为每小时800公里；(4)进入该区域的飞机在到达区域边缘时，与区域内飞机的距离应在60公里以上；(5)最多需考虑6架飞机；(6)不必考虑飞机离开此区域后的情况。请你对这个避免碰撞的飞行管理问题建立数学模型，列出计算步骤，对以下数据进行计算（方向角误
4、差不超过0.01度），要求飞机飞行方向角调整的幅度尽量小。设该区域4个顶点的坐标为（0，0），（160，0），（160，160），（0，160）。记录数据为：飞机编号横坐标纵坐标方向角（度）1150140243285852363150155220.54145501595130150230新进入0052注：方向角指飞行方向与x轴正向的夹角。二、问题分析此问题很容易想到以飞机调整的飞行角度平方和作为目标函数，而以每两架飞机之间的最小距离不超过8km，各飞机飞行角度调整的值不超过30°为约束条件。如此得出的是一
5、个非线性模型，在计算上可能会复杂些，但一目了然。三、符号说明表示表示时间；分别表示第架飞机的横纵坐标（问题中已给出）；表示第架飞机的飞行方向角（问题中已给出）；表示时刻第架飞机与第架飞机间的距离；表示飞机的飞行高度（）。四、模型的建立由题意可知，目标函数是约束条件为和其中利用微积分的知识可求出，由这里将代入即可求出。于是本问题的一个数学模型为引入记号：，，（g是由按构成的向量，在下面的程序中计算），则模型变为其中，。五、模型的求解调用Matlab命令fmincon求解，先写两个M函数airfun.m和ai
6、rfunco.m如下：M函数airfun.mfunctionf=airfun(delta)f=delta*delta';M函数airfunco.mfunction[c,ceq]=airfunco(delta)x0=[150851501451300];y0=[140,85,155,50,150,0];alpha0=[243236220.515923052]*pi/180;v=800;co=cos(alpha0+delta);si=sin(alpha0+delta);fori=2:6forj=1:i-1t(i
7、,j)=(x0(i)-x0(j))*(co(i)-co(j));t(i,j)=t(i,j)+(y0(i)-y0(j))*(si(i)-si(j));t(i,j)=-t(i,j)/v;t(i,j)=t(i,j)/((co(i)-co(j))^2+(si(i)-si(j))^2);ift(i,j)<0,d(i,j)=1000;elsed(i,j)=(x0(i)-x0(j)+v*t(i,j)*(co(i)-co(j)))^2;d(i,j)=d(i,j)+(y0(i)-y0(j)+v*t(i,j)*(si(i)-
8、si(j)))^2;endendendc=64-[d(2,1),d(3,1:2),d(4,1:3),d(5,1:4),d(6,1:5)];ceq=[];在Matlab命令窗口计算如下：>>deltaini=zeros(1,6);>>vlb=-pi*ones(1,6)/6;vub=pi*ones(1,6)/6;>>options=optimset('LargeScale','off');>>[dt,fval]=fmincon(@
