数学建模太阳影子定位.doc

word西 安 邮 电 大 学(理学院)数学建模报告题目： 太阳影子定位问题 班级：信息工程1403班 学号： 03144079 ： 侯思航 成绩：2016年6月30日 一、 摘要 本文针对太阳影子定位技术，通过太阳与地球相对运动的规律，建立杆长、影长、经纬度、时间、日期的关系，建立模型综合分析了不同地点，不同的时间，不同的季节时影子长度的形成规律与变化趋势，运用了软件进展分析，得出不同地区影子变化的模型最后将具体情况运用到建立的模型中，对实际问题进展可行性分析，根据条件的改变完善对模型的应用和实用性检验 第一问中，我们通过两种太阳高度角的表示方法建立等式关系，根据控制变量法，分析出影子长度分别与经、纬度、杆长、时间、日期的关系然后，根据时差计算关系，当时间在9:00-15:00时，天安门广场的时间，并应用建立的模型 第二问中，首先根据影子坐标求出影子的长度，拟合时间与影子长度的函数，找出影子长度的最低的点，从而根据时间求出当地经度，由于误差的存在，我们将经度、杆长、纬度给定一定围，根据第一问公式进展搜索，从而确定可能的地点。

关键字：〔宋体小四号 〕 真太阳时 平太阳时 赤纬角 太阳高度角 熵值法二、 问题提出如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法 1. 建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用你们建立的模型画出2015年10月22日时间9:00-15:00之间天安门广场〔北纬39度54分26秒,东经116度23分29秒〕3米高的直杆的太阳影子长度的变化曲线 2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点将你们的模型应用于附件1的影子顶点坐标数据，给出假如干个可能的地点3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期将你们的模型分别应用于附件2和附件3的影子顶点坐标数据，给出假如干个可能的地点与日期4．附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为2米请建立确定视频拍摄地点的数学模型，并应用你们的模型给出假如干个可能的拍摄地点三、 问题分析 第一问：根据物体在太照射下将产生影子的自然现象,研究物体影子的形成原理, 通过分析太线照射物体的角度的日变化和年变化，引起物体影子的长度和朝向有规律地变化来建立数学模型。

利用Matlab软件绘出影子长短随时间变化的图像将问题中所给参数带入，解决问题由于太线照射物体的角度的日变化和年变化，引起物体影子的长度和朝向有规律地变化第二问：通过对附件所给的影子坐标的数据，求出影子的长度，然后通过第一问的相关公式，对影长和时间的关系进展拟合，得到一个二次方程，得出影长的最低值，从而可知正午时间，再算出经度四、 建模过程 第一问 1． 模型假设 〔1〕：假设单一光源〔太〕照射 〔2〕：直杆严格垂直于水平地面 〔3〕：被照射直杆的形状不会影响影子的长度 〔4〕：将整个天空视为一个天体圆 〔5〕：不考虑大气折射 〔6〕：问题中给出的数据可靠 2．定义符号说明 符号 含义 符号 含义 aa 偏磁角〔赤纬角〕 e 时差 e:e 太阳高度角 er 入射角 b b 经度 gg 纬度 t 时间 A 太阳方位角 N 自1月1日算起的第几天 L 影子长度 ww 时角 :H 杆长3． 模型建立： 以杆影在下产生影子端点移动的轨迹，代替太阳运行轨迹。

运用相对运动原理，将地球自转与绕太阳公转的运动简化为地球不动，太阳绕地球转动〔1〕 计算磁偏角〔赤纬角〕 全年之中，每一天太阳和地球的运转与天体圆赤道之间所形成的夹角，也就是所谓的磁偏角a都不同，会在+23.45与-23.45之间变化，其计算公式为： a=23.45sin[2p(28+4N) ]/365pa°+= (1) 〔2〕由时间计算当地时间： 按太阳运行位置，世界采取了时差制度并且遵循此制度，各国时间历法都以此制度 为根底按太阳运行位置，划分时区，每个时区相差15°〔每个时区相差1个小时〕当地时间s的计算公式： S=t-(120-R)/ 60*15 b-=-´ (2) 当所得值为负数时，加上24小时 〔3〕计算时角 因为地球自转一周约为24小时，所以，太阳每小时大约自东向西移动15°〔即360/24 ° 〕， 故时角w的计算公式为： ()w=15(12-s)° (3) w为正表示偏东，w为负表示偏西。

注意：计算中将其划为弧度制 （4） 计算太阳高度角 太阳高度角简称太阳高度(其实是角度)太阳高度是决定地球外表获得太阳热能数量的最重要的因素，它在数值上等于太阳在地球地平坐标系中的地平高度太阳高度角e的计算公式为： ()r=arcnsin(sincoscoscosarcw) (4)〔5〕利用太阳高度角、杆长与影长列出函数式如下列图，由立竿见影的测量方式，得出影长L公式为：L=H/tan P4，模型求解：由有几何学原理，tan pr在0

g：北纬39度54分26秒H 3米 注意：计算中将经、纬度划弧度制 将参数带入Matlab中，绘出影子长度变化曲线，如二所示：第二问 1. 模型假设： 假设附件1所给顶点坐标数据符合事实 2. 定义符号说明：A：方位角 ：高度角 aa：赤纬角 bb：物体所在地理纬度 g：当地时间 NN：从1月1日起距当地日期的天数 n：表示24小时制的时间数 t：太阳某位置的方位时间 ww：时角 L0：影子的长度 (x.y)：影子顶点坐标 3,模型建立： 因为竹竿相对太阳的位置所对应于地球上所在点的相对位置, 由该点的地理纬度、日期和时间3个因素来决定一般通过地平坐标系与赤道坐标系来同时表示太阳的位置, 也就是以太阳高度角q方位角A 与 赤纬角a、时角w来表示赤纬角是指地球赤道平面与太阳和地球中心的连线之间的夹角，高度角指太的入射方向和地平面之间的夹角, 可知高度角的围是000~90，方位角指经过球心O与太阳位置点在地平圈上投影点的直线与地平圈正南向OS所夹的角。

定义方位角坐标以正南向S点为起始00逆时针方向为负, 分180顺时针方向为正, 亦分0180;正北向N点为正负180±因此根据太阳位置的变化可以绘制出坐标网图, 在坐标网图中用同心圆来代表太阳高度圈, 用圆周上的刻度角来表示太阳的方位角(自南向西为正值, 自南向东为负值) 故通过竿影轨迹点的坐标可求出影长，然后通过拟合影长 相关计算公式如下： 影长公式根据韦达定理得： 五、 模型评价与改良 第一问中，我们所用的F检验没有T检验适宜，但是结果偏差不大 第二问中，topsis法灵活简便，操作方便并用熵值法确定权重，具有合理性，但是缺少相应的模型检验 第四问中，通过经验设定综合指标进展求解，简化了相应的数学模型，只是缺少对综合指标设立的检验，依据性不强 六、 参考文献 [1] Topsis综合评价法2012年9月8日 [2]数学建模成绩的评价与预测，2012年9月8日 [3]光亭 裘哲勇 《数学建模》 高等教育 2010年2月 [4]王宏洲 《数学建模优秀论文》 清华大学 2011年9月 / 。