CN105553907B - 细频率偏差估计方法及使用其的频率偏差估计装置 - Google Patents

Abstract

本发明提供一种细频率偏差估计方法及使用其的频率偏差估计装置。此细频率偏差估计方法包括接收第k帧，第k帧包括第一帧头以及第二帧头，其中第二帧头具有N行M列的第一符号矩阵，且第k帧之调制方式决定关联于第一符号矩阵的沃尔什‑阿达玛序列；将第一符号矩阵每一行中的所有符号映像至同一象限，以产生N行M列的第二符号矩阵；累加第二符号矩阵的每一行中的所有符号，以获得具有N个累加值的累加向量；根据累加向量进行自相关运算，以获得对应第k帧的第k自相关值；依据第k自相关值计算出细频率偏差。

Description

细频率偏差估计方法及使用其的频率偏差估计装置

技术领域

本发明涉及一种细频率偏差估计方法，且特别是一种适用于任一种帧结构的细频率偏差估计方法以及使用其的频率偏差估计装置。

背景技术

DVB-S2X为欧洲电信标准协会(European Telecommunications StandardsInstitute，ETSI)所提出的新一代卫星数字电视广播传输标准。对此，相较于上一代DVB-S2的基础标准，DVB-S2X适用于在超低信噪比(Very Low Signal Noise Ratio，VL-SNR)下进行信号传输，其中超低信噪比约为-10dB。

DVB-S2X采用了更高阶数的PSK调制，最高可达256APSK，使得DVB-S2X的频谱利用率提升。为了在超低信噪比下进行信号传输，DVB-S2X在调制信号时还引入了新的结构。进一步说，DVB-S2标准下的调制信号仅包括物理层信令(Physical Layer Signaling，PLS)帧头及数据符号(Data symbol)，但在DVB-S2X标准下，调制信号还包括了超低信噪比帧头。

然而，超低信噪比帧头的符号数会随着不同的调制方式而改变。举例来说，使用QPSK进行调制的调制信号与使用π/2BPSK进行调制的调制信号彼此具有不同结构的超低信噪比帧头。也就是说，DVB-S2的频率偏差估计方法在调制方式未知的情形下，无法直接使用在DVB-S2X。因此， 需要一种细频率偏差估计方法之设计，能够在任何一种调制方式下对帧进行细频率偏差估计。

发明内容

本发明实施例提供一种细频率偏差估计方法。此细频率偏差估计方法适用于在超低信噪比下进行信号传输的接收器，其特征在于，细频率偏差估计方法包括步骤A：接收第k帧，第k帧包括第一帧头以及第二帧头，其中第二帧头具有N行M列的第一符号矩阵，M、N为大于0之整数，k为帧索引值，且第k帧之调制方式决定关联于第一符号矩阵的沃尔什-阿达玛序列。步骤B：将第一符号矩阵每一行中的所有符号映像至同一象限，以产生N行M列的第二符号矩阵。步骤C：累加第二符号矩阵的每一行中的所有符号，以获得具有N个累加值的累加向量。步骤D：根据累加向量进行自相关运算，以获得对应第k帧的第k自相关值。步骤E：依据第k自相关值计算出细频率偏差。

本发明实施例提供一种频率偏差估计装置。此频率偏差估计装置用以执行细频率偏差估计方法。频率偏差估计装置适用于在超低信噪比下进行信号传输的接收器，其特征在于，频率偏差估计装置包括自相关运算模块以及频率偏差运算单元。自相关运算模块用以接收第k帧，并根据第k帧进行自相关运算，k为帧索引值。频率偏差运算单元电性连接自相关运算模块，用以计算细频率偏差。细频率偏差估计方法包括步骤A：接收第k帧，第k帧包括第一帧头以及第二帧头，其中第二帧头具有N行M列的第一符号矩阵，M、N为大于0之整数，k为帧索引值，且第k帧之调制方式决定关联于第一符号矩阵的沃尔什-阿达玛序列。步骤B：将第一符号矩阵每一行中的所有符号映像至同一象限，以产生N行M列的第二符号矩阵。步骤C：累加第二符号矩阵的每一行中的所有符号，以获得具有N个累加值的累加向量。步骤D：根据累加向量进行自相关运算，以获得对应第k帧的第k自相关值。步骤E：依据第k自相关值计算出细频率偏差。

综上所述，本发明实施例所提供的细频率偏差估计方法及使用其的频率偏差估计装置可不受调制方式变动的影响并对帧进行细频率偏差估计。换句话说，不论DVB-S2X使用哪一种调制方式来调制帧，本发明实施例所提供的细频率偏差估计方法及使用其的频率偏差估计装置均可对帧进行细频率偏差估计。

为使能更进一步了解本发明之特征及技术内容，请参阅以下有关本发明之详细说明与附图，但是此等说明与所附附图仅系用来说明本发明，而非对本发明的权利范围作任何的限制。

附图说明

图1是本发明实施例之频率偏差估计装置的示意图。

图2是本发明实施例之第k帧的结构示意图。

图3是本发明实施例之第二符号矩阵中每一行的关连性示意图。

图4是本发明实施例之细频率偏差估计方法在不同调制方式下的仿真结果图。

图5是本发明实施例之细频率偏差估计方法在不同帧数下的仿真结果图。

图6是本发明实施例之细频率偏差估计方法的流程图。

图7是本发明实施例之计算细频率偏差的流程图。

具体实施方式

在下文将参看随附附图更充分地描述各种例示性实施例，在随附附图中展示一些例示性实施例。然而，本发明概念可能以许多不同形式来体现， 且不应解释为限于本文中所阐述之例示性实施例。确切而言，提供此等例示性实施例使得本发明将为详尽且完整，且将向本领域的技术人员充分传达本发明概念的范畴。在诸附图中，可为了清楚而夸大示层及区之大小及相对大小。类似数字始终指示类似组件。

应理解，虽然本文中可能使用术语第一、第二、第三等来描述各种组件或信号等，但此等组件或信号不应受此等术语限制。此等术语乃用以区分一组件与另一组件，或者一信号与另一信号。另外，如本文中所使用，术语「或」视实际情况可能包括相关联之列出项目中之任一者或者多者之所有组合。

请参阅图1，图1是本发明实施例之频率偏差估计装置的示意图。频率偏差估计装置1适用于在超低信噪比下进行信号传输的接收器，其中超低信噪比约为-10dB。频率偏差估计装置1包括了自相关运算模块10、极性判断单元11以及频率偏差运算单元12。自相关运算模块10电性连接于极性判断单元11以及频率偏差运算单元12。极性判断单元11电性连接于频率偏差运算单元12。

自相关运算模块10用以接收第k帧，并根据第k帧进行自相关运算以获得对应第k帧的第k自相关值，k为帧索引值。极性判断单元11用以依据相关于第k自相关值的第k帧自相关累加值的实数部分决定正负极性。频率偏差运算单元12依据自相关运算模块10计算的第k自相关值以及极性判断单元11得出的正负极性去计算细频率偏差。细频率偏差表示在一个特定中心频率下，允许频率偏差的值。

当频率偏差估计装置1接收到第k帧(例如第0帧)后会先对第k帧做粗频率偏差估计，以得出粗频率偏差。根据粗频率偏差，频率偏差估计装置1可得知频率偏差的估计范围。接着，根据此估计范围，频率偏差估计装置再对第k帧作细频率偏差估计。

在细频率偏差估计的第一步中，自相关运算模块10接收第k帧。请参阅图2，图2是本发明实施例之第k帧的结构示意图。第k帧的结构包括物理层信令帧头20、超低信噪比帧头21以及数据符号22，其中超低信噪比帧头21具有N行M列的第一符号矩阵h_k(n,m)，M、N为大于0之整数。

以本实施例来说，第一符号矩阵h_k(n,m)为16行56列的矩阵。超低信噪比帧头21除了第一符号矩阵h_k(n,m)中的896个符号外，于第一符号矩阵h_k(n,m)的首尾分别包括两个0，使得超低信噪比帧头21一共包括900个符号。于其他实施例中，超低信噪比帧头21亦可不为900个符号。或者，第一符号矩阵h_k(n,m)亦可不为16行56列的矩阵。总而言之，本发明实施例并不限制超低信噪比帧头21中包含的符号数以及超低信噪比帧头21的结构。

值得一提的是，第k帧的调制方式决定关联于第一符号矩阵h_k(n,m)的沃尔什-阿达玛(Walsh-Hadamard)序列。进一步说，当第k帧系使用QPSK进行调制时所对应的沃尔什-阿达玛序列与使用π/2BPSK进行调制时所对应的沃尔什-阿达玛序列将有所不同。

举例来说，沃尔什-阿达玛序列包括了N个元素p_i，元素p_i为1或-1，i为介于0至N-1的整数。至于p_i取值是1还是-1则是与第k帧的调制方式相关。将沃尔什-阿达玛序列映射到一个N行N列的矩阵P的主对角线上，矩阵P除了主对角线上的元素外，其余元素皆为0。以本实施例来说，矩阵P对应第一符号矩阵h_k(n,m)为16行16列的矩阵。矩阵P的表示如下：

方程式(1)

其中p_i∈{-1,1}，0≤i≤N-1。

接着，自相关运算模块10将第一符号矩阵h_k(n,m)每一行中的所有符号映像至同一象限，以产生N行M列的第二符号矩阵x_k(n,m)。详细地说，自相关运算模块10对第一符号矩阵h_k(n,m)乘上-2，并加上一个元素全为1的单位矩阵，以产生第三符号矩阵c_k(n,m)。自相关运算模块10将第三符号矩阵c_k(n,m)的偶数列与奇数列所有元素分别乘上1与e^-jπ/2，以产生第四符号矩阵s_k(n,m)。最后，自相关运算模块10将第四符号矩阵s_k(n,m)乘上一个解扰矩阵d_k(n,m)，便可产生第二符号矩阵x_k(n,m)，其中解扰矩阵d_k(n,m)∈{1,-1,j,-j}。因此，第四符号矩阵s_k(n,m)乘上解扰矩阵d_k(n,m)的动作可看成加减法运算。

请参阅图3，图3是本发明实施例之第二符号矩阵中每一行的关连性示意图。第二符号矩阵x_k(n,m)可看成由N的区块所组成，且每个区块中的符号属于同一个象限。基于沃尔什-阿达玛序列的特性，第二符号矩阵x_k(n,m)中奇偶相邻的两个区块中的符号相同或相反，例如相邻的两个区块中的符号全为1，或者其中一个区块中的符号全为1，而另一个区块中的符号全为-1。由于第二符号矩阵x_k(n,m)具有沃尔什-阿达玛序列的特性，故利用第二符号矩阵x_k(n,m)进行自相关运算可得出第k帧的最大自相关值。

自相关运算模块10累加第二符号矩阵x_k(n,m)的每一行中的所有符号，以获得具有N个累加值的累加向量r_k(n)。透过累加第二符号矩阵x_k(n,m)的每一行中的所有符号，可提升信号的能量并抑制噪声的能量，进而提升信噪比。详细方程式如下：

方程式(2)

自相关运算模块10根据累加向量r_k(n)进行自相关运算，以获得对应第k帧的第k自相关值R_k。详细方程式如下：

方程式(3)

重复执行上述步骤共N_f次，以分别获得第0帧至第N_f-1帧所分别对应的第0自相关值R₀至第N_f-1自相关值其中N_f为累加帧数。

在算出第0帧至第N_f-1帧各自的自相关值R₀至后，自相关运算模块10将第0帧对应的第0帧自相关累加值Γ₀初始化为第0自相关值R₀，并且依据第k’自相关值R_k'与第k’-1帧的第k’-1帧自相关累加值Γ_k'-1计算出第k’帧对应的第k’帧自相关累加值Γ_k'，其中k’为1至N_f-1的整数。换言之，自相关运算模块10由第0帧对应的第0帧自相关累加值Γ₀依序计算出第1帧对应的第1帧自相关累加值Γ₁至第N_f-1帧对应的第N_f-1帧自相关累加值

详细地说，若第k’-1帧自相关累加值Γ_k'-1与第k’自相关值R_k'相加后的绝对值大于等于第k’-1帧自相关累加值Γ_k'-1减去第k’自相关值的绝对值R_k'，则第k’帧自相关累加值Γ_k'为第k’-1帧自相关累加值Γ_k'-1与第k’自相关值的相加结果R_k'。反之，若第k’-1帧自相关累加值Γ_k'-1与第k’自相关值R_k'相加后的绝对值小于第k’-1帧自相关累加值Γ_k'-1减去第k’自相关值R_k'的绝对值，则第k’帧自相关累加值Γ_k'为第k’-1帧自相关累加值Γ_k'-1减去第k’自相关值R_k'的结果。详细方程式如下：

Γ₀＝R₀ 方程式(4)

方程式(5)

在自相关运算模块10计算出第k’帧自相关累加值Γ_k'(即第N_f-1帧对应的第N_f-1帧自相关累加值)后，极性判断单元11依据第N_f-1帧自相关累加值决定正负极性η。进一步说，极性判断单元11依据第N_f-1帧自相关累加值的实数部分的正负决定正负极性η。这么做的意义在于，由于细频率偏差估计系估计一、四象限中符号的频率偏差，故需要将二、三象限中的符号乘上其实数部分的正负极性η，使得二、三象限中的符号投射至一、四象限。详细方程式如下：

η＝sign{Re{Γ_k'}} 方程式(6)

最后，频率偏差运算单元12累加第0帧自相关累加值Γ₀至第N_f-1帧自相关累加值以产生加总数值Γ，对加总数值Γ乘上正负极性η后的结果取其幅角，并将幅角除以M及2π，即可获得细频率偏差详细方程式如下：

方程式(7)

方程式(8)

请参阅图4，图4是本发明实施例之细频率偏差估计方法在不同调制方式下的仿真结果图。图4之纵轴表示细频率偏差单位为均方根误差(Root-Mean-Square-Error,RMSE)，而横轴表示信噪比，单位为dB。图4之曲线S100、S200、S300、S400、S500分别代表调制方式为QPSK(码率2/9、数据符号22为正常长度)、π/2BPSK(码率1/5、数据符号22为中等长度)、π/2BPSK(码率11/45、数据符号22为中等长度)、π/2BPSK(码率1/5、数据符号22为短长度)、π/2BPSK(码率4/15、数据符号22为短长度)时进行细频率偏差估计时所得到的细频率偏差其中正常长度代表数据符号22为64800个符号、中等长度代表数据符号22为32400个符号、短长度代表数据符号22为16200个符号。

由图4可知，不论使用什么样的调制方式对帧进行调制，本发明实施例所提供的细频率偏差估计方法皆可准确地估算细频率偏差换句话说，本发明实施例所提供的细频率偏差估计方法的性能与调制方式无关。

请参阅图5，图5是本发明实施例之细频率偏差估计方法在不同帧数下的仿真结果图。图5之曲线S100’、S200’、S300’、S400’、S500’分别代表累计1个帧、2个帧、4个帧、8个帧、16个帧的自相关值进行细频率 偏差估计时所得到的细频率偏差由此可知，参与运算的帧越多的话，本发明实施例之细频率偏差估计方法所估算出来的细频率偏差越精准。综使信噪比在-12dB下，细频率偏差仍能控制在10^-4范围内。

本发明实施例系累加第0帧至第N_f-1帧各自的自相关值进行细频率偏差估计。于其他实施例中，频率偏差运算单元12亦可只接收一个帧(例如第0帧)便进行细频率偏差估计。然而，优选地，出于提高低信噪比下运算的性能，本发明实施例之细频率偏差估计方法会累计多个帧的自相关值进行细频率偏差估计。

请参阅图6，图6是本发明实施例之细频率偏差估计方法的流程图。于步骤S601，开始细频率偏差估计方法。于步骤S602，接收第k帧，k为0至N_f-1的整数，第k帧包括物理层信令帧头与超低信噪比帧头，其中超低信噪比帧头具有N行M列的第一符号矩阵。于步骤S603，将第一符号矩阵每一行中的所有符号映像至同一象限，以产生N行M列的第二符号矩阵。

于步骤S604，累加第二符号矩阵的每一行中的所有符号，以获得具有N个累加值的累加向量。于步骤S605，根据累加向量进行自相关运算，以获得对应第k帧的第k自相关值。于步骤S606，判断是否获得第0帧至第N_f-1帧所分别对应的第0自相关值至第N_f-1自相关值。于步骤S607，根据第0自相关值至第N_f-1自相关值计算细频率偏差。于步骤S608，结束细频率偏差估计方法。

请参阅图7，图7是本发明实施例之计算细频率偏差的流程图。于步骤S701，开始根据第0自相关值至第N_f-1自相关值计算细频率偏差。于步骤S702，将第0帧对应的第0帧自相关累加值初始化为第0自相关值。于步骤S703，依据第k’自相关值与第k’-1帧的第k’-1帧自相关累加值计算出第k’帧对应的第k’帧自相关累加值，其中k’为1至N_f-1的整数。

于步骤S704，依据第N_f-1帧对应的第N_f-1帧自相关累加值的实数部分的正负决定正负极性。于步骤S705，累加第0帧自相关累加值至第N_f-1帧自相关累加值，以产生加总数值。于步骤S706，对加总数值乘上正负极性后的结果取其幅角，并将幅角除以M及2π，以获得细频率偏差。于步骤S707，结束计算细频率偏差。

综上所述，本发明实施例所提供的细频率偏差估计方法及使用其的频率偏差估计装置可不受调制方式变动的影响并对帧进行细偏差估计。换句话说，不论DVB-S2X使用哪一种调制方式来调制帧，本发明实施例所提供的细频率偏差估计方法及使用其的频率偏差估计装置均可对帧进行细频率偏差估计。此外，本发明实施例所提供的细频率偏差估计方法的性能稳定，且只要变更参与运算的帧的数量即可调整细频率偏差估计方法的性能。

另一方面，本发明实施例所提供的细频率偏差估计方法的运算量适中。当参与运算的帧的数量为N_f时，执行复数乘法的次数为N_f(N-1)次。当参与运算的帧的数量为1时，执行复数乘法的次数仅为(N-1)次。由此可知，本发明实施例所提供的细频率偏差估计方法的运算量并不会太多，使用者亦可自行增加或减少参与运算的帧的数量，以调整细频率偏差估计方法整体的运算量。

以上所述，仅为本发明最优选之具体实施例，惟本发明之特征并不局限于此，本领域的技术人员在本发明之领域内，可轻易思及之变化或修饰，皆可涵盖在以下本案之专利范围。

【符号说明】

1：频率偏差估计装置

10：自相关运算模块

11：极性判断单元

12：频率偏差运算单元

20：物理层信令帧头

21：超低信噪比帧头

22：数据符号

S100、S200、S300、S400、S500：曲线

S100’、S200’、S300’、S400’、S500’：曲线

S601～S608：步骤流程

S701～S707：步骤流程。

Claims (14)

1.一种细频率偏差估计方法，适用于在超低信噪比下进行信号传输的一接收器，其特征在于，该细频率偏差估计方法包括：

步骤A：接收一第k帧，该第k帧包括一第一帧头以及一第二帧头，其中该第二帧头具有N行M列的一第一符号矩阵，M、N为大于0之整数，k为一帧索引值，且该第k帧之一调制方式决定关联于该第一符号矩阵的一沃尔什-阿达玛序列；

步骤B：将该第一符号矩阵每一行中的所有符号映像至同一象限，以产生N行M列的一第二符号矩阵；

步骤C：累加该第二符号矩阵的每一行中的所有符号，以获得具有N个累加值的一累加向量；

步骤D：根据该累加向量进行一自相关运算，以获得对应该第k帧的一第k自相关值；以及

步骤E：依据该第k自相关值计算出一细频率偏差；

其中该细频率偏差估计方法从k＝0至N_f-1，执行该步骤A～D共N_f次，以分别获得该第0帧至该第N_f-1帧所分别对应的该第0自相关值至该第N_f-1自相关值，其中N_f为一累加帧数；以及该步骤E是根据该第0自相关值至该第N_f-1自相关值计算出该细频率偏差；

其中该步骤E是依据相关于该第k自相关值的一第k帧自相关累加值的实数部分决定一正负极性，以及依据该第k自相关值及该正负极性计算该细频率偏差。

2.根据权利要求1所述的细频率偏差估计方法，其特征在于，其中该步骤E包括：

步骤E1：将该第0帧对应的一第0帧自相关累加值初始化为该第0自相关值，并且依据该第k’自相关值与该第k’-1帧的一第k’-1帧自相关累加值计算出该第k’帧对应的一第k’帧自相关累加值，其中k’为1至N_f-1的整数；

步骤E2：依据该第N_f-1帧对应的该第N_f-1帧自相关累加值决定一正负极性；以及

步骤E3：依据该正负极性与该第0帧自相关累加值至该第N_f-1帧自相关累加值计算出该细频率偏差。

3.根据权利要求2所述的细频率偏差估计方法，其特征在于，其中于该步骤E1中，若该第k’-1帧自相关累加值与该第k’自相关值相加后的绝对值大于等于该第k’-1帧自相关累加值减去该第k’自相关值的绝对值，则该第k’帧自相关累加值为该第k’-1帧自相关累加值与该第k’自相关值的相加结果；若该第k’-1帧自相关累加值与该第k’自相关值相加后的绝对值小于该第k’-1帧自相关累加值减去该第k’自相关值的绝对值，则该第k’帧自相关累加值为该第k’-1帧自相关累加值减去该第k’自相关值的结果。

4.根据权利要求2所述的细频率偏差估计方法，其特征在于，其中于该步骤E2中，依据该第N_f-1帧自相关累加值的实数部分的正负决定该正负极性。

5.根据权利要求2所述的细频率偏差估计方法，其特征在于，其中于该步骤E3中，累加该第0帧自相关累加值至该第N_f-1帧自相关累加值以产生一加总数值，对该加总数值乘上该正负极性后的结果取其幅角，并将该幅角除以M及2π，以获得该细频率偏差。

6.根据权利要求1所述的细频率偏差估计方法，其特征在于，其中于该步骤B中，对该第一符号矩阵乘上-2后，加上一个元素全为1的一单位矩阵后，产生一第三符号矩阵，将该第三符号矩阵的偶数列与奇数列所有元素分别乘上1与e^-jπ/2，以产生一第四符号矩阵，以及将该第四符号矩阵乘上一个解扰矩阵，以产生该第二符号矩阵。

7.根据权利要求1所述的细频率偏差估计方法，其中该细频率偏差估计方法适用于一延伸的第二代数字卫星广播系统中的该接收器。

8.一种频率偏差估计装置，用以执行一细频率偏差估计方法，该频率偏差估计装置适用于在超低信噪比下进行信号传输的一接收器，其特征在于，该频率偏差估计装置包括：

一自相关运算模块，用以接收一第k帧，并根据该第k帧进行一自相关运算，k为一帧索引值；以及

一频率偏差运算单元，电性连接该自相关运算模块，用以计算一细频率偏差；

其中，该细频率偏差估计方法包括：

步骤A：该自相关运算模块接收该第k帧，该第k帧包括一第一帧头以及一第二帧头，其中该第二帧头具有N行M列的一第一符号矩阵，M、N为大于0之整数，且该第k帧之一调制方式决定关联于该第一符号矩阵的一沃尔什-阿达玛序列；

步骤B：该自相关运算模块将该第一符号矩阵每一行中的所有符号映像至同一象限，以产生N行M列的一第二符号矩阵；

步骤C：该自相关运算模块累加该第二符号矩阵的每一行中的所有符号，以获得具有N个累加值的一累加向量；

步骤D：该自相关运算模块根据该累加向量进行该自相关运算，以获得对应该第k帧的一第k自相关值；以及

步骤E：该频率偏差运算单元依据该第k自相关值计算出该细频率偏差；

其中该自相关运算模块从k＝0至N_f-1，执行该步骤A～D共N_f次，以分别获得该第0帧至该第N_f-1帧所分别对应的该第0自相关值至该第N_f-1自相关值，其中N_f为一累加帧数；以及该步骤E是该频率偏差运算单元根据该第0自相关值至该第N_f-1自相关值计算出该细频率偏差；

其中该步骤E是依据相关于该第k自相关值的一第k帧自相关累加值的实数部分决定一正负极性，以及依据该第k自相关值及该正负极性计算该细频率偏差。

9.根据权利要求8所述的频率偏差估计装置，其特征在于，其中该步骤E包括：

步骤E1：该自相关运算模块将该第0帧对应的一第0帧自相关累加值初始化为该第0自相关值，并且依据该第k’自相关值与该第k’-1帧的一第k’-1帧自相关累加值计算出该第k’帧对应的一第k’帧自相关累加值，其中k’为1至N_f-1的整数；

步骤E2：该频率偏差估计装置之一极性判断单元依据该第N_f-1帧对应的该第N_f-1帧自相关累加值决定一正负极性；以及

步骤E3：该频率偏差运算单元依据该正负极性与该第0帧自相关累加值至该第N_f-1帧自相关累加值计算出该细频率偏差。

10.根据权利要求9所述的频率偏差估计装置，其特征在于，其中于该步骤E1中，若该第k’-1帧自相关累加值与该第k’自相关值相加后的绝对值大于等于该第k’-1帧自相关累加值减去该第k’自相关值的绝对值，则该第k’帧自相关累加值为该第k’-1帧自相关累加值与该第k’自相关值的相加结果；若该第k’-1帧自相关累加值与该第k’自相关值相加后的绝对值小于该第k’-1帧自相关累加值减去该第k’自相关值的绝对值，则该第k’帧自相关累加值为该第k’-1帧自相关累加值减去该第k’自相关值的结果。

11.根据权利要求9所述的频率偏差估计装置，其特征在于，其中于该步骤E2中，该极性判断单元依据该第N_f-1帧自相关累加值的实数部分的正负决定该正负极性。

12.根据权利要求9所述的频率偏差估计装置，其特征在于，其中于该步骤E3中，该频率偏差运算单元累加该第0帧自相关累加值至该第N_f-1帧自相关累加值以产生一加总数值，对该加总数值乘上该正负极性后的结果取其幅角，并将该幅角除以M及2π，以获得该细频率偏差。

13.根据权利要求8所述的频率偏差估计装置，其特征在于，其中于该步骤B中，该自相关运算模块对该第一符号矩阵乘上-2后，加上一个元素全为1的一单位矩阵后，产生一第三符号矩阵，该自相关运算模块将该第三符号矩阵的偶数列与奇数列所有元素分别乘上1与e^-jπ/2，以产生一第四符号矩阵，以及该自相关运算模块将该第四符号矩阵乘上一个解扰矩阵，以产生该第二符号矩阵。

14.根据权利要求8所述的频率偏差估计装置，其中该频率偏差估计装置适用于一延伸的第二代数字卫星广播系统中的该接收器。