CN105490705B - 基于ilsp算法的同步ds‑cdma信号盲解扩方法 - Google Patents
基于ilsp算法的同步ds‑cdma信号盲解扩方法 Download PDFInfo
- Publication number
- CN105490705B CN105490705B CN201610008345.0A CN201610008345A CN105490705B CN 105490705 B CN105490705 B CN 105490705B CN 201610008345 A CN201610008345 A CN 201610008345A CN 105490705 B CN105490705 B CN 105490705B
- Authority
- CN
- China
- Prior art keywords
- matrix
- signal
- code
- spreading
- sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7097—Interference-related aspects
- H04B1/7103—Interference-related aspects the interference being multiple access interference
- H04B1/7105—Joint detection techniques, e.g. linear detectors
- H04B1/71057—Joint detection techniques, e.g. linear detectors using maximum-likelihood sequence estimation [MLSE]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03057—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception with a recursive structure
- H04L25/0307—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception with a recursive structure using blind adaptation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明属于盲解扩技术领域,尤其涉及一种基于ILSP算法的同步DS‑CDMA信号盲解扩方法。本发明首先对接收到的同步短码DS‑CDMA信号的协方差矩阵进行特征分解,得到由各用户扩频码波形序列张成的信号子空间估计,然后根据扩频码的有限符号集特性,采用ILSP算法求得同步短码DS‑CDMA信号的扩频码序列矩阵估计,最后由估计得到的信号扩频序列矩阵解扩得到用户的信息码矩阵,完成盲解扩。本发明方法具有优良的盲解扩性能。
Description
技术领域
本发明属于盲解扩技术领域,尤其涉及一种基于ILSP算法的同步DS-CDMA信号盲解扩方法。
背景技术
直接序列扩频(Direct Sequence Spread Spectrum,DSSS)通信技术是现代通信技术中最常见技术之一,由于其具有抗干扰、低截获率等优点,被广泛的应用于民用通信、军用通信、卫星导航等领域,具有很高的研究价值。DSSS通信技术是在信号发送时,利用高速率的扩频序列与信息码序列相乘,使得信号的频谱被扩宽,谱密度降低,在接收方,再用与发送时相同的高速序列与接收到的信号相乘,使信号的频谱恢复完成解扩,而噪声和干扰的频谱却被扩展了,谱密度降低,这样就很容易将噪声和干扰从信号中滤除。
对于合作接收方,可以用已知的扩频序列进行DSSS信号解扩;而对于非合作方,由于扩频序列未知,因此对信号的解扩需要事先利用截获信号样本估计出扩频序列。对于直接序列码分多址(Direct Sequence—Code Division Multiple Access,DS-CDMA)信号的盲解扩问题,目前较为有效的方法为基于期望最大化(Expectation Maximization,EM)算法的盲解扩方法,该方法首先对样本协方差矩阵进行特征分解,得到由各用户扩频码波形序列张成的信号子空间估计,并将观测矩阵投影到信号子空间中。然后,再利用信息码元序列的有限符号集特性,将投影得到的信号样本建模为高斯混合分布模型,通过EM算法得到扩频矩阵的极大似然估计,进一步得到信号幅度、扩频码波形序列及信息码元序列的估计,完成盲解扩。基于EM算法的盲解扩方法主要包括样本协方差矩阵的计算与特征分解及高斯混合分布模型的极大似然估计两部分,而高斯混合分布模型的极大似然估计部分的计算复杂度会随用户个数的增加剧烈增加,这对于解决多用户问题是十分不利的。
发明内容
本发明的目的在于克服现有技术的不足,提出了基于ILSP算法的同步DS-CDMA信号盲解扩方法,该低信噪比条件下具有优良的盲解扩性能。
本发明的技术方案为:
首先对接收到的同步短码DS-CDMA信号的协方差矩阵进行特征分解,得到由各用户扩频码波形序列张成的信号子空间估计,然后根据扩频码的有限符号集特性,采用ILSP算法求得同步短码DS-CDMA信号的扩频码序列矩阵估计,最后由估计得到的信号扩频序列矩阵解扩得到用户的信息码矩阵,完成盲解扩。
基于ILSP算法的同步DS-CDMA信号盲解扩方法,包括如下步骤:
S1、非合作接收方截获到的DS-CDMA信号经码元速率采样后为:
其中,R为用户个数,r=1,2,3,...,R,Ar为第r个用户的信号幅度,M为信息码个数,br为第r个用户的采用BPSK调制的信息码序列,cr为第r个用户的扩频序列,L为扩频因子,v(n)为方差为σ2的高斯白噪声;
S2、求取截获的同步短码DS-CDAM信号的扩频码序列张成的信号子空间矩阵,具体为:
S21、同步短码DS-CDMA信号建模成矩阵模型Y=CABT+V,其中,L×R维矩阵C为扩频码矩阵,A=diag([A1,A2,...,AR])为信号幅度对角阵,M×R维矩阵B为BPSK调制的信息码矩阵,V为方差为σ2的高斯白噪声矩阵;
S22、由信号矩阵模型求得信号样本的协方差矩阵的估计和扩频序列张成的信号子空间估计矩阵其中,所述为信号协方差矩阵特征分解,最大的R个特征值对应的特征向量所组成的矩阵;
S3、由S2所得信号子空间估计矩阵根据扩频序列的有限符号集特性,采用ILSP算法来求信息码矩阵C的估计矩阵其中,R×R维矩阵为模糊矩阵,具体实现步骤为:
S31、令i=0,随机初始化R×R维矩阵
S32、令i=i+1,计算计算其中,sign()为取符号运算;
S33、重复S32直到方法收敛或达到最大迭代次数A,其中,所述方法收敛的条件为ε为收敛的门限值,最大迭代次数A为经验值;
S34、重复S31-S33H次,选择迭代效果最好的一次作为最终结果;
S4、由S3得到扩频码估计矩阵对截获信号解扩得到用户信息码矩阵估计,完成盲解扩。
进一步地,S33所述大迭代次数A为50次,ε=1×10-9。
进一步地,S34所述H=40。
本发明的有益效果是:
本发明首先对截获同步短码DS-CDMA信号的样本协方差矩阵进行特征分解,得到由各用户扩频码序列张成的信号子空间估计,然后根据扩频码的有限符号集特性,采用ILSP算法获得扩频码序列矩阵的估计值,最后利用扩频序列矩阵估计对截获信号样本进行解扩,得到各用户的信息码估计。计算机仿真表明本发明在低信噪比条件下具有优良的盲解扩性能。
附图说明
图1是本发明基于ILSP算法同步DS-CDMA信号盲解扩方法的具体实施方式流程图。
图2是本发明具体实施1,不同用户个数信号误码率随信噪比变化曲线示意图。
图3是本发明具体实施2,不同用户个数信号误码率随信息码个数变化曲线示意图。
具体实施方式
下面结合实施例和附图,详细说明本发明的技术方案。
具体实施1、
对不同用户个数的信息码矩阵误码率随SNR(Signal Noise Rate,信噪比)变化进行仿真,同时与相同情况下的合作通信进行对比。
如图1所示,本发明实现同步DS-CDMA信号盲解扩方法包括以下步骤:
S1、非合作接收方截获到的DS-CDMA信号经码元速率采样后为:
其中,R为用户个数,r=1,2,3,...,R,Ar为第r个用户的信号幅度,M为信息码个数,br为第r个用户的采用BPSK调制的信息码序列,cr为第r个用户的扩频序列,L为扩频因子,v(n)为方差为σ2的高斯白噪声;
S2、求取截获的同步短码DS-CDAM信号的扩频码序列张成的信号子空间矩阵,具体为:
S21、同步短码DS-CDMA信号建模成矩阵模型Y=CABT+V,其中,L×R维矩阵C为扩频码矩阵,A=diag([A1,A2,...,AR])为信号幅度对角阵,M×R维矩阵B为BPSK调制的信息码矩阵,V为方差为σ2的高斯白噪声矩阵;
S22、由信号矩阵模型求得信号样本的协方差矩阵的估计和扩频序列张成的信号子空间估计矩阵其中,所述为信号协方差矩阵特征分解,最大的R个特征值对应的特征向量所组成的矩阵;
S3、由S2所得信号子空间估计矩阵根据扩频序列的有限符号集特性,采用ILSP算法来求信息码矩阵C的估计矩阵其中,R×R维矩阵为模糊矩阵,ILSP算法每次迭代过程将前次估计值作为新的零点,直到最后收敛或达到最大迭代次数为止,具体实现步骤为:
S31、令i=0,随机初始化R×R维矩阵
S32、令i=i+1,计算计算其中,sign()为取符号运算;
S33、重复S32直到方法收敛或达到最大迭代次数A,其中,所述方法收敛的条件为ε为收敛的门限值,最大迭代次数A=50,ε=1×10-9;
S34、为了提高方法性能,重复S31-S33H次,选择迭代效果最好的一次作为最终结果,其中,H=40;
S4、由S3得到扩频码估计矩阵对截获信号解扩得到用户信息码矩阵估计,完成盲解扩。
进行1000次蒙特卡洛实验,最终得到的信息码矩阵误码率随SNR变化曲线图,并与合作情况(扩频码矩阵已知信息码矩阵可由B=sign(CHY)得到)对比,如图2所示。从图中可以看出本发明提出的盲解扩方法在低信噪比情况下性能良好,本方法的误码率随信噪比的增加逐渐逼近合作通信,证明本方法的有效性。
具体实施2、
对不同用户个数的信息码矩阵误码率随信息码个数变化进行仿真,并与合作情况(扩频码矩阵已知信息码矩阵可由B=sign(CHY)得到)对比。本次实施中信噪比固定为-10dB,信息码数量变化范围为10到390个,其他条件与实施1相同。通过1000次蒙特卡罗仿真实验得到信号误码率随信息码个数变化曲线图如图3所示。从图中可以看出随着信息码个数的增加,本方法的性能逐渐逼近合作通信,从而证明本方法的有效性。
Claims (3)
1.基于ILSP算法的同步DS-CDMA信号盲解扩方法,其特征在于,包括如下步骤:
S1、非合作接收方截获到的DS-CDMA信号经码元速率采样后为:
其中,r=1,2,3,...,R,Ar为第r个用户的信号幅度,M为信息码个数,br为第r个用户的采用BPSK调制的信息码序列,cr为第r个用户的扩频序列,L为扩频因子,v(n)为方差为σ2的高斯白噪声;
S2、求取截获的同步短码DS-CDAM信号的扩频码序列张成的信号子空间矩阵,具体为:
S21、同步短码DS-CDMA信号建模成矩阵模型Y=CABT+V,其中,L×R维矩阵C为扩频码矩阵,A=diag([A1,A2,...,AR])为信号幅度对角阵,M×R维矩阵B为BPSK调制的信息码矩阵,V为方差为σ2的高斯白噪声矩阵;
S22、由信号矩阵模型求得信号样本的协方差矩阵的估计和扩频序列张成的信号子空间估计矩阵其中,所述为信号协方差矩阵特征分解,最大的R个特征值对应的特征向量所组成的矩阵;
S3、由S2所得信号子空间估计矩阵根据扩频序列的有限符号集特性,采用ILSP算法来求信息码矩阵C的估计矩阵其中,R×R维矩阵为模糊矩阵,具体实现步骤为:
S31、令i=0,随机初始化R×R维矩阵
S32、令i=i+1,计算计算其中,sign()为取符号运算;
S33、重复S32直到方法收敛或达到最大迭代次数A,其中,所述方法收敛的条件为ε为收敛的门限值,最大迭代次数A为经验值;
S34、重复S31-S33H次,选择迭代效果最好的一次作为最终结果;
S4、由S3得到扩频码估计矩阵对截获信号解扩得到用户信息码矩阵估计,完成盲解扩。
2.根据权利要求1所述的基于ILSP算法的同步DS-CDMA信号盲解扩方法,其特征在于:S33所述最大迭代次数A为50次,ε=1×10-9。
3.根据权利要求1所述的基于ILSP算法的同步DS-CDMA信号盲解扩方法,其特征 在于:S34所述H=40。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008345.0A CN105490705B (zh) | 2016-01-07 | 2016-01-07 | 基于ilsp算法的同步ds‑cdma信号盲解扩方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008345.0A CN105490705B (zh) | 2016-01-07 | 2016-01-07 | 基于ilsp算法的同步ds‑cdma信号盲解扩方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105490705A CN105490705A (zh) | 2016-04-13 |
CN105490705B true CN105490705B (zh) | 2018-04-13 |
Family
ID=55677471
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610008345.0A Active CN105490705B (zh) | 2016-01-07 | 2016-01-07 | 基于ilsp算法的同步ds‑cdma信号盲解扩方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105490705B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107294560B (zh) * | 2017-07-03 | 2019-08-02 | 电子科技大学 | 异步短码ds-cdma信号盲解扩方法 |
CN107342786B (zh) * | 2017-07-03 | 2019-06-04 | 电子科技大学 | 基于ilsp算法的信号盲解扩方法 |
CN109150307B (zh) * | 2018-07-17 | 2020-07-31 | 北京邮电大学 | 一种可见光通信中的通信数据的获取方法、装置和设备 |
CN110808760A (zh) * | 2019-10-25 | 2020-02-18 | 重庆邮电大学 | 基于ilsp-cma的同步ds-cdma信号伪码序列和信息序列联合盲估计 |
CN116232809B (zh) * | 2023-01-12 | 2024-04-19 | 电子科技大学 | 一种同步长码ds-cdma信号盲估计方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005078949A1 (ja) * | 2004-02-12 | 2005-08-25 | Nec Corporation | 移動通信システム及びそれに用いる無線装置 |
CN101447819A (zh) * | 2009-01-04 | 2009-06-03 | 北京交通大学 | Mimo-cdma基站空间复用接收机系统 |
CN103414670A (zh) * | 2013-08-21 | 2013-11-27 | 电子科技大学 | 一种基于半定规划的长码dsss信号盲解扩方法 |
CN104917544B (zh) * | 2014-03-13 | 2017-07-14 | 大唐电信科技产业控股有限公司 | 一种扰码选择方法及装置 |
-
2016
- 2016-01-07 CN CN201610008345.0A patent/CN105490705B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105490705A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105490705B (zh) | 基于ilsp算法的同步ds‑cdma信号盲解扩方法 | |
CN105939188B (zh) | 一种underlay频谱共享下时频重叠MASK信号的载频估计方法 | |
CN105429673B (zh) | 基于张量低秩分解的长码ds‑cdma信号盲分离方法 | |
CN105634544B (zh) | 一种同步长码ds-cdma信号的盲解扩方法 | |
CN105680903B (zh) | 周期长短码直扩码分多址信号多伪码估计方法 | |
CN105791184A (zh) | 基于相似度的非周期长码直扩信号扩频序列盲估计方法 | |
CN103297160A (zh) | 基于归一化特征值的拟合优度检验的频谱感知方法及装置 | |
CN105871413A (zh) | 低信噪比直接序列扩频信号检测方法 | |
CN111953380A (zh) | 基于范数拟合的非周期长码直扩信号时延估计方法及系统 | |
Wei et al. | Blind estimation of PN sequence of DS-CDMA signal in multipath | |
CN108111189A (zh) | 基于Hebb规则优化的扩频码序列识别与同步方法 | |
CN104539312B (zh) | 基于序贯蒙特卡罗的扩频码及信息序列联合估计方法 | |
CN104022794B (zh) | 直接序列扩频系统脉冲干扰抑制装置及方法 | |
CN114070353B (zh) | 一种同步长码ds-cdma信号盲解扩方法 | |
CN109951406A (zh) | 联通异步wcdma实采信号用户数和扰码盲估计方法 | |
CN104506466B (zh) | 一种多载波cdma信号扩频码及信息序列盲估计方法 | |
CN107294560B (zh) | 异步短码ds-cdma信号盲解扩方法 | |
CN110808760A (zh) | 基于ilsp-cma的同步ds-cdma信号伪码序列和信息序列联合盲估计 | |
CN103986559B (zh) | 一种对压缩采样信号的五阶循环累积量估计算法 | |
CN105227504B (zh) | 一种mb-ofdm uwb系统的信道估计方法及系统 | |
CN107342786B (zh) | 基于ilsp算法的信号盲解扩方法 | |
CN107749771B (zh) | 一种基于似然反馈的并行干扰消除算法 | |
CN107733820A (zh) | 多径信道环境条件下的ds‑cdma信号盲估计方法 | |
Wu et al. | Blind channel estimation based on subspace for multicarrier CDMA | |
CN105530696B (zh) | 一种多用户干扰消除方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |