数学建模2015解析.doc

关于太阳影子定位的数学模型 摘要 利用球面天文学和空间解析几何等知识，并运用matlab、lingo等软件进行数值求解，解决了关于物体的太阳影子变化确定拍摄地点的经纬度及拍摄日期时间的问题。 针对问题一，通过引进太阳高度角、赤纬角等变量，得到太阳影子的变化与经度，纬度，杆长和时间的函数关系，求出直杆在不同时刻所对应的影长，利用matlab编程得到直杆影长的变化曲线，曲线大致呈对称形状，物体影长由长变短，在中午12点14分左右直杆的影长达到最短且为3.8米左右，12点之后逐渐变长。 针对问题二，运用最小二乘法将附件1所给的影子的顶点坐标进行拟合，通过误差分析取出最优的两个坐标点，并结合几何知识，建立模型求解出直杆的高度为2.04米。我们运用matlab对影子的长度进行二次拟合，求出一天中影子最短出现的时间，并求解出当地经度为。根据附件1，结合问题一中的数学模型，得到当地的纬度为，所对应的地区为海南。由于最优坐标的选取不同和三角函数的周期性，可以得到不同的地点。 针对问题三，对附件2的数据进行椭圆曲线拟合，应用地平式日晷原理，计算出当地的纬度为，经度为，数据记录的时间为2015年10月20日，所对应的地区为西藏自治区的阿里地区。同理，附件三对应的地区在越南附近，日期为2015年8月22日。 针对问题四，使用matlab将视频文件转化成帧数，再设定三分钟取出一个帧数，之后对这些图像进行灰度处理，根据它们的像素点从而来确定杆底、杆高及影子顶点的坐标，从而计算出影子的长度，用Excel统计出影长随时间变化的二次函数，根据问题二的模型求解出当地的坐标为（，），位置为内蒙古自治区。 关键词： 球面天文学 空间解析几何 拟合 地平式日晷 灰度处理 一、问题重述 1.1问题背景 所谓太阳影子定位技术，就是通过分析同一地点的同一物体在不同时间的太阳影子变化规律，从而确定出视频拍摄的地点和日期。由于各物体在不同时间或者不同地点的太阳高度角不同，因此各物体对应的影长也会有所不同，因此需要我们对研究太阳影子与时间，地点等参数的关系，建立出相应的数学模型，从而确定物体的拍摄地点及拍摄日期。 1.2 解决问题 建立数学模型，解决下列问题 1. 分析影子长度的变化规律，建立出在已知日期的直杆影子长度变化的数学模型，并应用此模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。 2.已知某固定直杆在水平地面上的太阳影子顶点坐标数据及测量日期，未知直杆的高度，建立相应数学模型从而确定直杆所处的地点。并应用于附件1数据，求解出若干个可能的地点。 3.根据固定直杆在水平地面上的坐标建立数学模型，确定直杆所处的地点和日期。并分别应用于附件2和附件3的数据，解出若干个可能的地点与日期。 4．根据附件4中已知拍摄的日期的一段视频，给出数学模型，估计出直杆的高度为2米。并给出若干个可能的拍摄地点。并且在拍摄日期未知的条件下，建立相关模型，确定出拍摄地点与日期。 模型假设 （1）假设收集数据的过程中天气良好，没有乌云或者其他遮挡物对太阳阴影的顶点坐标产生影响； （2）假设直杆在数据记录的过程中一直处于直立状态，没有受到风吹导致直杆倾斜等影响； （3）假设太阳从日出到日落一直在进行匀速运动，转过的角度为，且经过的时间为12个小时； （4）不考虑大气层的折射对太阳光照的影响。 （5）忽略对影长及其他数据测量时造成的测量误差。 （6）坐标系以直杆底端为原点O，水平地面为xOy面，Z轴垂直于xOy面，x轴，y轴，z轴满足右手法则。 三、名词解释与符号说明 根据球面天文学的原理，本文可以用到以下名词： （1）太阳高度角：通俗地讲，就是太阳光的入射方向和地平面之间所形成的夹角；而从专业的角度来说，则指某地太阳光线与通过该地与地心相连的地表切线的夹角，如图1中的。 （2）太阳时角：指的是日面中心的时角，换言之就是从观测点天球子午圈沿天赤道量至太阳所在时圈的角距离。 （3）真太阳时：太阳视圆面中心连续两次上中天的时间间隔叫做真太阳日。一个真太阳日又分为24真太阳时.这个时间系统称为真太阳时。 （4）赤纬角：赤纬角又称为太阳赤纬，它是指地球赤道平面与太阳和地球中心的连线之间的夹角，如图1中的。 本文还会用到如下的符号： ：当地纬度，如图1中的； ：直杆的高度； ：太阳影子的长度； ：直杆的顶点与其太阳影子顶点的夹角； ：太阳转过的角度； 利用软件soildworks建立相应的三维模型，则几个变量之间相对应的关系如下图1： 图1 球体、、、的关系图 问题分析 对于问题一，需要我们建立出影子长度的变化曲线，经资料差知，北京所对应的时区为东八区