太阳影子定位技术研究.doc

太阳影子定位技术的研究摘要影子在生活中非常常见，而太阳影子与物体本身、物体所在的地理位置、当地的日期和当天时间都有密切的关系，本文主要研究影长随物体本身的高度（）、物体所在的地理位置（与有关）、当地的日期（与有关）和当天时间（与t有关）的变化规律，从而根据影子变化的相关信息推测出物体所处的位置。首先，我们查找到了包含物体本身的高度、当地时间、季节和地理位置这四个因素的太阳高度角公式： ，然后采用控制变量法分别做出影长随单个因素的变化图像，进而分析影长与各个因素之间的联系。当物体本身的高度、季节、地理位置都已知时，则可以直接作出影子长度随时间的变化关系图。由图可知：2015年10月22日天安门广场3米高的直杆的太阳影子在9:00到15:00之间的变化大致呈对称分布，最小影长为：3.663米、最大影长为：6.616米。其次，考虑到影子顶点的位置随太阳方位角会变化，所以需要考虑太阳方位角（），对问题一的模型进行修正，消去杆长之后，利用附件1中每隔3分钟的影子长度变化差列出方程组，解得附件1的地点可能为：（北纬18度46分2秒,东经111度23分29秒），即海南附近。然后，我们根据附录二，绘制影长随时间的变化曲线，再对其进行拟合分析，得到当地经度为74度26分。由于纬度和赤纬不随时间变化，所以我们采用最小二乘法进行拟合，最终解得、N=275，即在新疆地区（北纬39度20分，东经74度26分），日期为10月2日。同时对于附录三，采用同样的处理办法，得到、N=325，即在湖北地区（北纬31度27分，东经108度44分），日期为11月23日。最后，以视屏中旗杆所在地为圆心建立坐标系，每隔3分钟对视屏进行截屏，得到一组和附件1类似的杆影的顶点坐标随时间变化的数据，把模型二对附件一的处理方法运用到这组数据中，确定出视频的拍摄位置在陕西境内。当拍摄日期未知时，采用模型三处理这组数据，也得到视频的拍摄地点在陕西境内。关键词：太阳高度角 太阳方位角 控制变量 最小二乘法拟合 19一问题重述视频数据分析在现代生活非常常见，而确定视频的拍摄地点和拍摄日期更是其重中之重，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。太阳影子的变化与物体本身的高度及所处的地理位置、时间和季节等因素有密切的关系，知道其中的几个因素便可以对物体所处地点进行定位。因而我们需要建立数学模型来确定这些因素之间的关系，对此进行以下几步工作：一、建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用所建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。二、根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。并将所建的模型应用于附件1的影子顶点坐标数据，得到若干个可能的地点。三、根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。并将所建的模型分别应用于附件2和附件3的影子顶点坐标数据，得到若干个可能的地点与日期。四、建立确定视频拍摄地点的数学模型，附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出了直杆的高度为2米。应用所建立的模型给出若干个可能的拍摄地点。二 问题分析太阳影子的变化与物体本身的高度、时间、季节和地理位置等因素都有重要关系，因此研究影子的变化规律有重要的意义。针对问题一，我们要查找包含物体本身的高度、当地时间、季节和地理位置的函数表达式，再对影长随单个因素的变化趋势进行分析。当物体本身的高度、季节、地理位置都是已知时，可以直接作出影子长度随时间的变化关系图。针对问题二，在问题一的基础上，考虑到影子顶点的位置随太阳方位角会变化，所以需要对问题一的模型进行修正，由于h和杆高、影长联系在一起的，而附件1中的杆高未知，所以要先把它代入之前的太阳高度角公式消去h, 得到矫正以后的太阳高度角公式，然会利用附件1中每隔3分钟的影子长度变化差，列出方程组，解得当地时间t, 用它和北京时间的时间差算出它与北京的经度差，进而算出当地的经度。同时，利用附件1中的影长、日期和当地的时间，可以算出当地的纬度。对于问题三：我们根据附录，绘制影长随时间的变化曲线，再对其进行拟合分析，得到函数关系，根据函数关系，我们预测太阳影长最短时的北京时间、计算时差，又由于纬度和赤纬不随时间的变化，所以我们采用最小二乘法进行拟合，确定参数和。对于问题四:由于附件4给了一个视频，所以我们先每隔3分钟对视频截屏，以旗杆与地面的交点为圆心建立坐标系，得到影的顶点坐标，把视屏数据转化为类似于附件1的数字数据，再运用模型二对这组数据列方程组，求解出视频的拍摄地点。三模型假设1.假设太阳是固定不动的。2.假设地球离太阳的距离无穷远，太阳光线是平行光。3.假设地球为球体。4.假设地球的运行轨道为圆。5.假设大气层均匀，太阳光线穿过大气层时没有折射。6.不考虑高山遮挡、海拔高度等因素的影响。四符号说明1.：太阳高度角。2.：地球纬度。3.：太阳赤纬。4.：时角。5.：太阳方位角。6.：影长。7.：物体高度。8.T：当地时间。9.x：影子的横坐标。10.y：影子的纵坐标。五 模型的建立与求解5.1建立模型一太阳影子的变化与物体的本身的高度及所处的地理位置、时间和季节等参数有密切的关系。 根据地理学知识，我们知道，物体高度越高，太阳影子越长，物体所处纬度越远离太阳直射点的纬度，太阳影子越长，查阅资料得： （1）其中，表示太阳高度角，表示地理纬度，表示太阳赤纬，表示时角。时角t是以正午十二点为零度开始算，每小时为15度。即十四点和十点分别为30度和-30度。太阳赤纬是地球赤道平面与太阳和地球中心的连线之间的夹角。因赤纬值日变化很小，一年内任何一天的赤纬角可用下式计算： （2） （3）式中N为日数，自每年1月1日开始计算。由于几何关系，太阳高度角与影长有如下关系（4）其中为物高，为影长，联立（1）、（4）可得影长与各个参数的关系，即为最终所建数学模型。（5）5.1.1模型一的求解通过对模型的分析，我们知道太阳影子的变化与物体的本身的高度及所处的地理位置、时间和季节等参数有密切的关系，因此我们采用控制变量法来分别分析影子长度与各个因素之间的关系：(1) 太阳影子和纬度的变化关系： 我们保持物体的本身的高度及时间、季节等参数不变，单独分析太阳影子和纬度之间的变化关系。对此我们作了以下图像：图1太阳影子随纬度变化通过对图1分析可知，对于同一日期，同一时刻，纬度由直射点向北或向南变化时，太阳影子的长度总是越来越长。(2)太阳影子和物高的变化关系我们保持物体所处的纬度及时间、季节等参数不变，单独分析太阳影子和物体高度的变化关系。对此我们作了以下图像： 图2 太阳影长随物高的变化通过图2分析可知，同一日期，同一时刻，同一纬度不同物体的太阳影子与物体的高度成正比。(3)太阳影子长度和季节的变化关系我们保持物体本身高度及所处的纬度和时间不变，单独分析太阳影子和季节的变化关系。太阳影子的变化与物体所处纬度的不同呈不同的规律:a.当纬度大于时，保持物体本身高度及所处的纬度和时间不变，太阳影子长度随日期的变化如下图所示：图3 南北回归线外影长随日期的变化由图知，物体所处纬度大于北回归线度数时，太阳影子长度自1月1日起随日期的增加先减小后增加。b.当纬度介于和之间时，保持物体本身高度及所处的纬度和时间不变，太阳影子长度随日期的变化如下图所示：图4 南北回归线内影长随日期的变化由图可知，物体所处纬度小于北回归线度数时，由于太阳直射点的移动，太阳影子长度自1月1日起随日期的增加先减小后增大再减小最后又增大。对于问题一，应用我们所建立的模型，对2015年10月22日北京9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化进行了如下求解：联立(2)、（3）式可知当N=295，把和 (39.97度)代入(5)得到太阳影长y随地方时t的变化关系：图5 影长在9点到12点的变化规律由图5可知，2015年10月22日天安门广场3米高的直杆的太阳影子在9:00到15:00之间的变化规律为：先减小后增大，呈对称分布，太阳影子长度分别在9：00和15：00点时出现最大值，最大影长约为：6.616米，在12：00左右时出现最小值，最小影长为：3.663米。5.2 建立模型二 经过对模型一的分析，我们知道，太阳影子的长短与太阳高度角有关，而太阳影子的朝向却与太阳方位角有关，因此，我们通过查阅资料，得到太阳方位角有如下公式 (6)其中为太阳方位角阴影长度的南北分量遮挡物高度太阳方位角阴影长度太阳光线h阴影长度的东西分量太阳高度角图6 太阳方位角示意图又根据图6的几何关系知： (7)其中，x为影子的横坐标。通过联立（1）（6）（7）可得模型二： （8）5.2.1 模型二的求解附件一提供了影子在一段时间内的坐标，利用上述模型，根据相同时间间隔内影长的增长量得到下列方程组： （9）其中由（7）式给出，由（2）式给出。运用mathematica8.0解方程组（9）得出（单位：弧度），则由公式： （10）得当地时间T=14.368时=14时22分。它与北京时间的时间差为20分钟，由时差可得该地的经度与北京的经度相差5度，所以该地的经度为111度23分29秒。联立（1）、（6）式得： （11）将问题二中已知的、值及求解出的值代入（11）式中，解得: =20度46分2秒，综上所述，附件1可能为：（北纬20度46分2秒,东经111度23分29秒），即海南附近。5.3模型三的建立首先，我们根据附录二，绘制影长随时间的变化曲线，再对其进行拟合分析得到函数关系，预测到当太阳影长最短时的北京时间，得到当地时间与北京时间的差值，再根据残差进行修正，得到当地的经度。而纬度和赤纬不随时间的变化，所以我们采用最小二乘法对公式（8）进行拟合，确定参数和，并运用mathematical8.0进行求解，确定的值。5.3.1模型三的求解对附件2进行曲线拟合得到图7：图7 附件2的曲线图由图7得到函数关系：，预测当太阳影子长度最短时北京时间是15:12，即当地时间与北京时间的差值为3.2小时，经过修正后，时间差值为2.73小时，即当地经度与北京经度的差值为40.95，得到当地的经度为74度26分。由最小二乘法对公式（8）进行拟合，最终确定、N=275。即在新疆地区（北纬39度20分，东经74度26分），日期为10月2日。同上，我们对于附录三采同样的处理方法，拟合曲线得到图8：图8 附件3的曲线图由图8得到函数关系: ,预测当太阳影子长度最短时北京时间是12:44，即当地时间与北京时间的差值为0.73个小时，经过修正后，时间差值为0.51个小时，即当地经度与北京经度的差值为7.65，得到当地的经度为108度44分。由最小二乘法对公式（8）进行拟合，最终确定，N=325，即在湖北地区（北纬31度27分，东经108度44分），日期为11月23日。5.4模型四的建立与求解：附件4只给了一个视频，所以我们先每隔3分钟对视频截屏，以旗杆与地面的交点为圆心建立坐标系，得到杆影的顶点坐标，把视屏数据转化为类似于附件1的数字数据：表1 视频中影长随时间的变化关系 时间：7月13号 北京时间X坐标Y坐标影长9:0320.160.188.100280.0523639:0619.760.186.352310.0523639:0919.420.1584.867930.0523619:1218.420.1582.246010.0523619：1518.410.1680.454740.0523619：1818.20.2179.539290.0523589：2117.830.2177.92250.0523579：2417.460.3276.313010.0523479：2717.040.3274.477930.0523469：3016.770.3273.298240.0523459：3316.560.3772.385260.052338要确定视频拍摄地的地理位置，就必须知道杆的高度、太阳影子的变化规律及当地的时间。而题中杆的高度、当地时间都是未知的，因此我们首先要消去杆高（）这个未知量，进而求得当地时间（t），对此我们进行如下处理：联立（1）、（5）、（6）并利用相同时间间隔内影长的增长量得到下列方程组： （12）N=194 （13） 联立（2）、（3）、（13）式解得 ： 0.383 （14）把、 、 、 代入（12）式解得：t=-0.883386（单位：弧度）则由公式（9）得:当地时间T=8：38 ,它与北京时间的时间差为25分钟，由时差可得该地的经度与北京的经度相差6.25度，所以该地的经度为110度8分29秒。把t=-0.5（单位：弧度）、=0.383、=0.999975代入（10）式，解得: =36度40分9秒综上所述，视频拍摄地大概位于（北纬36度40分9秒,东经110度8分29秒），在陕西省境内。当拍摄日期未知时，经分析发现，和问题三相似，所以我们采用模型三进行求解：首先，我们根据表1，绘制影长随时间的变化曲线，再对其进行拟合分析得到函数关系，预测到当太阳影长最短时的北京时间，得到当地时间与北京时间的差值，再根据残差进行修正，得到当地的经度。而纬度和赤纬不随时间的变化，所以我们采用最小二乘法对公式（8）进行拟合，确定参数和，并运用mathematical8.0进行求解，确定的值。对表2中的数据进行曲线拟合得到图9：图9 视频中影长随时间的变化关系图由图9得到函数关系：，预测当太阳影子长度最短时北京时间是12：25，即当地时间与北京时间的差值为0.47个小时，经过修正后，时间差值为0.41个小时，即当地经度与北京经度的差值为6.25，得到当地的经度为110度15分。由最小二乘法对公式（8）进行拟合，最终确定36度40分，0.385，N=193，即在西安地区（北纬36度40分，东经110度15分），日期为7月12日。利用模型三所求解的日期与视频所给出的日期在误差允许的范围内，说明我们的模型合理。参考文献1 百度百科，http/.www.baike.baidu,2015年9月11日星期一2 何银涛，张梅，黄华，垂直单轴跟踪光伏支架方针间距计算界面设计，太阳能，2015（7）3 武琳，基于太阳阴影轨迹的经纬度估计技术研究，天津大学，2010。 附录研究参数规律Private Sub Command1_Click()N = Text1.Textb = (2 * 3.1415926 * (N - 1) / 365)delta = (0.006918 - 0.399912 * Cos(b) + 0.070257 * Sin(b) - 0.006758 * Cos(2 * b) + 0.000907 * Sin(2 * b) - 0.002697 * Cos(3 * b) + 0.00148 * Sin(3 * b)x = 3t = 0phi = 10 * 3.1415926 / 180y = Sqr(x / (Sin(phi) * Sin(delta) + Cos(phi) * Cos(delta) * Cos(t) 2 - x 2)Print yEnd Sub求解方位角F=(Cost1-Cost1+3)/(Cost1+3-Cost1+6)*(Sqrt1-1.217*(Sint1+3)2-Sqrt1-1.208*(Sin6+t1)2)+Sqrt1-1.217(Sint1+3)2-Sqrt1-1.226*(Sint1)2;PlotF,t1,0,p/2ExportE:/a.txt,Tablet1,F,t1,0,p/2,0.001,Table求解时间与影长Private Sub Command1_Click()Dim Pi As DoublePi = 4 * Atn(1)phi = 39.906 * Pi / 180delta = -0.1880374t = (Text3.Text) * Pi / 180y = Sqr(3 / (Sin(phi) * Sin(delta) + Cos(phi) * Cos(delta) * Cos(t) 2 - 9)Print yEnd Sub问题二代码SolveSqrt 1 - (Cos0.1841*Sin1.01)/(1.0365/ Sqrt1.03652 + 0.49732)2 = SinPhi*Sin0.1841 + CosPhi*Cos0.1841*Cos1.01问题三代码fp = TablePrimex, x, 201.109142, 1.104962, 1.101012, 1.097295, 1.093757, 1.090422, 1.087214, 1.084188, 1.08132, 1.078567, 1.075945, 1.073407, 1.071014, 1.068697, 1.066489, 1.064357, 1.062298, 1.060331, 1.058438, 1.056612, 1.054833FindFitfp, (Tand*Sinf)/(x - Cosf), d, f, x目 录第一章 总论11.1项目名称与承办单位11.2研究工作的依据、内容及范围11.3编制原则31.4项目概况31.5技术经济指标51.6结论6第二章 项目背景及建设必要性82.1项目背景82.2建设的必要性9第三章 建设条件113.1项目区概况113.2建设地点选择错误！未定义书签。3.3项目建设条件优劣势分析错误！未定义书签。第四章 市场分析与销售方案134.1市场分析134.2营销策略、方案、模式14第五章 建设方案155.1建设规模和产品方案155.2建设规划和布局155.3运输185.4建设标准185.5公用工程205.6工艺技术方案215.7设备方案215.8节能减排措施24第六章 环境影响评价256.1环境影响256.2环境保护与治理措施266.3评价与审批28第七章 项目组织与管理297.1组织机构与职能划分297.2劳动定员297.3经营管理措施307.4技术培训30第八章 劳动、安全、卫生与消防318.1编制依据及采用的标准318.2安全卫生防护原则318.3自然灾害危害因素分析及防范措施328.4生产过程中产生的危害因素分析及防范措施328.5消防编制依据及采用的标准348.6消防设计原则358.7火灾隐患分析358.8总平面消防设计358.9消防给水设计368.10建筑防火368.11火灾检测报警系统378.12预期效果37第九章 项目实施进度389.1实施进度计划389.2项目实施建议38第十章 项目招投标方案4010.1招标原则4010.2项目招标范围4010.3投标、开标、评标和中标程序4010.4评标委员会的人员组成和资格要求42第十一章 投资估算和资金筹措4311.1投资估算4311.2资金筹措及使用计划45第十二章 财务评价4712.1费用与效益估算4712.2财务分析4812.3不确定性分析4912.5财务评价结论50第十三章 建设合理性分析5113.1产业政策符合性分析5113.2清洁生产符合性分析5113.3规划符合性分析5113.4项目建设环保政策符合性分析5113.5环境承载性分析5113.6结论52第十四章 结论与建议53