通信系统及其相位误差估计方法
【专利摘要】本发明提供一种通信系统及其相位误差估计方法，其中包含一接收模块与一相位误差估计模块。该接收模块接收一输入信号x。该输入信号x于极座标系统中具有一输入相位θ。利用针对函数f(x,θ)的自然对数施以偏微分运算，该相位误差估计模块产生该输入信号x的一相位误差估计值。f(x,θ)代表该接收模块接收到该输入信号x的一机率函数。
【专利说明】
通信系统及其相位误差估计方法
技术领域
[0001] 本发明与通信系统相关，并且尤其与通信系统中的相位误差估计技术相关。
【背景技术】
[0002] 随着电子相关技术的进步，各种类型的通信设备愈来愈普及。目前有许多采用同 步解调（coherent demodulation)技术的数位通信系统，例如第二代数位电视卫星广播 (digital video broadcasting-satellite, DVB-S2)系统。由于此类系统的接收端必须 要能准确估计接收到的信号存在多大的相位误差，因此以可加性白色高斯杂讯（additive white Gaussian noise，AWGN)模型为基础的相位误差估计近年来在信号处理领域中成为 一个被重视且广泛研究的议题。在现行技术中，于接收器追踪相位变化的相位误差侦测 器（phase error detector, PED)可分为无资料辅助（non-data aided, NDA)、决策导引 (decision directed, DD)、软性决策导引（soft decision directed, SDD)等几大类。
【发明内容】
[0003] 本发明提出一种新的通信系统及其相位误差估计方法，可满足克拉梅-劳下限 (Cramer-Rao bound, CRB)，亦即可最小化接收端的相位误差估计值和真实相位误差值间的 差异。本发明的概念可应用在多种通信系统。
[0004] 根据本发明的一具体实施例为一种通信系统，其中包含一接收模块与一相位误差 估计模块。该接收模块用以接收一输入信号X。该输入信号X于极座标系统中具有一输入 相位θ。f( x，Θ)代表该接收模块接收到该输入信号X的一机率函数。该相位误差估计模 块根据下列方程式产生该输入信号X的一相位误差估计值f :
[0005]
[0006] 根据本发明的另一具体实施例为一种通信系统，其中包含一接收模块与一相位误 差估计模块。该接收模块用以接收一输入信号。该输入信号对应于一传送端利用一 Μ相-相 位偏移调制（M-phase shift keying，M_PSK)所产生的一输出信号。数值Μ为二的任一正 整数次方。该输入信号于极座标系统中具有一输入信号半径r与一输入相位Θ。该相位误 差估计模块根据该输入信号半径r产生一第一谐波系数匕(r)与一第二谐波系数h 2 (r)，并 根据该输入相位Θ产生一第一正弦函数sin (Μ Θ )与一第二正弦函数sin (2M Θ )。随后，该 相位误差估计模块将该第一谐波系数h (r)乘以该第一正弦函数sin (Μ Θ )，以产生一第一 运算结果，将该第二谐波系数h2 (r)乘以该第二正弦函数sin (2M Θ )，以产生一第二运算结 果。根据该第一运算结果与该第二运算结果，该相位误差估计模块产生该输入信号的一相 位误差估计值§。
[0007] 根据本发明的另一具体实施例为一种应用于一通信系统的相位误差估计方法，用 以估计一输入信号X的一相位误差。首先，该相位误差估计方法执行一判断步骤，判断该输 入信号X于极座标系统中的一输入相位θ。f(X, Θ)代表接收到该输入信号X的一机率函 数。随后，该输入信号χ的一相位误差估计值_:i根据下列方程式被产生：
[0008]
[0009] 根据本发明的另一具体实施例为一种应用于一通信系统的相位误差估计方法，用 以估id输入信号X的一相位误差。该输入信号X对应于一传送端利用一 Μ相-相位偏移 调制（M-PSK)所产生的一输出信号。数值Μ为二的任一正整数次方。首先，该相位误差估 计方法执行一判断步骤，判断该输入信号于极座标系统中的一输入信号半径r与一输入相 位Θ。根据该输入信号半径 r，一第一谐波系数匕(r)与一第二谐波系数h2(r)被产生。根 据该输入相位Θ，一第一正弦函数sin (Μ Θ )与一第二正弦函数sin (2M Θ )被产生。随后， 该第一谐波系数h (r)被乘以该第一正弦函数sin (Μ Θ )，产生一第一运算结果；该第二谐 波系数h2(r)被乘以该第二正弦函数sin(2ΜΘ)，产生一第二运算结果。随后，根据该第一 运算结果与该第二运算结果，该输入信号的一相位误差估计值?被产生。
[0010] 关于本发明的优点与精神可以藉由以下发明详述及附图得到进一步的了解。
【附图说明】
[0011] 图1 (A)和图1 (B)为根据本发明的一实施例中的通信系统的功能方块图。
[0012] 图2为根据本发明的一实施例中的相位误差估计方法的流程图。
[0013] 符号说明
[0014] 100 :通信系统 12 :接收模块
[0015] 12A:自动增益控制单元 14:相位误差估计模块
[0016] S21~S25:流程步骤
【具体实施方式】
[0017] 本说明书中的数学表示式用以说明与本发明的实施例相关的原理和逻辑，除非有 特别指明的情况，否则不对本发明的范畴构成限制。本发明所属技术领域中具有通常知识 者可理解，有多种技术可实现这些数学式所对应的物理表现形式。
[0018] 以下首先说明根据本发明的通信系统及相位误差估计方法可以使相位误差估计 结果满足克拉梅-劳下限（Cramer-Rao bound, CRB)〇
[0019] 假设一通信系统接收端接收输入信号x，该输入信号x在利用极座标系统（polar coordinate system)表示时具有一输入相位Θ。定义可观测量S(x)做为输入相位Θ的 无偏差估计值（unbiased estimator)，则可观测量S(x)的期望值为：
[0020] E [S (X) ] = J" S (X) · f (X，Θ ) dx = Θ， （式 1)
[0021] 其中f(x，Θ )代表输入信号X的一机率函数。
[0022] 该通信系统接收端产生的相位误差估计值与真实相位误差的估计差值ε (X)为：
[0023] ε (χ) = S(x)- Θ cat,(式 2)
[0024] 其中代表输入信号χ的真实相位误差。根据本发明的通信系统与相位误差 估计方法的主要目标即在于令估计差值ε (χ)的变异值（variance)具有一特定下限。且 基于相位无偏差的假设，估计差值ε (χ)的期望值为零：
[0025] Ε[ ε (x)] = / [S(x)~ Θ act] · f (χ, Θ )dx = 0〇 (式 3)
[0026] 在正规情况（regularity condition)下
I"先，对期望值E [ ε (χ)] 施以偏微分可得到下列关系式：
[0027]
[0028]
[0029] = / ε (χ) · ν ( θ , χ) · f (χ, Θ ) dx- / f (x, Θ ) dx = 0, (式 4)
[0030] 其中 i" f(x，Θ )dx 等于 1，所以 i" ε (x) · v( Θ，x) · f(x，Θ )dx = 1
[0031] 且利用数学公式
[0032]
[0033] 及
[0034] E [v ( θ，χ) ] = 〇 (式 6)
[0035] 可以得到估计差值ε (χ)与ν( θ，χ)的共变异数（covariance)等于1 :
[0036] Ε [ ε (χ) · v ( θ，χ) ] = E [S (χ) · v ( θ，χ) ] - Θ actE [v ( θ，χ) ] = E [S (χ) · v ( θ，χ)] =1〇
[0037] ......(式 7)
[0038] 接着，利用1彡StDev( ε ) XStDev(v)的数学公式，当估计差值ε (χ)与ν( θ，χ) 间的相关系数< 1，可得出符合克拉梅-劳下限（CRB)定义的关系式如下：
[0039] Ε [ ε (χ) · v ( θ，χ) ] 2= 1 < Ε [ ε 2 (χ) ] · Ε [V2 ( θ，χ)]，（式 8)
[0040]
[0041] 由前述公式得知，只要估计差值ε (χ)与ν( θ，χ)彼此相关（correlation)，便可 令上述符合克拉梅-劳下限（CRB)的条件成立。故本发明的通信系统以及相位误差估计方 法利用ν(θ，χ)来估测估计差值ε (X)，藉此提供能满足克拉梅-劳下限（CRB)的相位误 差估计结果。
[0042] 根据本发明的一具体实施例为一种通信系统，其功能方块图绘示于图1(A)。通信 系统100包含一接收模块12与一相位误差估计模块14。接收模块12用以接收输入信号 X，输入信号χ可转换成极座标系统以输入半径r及输入相位Θ表示。f(x，Θ)代表接收模 块12接收到输入信号χ的一机率函数。相位误差估计模块14用以根据下列方程式产生输 入信号χ的相位误差估计值Θ :
[0043]
[0044] 须说明的是，本发明的范畴并不限于以某种特定组态或架构来实现通信装置100， 且电路设计者可自行根据实际信号和传输环境的特性决定机率函数f(x，Θ)的设定。实务 上，在决定机率函数f(x，Θ)之后，式10可被展开为其它型态的运算式，进而被适当简化。 另一方面，负责产生相位误差估计值f的相位误差估计模块14可被实现为固定式及/或可 编程数位逻辑电路，包含可编程逻辑门阵列、特定应用集成电路、微控制器、微处理器、数位 信号处理器，与其他必要电路。或者，相位误差估计模块14亦可被设计为透过执行记忆体 (未绘示）中所储存的处理器指令，来完成多种任务。
[0045] 以下段落以采用Μ相-相位偏移调制（M-phase shift keying, M-PSK)的通信系 统（数值Μ为二的任一正整数次方，例如2、4、8、16~等数值）为例。若暂不考虑信号在频 率和时域的偏移，接收模块12收到的信号可被表示为：
[0046] x[k] = a[k]eJ0+w[k]，k = 0,1，· · ·，（L0_l)，（式 11)
[0049] 其中符号〇代表一标准差。[0050] 式12可被改写为：「00511
[0047] 其中L。代表时间轴上的取样点数，a[k]代表传送端提供的理想信号，r[k]代表接 收模块12收到的信号，而w[k]代表可加性白色高斯杂讯（AWGN)。信号x[k]在极座标系统 中的角度Θ为待决定的未知数。对应于信号 x[0]的条件机率函数f(x[0] |a[0] ; Θ)可被 定义如下·
[0048]
[0052] ......(式 13)
[0053] 其中ak[0]的定义如下：
[0056] 其中符号Τ'代表：
[0057]
[0058] 对式13施以自然对数运算的结果为：
[0059]
[0060] ......(式 16)
[0061] 式16中的Inf (χ[0] ; Θ )即为式1〇中的Inf (X，Θ )。因此，对式16施以偏微分 可得到相位误差估计值f :
[0062]
[0063]
[0064]
[0065] 实务上，可利用Matlab?等模拟软体简化实际电路的复杂度，将式17进一步以傅 利叶级数表示，并将相位误差估计值I的计算结果贡献较低的谐波项次省略。于一较佳实 施例中，仅将第一谐波项次和第二谐波项次纳入考虑，式17经化简后成为：
[0066]
[0067] 其中第一谐波系数比(r)、第二谐波系数h2 (r)系与输入信号X于极座标系统中的 信号输入半径r相关，而第一正弦函数sin (Μ Θ )、第二正弦函数sin (2M Θ )与输入信号X 于极座标系统中的输入相位Θ相关。易言之，在判定输入信号X于极座标系统中的信号半 径r和输入相位Θ之后，相位误差估计模块14可根据输入信号半径r产生第一谐波系数 比(r)和第二谐波系数1! 2&)、根据输入相位Θ产生一第一正弦函数sin (ΜΘ)与一第二正 弦函数sin (2M Θ )。接着，相位误差估计模块14可将第一谐波系数匕(r)乘以第一正弦函数 sin(M Θ )，以产生一第一运算结果、将第二谐波系数h2(r)乘以第二正弦函数sin(2M Θ )，以 产生一第二运算结果。根据第一运算结果和第二运算结果，相位误差估计模块14即可产生 输入信号X的相位误差估计值#。
[0068] 如图1(B)所示，接收模块12可进一步包含一自动增益控制单元12A。于相位误差 估计模块14产生第一谐波系数h (r)与第二谐波系数1!2&)前，自动增益控制（AGC)单元 12A负责将该输入信号X调整为大致具有一信号强度α 2,其中α代表一标准化系数。于一 实施例中，相位误差估计模块14于产生第一谐波系数匕(r)与第二谐波系数1! 2&)时亦考 量该标准化系数α，且系根据下列方程式产生第一谐波系数匕(r)和第二谐波系数1!2&):
[0069]
[0070]
[0071] 其中山~七。各自代表一个常数。举例而言，当数值Μ等于8 (亦即通信装置100 系用以配合采用8PSK的通信规范），且于式17中代入符号能量/杂讯密度比（Es/No)等于 5. 6分贝，则会得到：
[0072]
[0073]
[0074] 须说明的是，第一谐波系数匕(r)和第二谐波系数h2 (r)的实际数值会与式17中 的各参数的大小相关，亦与式17被化简的方式相关，不以特定数值为限。第一谐波系数 比(r)和第二谐波系数1! 2&)的实际数值会与式17中的各参数的大小相关，亦与式17被化 简的方式相关，不以特定数值为限。自动增益控制技术和标准化系数α的选择方式为本发 明所属技术领域中具有通常知识者所知，于此不赘述。
[0075] 根据本发明的另一具体实施例为一种应用于一通信系统的相位误差估计方法，用 以估计一输入信号X的一相位误差。首先，该相位误差估计方法执行一判断步骤，判断该输 入信号X于极座标系统中的一输入相位Θ。f(X, Θ)代表接收到该输入信号X的一机率函 数。随后，该输入信号X的一相位误差估计值#根据下列方程式被产生：
[0076]
[0077] 根据本发明的另一具体实施例为一种应用于一通信系统的相位误差估计方法，用 以估计一输入信号X的一相位误差，其流程图绘示于图2。该输入信号X是对应于一传送端 利用一 Μ相-相位偏移调制（M-PSK，Μ为二的任一正整数次方）所产生的一输出信号。首 先，在步骤S21中，该相位误差估计方法执行一判断步骤，判断该输入信号于极座标系统中 的一输入信号半径r与一输入相位Θ。步骤S22为根据该输入信号半径r产生一第一谐波 系数h (r)与一第二谐波系数h2(r)。步骤S23为根据该输入相位Θ产生一第一正弦函数 sin(M Θ )与一第二正弦函数sin(2M Θ )。步骤S24为将该第一谐波系数hjr)乘以该第一 正弦函数sin (Μ Θ )，产生一第一运算结果，并将该第二谐波系数h2 (r)乘以该第二正弦函数 sin(2M Θ )，产生一第二运算结果。随后，步骤S25为根据该第一运算结果与该第二运算结 果产生该输入信号的一相位误差估计值旮。
[0078] 本发明所属技术领域中具有通常知识者可理解，在图2中，某些步骤的顺序可被 等效调换，并且不会影响该相位误差估计方法的整体效果。此外，先前在介绍通信系统1〇〇 时描述的各种操作变化亦可应用至图2中的相位误差估计方法，其细节不再赘述。
[0079] 须说明的是，本发明的附图包含呈现多种彼此关联的功能性模块的功能方块图。 这些图式并非细部电路图，且其中的连接线仅用以表示信号流。功能性元件及/或程序间 的多种互动关系不一定要透过直接的电性连结始能达成。此外，个别元件的功能不一定要 如附图中绘示的方式分配，且分散式的区块不一定要以分散式的电子元件实现。
[0080] 藉由以上较佳具体实施例的详述，希望能更加清楚描述本发明的特征与精神，而 并非以上述所揭示的较佳具体实施例来对本发明的范畴加以限制。相反地，其目的是希望 能涵盖各种改变及具相等性的安排于本发明所欲申请的专利范围的范畴内。
【主权项】
1. 一种通信系统，包含： 一接收模块，用W接收一输入信号X，该输入信号X包含一输入相位0 ; W及 一相位误差估计模块，用W根据下列方程式产生一相位误差估计値每?其中f(x，0 )代表该输入信号X的一机率函数。2. -种通信系统，包含： 一接收模块，用W接收一输入信号X，该输入信号X包含一输入信号半径r与一输入相 位0，其中该输入信号X是由一 M相-相位偏移调制所产生；W及 一相位误差估计模块，用W : 根据该输入信号半径r产生一第一谐波系数hi (r)与一第二谐波系数h2(r); 根据该输入相位0产生一第一正弦函数Sin(Me)与一第二正弦函数Sin(2M0); 将该第一谐波系数hi (r)乘W该第一正弦函数Sin (M 0 )，W产生一第一运算结果，将该 第二谐波系数h2 (r)乘W该第二正弦函数Sin (2M 0 )，W产生一第二运算结果；W及 根据该第一运算结果与该第二运算结果产生该输入信号的一相位误差估计值 5 其中该相位误差估计值易冉^n/(斬巧前傅利叶级数相关，f(x，0)代表该接收模 块接收到该输入信号X的一机率函数。3. 如权利要求2所述的通信系统，其特征在于，该接收模块进一步包含： 一自动增益控制单元，用于调整该输入信号X的一信号强度a 2,其中a代表一标准化 系数； 其中该第一谐波系数hi(r)与该第二谐波系数h2(r)与该标准化系数a有关。4. 如权利要求3所述的通信系统，其特征在于，该相位误差估计模块根据下列方程式 产生该第一谐波系数hi(r):其中di~d 7各自代表一个常数。5. 如权利要求3所述的通信系统，其特征在于，该相位误差估计模块根据下列方程式 产生该第二谐波系数h,(r):其中屯~d 1。各自代表一个常数。6. 如权利要求3所述的通信系统，其特征在于，该相位误差估计模块根据下列方程式 产生该相位误差估计值#:7. -种应用于一通信系统的相位误差估计方法，该相位误差估计方法包含： (a)产生该输入信号X的一输入相位0 ; W及 化）根据下列方程式产生一相位误差估计值疾乎 其中f(x，0 )代表接收到该输入信号X的一机率函数。8. -种应用于一通信系统的相位误差估计方法，用W估计一输入信号X的一相位误 差，该输入信号X对应于一传送端利用一 M相-相位偏移调制所产生的一输出信号，数值M 为二的任一正整数次方，该相位误差估计方法包含： (a)接收一输入信号X，该输入信号X包含一输入信号半径r与一输入相位0，其中该 输入信号X由一 M相-相位偏移调制所产生； 化）根据该输入信号半径r产生一第一谐波系数hi (r)与一第二谐波系数h2(r); (C)根据该输入相位0产生一第一正弦函数Sin(Me)与一第二正弦函数Sin (2M0); (d) 将该第一谐波系数hi (r)乘W该第一正弦函数Sin (M0)，W产生一第一运算结果， 将该第二谐波系数h2 (r)乘W该第二正弦函数Sin (2M 0 )，W产生一第二运算结果；W及 (e) 根据该第一运算结果与该第二运算结果产生该输入信号的一相位误差估计值召;： 其中该相位误差估计值贫与為化的傅利叶级数相关，f(x，0)代表该输入信 号X的一机率函数。9. 如权利要求8所述的相位误差估计方法，其特征在于，在步骤化）前进一步包含： 调整该输入信号X的一信号强度a 2,其中a代表一标准化系数； 其中该标准化系数a用于产生该第一谐波系数hi (r)与该第二谐波系数h2(r)。10. 如权利要求9所述的相位误差估计方法，其特征在于，步骤（b)包含根据下列方程 式产生该第一谐波系数hi (r):其中di~d 7各自代表一个常数。11. 如权利要求9所述的相位误差估计方法，其特征在于，步骤化）包含根据下列方程 式产生该第二谐波系数h2(r):其中屯~d 1。各自代表一个常数。12. 如权利要求9所述的相位误差估计方法，其特征在于，步骤（e)包含根据下列方程 式产生该相位误差估计值易：
【文档编号】H04L27/22GK105991494SQ201510052465
【公开日】2016年10月5日
【申请日】2015年2月2日
【发明人】周裕仁
【申请人】晨星半导体股份有限公司