一种短码直扩信号伪码序列盲估计方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN109150776A(43)申请公布日2019.01.04(21)申请号201810928344.7(22)申请日2018.08.15(71)申请人桂林电子科技大学地址541004广西壮族自治区桂林市七星区金鸡路1号(72)发明人纪元法范灼孙希延符强王守华严素清付文涛(74)专利代理机构北京中济纬天专利代理有限公司11429代理人石燕妮(51)Int.Cl.H04L25/02(2006.01)H04B1/707(2011.01)权利要求书2页说明书10页附图5页(54)发明名称一种短码直扩信号伪码序列盲估计方法(57)摘要本发明涉及一种短码直扩信号伪码序列盲估计方法，解决的是传统SVD算法在对实际无线信道下传输的直扩信号进行伪码序列盲估计时不能满足系统估计性能的要求、抗冲击噪声特性差的技术问题，通过采用方法包括：步骤一，按照两倍伪码周期长度的大小对接收信号进行分段；步骤二，构造接收信号的分数低阶观测矩阵；步骤三，将步骤二的分数低阶观测矩阵使用SVD分解算法求取最大奇异左向量；步骤四，使用改进失步点估计方法，由步骤三的最大奇异左向量估计伪码失步点位置，完成冲击噪声信道下对直扩信号的伪码序列盲估计的技术方案，较好的解决了该问题，可用于伪码序列盲估计中。CN109150776ACN109150776A权利要求书1/2页1.一种短码直扩信号伪码序列盲估计方法，其特征在于：所述方法包括：步骤一，按照两倍伪码周期长度的大小对接收信号进行分段；步骤二，构造接收信号的分数低阶观测矩阵；步骤三，将步骤二的分数低阶观测矩阵使用SVD分解算法求取最大奇异左向量；步骤四，使用改进失步点估计方法，由步骤三的最大奇异左向量估计伪码失步点位置，完成冲击噪声信道下对直扩信号的伪码序列盲估计。2.根据权利要求1所述的短码直扩信号伪码序列盲估计方法，其特征在于：所述步骤2采用分数低阶联合M估计的方法构造接收信号观测矩阵。3.根据权利要求2所述的短码直扩信号伪码序列盲估计方法，其特征在于：所述步骤二包括：步骤2.1，采用α稳定分布模型对冲击噪声进行建模，α稳定分布的特征函数为其中，α∈(0,2]是特征因子，β∈[-1,1]为对称参数，γ＞0为分散系数，-∞＜u＜∞为位置系数；步骤2.2，假设q段分段信号两两不相关，且统计独立，信号成分和噪声成分相互独立且均服从相同特征指数大小，及位置系数为零的SαS分布，计算出t时刻的接收信号矩阵k＝1,2,...,q代表接收信号的各分段数，ak(t)代表数据信息，sk(t)代表伪码序列信息，n(t)表示接收信号噪声成分；步骤2.3，假设w(t)＝ak(t)sk(t)，r(t)＝w(t)+n(t)，w(t)与n(t)相互独立，计算出接收信号的共变矩阵：[ri(t),rj(t)]α＝[wi(t)+ni(t),wj(t)+nj(t)]α＝[wi(t),wj(t)]α+[wi(t),nj(t)]α+[ni(t),wj(t)]α+[ni(t),nj(t)]α步骤2.4，计算出噪声分量的共变矩阵[ni(t),nj(t)]α＝γnδi,j；其中，γn＝[nk(t),nk(t)]α，δi,j为Kronecker函数；步骤2.5，根据噪声分量的共变矩阵，计算出接收信号的共变矩阵步骤2.6，计算出观测向量rk的共变矩阵＜α-1＞ΓR＝[r(t),r(t)]α＝AΓsA+γnI其中，ΓS＝diag(γs1,...,γsq)步骤2.7，计算出信号协方差矩阵p-2Ri,j＝E{ri(t)rj(t)|rj(t)|sign(ri(t)rj(t))},2CN109150776A权利要求书2/2页1＜p＜α；其中，E(·)表示期望。4.根据权利要求2或3所述的短码直扩信号伪码序列盲估计方法，其特征在于：步骤二还包括：步骤2.8，根据Huber类M估计加权函数，计算由分数低阶联合M估计构造的接收信号观测矩阵：p-2Ri,j＝E(ri(t)rj(t)|rj(t)|sign(ri(t)rj(t))p-2f(ri(t)rj(t)|rj(t)|sign(ri(t)rj(t))))，1＜p＜α；其中，Huber类M估计加权函数为5.根据权利要求1所述的短码直扩信号伪码序列盲估计方法，其特征在于：所述步骤四包括：步骤4.1，将SVD算法提取出来的包含完整伪码信息的最大奇异左向量r′按照两倍伪码周期长度进行移位分段，得到r′i＝r′(i:i+N)，i＝1,2,...,N；其中，i表示移位数，N为伪码信息的长度，r′(i:i+N)表示从信号r′移位值i处向后截取伪码长度N组成的信号成分；步骤4.2，计算元素补集，r′j＝[r′i(1:i-1),r′i(i+1:2N)]其中，[r′i(1:i-1),r′i