瞿孟虹，何晓霜，游 凌

(西南电子电信技术研究所 成都 610041)

同频混合信号参数联合最大似然递归估计

瞿孟虹，何晓霜，游 凌

(西南电子电信技术研究所 成都 610041)

为提高成对载波多址(PCMA)模式下参数估计精度，简化实现流程，提出一种混合信号条件下频偏、相偏、信道联合估计算法。该算法基于最大似然准则，利用数据辅助技术，推导出联合参数估计显性表达式，并给出递归求解方法。对该算法估计性能进行的定量分析表明，估计性能逼近混合信号条件下克拉美罗界(Cramer-Rao bound)，能适应频率选择性信道。计算机仿真结果验证了该算法的有效性。

同频干扰; 克拉美罗界; 联合估计; 最大似然

PCMA技术因能成倍地提高卫星频谱的利用效率，受到了广泛关注。由于在本地存储有自干扰信号的符号序列，因此只需对调制参数进行估计后就可重构自干扰信号波形并从接收混合信号中去除。在PCMA专利[1]中，参数估计是通过通信双方在正式通信前分别独自发送的一段训练序列进行的，需要进行握手且增加了系统时间开销。后续很多学者采用基于数据辅助的方法直接从混合信号中估计本地信号参数初始值，然后用迟早门进行跟踪[2]，其中，定时估计利用传统的相关估计[3]，频偏估计则应用Kay[4]、Fitz[5]、L&R[6]等经典算法。在经典算法基础上，考虑信道影响，文献[7]提出了联合信道频偏估计方法。文献[8]针对本地信号构建最大似然目标函数，提出了幅度、时延、相位的联合最大似然估计并进行三维寻优。上述方法的处理过程本质上仍是单信号模型下的求解，并未考虑混合信号模型。文献[9]给出了混合信号条件下联合定时估计算法，算法仅对定时偏差进行了联合估计且需预知混合信号功率比，实用性受到限制。

本文针对PCMA混合信号模型，构建了混合信号条件下最大似然目标函数，推导了混合信号频偏、相偏、信道联合估计表达式。利用单纯形法求频偏值，通过递归迭代求得其他估计量。算法估计性能逼近混合信号模型下CRB(Cramer-Rao)界，能适应频率选择性信道，计算机仿真验证了本文算法的有效性。

1 信号模型

设接收信号的模型为：

式中，当s(t)是常规调制信号时，有α1≠0，α2=0；Ti是第i路信号符号周期；Δωi是残余载波频偏；ϕi是相偏；τi是定时偏移，假定0≤τi

对于同符号速率合作方PCMA通信，有T1=T2=T，g1(t)=g2(t)=g(t )，经符号速率采样，可得：

式中，N为采样符号数；L为信道记忆长度；v1=Δω1T/2π，v2=Δω2T/2π。对于训练序列{a1,n;−L+1≤n≤N −1}、{a2,n;−L+1≤n≤N −1}，式(2)可描述为矩阵形式：

2 算法推导

接收向量r的最大似然函数为：

令：h=ejθ1αg，h=ejθ2αg，取对数后，最大似然函数可描述为：

对h1求导并赋值为0，有：

可得：

将式(7)代入式(5)，似然函数可表示为:

式中，C(v)=Γ(v)A(AHA)−1AHΓH(v)，将式(8)对h2求导并赋值为0，有：

可得：

式中，X(v1,v2)=[I-C(v1)]Γ(v2)B ，用X标识。

3 算法实现流程

2) 对式(11)取绝对值并取对数作为目标函数f(x)，设置单纯形初始坐标值分别为x1=反射参数α=1，扩展参数γ=1.5，收缩参数β=0.5，误差界ε=0.001；

3) 求最大值点xh、次大值点xg和最小值点xl；满足f(xh)=max{f(xi)}，f(xg)=max{f(xi)| xi≠xh}，f(xl)=min{f (xi)}，i=1,2,3。除xh外的形心为

5) 若f(x扩展)

6) 若f(x压缩)≤f(xh′)，则xh=x压缩，f(xh)= f(x压缩)，转到步骤7)；否则进一步收缩运算，xi=xi+0.5(xl−xi),i =1,2,3，f(xi)转到步骤7)；

4 性能分析

联合最大似然估计条件下，频偏1v、2v估计均方误差CRB推导如下：

利用式(13)对式(5)各分量分别求二阶偏导，有：

式中，

利用分块矩阵求逆公式，当D、R均可逆时，式(15)的逆矩阵可表示为[11]：

其中，

h1、h2估计均值及方差推导如下，将C(v)定义式代入式(10)有：

由式(17)和式(18)可知，混合信号条件下频偏估计CRB与信噪比Es/N0、信号功率比SINR、两信号频偏差有关。取Es/N0=16 dB，两信号频偏差为0，频偏估计CRB随SINR变化关系如图1所示，其中单信号模型下频偏估计CRB见文献[7]。

取SINR=8 dB，Es/N0=16，频偏估计CRB随两信号频偏差变化关系如图2所示。

由图1知，在两信号功率相当时，频偏估计CRB相同，随着SINR增大，小信号估计均方误差增大。混合信号条件下大信号频偏估计CRB大于同功率单信号的值。由图2知，混合信号条件下频偏估计CRB受频偏差影响不明显。

5 仿真分析

仿真条件：两路信号均为QPSK调制，数字下变频后3倍采样，帧头数据长度均为64符号，进行500次Monte-Carlo试验。

仿真一：大信号功率远大于小信号，定时误差为1/3T，本文算法与经典算法频偏估计均方误差随大信号信噪比变化关系如图3所示。

仿真二：大信号信噪比为16 dB，两信号功率比SINR在0～16 dB变化时，无定时误差，本文算法与经典算法频偏估计均方误差随功率比变化关系如图4所示。其中，仿真条件①指符号长度为64，仿真条件②指符号长度为128；一路信号全部符号信息已知，另一路信号仅帧头已知，未知符号设置为全零。

仿真三：仿真条件如仿真二，本文算法信道估计均方误差随SINR变化如图5所示。

仿真四：两信号功率比SINR为8 dB，大信号信噪比0～16 dB变化时，本文算法信道估计均方误差随信噪比变化如图6所示。

由仿真一可知，本文算法在大信号功率远大于小信号时，对大信号频偏估计性能优良，不受信道特性影响。由仿真二可知，本文算法对大信号频偏估计精度不受信号功率比SINR及信道特性影响，性能逼近混合信号频偏估计CRB。在利用一路信号全部符号信息条件下，通过对另一路信号未知符号填零，本文算法能具有更优良的性能。由仿真三四可知，本文算法信道估计精度逼近CRB且估计精度不受SINR影响。

6 计算量分析

本文算法最大的运算量消耗在频率联合估计式(11)目标函数寻优中。考虑到A(AHA)−1AH、B(BHB)−1BH均可预先计算直接调用即可。C(v)求1取需2N2次复乘，X(v1,v2)求取需(L+1)N2次复乘，N+(N−1)NL 次复加，X(v1,v2)HX(v1,v2)需NL2次复乘，(N−1)L2次复加，[X(v1,v2)HX(v1,v2)]−1为Hermite矩阵，考虑到分块矩阵求逆性质[12]需L3次复乘及L3次复加。一次复乘需要4次实乘及2次实加，一次复加需要2次实加，则式(11)共需O(N2L)次乘法及O(N2L)次加法运算。

7 结 束 语

本文针对PCMA混合信号模型，构建了混合信号条件下最大似然目标函数，推导了混合信号频偏、相偏、信道联合估计表达式。算法估计性能逼近混合信号模型CRB下界。

[1] DANKBERG M D, MILLER M J, MULLIGAN M G. Self-interference cancellation for two-party relayed communication: US, 5596439[P]. 1997-1-21.

[2] 廖灿辉. 卫星同频混合信号的单通道接收处理与应用技术研究[D]. 北京: 清华大学, 2010. LIAO Can-hui. Single-channel processing technology on satellite co-frequency signals and its application[D]. Beijing: TsingHua University, 2010.

[3] OERDER M, MEYR H. Digital filter and square timing recovery[J]. IEEE Trans on Commun, 1988, 36: 605-612.

[4] KAY S M. A fast and accurate single frequency estimator[J]. IEEE Trans on Acoustic, Speech, and Signal Processing, 1989, 37 (12):1987-1989.

[5] FITZ M P. Planar filtered techniques for burst mode carrier synchronization[C]//IEEE GLOBECOM'91. Phoenix, AZ, USA: IEEE Press, 1991.

[6] LUISE M, REGGIANNINI R. Carrier frequency recovery in all-digital modems for burst-mode transmissions[J]. IEEE Trans on Commun, 1995, 43(3): 1169-1178.

[7] MORELLI M, MENGALI U. Carrier-frequency estimation for transmissions over selective channels[J]. IEEE Trans on Commun, 2000, 48(9): 1580-1589.

[8] 万坚, 李明, 朱中梁. 一种新型卫星多址方式的分离新算法[J]. 系统仿真学报, 2006, 18(7): 1799-1801. WAN Jian, LI Ming, ZHU Zhong-liang. A novel separation algorithm of a new multiple access mode in satellite[J]. Journal of System Simulation, 2006, 18(7): 1799-1801.

[9] 廖灿辉, 周世东, 朱中梁. 基于最大似然的同频混合信号联合定时估计算法[J]. 系统工程与电子技术，2010, 32(6): 1121-1124. LIAO Can-hui, ZHOU Shi-dong, ZHU Zhong-liang. ML-based joint timing estimation algorithm for co-frequency signals[J]. Systems Engineering and Electronics, 2010, 32(6): 1121-1124.

[10] KAY S M. Fundamental of statistical signal processing[M]. Englewood Cliffs, NJ, USA: Prentice-Hall, 1993.

[11] KAILATH T, SAYED A H, HASSIBI B. Linear estimation[M]. Englewood Cliffs, NJ, USA: Prentice-Hall, 2000.

[12] 高飞, 王永良, 陈辉, 等. STAP中的矩阵求逆问题研究[J]. 雷达科学与技术, 2008, 6(3): 215-218.GAO Fei, WANG Yong-liang, CHEN Hui, et al. Study on matrix inversion for STAP[J]. Radar Science and Techno1ogy, 2008, 6(3): 215-218.

编 辑 张 俊

Joint ML Iterative Estimation of Carrier Frequency, Phase, Channel for Co-Frequency Signals

QU Meng-hong, HE Xiao-shuang, and YOU Ling
(Research Institute of Electronic and Telecommunication Technology of Southwest Chengdu 610041)

In order to improve the accuracy and simplify the produces of parameters estimation in paired carrier multiple access (PCMA), this paper proposes an algorithm of joint iterative estimation of carrier frequency, phase, and channel. Basing on maximum likelihood (ML) criterion, the joint estimator using data-aided technique is derived and the iterative resolution is given. The quantitative analysis of the joint estimator characteristics indicates that the accuracy the joint estimator almost achieves the Cramer-Rao bound (CRB), and therefore the algorithm fits for frequency selective channels. Computer simulation results show the validity of the algorithm.

co-channel interference; Cramer-Rao bound; joint estimation; maximum likelihood

.自动化技术.

TN95

A

10.3969/j.issn.1001-0548.2015.03.004

2014 − 02 − 20；

2015 − 01 − 13

国家部级预研基金

瞿孟虹(1983 − )，男，博士生，主要从事数字信号处理、阵列信号处理等方面的研究.