CN105119850A - 一种用于通信信号盲分离的快速不动点处理方法 - Google Patents
一种用于通信信号盲分离的快速不动点处理方法 Download PDFInfo
- Publication number
- CN105119850A CN105119850A CN201510377413.6A CN201510377413A CN105119850A CN 105119850 A CN105119850 A CN 105119850A CN 201510377413 A CN201510377413 A CN 201510377413A CN 105119850 A CN105119850 A CN 105119850A
- Authority
- CN
- China
- Prior art keywords
- separation
- prime
- sigma
- max
- matrix
- 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.)
- Pending
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明属于信号处理技术领域,尤其涉及一种高效的适用于通信信号盲分离的快速不动点处理方法。本发明针对现有的非圆快速不动点算法(nc-FastICA)不能适应选择最优非线性函数的缺陷,通过分析不同通信信号的特点,自适应地选择最优的非线性函数,提出了一种用于通信信号盲分离的快速不动点处理方法。本发明分离效率高,性能稳定,收敛速度快且在分离通信信号时较传统的nc-FastICA算法和ACMN算法有8~14dB左右的提升,可广泛用于通信信号盲分离。
Description
技术领域
本发明属于信号处理技术领域,尤其涉及一种高效的适用于通信信号盲分离的快速不动点处理方法。
背景技术
盲信号分离是信号处理中的一个重要问题,正是由于盲信号分离对源信号以及传输过程的先验知识并不要求,在语音、阵列、无线通信信号处理等领域有广阔的应用前景。独立分量分析(ICA)是解决盲信号分离的主要方法之一,其充分利用了源信号之间的独立性,通过寻找一组线性非正交变换,使得估计的源信号的统计独立性达到最大,以此恢复源信号。
ICA的模型为z=As,其中,观测信号 源信号 A是列满秩的混合矩阵,N表示信源的个数,M表示接收端传感器的数目。假设源信号相互独立、均值为零、并且功率是归一化的。ICA的目的就是找到一个分离矩阵W来恢复源信号,并且要恢复出来的源信号与真实的源信号只有幅度、相位或顺序的模糊。
分离效果的好坏可以用Amari指标来衡量,其定义为:
快速不动点算法(FastICA)作为独立分量分析(ICA)算法中的一种经典算法,其适用范围经历了从实数域到复数域以及非圆信号的过程。然而,传统的非圆快速不动点算法(nc-FastICA)在选择非线性函数的时候并没有考虑源信号的特殊形式,性能有很大的提升空间。自适应快速不动点算法(ACMN)虽然考虑到了信号的特殊形式并提出了自适应选择非线性函数的方法,但其性能只是在分离具有超高斯分布的广义高斯族信号时较传统的nc-FastICA有所提升,对于具有次高斯分布的通信信号并没有什么改善。迄今还没有人提出适用于通信信号盲分离的自适应快速不动点算法。所以,结合通信信号的特殊形式,提出一种高效能的、能自适应选择非线性函数的分离算法显得尤为重要。
发明内容
本发明为解决现有技术的不足,针对现有的非圆快速不动点算法(nc-FastICA)不能适应选择最优非线性函数的缺陷,通过分析不同通信信号的特点,自适应地选择最优的非线性函数,提出了一种用于通信信号盲分离的快速不动点处理方法。
一种用于通信信号盲分离的快速不动点处理方法,具体步骤如下:
S1、初始化参数,令n=1,i=1,nmax=50,ε=10-12,其中,
S2、白化接收数据,得到白化后的混合信号x,所述白化接收数据的具体过程如下:
S21、计算接收数据的自相关矩阵Rz=E{zzH},并做特征值分解Rz=UΛUH,其中,Λ=diag(λ1,λ2,…,λM),M表示接收端传感器的数目;
S22、计算白化矩阵其中,Λs=diag(λ1,λ2,…,λN),N表示信源的个数,Us为U的前N列组成的矩阵;
S23、根据S22所述白化矩阵V得到白化后的混合信号x=Vz;
S3、对混合矩阵进行预估计,具体步骤如下:
S31、置i←1;
S32、根据公式 对分离矩阵W的第i列wi进行更新,其中,G为非线性函数,G(z)=z2,G'为G的一阶导数,G”为G的二阶导数,n为更新次数;
S33、若i<N,则置i←i+1,返回步骤S32,若i≥N,则进入S4;
S4、对更新后的分离矩阵进行正交化处理:W←[WWH]-1/2W;
S5、对S4所述正交化处理的分离矩阵进行收敛判定,若收敛或者n=nmax,则输出预分离矩阵,根据公式 计算预分离后的信号并转到步骤S6,若不收敛或者n<nmax,则置n←n+1,返回步骤S3;
S6、初始化,令n=1,nmax=50,其中,nmax为最大迭代次数;
S7、对S5输出的预分离矩阵进行进一步地精确估计,输出最终的分离矩阵,具体为:
S71、置i←1;
S72、计算并自适应选择非线性函数;
S73、对S5输出的预分离矩阵的第i列wi按公式
S74、对S73所述更新后的分离矩阵的第i列进行如下的归一化处理wi=wi/||wi||2;
S75、判断S74所述更新后的分离矩阵的第i列是否收敛,如果收敛或者n=nmax,则输出分离矩阵的第i列,进入S76,如果不收敛且n<nmax,则置n←n+1返回S73;
S76、如果i<N,则置i←i+1,返回S72,否则输出最终的分离矩阵。
进一步地,S5所述收敛判定的方式为||W(n)HW(n+1)-IN||F<ε,其中,ε为常数,是经验值。
进一步地,S72所述非线性函数的选择方法为:如果ui>0.1,则G(z)=z20;如果ui≤0.1,则G(z)=z4。
进一步地,S75所述判断是否收敛的准则为
本发明的有益效果是:
本发明分离效率高,性能稳定,收敛速度快且在分离通信信号时较传统的nc-FastICA算法和ACMN算法有8~14dB左右的提升,可广泛用于通信信号盲分离。
附图说明
图1为本发明算法预分离部分流程图。
图2为本发明算法精确分离流程图。
图3为本发明算法性能随采样数变化曲线图(源信号为4PSK信号和16PSK信号)。
图4为本发明算法性能随采样数变化曲线图(源信号为4PAM信号和8PSK信号)。
具体实施方式
下面结合实施例和附图,详细说明本发明的技术方案。
如图1所示,一个主网络和认知网络共存的系统。
在第一个时隙内,CT传输数据给SR,同时CR2也收到该信息。为了尽可能的消除用户CR1和CR2之间的干扰,提高认知系统的吞吐量,假设若是CR2接收到CR1信号达到其可解码的SINR门限γth,则CR2可以对CT端传输给CR1的信号进行正确译码,正确译码的结果用于后面的干扰消除,否则不可以译码。
假设CT的最大发射功率为PC,max,CT以功率p1发送信号,则SR接收到的SINR表示为CT对PR的干扰功率为I1=gcpp1,CR2接收到的SINR为其中,gcr为CT与SR之间的链路增益,在后续内容中将节点之间的链路增益均表示为gij的形式,i和j表示为对应的节点,σ2为SR端的噪声功率,所有的噪声功率均归一化表示为σ2。
在第二时隙内,SR以自身的发射功率p2对接收到的CT的信号进行放大处理,并转发给CR1,则CR1接收到的SINR为
其中,为SR的功率放大增益,假设SR的最大发射功率为PR,max。CR1的吞吐量表示可以为
同时,CT以功率p3发送信息给CR2,由于SR发送信号给CR1,若在第一个时隙内CR2接收到的SINR满足则CR2无法正确译码并去掉关于CR1的干扰信号,CR2的SINR表示为
若在第一时隙内CR2接收到的SINR满足γ′2≥γth,则CR2可以利用模拟网络编码(ANC)技术将关于CR1的干扰信号β2gr2gcrp1去掉。此时CR2的SINR表示为
CR2的吞吐量可以表示为认知系统的总容量表示为C=c1+c2。
在第二时间片内CT和SR同时发送信号,共同对PR产生的干扰表示为
I2=p3gcp+p2grp(4)
功率优化问题可以写成:
通过化简,将(5)中关于p1、p2和p3的功率约束表示为如下形式:
在p2和p3固定的条件下令γ1对p1求导得到,
可知CR1的吞吐量C1是关于p1的增函数。
当CR2不可以解码CR1的信息时,即无法进行干扰消除时,得到CR2的吞吐量可以表示为如下形式:
此时CR2的吞吐量C2与p1无关。
同理,当CR2可以解码CR1的信息时,对CR2进行了干扰消除后,得到CR2的吞吐量可以表示为如下形式:
CR2的SINR对p1进行求导得到:
此时CR2的吞吐量C2是p1的增函数。认知系统总的吞吐量表达式为:
C=C1+C2(11)
综合可知总吞吐量表达式C(p1,p2,p3)始终是关于p1的增函数,并且在(5)的功率约束下p1的最优值是:
这里只是暂且认定p1的最优形式,如果后面的分析得到p2最优功率为0,从节省功率的角度考虑,重新令p1的最优功率也为0。
在确定了最优的p1之后,认知系统总吞吐量的表达式就变为关于p2和p3的函数,根据CR2的干扰的不同进一步讨论:
Case1:如果γ′2<γth,则CR2不能通过前一时隙的解码消除SR带来的CR1信号的干扰,总的吞吐量可以写为如下的表达式:
在假设p2固定的条件下,由(13)可知C′是关于p3的增函数,则可以得到此时最优的p3为也就是说任意给定一个p2,最优的p3形式就确定了。由于最优的p3表达式中包含变量p2,所以需要进行讨论得到最优的p2。对此,分情况讨论:
如果 那么最优的p3只能满足 的形式;如果 那么最优的p3需要从和这两种情况分类讨论后进行比较得到。由于在时的分析讨论求解最优的p2的过程与在时讨论的情况相同,为简略起见,只给出在下的分析过程。
当 即 时,最优的p3为 此时将关于p2的功率约束条件化简为:
令C′中log2()内的表达式γ′对p2进行求导得到:
下面通过详细的分析讨论得到最优的p2:
如果a1>0,且则C′是关于p2的增函数,有:
如果a1>0,且令求得其两个根表示为:
若θ2≤0,则C′在p2≥0范围内是关于p2的增函数,则有
若则C′在0≤p2≤θ2范围内是关于p2的减函数,在p2≥θ2范围内是增函数,有
若 则
若则C′在p2≤θ1范围内是关于p2的增函数,则有
若 则
若 则有:
如果a1<0,且则C′是关于p2的减函数,则
如果a1<0,且和前面的分析讨论过程类似,这里省略。
如果a1=0,令 得: 若a1=0,令 得:
若b1>0,则C′在p2≤θ3的范围内为减函数,在p2≥θ3范围内为增函数。
若θ3≤0,则
若 则
若 则
若b1<0,分析和b1>0时类似,这里省略。若b1=0,则若c1≥0,则 若c1<0,
在上面的分析过程中,当求解得到的时,从节省功率的角度考虑,令当时,由于CRBS在传输CR2的信息时对CR1不产生干扰,等效于CR1解码消除了CR2信息的干扰。得到的最优解,记为(p1,i,p2,i,p3,i)。
当 即 时,最优的p3为 将关于p2的功率约束条件化简为:
将带入(13)C′中,并对log2()内的表达式γ′对p2进行求导得到:
将得到的最优解记为(p1,ii,p2,ii,p3,ii)。
由于在时的分析讨论求解最优的p2的过程与在时讨论的情况类似,将最后得到的最优解记为(p1,iii,p2,iii,p3,iii)。
综合上述分析,全局的最优功率分配可以表示为:
Case2:如果γ′2≥γth,则CR2可以通过前一时隙解码消除SR带来的CR1信号的干扰,则有:
如果那么最优的p3需要从和这两种情况分类讨论得到。同样利用(Case1)中的分析方法进行分类讨论分析,得到:
当 即 时,最优的p3为 对于p2的功率约束如(14)形式,令C″中log2()内的表达式γ″对p2进行求导得到:
将得到的最优解记为(p1,iv,p2,iv,p3,iv)。
当 即 时,最优的p3为 对于p2的功率约束如(18)形式,令C″中log2()内的表达式γ″对p2进行求导得到:
将得到的最优解记为(p1,v,p2,v,p3,v)。
同样,如果那么最优的p3只能满足的形式,将最后得到的最优解记为(p1,vi,p2,vi,p3,vi)。
将上述两种解联合分析之后,全局的最优功率分配可以表示为:
本实施例采用MonteCarlo的仿真方法,将本发明中提到的一种中继辅助链路和直接链路协调传输的半双工认知中继系统中基于模拟网络编码干扰消除和没有干扰消除以及中继链路与直接链路独立传输情况下的功率优化方法进行比较。仿真参数为:所有的信道信息均服从均值为0,方差为1的复高斯分布;路径损耗指数为4;最大传输功率为10dB,噪声功率归一化为1;蒙特卡洛仿真次数为100000。系统中各节点初始位置如图1所示。
考虑以下三种方案的对比:
方案1:三个时隙独立传输
方案2:二个时隙未考虑干扰消除的联合传输
方案3:二个时隙考虑干扰消除的联合传输
图2显示了认知系统的和速率与不同干扰门限的变化关系,可以看出,在干扰门限数值是1的时候,方案3比方案2的性能高出了15%。所有方案的系统速率都随着可容忍的系统干扰门限的增大而增大,而不同方案的差异性也随之越来越明显的体现出来。
图3和图4分别比较了认知系统的和速率以及认知用户速率随不同中继放置位置的变化关系。从图中可以看出,所有方案的和速率随着中继距离基站的距离越来越远,所呈现的性能均是先提升后下降。在方案1中,CR2的速率与中继的位置无关,因此没有发生变化,系统的和速率有CR1决定。当中及距离基站端不是很远的时候,AF中继会出现最大的速率值,方案2也具有类似的趋势。然而,方案2中CR2的速率与方案1相比,有明显的提高。在方案3中,CR2处来自中继的干扰可以被消除掉,当中继很远的时候,干扰消除的效果也随之不明显,因此方案3的性能最终趋同于方案2。
Claims (4)
1.一种用于通信信号盲分离的快速不动点处理方法,其特征在于,包括如下步骤:
S1、初始化参数,令n=1,i=1,nmax=50,ε=10-12,其中,
S2、白化接收数据,得到白化后的混合信号x,所述白化接收数据的具体过程如下:
S21、计算接收数据的自相关矩阵Rz=E{zzH},并做特征值分解Rz=UΛUH,其中,Λ=diag(λ1,λ2,…,λM),M表示接收端传感器的数目;
S22、计算白化矩阵其中,Λs=diag(λ1,λ2,…,λN),N表示信源的个数,Us为U的前N列组成的矩阵;
S23、根据S22所述白化矩阵V得到白化后的混合信号x=Vz;
S3、对混合矩阵进行预估计,具体步骤如下:
S31、置i←1;
S32、根据公式 对分离矩阵W的第i列wi进行更新,其中,G为非线性函数,G(z)=z2,G'为G的一阶导数,G”为G的二阶导数,n为更新次数;
S33、若i<N,则置i←i+1,返回步骤S32,若i≥N,则进入S4;
S4、对更新后的分离矩阵进行正交化处理:W←[WWH]-1/2W;
S5、对S4所述正交化处理的分离矩阵进行收敛判定,若收敛或者n=nmax,则输出预分离矩阵,根据公式 计算预分离后的信号并转到步骤S6,若不收敛或者n<nmax,则置n←n+1,返回步骤S3;
S6、初始化,令n=1,nmax=50,其中,nmax为最大迭代次数;
S7、对S5输出的预分离矩阵进行进一步地精确估计,输出最终的分离矩阵,具体为:
S71、置i←1;
S72、计算并自适应选择非线性函数;
S73、对S5输出的预分离矩阵的第i列wi按公式
S74、对S73所述更新后的分离矩阵的第i列进行如下的归一化处理wi=wi/||wi||2;
S75、判断S74所述更新后的分离矩阵的第i列是否收敛,如果收敛或者n=nmax,则输出分离矩阵的第i列,进入S76,如果不收敛且n<nmax,则置n←n+1返回S73;
S76、如果i<N,则置i←i+1,返回S72,否则输出最终的分离矩阵。
2.根据权利要求1所述的.一种用于通信信号盲分离的快速不动点处理方法,其特征在于:S5所述收敛判定的方式为||W(n)HW(n+1)-IN||F<ε,其中,ε为常数,是经验值。
3.根据权利要求1所述的.一种用于通信信号盲分离的快速不动点处理方法,其特征在于:S72所述非线性函数的选择方法为:如果ui>0.1,则G(z)=z20;如果ui≤0.1,则G(z)=z4。
4.根据权利要求1所述的.一种用于通信信号盲分离的快速不动点处理方法,其特征在于:S75所述判断是否收敛的准则为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510377413.6A CN105119850A (zh) | 2015-06-30 | 2015-06-30 | 一种用于通信信号盲分离的快速不动点处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510377413.6A CN105119850A (zh) | 2015-06-30 | 2015-06-30 | 一种用于通信信号盲分离的快速不动点处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105119850A true CN105119850A (zh) | 2015-12-02 |
Family
ID=54667745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510377413.6A Pending CN105119850A (zh) | 2015-06-30 | 2015-06-30 | 一种用于通信信号盲分离的快速不动点处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105119850A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108390839A (zh) * | 2018-02-02 | 2018-08-10 | 东南大学 | 一种基于不动点的非线性信道均衡方法 |
CN110138478A (zh) * | 2019-05-30 | 2019-08-16 | 电子科技大学 | 一种用于非圆信号的多天线频谱感知方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103905352A (zh) * | 2012-12-27 | 2014-07-02 | 中国人民解放军理工大学 | 基于变异粒子群粒子滤波的单通道扰信盲分离方法 |
CN104485113A (zh) * | 2014-12-23 | 2015-04-01 | 长沙理工大学 | 一种多故障源声发射信号分离方法 |
-
2015
- 2015-06-30 CN CN201510377413.6A patent/CN105119850A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103905352A (zh) * | 2012-12-27 | 2014-07-02 | 中国人民解放军理工大学 | 基于变异粒子群粒子滤波的单通道扰信盲分离方法 |
CN104485113A (zh) * | 2014-12-23 | 2015-04-01 | 长沙理工大学 | 一种多故障源声发射信号分离方法 |
Non-Patent Citations (1)
Title |
---|
GUOBING QIAN,PING WEI,HONGSHULIAO: ""Efficient Variant of Noncircular Complex FastICA Algorithm for the Blind Source Separation of Digital Communication Signals"", 《SPRINGER SCIENCE+BUSINESS MEDIA NEW YORK 2015》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108390839A (zh) * | 2018-02-02 | 2018-08-10 | 东南大学 | 一种基于不动点的非线性信道均衡方法 |
CN110138478A (zh) * | 2019-05-30 | 2019-08-16 | 电子科技大学 | 一种用于非圆信号的多天线频谱感知方法 |
CN110138478B (zh) * | 2019-05-30 | 2021-03-30 | 电子科技大学 | 一种用于非圆信号的多天线频谱感知方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109347609B (zh) | 下行noma通信系统中基于动态swipt的协作传输方法 | |
CN105049383A (zh) | 一种fdd大规模mimo系统下行链路训练序列设计方法 | |
CN105827273A (zh) | 多小区大规模mimo系统用户双天线导频干扰消除方法 | |
CN105611633A (zh) | 基于swipt的波束赋形方法的接收机资源分配方法 | |
CN107318169A (zh) | 一种swipt‑idma系统中基于不完全信道状态信息的功率与时分因子联合分配方法 | |
CN104581780A (zh) | 一种基于预处理的分枝剪枝联合网络优化和波束成形方法 | |
CN103533613B (zh) | 一种过时信道信息下的放大转发中继选择方法 | |
CN103209051A (zh) | 一种协作多点联合传输系统在多用户场景下的两步预编码方法 | |
CN105450275A (zh) | 基于能效最优的多用户大规模天线中继系统天线选择方法 | |
CN105407535A (zh) | 一种基于约束马尔可夫决策过程的高能效资源优化方法 | |
CN103248461A (zh) | 基于波束成形的多小区干扰对齐迭代算法 | |
CN102547955B (zh) | 基于信漏噪比的分布式干扰消除方法 | |
CN105246158A (zh) | 基于高信噪比的能效最大化多天线中继系统功率分配方法 | |
CN112215335B (zh) | 一种基于深度学习的系统检测方法 | |
CN104754712B (zh) | 一种基于af认知中继协作传输的联合功率优化分配方法 | |
CN102347820A (zh) | 一种多小区协作无线通信系统联合编解码方法 | |
CN103209427B (zh) | 一种源用户基于用户信道质量的协作用户选择方法 | |
CN105119850A (zh) | 一种用于通信信号盲分离的快速不动点处理方法 | |
CN105636188A (zh) | 认知解码转发中继系统的功率分配方法 | |
CN103259585B (zh) | 基于收发机损耗的下行链路波束成形方法及其系统 | |
CN105162503A (zh) | 一种大规模mimo系统中多用户波束形成与天线选择的联合设计方法 | |
CN104506226A (zh) | 双层异构网络中基于协作毫微微蜂窝的干扰抑制预编码方法 | |
CN106850020A (zh) | 非理想信道状态下mimo系统中基于组合干扰对齐方法 | |
CN104954055B (zh) | 一种低复杂度的多用户信能同传系统能效优化方法 | |
CN110601736B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151202 |
|
WD01 | Invention patent application deemed withdrawn after publication |