全国大学生数学建模比赛A题一等奖论文Word文档格式.docx
《全国大学生数学建模比赛A题一等奖论文Word文档格式.docx》由会员分享，可在线阅读，更多相关《全国大学生数学建模比赛A题一等奖论文Word文档格式.docx（50页珍藏版）》请在冰豆网上搜索。
北纬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两地经度差
太线与直杆的夹角
直杆长度
米
直杆影长
地方时
时间
直杆所在地的经度
时
问题二、三符号系统
附件1中第一组坐标的y值
极径
极角
度
问题四符号系统
固定杆长度
实际长度与灰度值坐标下的转换比例
投影系统
4.问题一的建模与求解
4.1问题分析
在问题一中，为了描述直杆影子长度变化的动态过程，首先以直杆为z轴，建立空间三维坐标对直杆影子的变化进行数学抽象。
再将地球作为规则球体建立球面坐标系，利用空间解析几何与平面解析几何的知识，对两个坐标系中的相关向量与角度进行分析，分析出影响影子长度的参数，得到影子端点在坐标系中的位置表达式。
由此可以求出影子长度随各个参数的变化规律。
建模流程图如下所示：
图4.1问题模流程图
4.2模型准备
为了建模的方便，先给出一些地理名词的解释和一些数据的预处理方法。
4.2.1名词解释[6]
地方时：
以一个地方太阳升到最高的地方时间为正午12时，将连续两个正午12时之间等分为24个小时，所成的时间系统。
它是观测者所在的子午线的时间。
时间：
是中国采用时区的区时作为标所在的东八准时间。
时间并不是（东经116.4°
）地方的时间，而是东经120°
地方的地方时间。
太阳赤纬：
是地球赤道平面与太阳和地球中心的连线之间的夹角。
太阳直射点：
地球表面太射入角度（即太阳高度角）为90度的地点，它是地心与日心连线和地球球面的交点。
太阳高度角：
对于地球上的某个地点，太阳高度角是指太的入射方向和地平面之间的夹角；
专业上讲是指某地太线与通过该地与地心相连的地表切线的夹角。
4.2.2数据预处理
（1）经纬度转换
在问题一中，天安门广场的坐标是用经纬度（度分秒）的形式给出的。
为了下面建模求解的方便，将其统一转换成以“度”为单位。
换算方法为：
分位数除以60，秒位数除以3600。
所以，天安门广场的纬度可以转换为：
经度可以转换为：
（2）时间与地方时的转换[9]
问题中所给出的时刻为时间，而时间指的是东经120°
地方的地方时，并不是问题中地点的地方时。
所以先要将所给的时间转换成相应的地方时。
转换规则为：
东经度<
120度地区，每减少1度，减4分钟；
东经度>
120度地区，每增加1度，加4分钟。
所以有转化公式：
其中，E表示直杆所在地点的经度，是时间，是直杆所在地方的地方时。
用此公式对问题一中的时间进行操作，得到直杆所在地的地方时，如下表所示：
表4.1天安门的地方时与时间的转换
9:
00
30
10:
11:
12:
8:
75
25
13:
14:
15:
4.3模型的建立
要研究影子的变化，需要建立空间三维坐标对直杆影子的变化进行数学抽象。
通过对直杆和地球分别建立了两个空间直角坐标系，用空间解析几何和向量知识，可以确定两个坐标系上各点之间的位置和角度关系。
4.3.1建立直杆处空间三维坐标系
根据假设，视太线为平行光，以直杆所在地点的正向为轴，以正北方向为轴，以直杆直立即垂直于平面的方向为轴，建立空间直角坐标系，得到直杆在平面的投影与光线的位置关系，如下图所示：
图4.2直杆空间三维坐标系
其中，是与过A处的经线相切的方向向东的单位向量；
是A处地平面方向向北的单位向量。
AH是A处垂直于平面的直杆，AF是该直杆在平面的投影，HF是当天太线的照射方向，照射方向与直杆所成角度。
4.3.2建立直杆在地球上的宏观空间球面坐标系
根据假设，可视地球为规则球体，过直杆底端A处的经线与赤道交于D点，B点为某日的太阳直射点，过B点的经线与赤道交于C点。
以O为原点，以OD所在直线为x轴，以地轴ON所在直线为z轴建立空间直角坐标系，如图4.3所示：
图4.3直杆在地球上的空间三维坐标系
4.3.3确定各点之间的位置和角度关系
设地球半径为R，，则有：
（1）为直杆所处位置的纬度数，并且。
若A地在北半球，则，若A地在北半球，则。
（2）为太阳直射点B地的纬度，亦即上面提到的赤纬，并且。
（3）为A地与B地的经度差，t是地方时。
对于某日A地白昼t时刻：
。
（4），证明过程如下：
由假设可知，太线是一簇簇的平行线，所以，如图4.4，圆O’是过A,B两地的大圆，于是，证毕。
由以上分析可得：
图4.4过A、B两地，以地球中心为圆心的圆
升级会员 