根据影子判断地理位置Word文件下载.docx

1、通过提取视频特定帧，测量出杆的影子随时间变化的实际长度数据，与问题三类似，结合模型用matlab曲线拟合求解即可，或者取多组数据用lingo软件求解方程组。模型假设1.假设地球公转轨迹近似为圆。2.忽略太阳光线进入大气层时的折射误差。3.假设地面是水平的且直杆垂直地面。4.忽略太阳直射点纬度一天内的变化。5.假设所给数据准确可靠。 定义和符号说明H：杆长 L：杆的影长n：为从1月1号开始的天数N：为当地经度：为当地纬度t：时间（小时）时角区间均值解出现在某区间的频率模型建立与求解 问题一：在水平地面上，垂直竖杆与在地面上的影长是一个直角三角形的两条直角边，而太阳高度角是杆长H的对角，运用正切公

2、式可得：这样就把影长问题转化成了高度角问题。 一、构造高度角与各参数的关系图: 二、高度角与各参数关系模型构建 影响高度角的参数有太阳时角、赤纬角、杆的位置（经度、纬度）、时间。 时角，地球自转每小时15度，北京时间是以东经120度的为准的，故考虑时差问题，可得:，赤纬角是地球赤道平面与太阳和地球中心的连线之间的夹角，也就是太阳直射点的纬度。是因为地球公转形成的，因为地轴方向不变，所以赤纬角随地球在运行轨道上的不同点具有不同的数值。赤纬角以年为周期，在南北纬23.45度范围内移动，考虑其一年内运动规律，得赤纬角： 结合当地纬度，得太阳高度角与各参数间的关系模型（推算过程请见附件）：影子长度与各

3、参数关系模型： 联合公式可得：3、根据模型求解： 将问题一中的地点位置，日期，杆长代入式中，用matlab软件画出了直杆影子从9：00到15：00的长度变化曲线。（代码请见附件）。 影长曲线与时间的关系图如下图所示：问题二：根据影子顶点数据来求地点，由于直杆高度无法确定，所以很难根据影子长度来求太阳高度角问题，于是我们利用顶点坐标的轨迹来解决。利用matlab将题中给出杆的影子顶点运动轨迹拟合了出来，如图所示，其中坐标原点为杆的位置。可知影子发生了角度变化，以此来确定影子顶点与太阳的运动关系。1、立体问题转化为平面问题 我们通过将太阳与地面直杆的立体方位关系转化为平面关系来建立模型，把地球的自

4、转看成太阳在直射点轨迹上方绕地球转，并且短时间内可近似看作太阳是作由东向西的直线运动。并且存在两个关于直射点线（图中太阳运动轨迹）对称的位置，它们的影子变化规律也是对称的。以下模型建立是视y轴指向北，x轴指向东（若y轴以南为正方向，x轴以东为正方向，模型建立同理）。 如上图所示，此图可看作从太阳方向往地平面方向看去的俯视图，当太阳从A点运动到B点时，杆影顶点从a移动到b，O点为杆的位置。由于太阳、杆顶、影子顶点永远在一条直线上，据对顶角相等原理，可知太阳在时间段内移动过的经度差等于影子前后转动的角度差。从而可以根据角度关系来确定杆位置。2、模型建立 设北京时间t1t2时刻太阳从A到B，可知太阳

5、转过的经度为15（t2-t1）， 由于北京时间12点时太阳所在的经度为东经120度，所以t1时刻太阳所在A点经度为120-15（t1-12）。 太阳运动轨迹为直射点轨迹，根据问题一中的公式:可知太阳运动路线的纬度。 影子从a点（x1,y1）到b点（x2,y2） ，可知。 假设经度差为N，纬度差为 由三角关系解方程组： 又因杆位置纬度+, 经度N=15（t2-t1）+ 解得杆位置与参数间的关系模型：故若一直影子顶点坐标随时间变化的数据和日期，则可通过模型求出杆的位置。3、模型求解 我们把题中给出的影子顶点坐标数据和日期代入模型进行求解，每相邻两组坐标数据求一次解，利用matlab软件求得各组解的

6、坐标分布图，如图所示： 由于各解的分布为离散型，所以我们对解进行权重处理。我们把各解的坐标看为离散型随机变量，利用刻画离散型随机变量概率分布的方法来处理结果。把解出现在划分的区间的个数作为该区间的频数，频率作为所占的权重（在某区间出现的频率越大，则区间所占的权重越大），为该区间的均值，最后结果为 对经度处理： 区间 频数 频率 （110,111 1 0.05 （111,112 （112,113 3 0.15 （113,114 （114,115 （115,116 9 0.45最后计算得经度=114.15同理可求得纬度=27.27，此处不赘述。经查询，经纬度为114.15和27.27的地点位于江西

7、吉安市安福县。 若影子顶点坐标y轴是以南为正方向，x轴以东为正方向。则该杆位置与上述求得的位置关于直射点轨迹对称，计算得杆位置为 经度N=114.15 纬度=6.25。位于印度尼西亚。 四、模型检验 将问题中的问题三：问题三中已知的参数只有影子定点坐标和时间，所以无法利用问题二中的模型来求杆的具体位置。我们联合问题一和问题二建立的模型（以下分别称为模型一、模型二），然后用matlab对坐标数据进行曲线拟合，得出影长与时间的关系式，从而确定各参数值。1、模型结合 模型一： 模型二：因为已经知道影子顶点坐标数据，故利用模型二先进行预处理，可以得到杆位置经度N，由于日期未知，故模型二无法确定纬度。把

8、求得的经度代入模型一，可整理成函数关系式，结合以下公式 输入影长与时间关系的数据利用matlab，据最小二乘法原理来进行曲线拟合，并且关系式为给定的模型一，根据拟合结果确定参数、 则杆子位置的经纬度和日期就可以求得。2、模型求解 （一）附件二求解（1）、利用模型二求经度 求题方法于问题二中已说明，此处不详解。将题给的附件二中坐标数据代入模型二求解，整理得经度概率分布表：区间频数（80,853（85,906 0.3（90,954 0.2（95,1005 0.25（100,1052 0.1计算得经度N=91.75度。（2）利用模型一拟合曲线 首先，把坐标数据处理成影长随时间的关系，如下所示，为题中

9、给出的附件2影长与时间关系数据表格：时间/h12.6813.6814.6815.6816.6817.6818.6819.6820.6821.6822.68影长/m23.6824.68 把求得的经度N=91.75影长与时间的关系式L=，用matlab软件拟合曲线，进行参数估计，拟合结果如下图， 拟合结果得出了参数纬度=5.326度、天数n=357，从而就确定了杆位置为经度N=91.75,纬度=5.326。日期为12月22日。 与问题二同理，赤纬角=21.8，存在关于太阳直射点轨迹的对称点=38.3，N=91.75。（二）附件三求解 附件三数据求解过程与附件二同理，此处不赘述。（详情见附录）求解结

10、果为 杆位置经度N=115.8，纬度=5.6 日期为12月20日。问题四 首先，根据视频中可得在40分钟里直杆在太阳下影长变化规律。每隔两分钟截一张图可得20张图片。利用photoshop软件中的度量工具可得图片中直杆的高度与影子的长度可得一组数据。由于直杆的高度为2米，可得实际影长的长度数据如下：8:5557599:010305070911132.402.382.352.322.292.272.242.202.182.15151719212325272931332.122.102.072.052.011.991.971.941.911.89（1） 日期已知，为7月13日，将天数n代入，用ma

11、tlab软件对影长与时间的关系式进行曲线拟合，确定关系式 得曲线图：得参数经度N=-18.23度 纬度=21.86度 （2） 日期不确定，则有三个不确定的参数，三个未知量可用三个方程组求得解。用lingo软件输入三组数据，即可确定三个参数。 误差分析：1、不同时刻大气对光线的折射率不同；2、地球表面不是完全水平面；3、量取视频中影子长度时由于影子是斜的和视角问题，所量值非其实际值；4、数据计算取小数保留；5、曲线拟合存在误差；模型评价 模型优点：1、简洁性：模型原理简易，不需要太复杂的数学知识，易为大众所了解掌握。2、确定性：模型基本都是由公式组成，有较强的严谨性，计算出的都是确定的解。3、可

12、靠性，模型建立在坚实可靠的数学基础上。 模型局限性：1.计算公式过于繁琐，必要时得依赖计算机求解。2.模型只是建立的假设条件只在短时间内成立。3.应用曲线拟合只是大体上估计，存在较大偏差，只能拟合出大概解。 模型改进与推广： 模型建立条件是没有考虑大气层的折射问题，这可能在一定程度上影响模型准确度。所以应该在这方面改进，使其更接近于现实。参考文献：【1】方荣生，太阳能应用技术，北京：中国农业机械出版社，1985。【2】王松桂，概率论与数理统计，北京：科学出版社，2011。【3】杨长青. 正午太阳高度角模型”的设计、制作与运用. 内蒙古师范大学地理科学学院 2011-08-01 。【4】 王国安

13、. 太阳高度角和日出日落时刻太阳方位角一年变化范围的计算. 河南省气候中心 2007-09-15。【5】刘卫国，MATLAB程序设计教程，北京：中国水利水电出版社，2005。附录问题一相关代码：b=9:0.05:15;c=15*（b-12）-（120-116.39139）;a=23.45*sin（2*3.1415926*（284+295）/365）;i=39.907222;t=sind（i）*sind（a）+cosd（i）*cosd（a）*cosd（c）;x=asind（t）;d=3./tand（x）;plot（b,d）问题二相关代码： x=1.0365 1.0699 1.1038 1.138

14、3 1.1732 1.2087 1.2448 1.2815 1.3189 1.3568 1.3955 1.4349 1.4751 1.516 1.5577 1.6003 1.6438 1.6882 1.7337 1.7801 1.8277;y=0.4973 0.5029 0.5085 0.5142 0.5198 0.5255 0.5311 0.5368 0.5426 0.5483 0.5541 0.5598 0.5657 0.5715 0.5774 0.5833 0.5892 0.5952 0.6013 0.6074 0.6135t=14.7 14.75 14.8 14.85 14.9 14.

15、95 15 15.05 15.1 15.15 15.2 15.25 15.3 15.35 15.4 15.45 15.5 15.55 15.6 15.65 15.7for i=1:20e（i）=120-（t（i）-12）*15+15*（t（2）-t（1）*x（i）*y（i+1）/（y（i）*x（i+1）-x（i）*y（i+1）;n（i）=10.5109+15*（t（2）-t（1）*y（i）*y（i+1）/（y（i）*x（i+1）-x（i）*y（i+1）;endscatter（e,n,r*）figurescatter（n,e,b*杆位置离散图：频率10.050.1590.45问题三中附件二杆位置

16、离散图：相关代码：x=1.2352 1.2081 1.1813 1.1546 1.1281 1.1018 1.0756 1.0496 1.0237 0.998 0.9724 0.947 0.9217 0.8965 0.8714 0.8464 0.8215 0.7967 0.7719 0.7473 0.7227;y=0.173 0.189 0.2048 0.2203 0.2356 0.2505 0.2653 0.2798 0.294 0.308 0.3218 0.3354 0.3488 0.3619 0.3748 0.3876 0.4001 0.4124 0.4246 0.4366 0.4484

17、;t=12.68333333 12.73333333 12.78333333 12.83333333 12.88333333 12.93333333 12.98333333 13.03333333 13.08333333 13.13333333 13.18333333 13.23333333 13.28333333 13.33333333 13.38333333 13.43333333 13.48333333 13.53333333 13.58333333 13.63333333 13.68333333; e（i）=120-（t（i）-12）*15-15*（t（2）-t（1）-15*（t（2）

18、-t（1）*x（i+1）*y（i）/（y（i+1）*x（i）-x（i+1）*y（i）; n（i）=15*（t（2）-t（1）*y（i）*y（i+1）/（y（i+1）*x（i）-x（i+1）*y（i）;stem（e,n,stem（n,e,y=1.247256205 1.22279459 1.198921486 1.175428964 1.152439573 1.12991747 1.10783548 1.086254206 1.065081072 1.044446265 1.024264126 1.004640314 0.985490908 0.966790494 0.948584735 0.9

19、30927881 0.91375175 0.897109051 0.880973762 0.865492259 0.850504468x=12.68333333 12.73333333 12.78333333 12.83333333 12.88333333 12.93333333 12.98333333 13.03333333 13.08333333 13.13333333 13.18333333 13.23333333 13.28333333 13.33333333 13.38333333 13.43333333 13.48333333 13.53333333 13.58333333 13.

20、63333333 13.68333333;y=h*（1-（sind（a）*sind（23.45*sin（2*3.1415926*（284+n）/365）+cosd（a）*cosd（23.45*sin（2*3.1415926*（284+n）/365）*cosd（15*（x-12）-120+91.75）2）（1/2）/（sind（a）*sind（23.45*sin（2*3.1415926*（284+n）/365）+cosd（a）*cosd（23.45*sin（2*3.1415926*（284+n）/365）*cosd（15*（x-12）-120+91.75）附件三中位置离散图：附件三相关代码：x=

21、1.1637 1.2212 1.2791 1.3373 1.396 1.4552 1.5148 1.575 1.6357 1.697 1.7589 1.8215 1.8848 1.9488 2.0136 2.0792 2.1457 2.2131 2.2815 2.3508 2.4213;y=3.336 3.3299 3.3242 3.3188 3.3137 3.3091 3.3048 3.3007 3.2971 3.2937 3.2907 3.2881 3.2859 3.284 3.2824 3.2813 3.2805 3.2801 3.2801 3.2804 3.2812;t=13.15 1

22、3.2 13.25 13.3 13.35 13.4 13.45 13.5 13.55 13.6 13.65 13.7 13.75 13.8 13.85 13.9 13.95 14 14.05 14.1 14.15; e（i）=120-（t（i）-12）*15+15*（t（2）-t（1）*x（i）*y（i+1）/（y（i）*x（i+1）-x（i）*y（i+1）; n（i）=15*（t（2）-t（1）*y（i）*y（i+1）/（y（i）*x（i+1）-x（i）*y（i+1）;k*y=3.533142184 3.546768029 3.561797643 3.578100715 3.59575078

23、3 3.61493428 3.635425983 3.657218272 3.680541115 3.705167836 3.731278025 3.758917911 3.788087888 3.818701015 3.850809619 3.88458522 3.919911828 3.956875992 3.99553479 4.035750835 4.077863059;x=13.15 13.2 13.25 13.3 13.35 13.4 13.45 13.5 13.55 13.6 13.65 13.7 13.75 13.8 13.85 13.9 13.95 14 14.05 14.1

24、 14.15;y=h*（1-（sind（a）*sind（23.45*sin（2*3.1415926*（284+n）/365）+cosd（a）*cosd（23.45*sin（2*3.1415926*（284+n）/365）*cosd（15*（x-12）-120+115.8）2）（1/2）/（sind（a）*sind（23.45*sin（2*3.1415926*（284+n）/365）+cosd（a）*cosd（23.45*sin（2*3.1415926*（284+n）/365）*cosd（15*（x-12）-120+115.8）问题四相关代码：x=8.916666667 8.95 8.983333333 9.016666667 9.05 9.083333333 9.116666667 9.15 9.183333333 9.216666667 9.25 9.283333333 9.316666667 9.3