CN105262506B - 一种时频重叠高斯调幅通信信号分离方法 - Google Patents
一种时频重叠高斯调幅通信信号分离方法 Download PDFInfo
- Publication number
- CN105262506B CN105262506B CN201510570687.7A CN201510570687A CN105262506B CN 105262506 B CN105262506 B CN 105262506B CN 201510570687 A CN201510570687 A CN 201510570687A CN 105262506 B CN105262506 B CN 105262506B
- Authority
- CN
- China
- Prior art keywords
- value
- signal
- variable parameter
- amplitude modulation
- multidimensional variable
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000004891 communication Methods 0.000 title claims abstract description 24
- 238000000926 separation method Methods 0.000 claims abstract description 21
- 230000002068 genetic effect Effects 0.000 claims abstract description 9
- 230000006835 compression Effects 0.000 claims description 12
- 238000007906 compression Methods 0.000 claims description 12
- 230000035772 mutation Effects 0.000 claims description 6
- 230000010076 replication Effects 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 8
- 239000002245 particle Substances 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000000739 chaotic effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000010845 search algorithm Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000010353 genetic engineering Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明提供一种时频重叠高斯调幅通信信号分离方法，包括以下步骤：步骤一，建立信号模型，将接收到的混合信号转换成由若干个高斯调幅源信号组成，将混合信号的求解过程转化成求解多维变量参数的过程；步骤二，采用遗传算法计算出多维变量参数的初始值；步骤三，采用最小值搜索法计算出多维变量参数的最优值；步骤四，根据多维变量参数的最优值计算出各个源信号。能有效地实现单通道时频重叠高斯调幅通信信号的分离，该方法是先建立高斯调幅通信信号模型，求出所需要估计参数的初始值，利用该方法搜索到各个所需估计参数的最优值，根据得到的最优值恢复出各个高斯调幅源信号。
Description
技术领域
本发明涉及超宽带无线通信信号共存于信道中的信号处理技术领域，具体涉及一种单通道时频重叠高斯调幅通信信号分离的方法。
背景技术
超宽带(ultrawideband,UWB)通信是无线通信领域的一大突破,它基于高带宽实现高速率数据传输。相比其他无线通信系统,它具有高空间频谱效率、高测距精度、低截获概率、高抗多径衰落能力、低功耗、低成本和小体积等诸多优点和潜力。在实际的超宽带通信信号接收中，伴随着无线频谱资源的过度使用，信号会更为密集，多个信号分量往往以时频域重叠的方式存在于共信道中，而大多数情况下都是单天线信号接收的，这给后续信号处理带来了极大的困难，如果需要提取某个信号或者研究单个信号的特性，就需要对这类单通道时频重叠信号进行盲分离。
从国内外大量的文献可知,多项式拟合方法[Barbarossa S,Scaglione A,Giannakis G B.Product high-order ambiguity function for multicomponentpolynomial-phase signal modeling,IEEE Trans Signal Process,vol.46,Issue 3,1998,p.691-708]、利用能量算子计算各个信号的幅度和相位的方法[Cai X W,Wei P,XiaoX C.The single channel of time-frequency overlapping signal blind sourceseparation method based on energy operator,China science letter E:informationscience,2008,38:607-619]、利用各种不同的时频分布估计参数方法[Li M Z,Zhao HC.Research on parameters extraction of pseudo code phase modulation-carryfrequency modulation combined fuse signal based on the adaptive window lengthof improved B distribution,Acta Armamentarii,2011,32:543-547]和各种虚拟通道法等都可以对单通道时频重叠高斯调幅通信信号进行分离。近年来常用的方法是虚拟通道法，包括wavelet-ICA和EMD-ICA方法。Hong Hoonbin等对接收到的混合信号先进行小波分解，然后对于得到的分量进行ICA处理，将此方法用于机械故障诊断中去，虽然这种方法不受信号类型的限制，对于单频率的信号进行分离的效果较好，但是对于带宽信号不能彻底分离，仍保留了混合信号的信息，不能很好地分离出各个源信号，分离效果不是很理想[Hong H,Liang M.Separation of fault features from a single-channel mechanicalsignal mixture using wavelet decomposition,Mechanical Systems and SignalProcessing,2007,21:2025-2040]。毋文峰等提出利用总体经验模态分解EMD将单路混合信号分解为多个本征模态函数，将本征模态函数分量进行ICA分离恢复,将经验模态分解EMD-ICA的单通道盲源分离方法用于轴承和齿轮的仿真研究，正确分离出轴承和齿轮源信号，此法相对于时空法、小波分解法等分离效果明显提高，然而，对于频谱重叠过多的信号，每个本征模函数分量覆盖了一个宽的频率范围，致使分离效果不佳[Wu W F,Chen X H，Su XJ.Blind Source Separation of Single-channel Mechanical Signal Based onEmpirical Mode Decomposition,Chinese Journal of Mechanical Engineering,2011,47:12-16]。在实际中，对于特定模型的信号，可以采用模型估计的方法进行盲分离，即建立信号的模型，利用一些估计方法将模型中信号的参数的估计出来。王世元等采用参数估计的方法解决了原始混沌信号的盲分离问题，采用一种容积准则近似该映射的加权积分函数，基于由状态空间模型建模的参数，提出了一种新的参数估计，有效实现了混沌信号的重构。由于一定要建立合适的状态空间模型的缺点，导致该方法不适用于单通道时频重叠高斯调幅信号的盲分离[Wang S Y,Feng J C.A novel method of estimating parameterand its application to blind separation of chaotic signals,Acta Phys.Sin.,2012,61:170508]。张淑宁等提出了基于粒子滤波的单通道正弦调频混合信号的分离和参数估计,针对正弦调频混合信号频率无跳变的特征,提出了一种基于粒子滤波的相位差解混叠算法,并通过源信号相位差解决了算法中粒子滤波高维状态空间降维问题,提出了一种适合高维状态空间的似然函数模型,比较固定长度粒子估计值和真实值误差,进而准确衡量粒子权重.通过在重采样后引入MCMC转移,解决了静止参数下粒子多样性降低问题,有效提高粒子滤波迭代收敛速度,完成对单通道正弦调频混合信号的参数提取,并通过重构信号完成正弦调频混合信号分离[Zhang S N,Zhao H C,Xiong G.Separation andparameter estimation of single channel sinusoidal frequency modulated signalmixture sources based on particle filtering.Acta Phys.Sin.，2014,63:158401]。
发明内容
本发明为了解决现有技术中存在的上述缺陷和不足，提供了一种时频重叠高斯调幅通信信号分离方法，能有效地实现单通道时频重叠高斯调幅通信信号的分离，该方法是先建立高斯调幅通信信号模型，求出所需要估计参数的初始值，利用该方法搜索到各个所需估计参数的最优值，根据得到的最优值恢复出各个高斯调幅源信号。
为解决上述技术问题，本发明提供一种时频重叠高斯调幅通信信号分离方法，包括以下步骤：
步骤一，建立信号模型，将接收到的混合信号转换成由若干个高斯调幅源信号组成，将混合信号的求解过程转化成求解多维变量参数的过程；
步骤二，采用遗传算法计算出多维变量参数的初始值；
步骤三，采用最小值搜索法计算出多维变量参数的最优值；
步骤四，根据多维变量参数的最优值计算出各个源信号。
所述步骤一中，将接收到的混合信号x(t)转换成由M个高斯调幅源信号组成，表达式如下：
其中，表示各个源信号分量的幅度调制系数，t0表示各个源信号分量的时域中心，ωi表示各个源信号分量的时宽，P为各个源信号分量时宽的调制参数，为各个源信号分量的载频，θi表示各个源信号分量的初相位。
所述步骤二中，设定多维变量参数λ的表达式如下：
采用完全随机生成m组定长的二进制串作为群体规模为m的初始种群。参数λ中每个变量的变化范围为[λminλmax]，用q位二进制数p来表示，则有：
将初始种群的数据代入所设的模型：
根据个体适应值的评价函数来判断多维变量参数参数λ的估计值与真实值之间的接近程度，其中xz为个体，x为混合信号；当最大迭代次数大于等于50时，输出此时的多维变量参数λ初始值X0；当最大迭代次数小于50时，个体进行遗传操作，得到新的种群后，从而得到一个新的种群，重复所述步骤二。
所述步骤三中，将多维变量参数λ的初始值X0作为搜索的初始点，利用给定的单纯形的顶点函数值大小，确定单纯形函数的最高点和最低点，通过的反射、扩展及压缩操作构成新的单纯形，当估计值与混合信号的绝对误差E小于1e-6时，输出的解X1为多维变量λ的最优值；
绝对误差计算表达式为E＝|xz-x|，其中，xz为个体，x为混合信号；
反射、压缩和扩展系数为a，b，c均为常数，则反射、压缩和扩展操作分别为：
其中，n为变量的维数，k为迭代次数，表示单纯形形心的函数值，分别表示经过反射、压缩、扩展之后的函数值,x(k) h表示反射之后的函数值与第(n+1)点处函数值中的最小值。
将多维变量参数λ的最优值X1中的参数代入源信号表达式：
即可求出各个高斯调幅通信信号。
优选地，所述遗传操作包括选择复制、交叉和变异操作，其中，进行选择复制和交叉操作的个体的概率大于进行变异操作个体的概率。
本发明所达到的有益技术效果：本发明方法可以有效对单通道时频重叠高斯调幅通信信号进行盲分离，不受初值选取的影响，估计出的参数值精度高，参数的均方根误差降低了8个数量级。这说明该方法具有更高的估计精度且不受初值选取的影响，更适合于多变量参数的搜索。
附图说明
图1本发明流程示意图；
图2本发明之计算多维变量参数初始值和最优值流程示意图；
图3本发明具体实施例中测试用的原始高斯调幅通信信号；
图4本发明具体实施例中由三个源信号混合而成的混合信号图；
图5采用本发明的方法分离后的源信号图。
具体实施方式
为了能更好的了解本发明的技术特征、技术内容及其达到的技术效果，现将本发明的附图结合实施例进行更详细的说明。
下面结合附图和实施例对本发明专利进一步说明。
如图1-2所示，本发明提供一种时频重叠高斯调幅通信信号分离方法，包括以下步骤：
步骤一，建立信号模型，将混合信号的求解过程转化成求解多维变量参数的过程；接收到的混合信号x(t)转换成由M(M≥2)个高斯调幅源信号组成，表达式如下：
其中，表示各个源信号分量的幅度调制系数，t0表示各个源信号分量的时域中心，ωi表示各个源信号分量的时宽，P为各个源信号分量时宽的调制参数，fci为各个源信号分量的载频，θi表示各个源信号分量的初相位。
步骤二，采用遗传算法计算出多维变量参数的初始值：
设定多维变量参数λ的表达式如下：
采用完全随机生成m组定长的二进制串作为群体规模为m的初始种群。参数λ中每个变量的变化范围为[λminλmax]，用q位二进制数p来表示，则有：
将初始种群的数据代入所设的模型：
根据个体适应值的评价函数来判断多维变量参数参数λ的估计值与真实值之间的接近程度，其中xz为个体，x为混合信号；当最大迭代次数大于等于50时，输出此时的多维变量参数λ初始值X0；当最大迭代次数小于50时，个体进行遗传操作，所述遗传操作包括选择复制、交叉和变异操作，其中，进行选择复制和交叉操作的个体的概率大于进行变异操作个体的概率，得到新的种群后，从而得到一个新的种群，重复所述步骤二。
步骤三，采用最小值搜索法计算出多维变量参数的最优值：
将多维变量参数λ的初始值X0作为搜索的初始点，利用给定的单纯形的顶点函数值大小，确定单纯形函数的最高点和最低点，通过的反射、扩展及压缩操作构成新的单纯形，当估计值与混合信号的绝对误差E小于1e-6时，输出的解X1为多维变量λ的最优值；
绝对误差计算表达式为E＝|xz-x|，其中，xz为个体，x为混合信号；
反射、压缩和扩展系数为a，b，c均为常数，则反射、压缩和扩展操作分别为：
其中，n为变量的维数，k为迭代次数，表示单纯形形心的函数值，分别表示经过反射、压缩、扩展之后的函数值,x(k) h表示反射之后的函数值与第(n+1)点处函数值中的最小值。
步骤四，根据多维变量参数的最优值计算出各个源信号：
将多维变量参数λ的最优值X1中的参数代入源信号表达式：
即可求出各个高斯调幅通信信号。
为了验证本发明方法的效果，在matlab软件平台上实现了本发明方法，并与初值控制在远离真实值5％以内的最小值搜索法进行了对比，对比结果如表1所示。
根据表1的数据可以得出，本发明的混合搜索算法相较于最小值搜索算法，均方根误差明显减小，其估计精度更高，且不受初值选取的影响。
图3为中心频率分别为6MHz、6.5MHz、7MHz的单个高斯调幅通信信号。其接收到的混合信号x(t)如图4所示。分离恢复出各个源信号如图5所示。将图5与图3对比可以看出，分离效果很好。
表1分离结果参数对比
以上已以较佳实施例公布了本发明，然其并非用以限制本发明，凡采取等同替换或等效变换的方案所获得的技术方案，均落在本发明的保护范围内。
Claims (2)
1.一种时频重叠高斯调幅通信信号分离方法，其特征在于：包括以下步骤：
步骤一，建立信号模型，将接收到的混合信号的理论表达式转换成由若干个高斯调幅源信号组成，将混合信号的求解过程转化成求解多维变量参数的过程：
将接收到的混合信号的理论表达式x(t)转换成由M个高斯调幅源信号组成，表达式如下：
其中，表示各个源信号分量的幅度调制系数，t0表示各个源信号分量的时域中心，ωi表示各个源信号分量的时宽，为各个源信号分量的载频，θi表示各个源信号分量的初相位；
步骤二，采用遗传算法计算出多维变量参数的初始值：
设定多维变量参数λ的表达式如下：
采用完全随机生成m组定长的二进制串作为群体规模为m的初始种群，参数λ中每个变量的变化范围为[λminλmax]，用q位二进制数p来表示，则有：
将初始种群的数据代入所设的模型：
根据个体适应值的评价函数来判断多维变量参数参数λ的估计值与真实值之间的接近程度，其中xz为个体，x为实际接收到的混合信号；当最大迭代次数大于等于50时，输出此时的多维变量参数λ初始值X0；当最大迭代次数小于50时，个体进行遗传操作，得到新的种群后，从而得到一个新的种群，重复所述步骤二；
步骤三，采用最小值搜索法计算出多维变量参数的最优值：
将多维变量参数λ的初始值X0作为搜索的初始点，利用给定的单纯形的顶点函数值大小，确定单纯形函数的最高点和最低点，通过的反射、扩展及压缩操作构成新的单纯形，当估计值与混合信号的绝对误差E小于1e-6时，输出的解X1为多维变量λ的最优值；
绝对误差计算表达式为E＝|xz-x|，其中，xz为个体，x为实际接收到的混合信号；
反射、压缩和扩展系数为a，b，c均为常数，则反射、压缩和扩展操作分别为：
其中，n为变量的维数，k为迭代次数，表示单纯形形心的函 数值，分别表示经过反射、压缩、扩展之后的函数值,x(k) h表示反射之后的函数值与第(n+1)点处函数值中的最小值；
步骤四，根据多维变量参数的最优值计算出各个源信号：
将多维变量参数λ的最优值X1中的参数代入源信号表达式：
即可求出各个高斯调幅通信信号。
2.根据权利要求1所述的时频重叠高斯调幅通信信号分离方法，其特征在于：所述遗传操作包括选择复制、交叉和变异操作，其中，进行选择复制和交叉操作的个体的概率大于进行变异操作个体的概率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510570687.7A CN105262506B (zh) | 2015-09-09 | 2015-09-09 | 一种时频重叠高斯调幅通信信号分离方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510570687.7A CN105262506B (zh) | 2015-09-09 | 2015-09-09 | 一种时频重叠高斯调幅通信信号分离方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105262506A CN105262506A (zh) | 2016-01-20 |
| CN105262506B true CN105262506B (zh) | 2017-12-19 |
Family
ID=55102040
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510570687.7A Active CN105262506B (zh) | 2015-09-09 | 2015-09-09 | 一种时频重叠高斯调幅通信信号分离方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105262506B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105928877B (zh) * | 2016-04-27 | 2017-03-15 | 广东工业大学 | 强脉冲噪声环境下基于方差比盲分离的复杂结构材料高精度透视检测方法 |
| CN108449291B (zh) * | 2018-01-30 | 2021-09-10 | 西安电子科技大学 | 一种成对载波系统中参数的联合估计方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7542526B2 (en) * | 2004-06-04 | 2009-06-02 | Xg Technology, Inc. | Integer cycle event detection using wavelet pass filter system and method |
| US9088348B2 (en) * | 2007-06-18 | 2015-07-21 | Digi International Inc. | System and method for obtaining frequency and time synchronization in a wideband communication system |
| CN103117819B (zh) * | 2013-01-18 | 2015-08-05 | 宁波大学 | 一种基于压缩感知的脉冲超宽带信号检测方法 |
-
2015
- 2015-09-09 CN CN201510570687.7A patent/CN105262506B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105262506A (zh) | 2016-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | An efficient deep learning model for automatic modulation recognition based on parameter estimation and transformation | |
| Guo et al. | A comparison study of radar emitter identification based on signal transients | |
| CN103941091A (zh) | 基于改进emd端点效应的电力系统hht谐波检测方法 | |
| CN105162740B (zh) | 一种单通道时频重叠信号盲分离方法 | |
| CN104158633A (zh) | 一种基于混合高斯模型的最大似然调制识别方法 | |
| Carlen et al. | Inequalities for quantum divergences and the Audenaert–Datta conjecture | |
| Valipour et al. | Automatic digital modulation recognition in presence of noise using SVM and PSO | |
| CN105262506B (zh) | 一种时频重叠高斯调幅通信信号分离方法 | |
| Yang et al. | One-dimensional deep attention convolution network (ODACN) for signals classification | |
| Ghauri | KNN based classification of digital modulated signals | |
| Yang et al. | Modulation recognition using wavelet transform based on AlexNet | |
| Liu et al. | A novel signal separation algorithm for wideband spectrum sensing in cognitive networks | |
| CN105717490B (zh) | 基于时频分析的lfm信号分离及参数估计方法 | |
| CN102087313A (zh) | 一种卫星搜救信号的频率估计方法 | |
| Bradonjić et al. | The phase transition in inhomogeneous random intersection graphs | |
| Zhu et al. | Modulation classification in cognitive radios for satellite and terrestrial systems | |
| CN105763499A (zh) | 一种基于信号瞬时频率滑动平均的cpfsk调制识别方法 | |
| CN102868653B (zh) | 一种基于双谱和稀疏矩阵的数字调制信号分类方法 | |
| Zhang et al. | DFLNet: Deep federated learning network with privacy preserving for vehicular LoRa nodes fingerprinting | |
| Yin et al. | Modulation classification of MQAM signals based on density spectrum of the constellations | |
| Bari et al. | Identification of L-ary CPFSK in a fading channel using approximate entropy | |
| CN110636053B (zh) | 一种基于局部均值分解和支持向量机的网络攻击检测方法 | |
| Bari et al. | Order recognition of continuous-phase FSK | |
| CN113242201A (zh) | 基于生成分类网络的无线信号增强解调方法及系统 | |
| CN116484180B (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 |