a)開始于格林尼治時間-4712年1月1日12:00:00,它是 通過格林尼治時間而得到的。
[0055] JD = INT (365. 25 X (Y+4716)) +INT (30. 6001 X (M+l) +D+B-1524. 5 (1)
[0056] 其中：INT是計算項的整數(shù)部分，Y表示年，M表示月，D表示天數(shù)需要用十進制計 算，而B可由下式計算：
[0058] 進而，計算儒略歷書日JDE，儒略世紀JC，儒略歷書世紀JCE
[0060] 其中 Δ T 為世界時（Universal Time)與協(xié)調(diào)世界時（Coordinated Universal Time)之間時間差值。
[0061] 儒略歷書千年JME可由下式計算：
[0062] JME = JCE/10 (4)
[0063] 步驟S2,計算日心坐標地球經(jīng)度、煒度和半徑矢量。
[0064] 日心坐標地球經(jīng)度L，日心坐標地球煒度B，半徑矢量R可由如下公式獲取。
[0066] 上式中，Ap C1的具體數(shù)值可在星歷矩陣中通過查表法獲取，η是LO在表中的 行數(shù)，用類似的方法可以計算LU L2、L3、L4、L5,日心坐標地球經(jīng)度L由下式計算可得：
[0068] 類似上述的計算過程，可用同樣的方式查表算出日心坐標地球煒度B、半徑矢量 R0
[0069] 步驟S3,計算地心坐標地球經(jīng)度和煒度
[0070] 根據(jù)下式計算地心經(jīng)度Θ和地心煒度β :
[0072] 步驟S4,計算經(jīng)度和傾斜角章動。
[0073] 結(jié)合星歷表中的數(shù)據(jù)A1J1X1J 1和月球與太陽的平均距角X。、太陽的平均平近點 角X1、月球的平均平近點角X2、月球煒度參數(shù)X 3、春秋分時月球旋轉(zhuǎn)軌道與黃道偏角X4，利用 下式可計算出每一行中的經(jīng)度A φ JP傾斜角Δ ε 1<3
[0075] 則經(jīng)度的章動Δ φ和傾斜角的章動Δ ε的計算公式如下：
[0077] 步驟S5,計算黃道實際傾斜角。
[0078] 黃赤交角ε可以通過黃道平均傾斜角ε。表示如下：
[0079] ε = ε 0/3600+Δ ε (10)
[0080] 步驟S6,計算校正偏差。
[0081] 采用下式校正偏差Λ τ :
[0082] Δ τ = -20. 4898/(3600 · R) (11)
[0083] 步驟S7,計算太陽經(jīng)度。
[0084] 計算太陽視黃經(jīng)λ為：
[0085] λ=θ+Δφ + Δτ (12)
[0086] 步驟S8,計算給定格林尼治時間平均恒星時。
[0087] 格林尼治視太陽時V使用平均恒星時V。表示如下：
[0088] V = V0+ Δ φ X cos ( ε ) (13)
[0089] 步驟S9,計算太陽赤經(jīng)。
[0090] 將太陽赤經(jīng)弧度α制限制在〇度到360度之間，計算公式如下：
[0092] 步驟S10,計算太陽赤煒。
[0093] 計算太陽赤煒δ表示如下：
[0094] δ = arc sin (sin β X cos ε+cos β X sin ε X sin λ ) (15)
[0095] 太陽在天球赤道北面時δ為正，太陽在天球赤道南面時δ為負，并將δ以度為 單位表示。
[0096] 步驟Sll計算觀測點當?shù)貢r間角。
[0097] 計算觀測點的局部時間角H以格林尼治視太陽時V、觀測點地理經(jīng)度〇和太陽赤 經(jīng)弧度α表示為：
[0098] H= V + σ -α (16)
[0099] 步驟S12,計算視日赤經(jīng)。
[0100] 視日赤經(jīng)差Δ α可以用赤道平面視差ξ和觀測器的局部時間角H計算如下：
[0101]
[0102] 上式中妒為觀測點煒度，E為觀測點的海拔高度。
[0103] 那么視日赤經(jīng)α '和視日偏差δ '可用太陽赤煒δ和太陽赤經(jīng)α如下表示：
[0105] 步驟S13，計算視日時間角。
[0106] 視日時間角H'可通過觀測點的局部時間角H與視日赤經(jīng)差Δ α的差值獲取。
[0107] 步驟S14,計算太陽天頂角和高度角。
[0108] 太陽的高度角Θ為理想太陽高度角e。和考慮空氣折射偏差Ae之和。
[0110] 上式中P為觀測點年平均氣壓，T為觀測點年平均溫度。
[0111] 太陽天頂角與太陽高度角符合互余的關(guān)系。
[0112] 步驟S15,計算太陽方位角。
[0113] 太陽方位角?？捎上率奖硎荆?br>[0115] 至此北半球上某一確定地點的太陽入射光線的高度角和方位角可通過上述步驟 實時獲取，
[0116] 步驟102,根據(jù)當?shù)谿PS時鐘（具體年月日時分秒）、時區(qū)、經(jīng)煒度、氣壓、氣溫、海 拔高度等氣象信息利用太陽運動算法計算出該地點時間的太陽光入射實時方向角與高度 角，如圖4所示為本發(fā)明中太陽運動算法輸入輸出接口模塊示意圖，進而太陽光線信息作 為本發(fā)明定日鏡姿態(tài)角度算法的輸入?yún)?shù)數(shù)據(jù)
[0117] 步驟103,根據(jù)太陽入射光線的實時方位角和高度角、定日鏡和塔的位置排布，建 立太陽光線與定日鏡和塔的反射模型。
[0118] 對太陽能熱發(fā)電而言，將太陽光線準確高效的收集至關(guān)重要，因此有必要將太陽 光、定日鏡、接收塔綜合建模，分析定日鏡的控制策略。如圖5為本發(fā)明中建立的定日鏡自 動跟蹤系統(tǒng)空間建模示意圖，以地平面建立定日鏡理論模型，將正東方向定義為坐標系X 軸正方向，將正北方向定義為坐標系Y軸正方向，將天頂方向定義Z軸正方向（符合坐標系 定義右手定則），XYZ坐標系原點0點為接收塔與X-Y水平面交點，假設接收塔塔頂高度為 h，定日鏡中心點坐標為（x，y，0)，可以利用向量法獲取定日鏡方位角辦與高度角Θ。
[0119] 在模型建立時需要做以下假設：1)由于在實際的太陽光入射過程中，太陽光線并 非完全平行，存在一個約32'微小的角度差，在距離較遠時可能會出現(xiàn)一定的誤差，建模時 需假定投射到定日鏡上的所有太陽入射光線均為平行光；2)假設鏡面平整光滑，鏡面厚度 不計，鏡面中心點與鏡面的固定點之間兩點重合，不存在機械偏差偏差。
[0120] 本發(fā)明中定日鏡采用高度角-方位角的雙軸跟蹤方式，根據(jù)定日鏡高度軸和方向 軸這兩個軸，可以確定定日鏡當前此刻的姿態(tài)，如圖6所示為定日鏡運動過程分析圖，為了 簡述方便，對其中參數(shù)定義如下：
[0121] ?:太陽入射光線法向量；
[0122] h太陽光線經(jīng)定日鏡反射到吸熱器塔頂目標點向量；
[0123] 定日鏡鏡面的法向量，垂直于定日鏡鏡面；
[0124] Θ :定日鏡的高度角，定日鏡的鏡面與水平面之間的夾角；
[0125] 夢:定日鏡的方位角，定日鏡鏡面法向量與Y軸（正北方向）的夾角為定日鏡的方 位角。
[0126] 步驟104,采用時控的方法計算定日鏡方位角和高度角，實時調(diào)節(jié)鏡面的角度。
[0127] 定日鏡的姿態(tài)由高度角Θ和方位角供唯一確定，根據(jù)光學原理，為了保證定日鏡 能將反射光線準確的投射到吸熱器塔頂目標點上，根據(jù)計算好的太陽運動軌跡和定日鏡當 前位置，只需確定定日鏡的鏡面法向量?即可完成追日。
[0128] 根據(jù)坐標系定義，由太陽照射到定日鏡中心點B點，太陽光線入射單位向量;？.為：
[0130] 塔頂?shù)淖鴺藶椋?,0，h)，定日鏡中心B點坐標為（x，y，0)，那么從定日鏡中心點反 射到吸熱塔頂接收靶的方向單位向量F為：
[0132] 因為太陽入射光線經(jīng)過定日鏡實時反射到吸熱塔接收靶上，那么定日鏡的法線單 位向量《平分太陽入射光線單位向量和反射光線單位向量產(chǎn)，即：
[0134]