2015高教社杯全国大学生数学建模竞赛题目影子定位.docx

《2015高教社杯全国大学生数学建模竞赛题目影子定位.docx》由会员分享，可在线阅读，更多相关《2015高教社杯全国大学生数学建模竞赛题目影子定位.docx（16页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、2015高教社杯全国高校生数学建模竞赛题目（请先阅读“全国高校生数学建模竞赛论文格式标准”）A题 太阳影子定位如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子改变，确定视频拍摄的地点和日期的一种方法。1.建立影子长度改变的数学模型，分析影子长度关于各个参数的改变规律，并应用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的改变曲线。2.依据某固定直杆在程度地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模

2、型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。3. 依据某固定直杆在程度地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将你们的模型分别应用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。4附件4为一根直杆在太阳下的影子改变的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。假如拍摄日期未知，你能否依据视频确定出拍摄地点与日期？太阳影子定位摘要 本文通过分析物体的太阳影子改变，利用太阳影子定位技术建立确定视频拍摄的地点和日期的模型。针对问题一，首先通过分析知影子长度的改变主要

3、影响参数为：当地的经度 、纬度、时刻、直杆长度、季节（日期)等，引入地理学参数：太阳赤纬、时角及太阳高度角，建立一个可以刻画影子长度改变和各个参数间关系的模型：；其次以实例对模型进展检验，在误差可允许的范围内，认为模型正确；进而对模型采纳限制变量法分析影子长度关于各个参数的改变规律；然后求解出满意条件影子长度12时15分是最短，大约3.674米（表3）。影子长度的改变曲线（图5），9时至12时15分影子长度呈现下降趋势，12时15分之15时影子长度呈现上升趋势；最终考虑太阳照耀中发生折射现象的推广。 针对问题二， 关键词一、 问题重述：如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，

4、太阳影子定位技术就是通过分析视频中物体的太阳影子改变，确定视频拍摄的地点和日期的一种方法。1.建立影子长度改变的数学模型，分析影子长度关于各个参数的改变规律，并应用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的改变曲线。2.依据某固定直杆在程度地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。3. 依据某固定直杆在程度地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将你们的模型分别应

5、用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。4附件4为一根直杆在太阳下的影子改变的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。假如拍摄日期未知，你能否依据视频确定出拍摄地点与日期？二、 问题分析：针对问题一，首先通过分析知影子长度的改变主要影响参数为：当地的经度 、纬度、时刻、直杆长度、季节（日期)等，引入地理学参数：太阳赤纬、时角及太阳高度角，建立一个可以刻画影子长度改变和各个参数间关系的模型；其次以实例对模型进展检验，在误差可允许的范围内，认为模型正确；进而依据所建模型分析影子长度关于各个

6、参数的改变规律；然后做出影子长度的改变曲线；最终考虑太阳照耀中发生折射现象的推广。针对问题二，三、模型假设：1、求解此问题时忽视地球的自转2、不考虑太阳光线在穿过大气层时的折射、太阳的视面角、高山阻挡、海拔高度等因素。3、认为照耀到地球上的太阳光可以看成是平行光线，地球上某地的程度地面是地球球面上过该地的切面。四、符号说明：太阳赤纬：太阳时角：表示某地的地理纬度：表示某地的地理经度：太阳高度角五、问题一的模型建立与求解5.1影响影子长度参数确实定1、 太阳赤纬 太阳赤纬是地球赤道平面与太阳和地球中心的连线之间的夹角且以年为周期，在周年运动中任何时刻的赤纬值都是严格已知的，可用下式计算：式中称为

7、日角，即，这里,式中为积日，就是日期在年内的依次号（例如1月1日其积日为1，平年1212月31日的积日为365，闰年则为366等）。（式中表示取整数局部）故假如已知某日期的年、月、日，代入式（1），即可求得此日期的太阳赤纬值。2、 时角时角表示一天体是否通过了当地的子午圈，其数值表示该天体与当地子午圈的角间隔 ，并借用时间的单位以小时来计量，其中当地时间12点时的时角为零，令上午的时角为正，下午为负。某地的时角计算式如下：故若给定某地时刻的值，代入式（2），即可求出此时的时角值3、 太阳高度角 太阳高度角,是指太阳光线与地平面的夹角。应用球面三角形余弦公式，结合图形,可以推出随意时刻太阳高度角

8、的计算公式为： 进一步可以得到：式中，表示当地的地理纬度,表示太阳赤纬，表示太阳时角、表示太阳高度角。、的取值为北正南负。故若已知某地的地理纬度、太阳赤纬、太阳时角，代入式（4），即可求得太阳高度角。4、 直杆的长度直杆的影子始终在物体背着光源的一面，光从物体顶端照耀到地面形成影子，直杆的长度干脆影响着影子长度的改变。5.2模型一的建立 首先分析影子长度改变干脆受直杆长度和太阳高度角的影响，其关系式为： 其中表示影子的长度。 然后以太阳高度角与太阳赤纬、时角、地理纬度的关系为根底，即式（4），引进参数当地的经度，构建影子长度改变模型如下： 式中、分别代表当地的经度、纬度、太阳赤纬、太阳高度角和

9、时刻（北京时间）。编写程序计算时留意到，当时，反余弦函数值才存在，应采纳推断，当其大于时则进展下一个的计算。 5.3模型一的检验 以实际某地状况为例，已知当地的经度为，纬度为， 5.4影子长度关于各个参数的改变规律 依据影子长度改变模型，以影子长度为因变量，依次选择各个参数为自变量，其余参数看为固定值，进而描绘出影子长度关于各个参数改变规律。 1、影子长度关于直杆长度的改变规律 以直杆长度为自变量，以影子长度为因变量，太阳高度角为固定值，取直杆长度为米，间隔为0.5米，依据公式（5），得到与之间的关系，画出关系图像如图1：图1：影子长度与直杆长度的关系由图1我们可以很直观的的看出，在其他参数不

10、变时，直杆长度越长，影子长度越长，且两者的比值是不变的。 2、影子长度的日改变规律以每日的时间作为自变量，影子长度为因变量，其它参数为固定值，取为2015年10月22日北京时间8:00-16:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆，依据公式（6），得到影子长度的日改变规律。图2：影子长度和日的关系由于地球是个固体球,且自西向东自转,因此太阳高度也呈现出与之对应的日改变规律:地球上视察太阳为东升西落,早晨太阳从东方地平线上升起,晨昏线上太阳高度为，随着太阳的渐渐上升,太阳高度是渐渐增大的，影子的长度渐渐减小。当某地经线正对太阳光时,地方时为正午12点，

11、即为北京时间12时15分,此时太阳高度角到达一天中的最大，影子长度最小。之后太阳渐渐西落,太阳高度也渐渐变小，影子渐渐变长,到西方地平线落下时,没有影子。 3、影子长度随纬度改变规律 以纬度作为自变量，影子长度为因变量。其它参数为固定值，由于正午太阳高度角最能反映太阳辐射的强弱改变,故取2015年10月22日北京时间12时东经116度23分29秒3米高的直杆，依据公式（6），得到影子长度的随纬度的改变规律。图3：影子长度和纬度的关系正午太阳高度由直射点向南北两侧递减，故影子长度由直射点向南北两侧递增。夏至:，影子长度从北回来线向南北两侧递增;同样冬至:，影子长度从南回来线向南北两侧递增;春秋分

12、:，影子长度从赤道向南北两侧递增,且离直射点间隔 越近,与直射点纬度差越小,影子长度就越小。4、 影子长度随季节的改变规律 以季节作为自变量，影子长度为因变量。由于正午太阳高度角最能反映太阳辐射的强弱改变,故取2015年北京时间12时天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆，依据公式（6），得到影子长度的随季节的改变规律。图4：影子长度和季节的关系 北回来线及以北地区,6月22日正午太阳高度达一年中最大值，此时影子长度最小,12月22日达一年中最大值;南回来线以南地区,12月22日达一年中最小值,6月22日达一年中最大值;南北回来线之间,一年中因有两次太阳直射

13、时机,赤道至北回来线之间12月22日的正午影子长度达一年中最大值,而赤道至南回来线之间6月22日达一年中最大值。 5.5问题一的求解求解步骤：1）将2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆符号化为，；2） 计算太阳赤纬： 已知年份为2015年，积日为295，代入公式得： 得到太阳的赤纬角为3） 计算太阳时角： 需将当地时间转化为北京时间，对应北京时间，将代，得到2015年10月22日北京时间天安门广场时角值（如表1），其中。表1：2015年10月22日北京时间天安门广场时角统计表时刻9:009:3010:0

14、010:3011:0011:1211:24时角/度-48.61-41.11-33.61-26.11-18.61-15.61-12.61时刻11:3611:4812:0012:1212:1512:2412:36时角/度-9.61-6.61-3.61-0.610.142.395.39时刻12:4813:0013:3014:0014:3015:00时角/度8.3911.3918.8926.3933.8941.394） 计算太阳高度角：将所得到的、分别代入公式（4），得到2015年10月22日北京时间天安门广场太阳高度角值（如表2）表2：2015年10月22日北京时间天安门广场太阳高度角统计表时刻9:

15、009:3010:0010:3011:0011:1211:24太阳高度角/度22.1626.5330.4333.7536.3737.2037.89时刻11:3611:4812:0012:1212:1512:2412:36太阳高度角/度38.4538.8639.1239.2239.2339.1838.98时刻12:4813:0013:3014:0014:3015:00太阳高度角/度38.6338.1436.2933.6333.8926.375）求解影子长度改变： 应用建立的模型 将就得各个参数值代入模型，可以得到2015年10月22日北京时间9:00-15:00天安门广场3米高直杆的影子长度（见

16、表3），通过Matlab软件画出直杆的太阳影子长度的改变曲线（如图5）表3：2015年10月22日北京时间天安门广场3米高直杆的影子长度统计表时刻9:009:3010:0010:3011:0011:1211:24影长/米7.3676.0095.1074.4904.0733.9523.854时刻11:3611:4812:0012:1212:1512:2412:36影长/米3.7783.7243.6893.6753.6743.8543.778时刻12:4813:0013:3014:0014:3015:00影长/米3.7543.8214.0864.5095.1356.050图5：2015年10月22

17、日北京时间天安门3米高直杆的影子长度改变曲线结合图5和表3，可以看出2015年10月22日北京时间12点15分天安门3米高直杆的影子长度是最短的，大约3.674米。9时至12时15分影子长度呈现下降趋势，12时15分至15时影子长度呈现上升趋势。六、 问题二的模型建立与求解6.1直杆影长的计算已知某固定直杆在程度地面上的太阳影子顶点坐标，依据勾股定理表示出影长：6.3问题二的求解：依据附件1中北京时间14：42-15:42之间影长的横、纵坐标，通过勾股定理导出：计算出21个时刻对应的影长（见表4）表4：中北京时间当地该直杆影子长度统计表时刻14:4214:4514:4814:5114:5414

18、:5715:00影长/米1.1501.1821.2151.2491.2831.3181.353时刻15:0315:0615:0915:1215:1515:1815:21影长/米1.3891.4261.4631.5011.5401.5801.620时刻15:2415:2715:3015:3315:3615:3915:42影长/米1.6611.7031.7461.7901.8351.8811.9286.2经度确定模型6.2.1确定地区经度与时差的关系地球从经线划分，向东为：,向西为：，地球一周就是：。地球自转一周的时间为，所以可得地球每小时转过的角度为：。6.2.2建立二次拟合模型依据附件一算出的

19、影长，画出时刻-影长的改变曲线，如图6：图6：附件一所求影长与时刻关系图 以影长为因变量，时刻为自变量，用函数进展二次拟合，得到的二次函数为：。式中为影长，为时刻，去自变量画图，得到影长与时刻的二次拟合曲线，图7：影长与时刻的二次拟合曲线6.2.3确定当地经度因为二次拟合的曲线为抛物线，依据抛物线函数的对称轴公式：求解得到，可以断定该地在北京的西边，因为北京所在经度为，所以确定该地所在经度为。6.3纬度确定模型6.3.1影长比例不变原理图8：同一地区不同杆在不同时间下的影子示意图意图PQAMNBB杆A杆 表示1时刻的太阳光线 表示2时刻的太阳光线因为太阳光是平行的，所以,又因为均为直角三角形，

20、所以，得到，同理得到，进而得到。6.3.2遍历模型依据附件中的得到的21个影长，用后一时刻的影长比上前一时刻的影长得到20个比值。附件1中各相邻时刻影长比例统计表标号1234567比值1.028331.028001.027771.027341.027121.026841.02662标号891011121314比值1.026461.026121.026021.025811.025721.025501.02538标号151617181920比值1.025291.025201.025111.025121.024991.02501对于问题二给定杆长为定值，在随意给定一个纬度值，依据公式（6）求得21个

21、影长，并求得20个影长比值。依据之前的比例不变原理，假如，认为两个杆处在同一位置，但由于附件一的影子顶点坐标的测量出来的，所以影长比值不行能完全相等，可能会存在肯定误差，所以做修正,定义如下：对全部纬度进展求解，每一个纬度值都对应一个,求全部中的最小值，得到纬度值，即该地所处纬度。求解结果，的最小值为：，所处纬度为：。6.3其余可能点确实定已知2015年春分日是3月21日，此时太阳直射赤道，即纬线，夏至日是6月22日，太阳直射北回来线，即，这两个节气的天数差为93天，大致算出每天太阳直射纬度的改变为：。附件一中所给时间为2015年4月18日，与春分日日差28天，估计此时太阳直射，令之前求出的纬

22、度线以纬度线为对称轴，求出对称的纬度线为：。假如时刻定在2015年4月18日的15时，依据之前的地区经度与时差的关系，可以得到此时刻太阳直射的经度线为：,令之前求出的某地的经度线以经度线为对称轴，求出对称的纬度线为：。依据附件一的数据可得：状况一：以2015年4月18日14:42-15:42这个时段为基准，直杆可能在的位置有两个：（）和（）。状况二：以2015年4月18日15:00这个时刻为基准，直杆可能在的位置有四个，分别是：（）、（）、（）和（）。以上我们得到了考虑大气折射率影响的太阳高度角，从肯定程度上减小了误差，使得结果更具合理性。参考文献：1但尚铭，太阳高度角计算程序介绍，环境科学，

23、10卷4期：76-78,1988年2赖月喜，太阳高度改变的一般规律及图示分析，中学政史地(高三)，04期，48-52,2007年3百度百科，太阳高度角， :/baike.baidu /linkurl=c6pown3qPkE7Bht77RfuvtRmar03oJRxL15Fw24a7OrMySjPIDp1W7vYNFY8M7R0On8a1BsiI9qAfLcr79sPi_#ref_1_86609，2015/9/114百度百科，时角， :/baike.baidu /view/823827.htm，2015/9/115百度百科，太阳赤纬, :/baike.baidu /view/86281，2015/9/116