全国大学生数学建模2015年国二a题

1、2015高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍

2、、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员（打印并签名）：1.指导教师或指导教师组负责人（打印并签名）：日期:年月日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）： 太阳影子定位摘要本文研究了太阳影子定位问题，基于天球坐标系相关知识、球面几何理论以及相似度理论，对不同情况下的数据，建

3、立了相应的数学模型并得到了最优的匹配地点与日期。问题1中，利用球面三角形余弦定理给出了太阳高度角公式，并建立了影子长度变化的数学模型，定性的分析了影子长度关于时角、当地纬度以及赤纬角的变化规律：(1).时角的绝对值越大，影子长度越大；(2).在同一经度上(即时角一定)，当地纬度与此时的太阳赤纬之差越大，影子长度越大；(3).在同一纬度不同经度上，当地经度和此时太阳直射点所在的经度之差越大，影子长度越大。用所建的模型，得到了2015年10月22日北京时间9:00-15:00之间天安门广场3米高的直杆的太阳影子长度的变化曲线。问题2中，由太阳高度角、方位角、时角、赤纬角之间的关系，建立了杆的高度、

4、太阳高度角、太阳方位角等变量间的方程组，推导出影子的长度与太阳高度角满足的方程，该方程不依赖于杆的高度、坐标系的选取。利用附件1的数据并结合遍历相关知识，给出了最优匹配地点：最优地点为北纬24度，东经100度，(云南省沧市临翔区)；次优地点为北纬18度，东经110度(海南省三亚市)。问题3中，在问题2模型的基础上，对未知日期引入了日期参数，导出了相应的优化模型。由附件2的数据得出了最优匹配日期和地点：7月10日，北纬41度，东经80度，新疆维吾尔自治区阿克苏地区乌什县。由附件3的数据得出了最优匹配日期和地点：3月02日，北纬37度，东经110度，陕西省延安市延川县。问题4中，给出了杆的高度的近

5、似值，我们修正了问题3中的模型，在优化模型中，引入了变量h。对视频进行逐针读取，再灰度化处理，利用PhotOShop软件得到了影子坐标的和杆的坐标的像素，由它们之间的比例关系，得出了相应的坐标。利用坐标数据和杆的近似高度，通过优化模型，经过Matlab遍历，得到了最优的拍摄地点：最优的为北纬40度，东经113度，山西省大同市南郊区；其它次优的为北纬44度，东经110度，蒙古；北纬39度，东经110度，陕西省榆林市神木县；北纬39度，东经111度，陕西省榆林市府谷县；北纬39度，东经112度，山西省忻州市神池县；北纬40度，东经114度，河北省张家口市阳原县；北纬41度，东经115度，河北省张家

6、口市崇礼县。如果拍摄日期未知，最优的拍摄地点为北纬40，东经115，河北省张家口市蔚县，日期为6月22日。同时，我们也提出了分光理论，并建立了相应的模型和数学推导。由于时间的限制，不能完善该模型。针对问题二、三、四，提出了影子弧线的相似度模型。由于提取数据的误差，可能使得我们的计算结果与实际地点有稍许偏差。关键字:太阳影子定位，天球坐标系，遍历，时角，赤纬角，太阳高度角，太阳方位角问题重述太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。问题一有两部分组成：（1）：确定影子长度变化数学模型，分析影子长度关于各个参数的变化规律。（2）：用已建立的模型画出2

7、015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。问题二有两部分组成：（1）：根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。（2）：将模型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。问题三有两部分组成：（1）：根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。（2）：将所建的模型分别应用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。问题有三部分组成：（1）：根据附件4一根直杆在太阳下的影子变化

8、视频，建立确定视频拍摄地点的数学模型。（2）：应用已建立的模型给出若干个可能的拍摄地点。（3）：如果拍摄日期未知，根据视频确定出拍摄地点与日期。模型假设（1）：地球是球体；（2）：直杆影子长度为直线段长度而非弧长；（3）：到达地球的太阳光线为平行光线；（4）：太阳高度角为太阳光的入射方向与入射点切平面（地平面）的夹角；（5）：假设地面无障碍物，直杆影子都投影在一个平面上；（6）：假设附件4中拍摄角度对所提取的影长无影响；（7）：假设一天当中太阳赤纬角不变；符号说明表3.1符号说明表h直杆长度r直杆影子长度a太阳高度角t当地太阳时时角5赤纬角A光线与直杆的夹角当地维度当地经度注：规定北纬为正，南

9、纬为负模型建立及求解问题一的分析影子的长度主要依赖于两个因素，即直杆的长度和当地的太阳高度角。由推导可得，太阳高度角又依赖于当时的太阳赤纬，时角和当地的纬度。易知当地的纬度和时角；又可知赤纬角由日期确定，注意我们假设一天当中的赤纬角是不变的。问题一模型的建立及求解ana直杆高度、太阳高度角及影子长度三者有如下关系：(4.1)其中，r为影子长度，h为直杆高度，a为太阳高度角，其关系如图4.1所示。 图4.1直杆高度、太阳高度角及影子长度关系示意图在天文三角形中，已知两边：测者余纬，天体极距pnb=?-，它们的夹角，即天体半圆时角，应用边的余弦公式可以求得第三边：天顶距z=,2a,如图4.2所示。

10、E太阳條纬地平线5fl一W图4.2天球坐标系中的相关参数应用球面三角形边的余弦定理I，2】,有(4.2)cosa=cosbcosc+sinbsinccos,又a+a=兀/2,b+6=兀/2,c+Q=兀/2，所以cos（90o-a）=cos（90o）cos（90。Q）（4.3）+sin（90o-5）sin（90。-Q）sin（90。-w）即有TOC o 1-5 h z HYPERLINK l bookmark30 o Current Document sina=sinQsin5+cosQcos5cosw.（4.4）带入（4.1）式，可得（4.4）sinQsin5+cosQcos5coswh1-（

11、sinQsin5+cosQcos5cosw）2其中，时角w的计算公式为3: HYPERLINK l bookmark34 o Current Document w=15（t12）（度）,（4.5）s太阳赤纬角5（即太阳直射点的维度）：e（2兀（284n）（宀） HYPERLINK l bookmark38 o Current Document 5=23.45s_（度丿（4.6）（365丿其中n为日期序号，例如，1月1日为n=1,2月25日为n=56。下面我们将定性分析影子长度r关于各个参数（时角w、纬度Q及赤纬角5）的变化规律。首先，分析在一天之中时角对影子长度的影响。由假设可知赤纬角5在一天

12、之中是恒定的，又r=htana可知，太阳高度角a越大，影子长度越小；又由式（4.4）可得，时角w的绝对值w越大，a越小。所以可知，w越大，影子长度r越大。其次，分析纬度Q对影子长度的影响。我们在同一经度上（即时角w定）考虑此问题，当地纬度Q与此时的太阳赤纬之差越大，影子长度r越大。最后，分析赤纬角5对影子长度的影响。我们在同一纬度不同经度上考虑此问题，当地经度和此时太阳直射点所在的经度之差越大，影子长度r越大。应用我们所建立的模型（4.4），基于Matlabplot函数给出9:0015:00的影子长度，见图（4.3）。其中整点时刻的影子长度见表4.1.表4.19:0015:00间整点时刻直杆影

13、子的长度北京时间9101112131415影长（米）7.81055.36524.26753.85853.98434.68746.2862351i1119101112131415时间（单位：t）图4.39:0015:00天安门广场直杆影子长度的变化曲线问题二的分析首先建立一个坐标系（标准坐标系），附件数据所在的坐标系为旧坐标系，标准坐标系与旧坐标系之间有一个旋转角。通过推导发现，影子长度与相应时刻太阳高度角的余切之比为定值，故可以直接利用旧坐标系下的影长数据，而避开了未知杆长、旋转角及标准坐标系下的坐标分量，从而简化了问题的复杂度和求解难度。问题二模型的建立及求解设a为i时刻的太阳高度角，卩为i

14、时刻太阳方位角（太阳方位角为正北方向按ii顺时针旋转到太阳投影点所旋转的角度）团,r.为i时刻影子长度，（x，y）为旧坐标i系下的杆顶点投影的坐标（旧坐标系为附件中数据所在的坐标系），建立标准坐标系，以杆低端为原点，x轴朝向正南，y轴朝向正东。（X,y）为标准坐标系下杆顶点投影的坐标。在标准坐标系下，有(4.7)hh二二tana丫X2+y2所以(4.8)h-cota-,:x2+y2又0为i时刻太阳方位角，则兀-卩为旧坐标系和标准坐标系之间的旋转角度,ii故(4.9)hcotasin(兀-0)-yTOC o 1-5 h z HYPERLINK l bookmark18 o Current Doc

15、ument viiihcotacos(兀-0)-xViiiy_sin(兀一0)_sin0ilxcos(兀一0)一cos0iii一tan0i(4.10)也就是直杆影子顶点横纵坐标的比值与杆子的高度无关，仅仅与太阳的方位角有关。hcota因为(4.11)所以h+y2i-cota恒定(4.12)由上面分析可知，新旧坐标系下的坐标，杆的高度，太阳高度角，太阳方位角，满足如下方程组：hcotacos(兀-0)xiiihcotasin(兀-0)yx_y_icos0一sin0iisin0 xcos0y(4.13)hcotai上述方程组为含多个未知量的非线性的方程，求解非常困难。但是由方程组可以得到如下关系式

16、：由此得到hcotah(4.14)cotaicotajir就2+y2jJ(4.15)若所求点恰好为附件所给数据的拍摄地点，必满足工(cota-r一cota-r=0(4.16)ijji从而建立目标函数如下：minf-S(cota-r一cota-r)(4.17)ijjii,j此目标函数和杆的高度、两个坐标系均无关，大大简化了模型的求解难度，加快了求解速度。为了找出最佳匹配的点，我们需要在球面上，计算出一些点的某段时刻的高度角，带入目标函数，求出最佳匹配点。在Matlab环境下，基于遍历思想，求出最佳匹配点如下：最优：北纬24度东经100度(云南省沧市临翔区)次优：北纬18度东经110度(海南省三亚

17、市)基于目标函数(4.17)，给出影子长度随时间变化的规律，见图(4.4)，其中不同时刻影子的长度见表4.2.北京时间/min图4.4附件1中影子长度的变化曲线表4.2不同时刻影子的长度北京时间12:4112:4412:4712:5012:5312:5612:5913:0213:0513:0813:11影子1.2471.2221.1981.1751.1521.1291.1071.0861.0651.0441.024长度25679592142944917835254081446264北京时间13:1413:1713:2013:2313:2613:2913:3213:3513:3813:41影子1

18、.0040.9850.9660.9480.9300.9130.8970.8800.8650.850长度6449179585928752109974492504问题三的分析在问题二中的目标函数里面，由于日期是已知的，即赤纬角是已知的；而问题三中，日期是未知的，即赤纬角是未知的。故太阳高度角和当地经纬度、日期均有关，这就导致所建的目标函数还是日期的函数，故需要在问题二的基础上，还须对日期进行遍历。问题三模型的建立及求解建立目标函数如下：minf=工(cota-r-cota-r)4ijijji(4.18)i,jsina=sinsin8+coscos8cos在Matlab平台下，基于遍历思想，求出最佳

19、匹配点如下：附件二的最优解为：日期1917月10日北纬41度东经80度新疆维吾尔自治区阿克苏地区乌什县；基于目标函数(4.18)，给出影子长度随时间变化的规律，见图(4.5)，其中不同时刻影子的长度见表4.3北京时间/min图4.5附件2中影子长度的变化曲线北京时间/min表4.3不同时刻影子的长度北京时间12:4112:4412:4712:5012:5312:5612:5913:0213:0513:0813:11影子1.2471.2221.1981.1751.1521.1291.1071.0861.0651.0441.024长度25679592142944917835254081446264

20、北京13:1413:1713:2013:2313:2613:2913:3213:3513:3813:41时间影子1.0040.9850.9660.9480.9300.9130.8970.8800.8650.850长度6449179585928752109974492504附件三的最优解为：日期613月02日北纬37度东经110度陕西省延安市延川县基于目标函数(4.18)，给出影子长度随时间变化的规律，见图(4.6)，其中不同时刻影子的长度见表4.4图4.6附件2中影子长度的变化曲线表4.4不同时刻影子的长度北13:0913:1213:1513:1813:2113:2413:2713:3013:

21、3313:3613:39京时间影3.5333.5463.5613.5783.5953.6143.6353.6573.6803.7053.731子长度142768798101751934426218541168278北京时间13:4213:4513:4813:5113:5413:5714:0014:0314:0614:09影子长度3.7589183.7880883.8187013.850813.8845853.9199123.9568763.9955354.0357514.077863问题四的分析由于影子长度的数据未知，所以首先需要对视频进行逐针读取，再灰度化处理，利用Photoshop软件提取

22、影子长度数据。在该问题中，由于给出了直杆高度的估计值，我们便可以进一步优化目标函数，找出最佳匹配点。问题四模型的建立及求解问题四给出了直杆高度的估计值，基于问题二和问题三的分析，我们可以进一步优化目标函数为：minf=Y(cota-r-h)2(4.19)i求出最佳匹配点如下：山西省大同市南郊区蒙古陕西省榆林市神木县陕西省榆林市府谷县山西省忻州市神池县河北省张家口市阳原县河北省张家口市崇礼县sina=sinsin8+coscos8cos在Matlab环境下，基于遍历思想，最优：北纬40度东经113度次优：北纬44度东经110度北纬39度东经110度北纬39度东经111度北纬39度东经112度北纬

23、40度东经114度北纬41度东经115度若拍摄日期未知，可以得到最优解，即可以确定拍摄地点和日期，日期和地点为：日期2076月22日北纬40度东经115度河北省张家口市蔚县模型的改进与评价问题一模型的改进与评价原问题一中所建的模型的求解比较复杂，把地球看成了一个圆球体，而地球是一个椭球体，所以该模型舍去了一些高阶微小量，只是给出了一个近似模型。我们需要注意，可以求得任何一地点在任何时刻的太阳直射点，一般情况下，该直射点和直杆处的经度和纬度是不同的，而我们假设太阳光线为平行光，基于此，可以把直杆处的光线沿着经线的切线方向和纬线的切线方向进行分解，分别得到两个方向的单位化向量，然后将经线切线的单位

24、向量除以直杆处纬度和此时赤纬角的差绝对值的余弦，得到一向量；再将纬线切线的单位向量除以直杆处经度和直射点处的经度的差绝对值的余弦，得到另外一向量，分别将这两个向量加权合成光线向量，进一步可得影子的长度。假设A处为太阳直射点，B处为直杆所在点，把直杆处的光线沿着经线的切线方向和纬线的切线方向进行分解,即分别在同一经度不同纬度下，同一纬度不同经度下考虑，如图5.1所示，假设北纬为正，南纬为负。首先考虑，在同一经度不同纬度下的情况，如图5.2所示。直杆处纬度和赤纬之差的绝对值为Y=|0-5假设在B点建立空间直角坐标系如下：如图5.3所示。再考虑同一纬度不同经度的情况，如图5.4所示。假设在D点建立空

25、间直角坐标系如下：如图5.5所示。设V为单位向量，因为太阳光朝向Z轴负半轴，所以V=(0,sina,-cosa)11设V为单位向量，因为太阳光朝向Z轴负半轴，所以V=(sin,0,-coso).22图5.2直杆处光线经线切线分解示意图(九2sin,九sina,一九cosa一11九coso)2tanA=cosa一九coso2九2sin2o+九2sin2a21其中，A为V与Z轴夹角，r=htanA.合成光线向量如图5.6所示。V=(cosO-sino,sin0-sin,cos)=(sino,0,coso)图5.6C处空间直角坐标系5.2问题二、三、四模型的改进与评价缺点：原模型由于所给数据为一小时内间隔的数据，所以模型不能很好的解决问题，只能得出最优解，