基于遍歷搜索算法的太陽影子定位技術研究

王一凡+張樂權+尹澤文+白昊東+陶玲玲

摘 要：該文對太陽影子定位技術進行了研究，在合理假設的基礎上，結合Cooper經驗公式，構建了消除物體長度的單目標優(yōu)化模型。通過MATLAB編寫遍歷搜索算法進行求解，最終得到物體所在的經緯度，并且誤差較小。由于該模型不需要物體的長度，適用范圍較廣，可以進一步應用到視頻數據分析定位中。

關鍵詞：太陽影子定位 Cooper經驗公式 單目標優(yōu)化模型 MATLAB遍歷搜索算法

中圖分類號：P128 文獻標識碼：A 文章編號：1674-098X（2017）04（c）-0237-02

伴隨著21世紀大數據時代的到來，數據的利用顯得越來越重要。對于生活中大量的視頻數據，如何充分利用這些數據進行定位，在軍事上和政治上顯得尤為重要。在某些特殊的情況下，可以通過視頻中的影子變化來定位；但是其核心還是太陽影子的定位問題。

太陽影子定位技術是通過在一段時間內，太陽影子隨著時間的變化來對所處位置進行定位。

對于太陽影子定位技術主要考慮以下3個問題。

（1）建立太陽影子長度與各影響因素的數學模型。

（2）考慮到實際情況，該文著重研究在物體長度未知的情況下，如何通過太陽影子的變化確定物體所在地點和拍攝日期。

（3）通過對模型的改進，來消除物體高度未知對模型的影響。

1 數據來源及模型假設

該文中物體太陽影子的頂點坐標數據來自2015全國大學生數學競賽A題附件[2]。

為了合理地解決問題提出以下假設：（1）不考慮大氣對太陽光的折射。（3）太陽的運行速度是勻速運動。（4）太陽光為平行光。

2 定位模型的建立

2.1 研究思路

首先分析影子長度關于太陽高度角、桿長、時角的關系式；其次確定方位角關于太陽高度角、赤緯角的表達式。聯(lián)立方程組得到最終的影長變化公式，再以影長變化平方最小建立單目標優(yōu)化模型。

2.2 基礎模型

2.2.1 影子長度

影子長度由太陽高度角hs和桿長直接決定，其中太陽高度角是太陽光射在地上與地面形成的夾角，可以得到初步的影長公式：

hs

2.2.2 太陽高度角hs

太陽高度角是太陽光的入射方向和地平面之間的夾角，更準確地說是太陽光的入射方向與地表切線形成的夾角，隨太陽赤緯角δ、太陽方位角的變化而變化，表達式為：

式中：φ表示緯度，ω表示太陽時角。

2.2.3 太陽赤緯角δ

太陽赤緯角是太陽直射點的緯度，其變化造成了南北半球晝夜的長短變化，直射點在北半球晝長夜短，直射點在南半球晝短夜長（以北半球為例），直射點緯度范圍在-N。由Cooper經驗方程得：

式中：n表示某時期距1月1日的天數，1月1日時n。

2.2.4 太陽時角ω

太陽時角是某時太陽光線與當地正午時間太陽光線所成的角度，每隔1個小時變化15°，當地正午時刻太陽時角為。查閱相關文獻[2]，可知：

式中：表示子午線當地標準時間， 表示當地經度。

2.3 目標優(yōu)化模型

由于特定時期、特定經緯度的影子太陽高度角和太陽方位角都可以求出，且一天中的不同時刻桿的影長僅與太陽高度角有關，故特定時期、特定經緯度的影長與太陽方位角都可求出。

常見的定位問題是已知影長和桿長尋找目標地點、拍攝日期，若是不知道桿的長度，那么目標模型需要去除桿長的影響。

對于任意，時刻，分別為對應時刻的太陽高度角，由影長公式hs，兩時刻理論影長比值為，從而就可以消除桿長的影響。

模型將實際影長比值與理論影長比值差值的平方和最小為目標函數：

注：由文獻附件提供的數據有21組，為使數據有明顯差異性，以間隔10個點為分組，即第1個數據和第11個數據為一組，第2個數據和第12個數據為一組，以此類推。

3 算法求解

對于定位問題，有擬牛頓迭代法、多重搜索算法、人工魚群算法等多種算法。其中遍歷算法是一種較為可行的、操作性強的方法。算法流程：

（1）設定初始狀態(tài)、循環(huán)條件、候選解確定準則。

（2）初步確定范圍。設定精度ε1、幅度，滿足判定準則，輸出候選解，繼續(xù)循環(huán)，直至循環(huán)結束。

（3）確定估計值。設定精度、幅度，以候選解為中心確定搜索范圍，分部式循環(huán)搜索。

在筆者模型求解過程中設定精度為，幅度為，以參考文獻[2]附件1數據為例，定位地點為（109.9E，19.4N），參考標準給的（110E，18N）。經緯度相差不超過，可以進行精確的定位。

誤差主要來源：模型中部分公式為經驗公式，存在誤差；幅度過大，可進行三重搜索。

4 結語

該文通過建立單目標優(yōu)化模型來研究太陽影子定位技術，研究了在物體長度未知的情況下，如何通過太陽影子的變化確定物體所在地點和拍攝日期。

創(chuàng)新點在于通過以影長變化平方和最小建立單目標優(yōu)化模型，消除了物體長度未知的影響。最后結合Matlab使用遍歷算法進行求解，解出較為精確的拍攝地址。由于沒有考慮大氣折射的影響，因此在后續(xù)的研究中可以將大氣折射等影響因素加入，從而進一步定位的精確度。

參考文獻

[1] 溫欣妍.Matlab R2016a從入門到精通[M].北京：清華大學出版社，2017.

[2] 2015年高教社杯全國大學生數學建模競賽賽題[EB/OL].（2016-09-22）.http：//special.univs.cn/service/jianmo/sxjmtmhb/2015/0910/1107276.shtml.

[3] Position of the Sun[EB/OL].（2016-11-13）.https：//en.wikipedia.org/wiki/Position_of_the_Sun.

[4] SundialTime Correction - Equation of Time[EB/OL].（2016-11-03）.http：//www.mb-soft.com/public3/equatime.html.

[5] 蔡志杰.太陽影子定位[J].數學建模及其應用，2015（4）：25-33.