太阳影子定位技术的研究毕业论文

1、 太阳影子定位技术的研究摘要影子在生活中非常常见，而太阳影子与物体本身、物体所在的地理位置、当地的日期和当天时间都有密切的关系，本文主要研究影长随物体本身的高度（丨）、物体所在的地理位置 （与有关）、当地的日期（与:.有关）和当天时间（与 t有关）的变化规律，从而根据影子变化的相关信息推测 出物体所处的位置。首先，我们查找到了包含物体本身的高度、当地时间、季节和地理位置这四个因素的太阳高度 角公式：sinh二sin sin ； cos co cost，然后采用控制变量法分别做出影长随单个因素的 变化图像，进而分析影长与各个因素之间的联系。当物体本身的高度、季节、地理位置都已知时， 则可以直接作

2、出影子长度随时间的变化关系图。由图可知：2015年10月22日天安门广场3米高的直杆的太阳影子在 9:00到15:00之间的变化大致呈对称分布，最小影长为：3.663米、最大影长为：6.616 米。sco柘cost其次，考虑到影子顶点的位置随太阳方位角会变化，所以需要考虑太阳方位角 （cosh ）， 对问题一的模型进行修正，消去杆长之后，利用附件1中每隔3分钟的影子长度变化差列出方程组，解得附件1的地点可能为：（北纬18度46分2秒冻经111度23分29秒），即海南附近。然后，我们根据附录二，绘制影长随时间的变化曲线，再对其进行拟合分析，得到当地经度为74度26分。由于纬度和赤纬不随时间变化，

3、所以我们采用最小二乘法进行拟合，最终解得 =39.34、=-0.056795、N=275，即在新疆地区（北纬 39度20分，东经74度26分），日期为10月2 日。同时对于附录三，采用同样的处理办法，得到=31.46 、- -0.3449、N=325，即在湖北地区（北纬31度27分，东经108度44分），日期为11月23日。最后，以视屏中旗杆所在地为圆心建立坐标系，每隔3分钟对视屏进行截屏，得到一组和附件1类似的杆影的顶点坐标随时间变化的数据，把模型二对附件一的处理方法运用到这组数据中，确定出视频的拍摄位置在陕西境内。当拍摄日期未知时，采用模型三处理这组数据，也得到视频的拍 摄地点在陕西境内。

4、关键词：太阳高度角太阳方位角控制变量最小二乘法拟合一问题重述视频数据分析在现代生活非常常见，而确定视频的拍摄地点和拍摄日期更是其重中之重，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。太阳影子的变化与物体本身的高度及所处的地理位置、时间和季节等因素有密切的关系，知道其中的几个因素便可以对物体所处地点进行定位。因而我们需要建立数学模型来确定这些因素之间的关系，对此进行以下几步工作：一、建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用所建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分2

5、6秒冻经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。二、根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。并将所建的模型应用于附件1的影子顶点坐标数据，得到若干个可能的地点。三、根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。并将所建的模型分别应用于附件2和附件3的影子顶点坐标数据，得到若干个可能的地点与日期。四、建立确定视频拍摄地点的数学模型，附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出了直杆的高度为2米。应用所建立的模型给出若干个可能的拍摄地点。二问题分析太阳影子的变化与物体本身的

6、高度、时间、季节和地理位置等因素都有重要关系，因此研究影子的变化规律有重要的意义。针对问题一，我们要查找包含物体本身的高度、当地时间、季节和地理位置的函数表达式，再对影长随单个因素的变化趋势进行分析。当物体本身的高度、季节、地理位置都是已知时，可以直 接作出影子长度随时间的变化关系图。针对问题二，在问题一的基础上，考虑到影子顶点的位置随太阳方位角会变化，所以需要对问题一的模型进行修正，由于h和杆高、影长联系在一起的，而附件1中的杆高未知，所以要先把它代入之前的太阳高度角公式消去h,得到矫正以后的太阳高度角公式，然会利用附件1中每隔3分钟的影子长度变化差，列出方程组，解得当地时间t,用它和北京时

7、间的时间差算出它与北京的经度差，进而算出当地的经度。同时，利用附件1中的影长、日期和当地的时间， 可以算出当地的纬度。对于问题三：我们根据附录，绘制影长随时间的变化曲线，再对其进行拟合分析，得到函数关系，根据函数关系，我们预测太阳影长最短时的北京时间、计算时差，又由于纬度和赤纬不随时间的变化，所以我们采用最小二乘法进行拟合，确定参数和。对于问题四：由于附件4给了一个视频，所以我们先每隔3分钟对视频截屏，以旗杆与地面的交点为圆心建立坐标系，得到影的顶点坐标，把视屏数据转化为类似于附件1的数字数据，再运用模型二对这组数据列方程组，求解出视频的拍摄地点。三模型假设假设太阳是固定不动的。假设地球离太阳

8、的距离无穷远，太阳光线是平行光。假设地球为球体。假设地球的运行轨道为圆。假设大气层均匀，太阳光线穿过大气层时没有折射。不考虑高山遮挡、海拔高度等因素的影响。四符号说明h :太阳高度角。:地球纬度。:太阳赤纬。t:时角。:太阳方位角。11 :影长。l :物体高度。T :当地时间。x :影子的横坐标。y :影子的纵坐标。五模型的建立与求解5.1建立模型一太阳影子的变化与物体的本身的高度及所处的地理位置、时间和季节等参数有密切的关系。根据地理学知识，我们知道，物体高度越高，太阳影子越长，物体所处纬度越远离太阳直射点的纬度， 太阳影子越长，查阅资料得：sinh 二 sin sin 心亠 cos cos

9、、cost （1）其中，h表示太阳高度角，表示地理纬度，表示太阳赤纬，t表示时角。时角t是以正午十二点为零度开始算，每小时为15度。即十四点和十点分别为 30度和-30度。太阳赤纬是地球赤道平面与太阳和地球中心的连线之间的夹角。因赤纬值日变化很小，一年内 任何一天的赤纬角 S可用下式计算:(2)(3)=0.006918 _0.399912cosb 0.070257sin b J0.006758cos(2 b) 0.000907sin(2 b) _0.002697cos(3b) 0.00148sin(3 b)b =2二(N -1)/365 式中N为日数，自每年1月1日开始计算。由于几何关系，太阳

10、高度角与影长有如下关系其中1为物咼，l1为影长，联立(1)、(4)(4)可得影长与各个参数的关系，即为最终所建数学模型l1sin sin、 cos : cos、cost-l2(5)5.1.1模型一的求解通过对模型的分析，我们知道太阳影子的变化与物体的本身的高度及所处的地理位置、时间和季节等参数有密切的关系，因此我们采用控制变量法来分别分析影子长度与各个因素之间的关系：(1)太阳影子和纬度的变化关系:我们保持物体的本身的高度及时间、季节等参数不变，单独分析太阳影子和纬度之间的变化关系。对此我们作了以下图像：图1太阳影子随纬度变化通过对图1分析可知，对于同一日期，同一时刻，纬度由直射点向北或向南变

11、化时，太阳影子 的长度总是越来越长。太阳影子和物高的变化关系我们保持物体所处的纬度及时间、季节等参数不变，单独分析太阳影子和物体高度的变化关系。对此我们作了以下图像：影长图2太阳影长随物高的变化通过图2分析可知，同一日期，同一时刻，同一纬度不同物体的太阳影子与物体的高度成正比。太阳影子长度和季节的变化关系我们保持物体本身高度及所处的纬度和时间不变，单独分析太阳影子和季节的变化关系。太阳影子的变化与物体所处纬度的不同呈不同的规律：a当纬度大于23 26时，保持物体本身高度及所处的纬度和时间不变，太阳影子长度随日期的 变化如下图所示：I 影子3530子15.10 -5050100150200250

12、300350400日期图3南北回归线外影长随日期的变化由图知，物体所处纬度大于北回归线度数时，太阳影子长度自1月1日起随日期的增加先减小后增加。b当纬度介于0和23 26之间时，保持物体本身高度及所处的纬度和时间不变，太阳影子长度随日期的变化如下图所示：图4南北回归线内影长随日期的变化由图可知，物体所处纬度小于北回归线度数时，由于太阳直射点的移动，太阳影子长度自1日起随日期的增加先减小后增大再减小最后又增大。对于问题一，应用我们所建立的模型，对2015年10月22日北京9:00-15:00之间天安门广场（北纬39度54分26秒，东经116度23分29秒）3米高的直杆的太阳影子长度的变化进行了如

13、下求解：联立、（3）式可知当N=295，、： - -0.181把、：和:（39.97度）代入得到太阳影长y随地方时t的变化关系:影长图5影长在9点到12点的变化规律 (8) 由图5可知，2015年10月22日天安门广场 3米高的直杆的太阳影子在9:00到15:00之间的变化规律为：先减小后增大，呈对称分布，太阳影子长度分别在 9 : 00和15 : 00点时出现最大值， 最大影长约为：6.616米，在12 : 00左右时出现最小值，最小影长为：3.663米。5.2建立模型二经过对模型一的分析，我们知道，太阳影子的长短与太阳高度角有关，而太阳影子的朝向却与太阳方位角有关，因此，我们通过查阅资料，

14、得到太阳方位角有如下公式.+ COS6 cost 2 sincosh其中为太阳方位角阴影长度的南北分量图6太阳方位角示意图又根据图6的几何关系知:I1sin其中，x为影子的横坐标。通过联立（1）（ 6）（ 7）可得模型二：l _ x（s in sin 6 + cos cos 6 cost）1COS cost5.2.1模型二的求解附件一提供了影子在一段时间内的坐标，利用上述模型，根据相同时间间隔内影长的增长量得到下列方程组:t2 -ti =0.1047cos?目sin211 屮 sin2寿cos2、sin212sin2 2cos2sin212sin2 22盂2丄cos sin t32 I-sin

15、 3cost - cost2 cost2 -cost3t _t2= 0.1047其中sin由（7）式给出，J.由（2）式给出。运用mathematica8.0 解方程组（9）得出I）= 0.62 （单位：弧度），则由公式:(10)T =12 (t 180 -门)15得当地时间T=14.368时=14时22分。它与北京时间的时间差为 20分钟，由时差可得该地的经度与北京的经度相差 5度，所以该地的经度为 111度23分29秒。联立（1）、（ 6）式得：cos2、sin2112I：-sin 1=sin sin 心亠 coscos、cost1(11)将问题二中已知的 二、1值及求解出的t1值代入（1

16、1、式中，解得：打=20度46分2秒，综 上所述，附件1可能为：（北纬20度46分2秒，东经111度23分29秒），即海南附近。5.3模型三的建立首先，我们根据附录二，绘制影长随时间的变化曲线，再对其进行拟合分析得到函数关系，预测到当太阳影长最短时的北京时间，得到当地时间与北京时间的差值，再根据残差进行修正，得到当地的经度。而纬度和赤纬不随时间的变化，所以我们采用最小二乘法对公式（8、进行拟合，确定参数和-：,并运用mathematical.。进行求解，确定：、：、N的值。5.3.1模型三的求解对附件2进行曲线拟合得到图 7:1.3图7附件2的曲线图由图7得到函数关系：h =23.3-2.9&

17、 0.098t2，预测当太阳影子长度最短时北京时间是15:12，即当地时间与北京时间的差值为3.2小时，经过修正后，时间差值为2.73小时，即当地经度与北京经度的差值为 40.95 ，得到当地的经度为 74度26分。由最小二乘法对公式（8）进行拟合，最终确定 即=39.34、： - -0.056795、N=275。即在新疆地区（北纬 39度20分，东经74 度26分），日期为10月2日。同上，我们对于附录三采同样的处理方法，拟合曲线得到图&图8附件3的曲线图由图8得到函数关系 九=51.56-7.55t0.296t2,预测当太阳影子长度最短时北京时间是12:44，即当地时间与北京时间的差值为0

18、.73个小时，经过修正后，时间差值为 0.51个小时，即当地经度与北京经度的差值为7.65 ，得到当地的经度为108度44分。由最小二乘法对公式（8）进行拟合，最终确定=31.46，、： - -0.3449 , N=325，即在湖北地区（北纬 31度27分，东经108度44分），日期为11月23日。5.4模型四的建立与求解：附件4只给了一个视频，所以我们先每隔3分钟对视频截屏，以旗杆与地面的交点为圆心建立坐标系，得到杆影的顶点坐标，把视屏数据转化为类似于附件1的数字数据：表1视频中影长随时间的变化关系时间：7月13号北京时间X坐标Y坐标影长2 v sin9:030028

19、0.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.0523

20、38要确定视频拍摄地的地理位置，就必须知道杆的高度、太阳影子的变化规律及当地的时间。而题中杆的高度、当地时间都是未知的，因此我们首先要消去杆高（丨）这个未知量，进而求得当地时间（t）,对此我们进行如下处理:联立（1 ）、（ 5）、（ 6）并利用相同时间间隔内影长的增长量得到下列方程组:t2 7 =0.0131cos2、sin211cos2sin212sin2 1sin2 2cos2、sin212cos2sin 2t3cost -cost2 cost2 - cost3(12)sin2 2sin2 3t3 -t2 =0.0131N=194联立（2）、（ 3）、（ 13）式解得：6 = 0.383(

21、13)(14)把、sin2 1、 sin2 2、Sih 3代入（12）式解得：=-0.883386 （单位：弧度）则由公式（9）得：当地时间T=8 : 38,它与北京时间的时间差为25分钟，由时差可得该地的经度与北京的经度相差 6.25度，所以该地的经度为110度8分29秒。把 t1=-0.5 （单位：弧度）、=0.383、sin2 讦0.999975 代入（10 ）式，解得： =36度40分9秒综上所述，视频拍摄地大概位于（北纬 36度40分9秒，东经110度8分29秒），在陕西省境 内。当拍摄日期未知时，经分析发现，和问题三相似，所以我们采用模型三进行求解：首先，我们根据表1,绘制影长随时

22、间的变化曲线，再对其进行拟合分析得到函数关系，预测到当太阳影长最短时的北京时间，得到当地时间与北京时间的差值，再根据残差进行修正，得到当地的经度。而纬度和赤纬不随时间的变化，所以我们采用最小二乘法对公式（8 ）进行拟合，确定参数和、：，并运用 mathematical.。进行求解，确定 、：、N的值。对表2中的数据进行曲线拟合得到图9 :0.900.880.860.840.820.780.760.740.729.009.019.029.039.049.059.06时间长080Equationy - Intercept + B1*xA1 + B2*xA2WeightNo WeightingRes

23、idual Sum ofSquares8.66643E-5Adj. R-Square0.9962ValueStandard ErrorBIntercept1310.76869366.49338BB1-286.9321781.17257BB215.710174.4946图9视频中影长随时间的变化关系图2由图9得到函数关系：|1 = 1310.76869-286.9321A 15.71017t，预测当太阳影子长度最 短时北京时间是12 : 25，即当地时间与北京时间的差值为 0.47个小时，经过修正后，时间差值为 0.41个小时，即当地经度与北京经度的差值为 6.25 ，得到当地的经度为 110度

24、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 Comma nd1_Click()N = Tex

25、tl.Text b = (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) phi = 10 * 3.1415926 / 180 y = Sqr(x / (Si n( phi) * Si n(delta) + Cos(phi) * Cos(delta) * Cos(t) A 2 - x A 2) Pr