技術特征：

1.一種基于塊狀導頻的迭代相關符號定時估計方法，包括如下步驟：

1)在接收端的每一幀數據的每個子幀中，用rⁱ表示第i個帶循環(huán)前綴的時域導頻序列，其中i＝1,2，去除其循環(huán)前綴得到時域導頻序列該是一個N點序列，N＝128,256,512,1024,1536,2048；

2)對時域導頻序列進行快速傅里葉變換FFT，得到N點頻域導頻序列Rⁱ，對該Rⁱ進行M/2點循環(huán)移位得到頻域導頻移位序列Fⁱ，再對該Fⁱ進行集中式子載波解映射得到M點接收頻域導頻序列Sⁱ(m)，并在接收端提取Sⁱ(m)對應的M點本地頻域導頻序列Cⁱ(m)，其中M是對應于N的，表示子載波的數量，M＝72,180,300,600,900,1200，m＝0,1,…,M-1；

3)利用接收頻域導頻序列Sⁱ(m)計算各個頻域導頻序列處信道二范數的二次冪

4)求出信道二范數的二次冪的最大值對應的頻域導頻編號取出i_max對應的接收頻域導頻序列和本地頻域導頻序列

5)將和分別分割成長度均為p的若干段，將分割的接收頻域導頻序列的相鄰兩段表示為S_a(k)和S_a+1(k+p)，將分割的本地頻域導頻序列的相鄰兩段表示為C_a(k)和C_a+1(k+p)，其中τ_max為系統(tǒng)最大符號定時誤差，a＝0,1,2,…,d₁-2，d₁＝M/p，a·p≤k≤(a+1)·p-1；

6)根據5)得到的S_a(k)、S_a+1(k+p)、C_a(k)、C_a+1(k+p)，計算和的分段相關值α_a：

其中上標*表示共軛轉置；

再利用反三角函數求出α_a對應的角度∠α_a＝angle(α_a)；

7)根據∠α_a計算各個分段相關值α_a對應的定時估計值：對該τ_a求算數平均值并四舍五入，得到第一次定時誤差的估計值其中表示向下取整；

8)根據估計值τ¹和2)中的子載波解映射方式構造M維修正對角矩陣Λ：

$\mathrm{\Λ}=diag(e^{j2\mathrm{\π}\·(-\frac{M}{2})\·{\mathrm{\τ}}^1/N},e^{j2\mathrm{\π}\·(-\frac{M}{2}+1)\·{\mathrm{\τ}}^1/N},...,e^{j2\mathrm{\π}\·(\frac{M}{2}-1)\·{\mathrm{\τ}}^1/N}),$

并利用該修正對角矩陣Λ計算接收頻域導頻修正序列：

其中Λ(m,m)表示矩陣Λ的第m行第m列元素；

9)將接收頻域導頻修正序列和本地頻域導頻序列分別分割成長度均為q的若干段，將分割的接收頻域導頻修正序列的相鄰兩項表示為Y_b(t)和Y_b+1(t+q)，將分割的本地頻域導頻序列的相鄰兩項表示為C_b(t)和C_b+1(t+q)，其中q>p，b＝0,1,2,…,d₂-2，d₂＝M/q，b·q≤t≤(b+1)·q-1；

10)根據9)得到的Y_b(t)、Y_b+1(t+q)、C_b(t)、C_b+1(t+q)，計算和的分段相關值β_b：

其中上標*表示共軛轉置；

再利用反三角函數求出β_b對應的角度∠β_b＝angle(β_b)；

11)根據∠β_b計算各個分段相關值β_b對應的定時估計值對該ν_b求算數平均值并四舍五入，得到第二次定時誤差的估計值

12)根據以上兩次定時估計，得到最終的接收子幀的符號定時誤差估計值τ＝τ¹+τ²。

2.根據權利要求1所述的方法，其中步驟1)中第i個帶循環(huán)前綴的時域導頻序列rⁱ，表示如下：

$r^i={\[r_{CP}^i,r_p^i\]}^T,$

其中，循環(huán)前綴表示為：

時域導頻序列表示為：

N＝128,256,512,1024,1536,2048，表示子載波個數，L＝18,20，表示循環(huán)前綴的長度，其中上標T表示轉置。

3.根據權利要求1所述的方法，其中步驟2)中的快速傅里葉變換，表示如下：

Rⁱ(n)＝FFT{rⁱ(k)}，

其中k＝0,1,…,N-1，n＝0,1,…,N-1。

4.根據權利要求1所述的方法，其中步驟2)中的N點頻域導頻序列Rⁱ，表示如下：

Rⁱ＝[Rⁱ(0),Rⁱ(1),…,Rⁱ(N-1)]^T。

5.根據權利要求1所述的方法，其中步驟2)中得到的頻域導頻移位序列Fⁱ，表示如下：

Fⁱ＝[Rⁱ(N-M/2),Rⁱ(N-M/2+1),…,Rⁱ(N-1),Rⁱ(0),Rⁱ(1),…,Rⁱ(M/2-1)]^T，

其中Rⁱ(·)括號中的值先從N-M/2依次增加至N-1，再從0依次增加至M/2-1。

6.根據權利要求1所述的方法，其中步驟2)中對頻域導頻移位序列Fⁱ進行集中式子載波解映射，是通過將Fⁱ與一個集中式子載波解映射矩陣D相乘得到：

D·Fⁱ＝[Sⁱ(0),Sⁱ(1),…,Sⁱ(m),…,Sⁱ(M-1)]^T，

其中D＝[I_M 0_M×(N-M)]，Sⁱ(m)是M點接收頻域導頻序列的第m項，m＝0,1,…,M-1。