4.如權利要求1所述的基于進退法的橫波分裂分析方法，其特征在于，所述步驟S4包括以下子步驟：

(1)對Δt方向零點的搜索；

第一步、確定范圍，范圍的上下界由初始點(θ_k,Δt_k)和比較點(θ_k,Δt_k+1)確定；Δt_k+1等于Δt_k加上步長h；如果COV'_Δt(θ_k,Δt_k)×COV'_Δt(θ_k,Δt_k+1)小于零，得到范圍[Δt_k,Δt_k+1]，進行下一步；如果COV'_Δt(θ_k,Δt_k)×COV'_Δt(θ_k,Δt_k+1)大于零，得到新的Δt_k+1，新Δt_k+1等于新的Δt_k加上新的h，新的Δt_k等于上一步的Δt_k+1，新的h是上一步h的兩倍；然后，再次判斷COV'_Δt(θ_k,Δt_k)×COV'_Δt(θ_k,Δt_k+1)的值是大于零還是小于零，確實是執行下一步或者繼續更新這兩個點；

第二步、確定另外兩個比較點(θ_k,l_t)和(θ_k,u_t)縮小范圍；其中，l_t和u_t的計算公式如下：

l_t＝Δt_k+0.382×(Δt_k+1-Δt_k),

u_t＝Δt_k+0.618×(Δt_k+1-Δt_k).

通過比較|COV'_Δt(θ_k,l_t)|與|COV'_Δt(θ_k,u_t)|的大小，更新l_t或者u_t：

如果|COV'_Δt(θ_k,l_t)||COV'_Δt(θ_k,u_t)|，則

Δt_k＝l_t,

l_t+1＝Δt_k+0.382×(Δt_k+1-Δt_k)

如果|COV'_Δt(θ_k,l_t)||COV'_Δt(θ_k,u_t)|，則

Δt_k+1＝u_t,

u_t+1＝Δt_k+0.618×(Δt_k+1-Δt_k)；

繼續迭代，直到|COV'_Δt(θ_k,l_t)|與|COV'_Δt(θ_k,u_t)|之間的差值小于精度ε_Δt，ε_Δt＝10^-6，即得到Δt方向的零點，Δt＝(l_t+u_t)/2；

(2)對θ方向零點的搜索；

第一步、由初始點(θ_k,Δt)得到比較點(θ_k+1,Δt)，θ_k+1等于θ_k加上步長h；然后，計算COV'_θ(θ_k,Δt)×COV'_θ(θ_k+1,Δt)的值，如果COV'_θ(θ_k,Δt)×COV'_θ(θ_k+1,Δt)0，則范圍為[θ_k,θ_k+1]，執行下一步；如果COV'_θ(θ_k,Δt)×COV'_θ(θ_k+1,Δt)0，繼續更新，得到新的θ_k+1，新的θ_k+1等于新的θ_k加上新的h，新的θ_k,等于上一個θ_k+1，新的h是上一個h的兩倍；然后，再次判斷COV'_θ(θ_k,Δt)×COV'_θ(θ_k+1,Δt)的值是大于零還是小于零，確實是執行下一步或者繼續更新這兩個點；

第二步、確定另外兩個比較點(l_θ,Δt)和(u_θ,Δt)縮小范圍；其中，l_θ和u_θ的計算公式如下：

l_θ＝θ_k+0.382×(θ_k+1-θ_k)

u_θ＝θ_k+0.618×(θ_k+1-θ_k)

然后，通過比較|COV'_θ(l_θ,Δt)|和|COV'_θ(u_θ,Δt)|的大小，更新l_t或者u_t：

如果|COV'_θ(l_θ,Δt)||COV'_θ(u_θ,Δt)|，則

θ_k＝l_θ

l_θ+1＝θ_k+0.382×(θ_k+1-θ_k)

如果|COV'_θ(l_θ,Δt)||COV'_θ(u_θ,Δt)|，則

θ_k+1＝u_θ

u_θ+1＝θ_k+0.618×(θ_k+1-θ_k)

繼續迭代，直到|COV'_θ(l_θ,Δt)|和|COV'_θ(u_θ,Δt)|的差值小于精度ε_θ，ε_θ＝10^-6，即得到θ方向的零點，θ＝(l_θ+u_θ)/2；

(3)計算點(θ,Δt)的COV'，如果COV'的值小于精度ε₃，ε₃＝10^-8，停止循環，將該點作為結果；反之，則繼續重復步驟(1)和(2)交替搜索方向，尋找新的Δt和θ方向的零點，直到結果滿足條件。