2015年全国大学生数学建模比赛a题一等奖论文.docx

全国大学生数学建模竞赛一等奖论文太阳影子定位问题摘要目前，如何确定视频的拍摄地点和拍摄日期是计算机视觉的热点研究问题，是视频数据分析的重要方面，有重要的研究意义本文通过建立数学模型，给出了通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的方法对于问题一，建立空间三维直角坐标系和球面坐标系对直杆投影和地球进行数学抽象，引入地方时、北京时间、太阳赤纬、杆长、太阳高度角等五个参数，建立了太阳光下物体影子的长度变化综合模型求解过程中，利用问题所给的数据，得到太阳赤纬等变量，将太阳赤纬等参量代入模型，求得了北京地区的 9:00 至 15:00 的影子长度变化曲线，当 12：09 时，影子长度最短；并分析出影长随这些参数的变化规律，利用控制变量法思想，总结了五个参数与影子长度的关系最后进行模型检验，将该模型运用于东京、西藏两地，得到了这两座城市的影长变化规律曲线，发现变化规律符合实际两地实际情况对于问题二，为了消除不同直角坐标系带来的影响，将实际坐标转换为二次曲线的极坐标，建立了极坐标下基于多层优化搜索算法的空间匹配优化模型求解时，先将未知点的直角坐标系的点转换为极坐标，然后设计了多层优化搜索算法，通过多次不同精度的搜索，最后得出实际观测点的经纬度为东经 E 北纬 N 。

同时对模型152进行验证，实地测量了现居住地的某个时间段的值，通过模型二来求解出现居住地的经纬度，分析了误差产生的原因：大气层的折射和拟合误差对于问题三，将极坐标转换后的基本模型转换为优化模型，建立了基于遗传算法的时空匹配优化模型将目标函数作为个体的适应度函数，将经度纬度及日期作为待求解变量，用遗传算法进行求解，得到可能的经度纬度及其日期：北纬 20 度，东经114 度，5 月 21 日；北纬 20 度，东经 114 度，7 月 24 日；东经 94.5 度，北纬 33.8 度，6 月 19 日最后，将遗传算法与多层优化搜索算法进行对比分析，得出遗传算法的求解效率和求解精度均优于多层次搜索算法对于问题四，首先将视频材料以 1min 为间隔进行采样得到 41 帧（静态图片） ，将这些静止图片先利用 matlab 进行处理，后进行阀值归一化处理，得到这些帧的灰度值矩阵在图片上建立参考模型，获得影子端点的参考位置利用投影系统和模型二，建立了基于图形处理的视频拍摄地点搜索模型利用模型二中多层搜索算法，求得满足精度的最优地点最优的地点是：东经 119,北纬 48.7，在内蒙古的呼伦贝尔市同时假设日期是未知量，将模型四与模型三相结合，得到了可能的地点和时间，并分析了可能出现误差的原因，最后回答了当视频日期未知，也可以确定其位置和日期。

最后，给出了模型的优缺点和改进方案关键词：极坐标化，多层优化搜索算法，遗传算法，图像处理，MATLAB全国大学生数学建模竞赛一等奖论文1. 问题重述1.1 问题背景随着现代科技的发展，日常生活中摄像机的应用越来越普遍无论是个人家庭还是组织单位，都通过摄像机来录制各种视频以分享信息，例如实时视频监控、记录自然景观、观测气象信息等而通过视频来确定拍摄地点的地理位置信息是目前计算机视觉领域的热点研究问题之一一个视频的地理位置能够提供当地气候、平均温度、平均降雨量、植物索引、地表概况、海拔高度和人口密度等大量背景信息[1]因此从视频中确定地理位置是一项有很大潜力应用空间的技术1.2 问题描述视频数据分析是视频处理过程中的重要环节，而如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法试建立数学模型讨论下列问题：1.建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用所建立的模型画出 2015 年 10 月 22 日北京时间 9:00-15:00 之间天安门广场 3 米高的直杆的太阳影子长度的变化曲线。

2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点将模型应用于附件 1 的影子顶点坐标数据，给出若干个可能的地点3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期将模型分别应用于附件 2 和附件 3 的影子顶点坐标数据，给出若干个可能的地点与日期4．附件 4 为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为 2 米请建立确定视频拍摄地点的数学模型，并应用此模型给出若干个可能的拍摄地点如果拍摄日期未知，能否根据视频确定出拍摄地点与日期？2. 问题分析2.1 问题一分析问题一要求分析投影长度随各参数的变化规律，建立影子长度变化的数学模型首先对直杆建立空间三维坐标系，将地球简化成规则球体建立球面坐标系在这两个坐标系中，通过几何证明，运用向量知识可分析出影响影子长度的各种参数，得出地球上某日白天某时刻影子顶端在地平面上的具体位置，由此可以给出影子长度的变化规律2.2 问题二分析问题二要求根据某固定直杆在水平地面上的太阳影子顶点坐标数据及日期数据，建立数学模型确定直杆所处的地点与第一问有相似之处，但分析附件所给数据，发全国大学生数学建模竞赛一等奖论文现附件中只给出 x、y 坐标值，而并没有给出 xy 轴的准确方向，所以考虑将直角坐标转换成极坐标，来消除由于不同坐标系选取所造成的影响。

2.3 问题三分析问题三与问题二有相似处，区别是第三问附件没有提供日期，需要根据直杆影子端点坐标确定直杆所在地点的经纬度和日期具体的日期可以由太阳直射点纬度来确定，而根据问题二中的模型，xy 坐标与太阳直射点纬度有关如果继续用第二问的模型来求解，需要不断改变太阳直射点纬度来拟合极坐标方程，这样做算法复杂度会很大所以考虑对问题二模型进行修改，不采用拟合，而直接建立与待求点经纬度以及日期有关的目标函数，通过约束经纬度范围来缩小待求点的可行域，从而简化算法复杂度 2.4 问题四分析问题四中，直接以视频的方式给出了固定杆长的距离变化规律将图片形式的影长变化规律以坐标的形式进行转换，转换为现实的坐标形式这样就可以利用问题二的模型，整合现有的算法，求出拍摄地点3. 模型假设与符号系统3.1 模型的假设（1）假设地球为一个规则的球体2）由于日地距离远大于地球半径，所以假设太阳光线为平行光3）假设地球上某地的水平地面是地球球面上过该地的切面4）假设不考虑太阳光线穿过大气层时所发生的折射5）假设一天中太阳直射点的纬度不变6）假设不考虑太阳的视面角、高山阻挡、海拔高度等因素的影响7）假设不考虑阴天没有阳光的情况。

3.2 符号系统问题一符号系统符号 意义 单位直杆所在地纬度值 度太阳直射点的纬度 度A、B 两地经度差 度太阳光线与直杆的夹角 度h直杆长度 米L直杆影长 米t 地方时 时全国大学生数学建模竞赛一等奖论文0tE北京时间直杆所在地的经度时度问题二、三符号系统符号 意义 单位直杆所在地纬度值 度太阳直射点的纬度 度1y附件 1 中第一组坐标的 y 值 米极径h直杆长度 米极角 度问题四符号系统符号 意义 单位L固定杆长度 米k实际长度与灰度值坐标下的转换比例 P投影系统 4. 问题一的建模与求解4.1 问题分析在问题一中，为了描述直杆影子长度变化的动态过程，首先以直杆为 z 轴，建立空间三维坐标对直杆影子的变化进行数学抽象再将地球作为规则球体建立球面坐标系，利用空间解析几何与平面解析几何的知识，对两个坐标系中的相关向量与角度进行分析，分析出影响影子长度的参数，得到影子端点在坐标系中的位置表达式由此可以求出影子长度随各个参数的变化规律建模流程图如下所示：全国大学生数学建模竞赛一等奖论文图 4.1 问题一建模流程图4.2 模型准备为了建模的方便，先给出一些地理名词的解释和一些数据的预处理方法。

4.2.1 名词解释[6]地方时：以一个地方太阳升到最高的地方时间为正午 12 时，将连续两个正午 12时之间等分为 24 个小时，所成的时间系统它是观测者所在的子午线的时间北京时间：是中国采用北京时区的区时作为标所在的东八准时间北京时间并不是北京（东经 116.4°）地方的时间，而是东经 120°地方的地方时间HF太阳赤纬：是地球赤道平面与太阳和地球中心的连线之间的夹角太阳直射点：地球表面太阳光射入角度(即太阳高度角 )为 90 度的地点，它是地心与日心连线和地球球面的交点 太阳高度角：对于地球上的某个地点，太阳高度角是指太阳光的入射方向和地平面之间的夹角；专业上讲是指某地太阳光线与通过该地与地心相连的地表切线的夹角4.2.2 数据预处理（1）经纬度转换在问题一中，天安门广场的坐标是用经纬度（度分秒）的形式给出的为了下面建模求解的方便，将其统一转换成以“度”为单位换算方法为：分位数除以 60，秒位数除以 3600所以，天安门广场的纬度可以转换为：395426=+546023=9.07经度可以转换为：1116（2）北京时间与地方时的转换[9]全国大学生数学建模竞赛一等奖论文问题中所给出的时刻为北京时间，而北京时间指的是东经 120°地方的地方时，并不是问题中地点的地方时。

所以先要将所给的北京时间转换成相应的地方时转换规则为：东经度120 度地区，每增加 1 度，加 4 分钟所以有转化公式： 0(E2)*,0=2t其中，E 表示直杆所在地点的经度， 是北京时间， 是直杆所在地方的地方时 0tt用此公式对问题一中的北京时间进行操作，得到直杆所在地的地方时，如下表所示：表 4.1 天安门的地方时与北京时间的转换北京时间 9:00 9:30 10:00 10:30 11:00 11:30 12:00地方时 8:75 9:25 9:75 10:25 10:75 11:25 11:75北京时间 12:30 13:00 13:30 14:00 14:30 15:00地方时 12:25 12:75 13:25 13:75 14:25 14:754.3 模型的建立要研究影子的变化，需要建立空间三维坐标对直杆影子的变化进行数学抽象通过对直杆和地球分别建立了两个空间直角坐标系，用空间解析几何和向量知识，可以确定两个坐标系上各点之间的位置和角度关系4.3.1 建立直杆处空间三维坐标系根据假设，视太阳光线为平行光，以直杆所在地点的正东方向为 轴，以正北方x向为 轴，以直杆直立即垂直于 平面的方向为 轴，建立空间直角坐标系，得到直y0xyz杆在 平面的投影与光线的位置关系，如下图所示：0x全国大学生数学建模竞赛一等奖论文图 4.2 直杆空间三维坐标系其中， 是与过 A 处的经线相切的方向向东的单位向量； 是 A 处地平面内E K方向向北的单位向量。

AH 是 A 处垂直于 平面的直杆，AF 是该直杆在 平面内xOyxOy的投影，HF 是当天太阳光线的照射方向，照射方向与直杆所成角度 FH4.3.2 建立直杆在地球上的宏观空间球面坐标系根据假设，可视地球为规则球体 ，过直杆底端 A 处的经线与赤道交于 D 点，B1O点为某日的太阳直射点，过 B 点的经线与赤道交于 C 点以 O 为原点，以 OD 所在直线为 x 轴，以地轴 ON 所在直线为 z 轴建立空间直角坐标系 ，如图 4.3 所示：xyz图 4.3 直杆在地球上的空间三维坐标系4.3.3 确定各点之间的位置和角度关系设地球半径为 R， ，则有：1,,,AODBCDOAB（1） 为直杆所处位置的纬度数，并且 90若 A 地在北半球，则 ，若 A 地在北半球，则 0（2） 为太阳直射点 B 地的纬度，亦即上面提到的赤纬，并且36''（3） 为 A 地与 B 地的经度差，t 是地方时对于某日 A 地白昼 t 时刻： 12)5(024)。