CN107332632B - MSK调制方式下Alpha稳定分布噪声参数的估计方法 - Google Patents
MSK调制方式下Alpha稳定分布噪声参数的估计方法 Download PDFInfo
- Publication number
- CN107332632B CN107332632B CN201710489292.3A CN201710489292A CN107332632B CN 107332632 B CN107332632 B CN 107332632B CN 201710489292 A CN201710489292 A CN 201710489292A CN 107332632 B CN107332632 B CN 107332632B
- Authority
- CN
- China
- Prior art keywords
- noise
- characteristic function
- characteristic
- alpha
- alpha stable
- 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
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/345—Interference values
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/20—Arrangements for detecting or preventing errors in the information received using signal quality detector
- H04L1/206—Arrangements for detecting or preventing errors in the information received using signal quality detector for modulated signals
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Electromagnetism (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Noise Elimination (AREA)
Abstract
本发明属于噪声参数估计技术领域,公开了一种MSK调制信号通过Alpha稳定分布衰落信道的噪声参数估计方法,包括:计算接收信号的经验特征函数;然后分别求出MSK信号和Alpha稳定分布噪声的特征函数;接着利用信号和噪声的独立性求出接收信号的特征函数;在任意两个相差比较小的时间点t1、t2上令接收信号的经验特征函数等于其特征函数联立方程组,通过对该方程组的求解,得到Alpha稳定分布噪声的特征参数和尺度参数的估计值。本发明在MSK调制信号通过Alpha稳定分布衰落信道时,对Alpha稳定分布噪声的特征参数和尺度参数具有良好的估计性能。
Description
技术领域
本发明属于噪声参数估计技术领域,尤其涉及一种MSK调制方式下Alpha稳定分布噪声参数的估计方法。
背景技术
传统的高斯分布不适用于深空和航天通信环境下噪声的建模,因为它们具有冲击性,即比高斯噪声更为频繁地出现大幅度的数据突变。而Alpha稳定分布却是理想的冲击噪声模型,有良好的理论基础,能够很好地刻画噪声的冲击特性,并且将高斯分布作为一个特例包括在内。由于对Alpha稳定分布噪声中的信号处理方法的研究尚处在起步阶段,所以,发展Alpha稳定分布噪声背景下的信号处理方法具有重要的理论和实际意义。由于MSK调制的相位连续、包络恒定、频偏较小,高带宽利用率和有记忆性,使得这种调制方式广泛应用于存在Alpha稳定分布噪声的通信中。对Alpha噪声参数进行估计,能够为后续的信号解调处理奠定良好的基础。目前针对观察信号具有Alpha稳定分布的情形下估计Alpha稳定分布参数的问题已有广泛的研究,如假设H0:r(t)=n(t),在这种情形下,提出了基于样本分位数、极值理论和ECFs等的算法。然而,针对观察信号是Alpha稳定分布信号和非Alpha稳定分布信号之和的情形下估计Alpha稳定分布参数的问题尚未解决,如假设H1:r(t)=s(t)+n(t)。V.Gautham Chavali等人针对第二种情形提出了新的基于ECF的算法,来解决数字幅度相位调制信号经过Alpha稳定分布噪声衰落信道后噪声参数的估计问题(V.GauthamChavali and Claudio R.C.M.da Silva Detectionof Digital Amplitude-PhaseModulated Signals in Symmetric Alpha-Stable Noise[J].IEEE TransactionsonCommunications,2012,60(11):3365-3375.)。Zhijiang Xu等人推导了MPSK和MQAM叠加加性Alpha稳定分布噪声后噪声参数的估计公式(Zhijiang Xu,Kang Wang and LinminMeng A Novel Parameters Estimation Scheme for Multivariate Symmetric AlphaStable Noise in High-Order Modulation System[C]//International Conference onWireless Communications,NETWORKING and Mobile Computing.IEEE,2008:1-4.)。由于其推导过程只基于实数域,所以理论上估计误差较大,此外,不同于无记忆的MQAM、MPSK和数字幅度相位调制信号,对于有记忆的MSK信号,如何求解其特征函数也是需要解决的技术难题之一。
综上所述,现有技术存在的问题是:目前的Alpha稳定分布噪声参数的估计存在误差较大。
发明内容
针对现有技术存在的问题,本发明提供了一种MSK调制方式下Alpha稳定分布噪声参数的估计方法。
本发明是这样实现的,一种MSK调制方式下Alpha稳定分布噪声参数的估计方法,所述MSK调制方式下Alpha稳定分布噪声参数的估计方法包括:
首先计算接收信号的经验特征函数;然后分别求出MSK信号和Alpha稳定分布噪声的特征函数;
然后利用信号和噪声的独立性求出接收信号的特征函数;
最后在两个时间点上令接收信号的经验特征函数等于其特征函数联立方程组,通过对该方程组的求解,可以得到Alpha稳定分布噪声的特征参数和尺度参数的估计值。
进一步,所述MSK调制方式下Alpha稳定分布噪声参数的估计方法包括:以下步骤:
步骤二,分别求出MSK信号和Alpha稳定分布噪声的特征函数;
步骤三,利用信号和噪声的独立性求出接收信号的特征函数;
步骤四,在两个时间点t1和t2上,根据ECF算法原理,经验特征函数是特征函数的无偏估计,令接收信号的经验特征函数等于其特征函数联立方程组,通过对该方程组的求解,可得到Alpha稳定分布噪声的特征参数和尺度参数的估计值。
进一步,所述步骤一中的信号模型为:
MSK调制信号经过对称Alpha稳定分布噪声衰落信道后到达接收端的信号为:
进一步,所述步骤二中的MSK信号的特征函数求解方法为:
设P(In=1)=1-P(In=-1)=q,q为区间[0,1]上的任意实数,即信源符号不一定是等概的,则状态转移矩阵为:
达到稳态时的概率p0、p1、p2、p3由下式方程组的解给出:
MSK调制信号的特征函数为:
其中,pr(sm)表示传输符号sm的概率,t→0且t>0,上式可泰勒展开到四阶得:
进一步,所述步骤二中Alpha噪声的特征函数求解方法为:
复Alpha稳定分布噪声为n(t)=nR(t)+jnI(t),nR(t)与nI(t)是相互独立且都服从参数为α,γ的对称Alpha稳定分布:
即n(t)服从参数为α,2γ的对称Alpha稳定分布。
进一步,所述步骤三中特征函数求解方法为:
接收信号的特征函数由以下公式求得:
其中,s为MSK信号,n为Alpha稳定分布噪声,φs(t)和φn(t)分别为信号和噪声的特征函数。
进一步,所述步骤四中特征参数和尺度参数的估计方法为:
求解方程组:
可得:
本发明的优点及积极效果为:MSK调制自身的优点以及本发明的所有推导基于复数域,使得本发明的估计性能优于同等条件下,即相同的传输速率和功率,8PSK调制方式下的估计性能:当样本点数为1000时,本发明的特征参数估计均方误差约为1.7×10-3,尺度参数估计均方误差约为0.008,而8PSK调制方式下分别约为4.2×10-3和0.016;当样本点数为10000时,本发明的特征参数估计均方误差约为2×10-4,尺度参数估计均方误差约为0.0008而8PSK调制方式下分别约为4×10-4和0.0018。当样本点取10000时,本发明在不同的特征参数和尺度参数取值下,特征参数的估计均方误差能够达到10-4数量级,在高斯条件(α=2)下特征参数的估计均方误差甚至能够达到10-7数量级,尺度参数的估计均方误差能够达到10-3数量级,说明本发明在Alpha稳定分布噪声衰落信道下是有效可行的。
附图说明
图1是本发明实施例提供的Alpha稳定分布噪声衰落信道后噪声参数的估计方法流程图。
图2是本发明实施例提供的状态转移图。
图3是本发明实施例提供的星座图。
图4是本发明实施例提供的在不同样本点数下特征参数的估计均方误差以及与现有8PSK调制方式下的性能对比,特征参数α=0.5。
图5是本发明实施例提供的在不同样本点数下尺度参数的估计均方误差以及与现有8PSK调制方式下的性能对比,尺度参数γ=1。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明所述的MSK调制方式下Alpha稳定分布噪声参数的估计方法,推导了MSK调制方式下特征参数和尺度参数的估计表达式。目前所查阅到的文献中有MPSK、MQAM和数字幅度相位调制方式下对噪声参数估计的推导,但是MSK调制方式下的估计问题尚未解决。
MSK信号作为CPM信号的一个特例,被广泛应用于深空和航天通信中,由于这类环境下的噪声具有冲击性,所以传统的高斯分布不适用于噪声建模,取而代之的是Alpha稳定分布。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示,本发明实施例提供的Alpha稳定分布噪声衰落信道后噪声参数的估计方法包括以下步骤:
S101:计算接收信号的经验特征函数;然后分别求出MSK信号和Alpha稳定分布噪声的特征函数;
S102:利用信号和噪声的独立性求出接收信号的特征函数;
S103:在两个时间点上令接收信号的经验特征函数等于其特征函数联立方程组,通过对该方程组的求解,可以得到Alpha稳定分布噪声的特征参数和尺度参数的估计值。
下面结合附图对本发明的应用原理作进一步的描述。
如图1所示,本发明为一种MSK信号经过Alpha稳定分布噪声衰落信道后噪声参数的估计方法,所述方法包括以下步骤:
步骤S1中的信号模型为:
MSK调制信号经过对称Alpha稳定分布噪声衰落信道后到达接收端的信号为:
其中,对于二进制MSK来说,是平坦衰落系数,假设其在观察间隔内为已知常数,为了简化推导过程,令其值为1,n(t)是服从对称Alpha稳定分布的噪声。经频率为fs的采样后,有ri=si+ni,i=1,2,…,L。
S2分别求出MSK信号和Alpha稳定分布噪声的特征函数。
步骤S2中MSK信号的特征函数求解方法为:
假设P(In=1)=1-P(In=-1)=q,则状态转移矩阵为:
达到稳态时的概率p0、p1、p2、p3可由下式方程组的解给出:
其中,pr(sm)表示传输符号sm的概率,由于ECF与CF只有在t→0时十分接近且估计误差的方差很小。不失一般性,假设t→0且t>0,上式可泰勒展开到四阶得:
步骤S2中Alpha噪声的特征函数求解方法为:
复Alpha稳定分布噪声为n(t)=nR(t)+jnI(t),nR(t)与nI(t)是相互独立且都服从参数为α,γ的对称Alpha稳定分布,即:
即n(t)服从参数为α,2γ的对称Alpha稳定分布。
S3利用信号和噪声的独立性求出接收信号的特征函数。
步骤S3中特征函数求解方法为:
由于MSK调制信号与Alpha稳定分布噪声是相互独立的,所以接收信号的特征函数可由以下公式求得:
其中,s为MSK信号,n为Alpha稳定分布噪声,φs(t)和φn(t)分别为信号和噪声的特征函数。
S4在两个时间点t1和t2上,根据ECF算法原理,即经验特征函数是特征函数的无偏估计,令接收信号的经验特征函数等于其特征函数联立方程组,通过对该方程组的求解,即可得到Alpha稳定分布噪声的特征参数和尺度参数的估计值。
步骤S4中两个参数的估计方法为:
求解方程组:
可得:
下面结合测试对本发明的应用效果作详细的描述。
为了测试本发明的检验统计量的性能,每个数据由2000次蒙特卡洛实验得到,参数设置如下:码元速率RB=2000波特,采样频率为fs=20000Hz,t1=0.1,t2=0.2。
表1表示在L=10000,α分别取0.5,1,1.5和2,γ分别取0.1、0.6、1.1、1.6、2.1和2.6时,特征参数α的估计均方误差,显然,特征参数的估计均方误差能够达到10-4数量级,在高斯条件(α=2)下特征参数的估计均方误差甚至能够达到10-7数量级,表2表示在L=10000,α分别取0.5,1,1.5和2,γ分别取0.1、0.6、1.1、1.6、2.1和2.6时,γ的估计均方误差,显然,特征参数的估计均方误差能够达到10-3数量级。图4为特征参数α=0.5时,在不同样本点数下特征参数的估计均方误差以及与现有的8PSK在实数域的估计方法的性能对比,图5为尺度参数γ=1,在不同样本点数下尺度参数的估计均方误差以及与现有的8PSK在实数域的估计方法的性能对比,显然,本发明两个参数的估计均方误差更小,性能更优。本发明在Alpha稳定分布噪声衰落信道下,对特征参数和尺度参数的估计有较好的性能。
表1在不同的特征参数和尺度参数取值下,特征参数的估计均方误差,样本点数L=10000。
表2在不同特征参数和尺度参数取值下,尺度参数的估计均方误差,样本点数L=10000。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种MSK调制方式下Alpha稳定分布噪声参数的估计方法,其特征在于,所述的MSK调制方式下Alpha稳定分布噪声参数的估计方法具体包括以下步骤:
步骤二,分别求出MSK信号和Alpha稳定分布噪声的特征函数;
步骤三,利用信号和噪声的独立性求出接收信号的特征函数;
步骤四,在任意两个相差比较小的时间点t1和t2上,根据ECF算法原理,经验特征函数是特征函数的无偏估计,令接收信号的经验特征函数等于其特征函数联立方程组,通过对该方程组的求解,可得到Alpha稳定分布噪声的特征参数和尺度参数的估计值;
所述步骤二中,不同于无记忆的MPSK、MQAM和数字幅度相位调制方式,对于有记忆的MSK信号,利用马尔可夫链建立状态转移矩阵来求解星座图上四个点的概率,具体求解方法为:
设P(In=1)=1-P(In=-1)=q,q为区间[0,1]上的任意实数,即信源符号不一定是等概的,则状态转移矩阵为:
达到稳态时的概率p0、p1、p2、p3由下式方程组的解给出:
所述步骤二中Alpha噪声的特征函数求解方法为:
复Alpha稳定分布噪声为n(t)=nR(t)+jnI(t),nR(t)与nI(t)是相互独立且都服从参数为α,γ的对称Alpha稳定分布:
即n(t)服从参数为α,2γ的对称Alpha稳定分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710489292.3A CN107332632B (zh) | 2017-06-24 | 2017-06-24 | MSK调制方式下Alpha稳定分布噪声参数的估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710489292.3A CN107332632B (zh) | 2017-06-24 | 2017-06-24 | MSK调制方式下Alpha稳定分布噪声参数的估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107332632A CN107332632A (zh) | 2017-11-07 |
CN107332632B true CN107332632B (zh) | 2020-10-09 |
Family
ID=60194235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710489292.3A Active CN107332632B (zh) | 2017-06-24 | 2017-06-24 | MSK调制方式下Alpha稳定分布噪声参数的估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107332632B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108563156A (zh) * | 2018-01-09 | 2018-09-21 | 四川文理学院 | 一种设备模块易于通信的接口及扩展接口系统 |
CN110224769B (zh) * | 2019-05-29 | 2021-05-04 | 浙江科技学院 | 一种通信系统中联合振幅与噪声方差的估计方法 |
CN112733086B (zh) * | 2021-01-22 | 2021-09-07 | 大连交通大学 | 一种基于q阶矩的时变相关指数估计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102882819B (zh) * | 2012-09-04 | 2015-02-04 | 西安电子科技大学 | 非高斯噪声下数字调制信号识别方法 |
CN103997475B (zh) * | 2014-05-29 | 2017-08-29 | 西安电子科技大学 | 一种识别Alpha稳定分布噪声下数字调制信号的方法 |
CN104052702B (zh) * | 2014-06-20 | 2017-12-08 | 西安电子科技大学 | 一种复杂噪声下数字调制信号的识别方法 |
-
2017
- 2017-06-24 CN CN201710489292.3A patent/CN107332632B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107332632A (zh) | 2017-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | Automatic modulation identification based on the probability density function of signal phase | |
CN107332632B (zh) | MSK调制方式下Alpha稳定分布噪声参数的估计方法 | |
CN106911374B (zh) | 一种低复杂度软输出空间调制检测方法 | |
JP6438045B2 (ja) | カバレージ改善のための低papr変調 | |
CN103023831B (zh) | 一种适用于突发波形的载波频偏估计方法 | |
CN103117964A (zh) | 一种60GHz毫米波通信系统下的信号检测方法与装置 | |
CN112350970B (zh) | 一种多相位频移键控调制、解调方法及设备 | |
CN111901269B (zh) | 可变调制指数的高斯频移键控调制方法、装置及系统 | |
Su et al. | A multistage CPE scheme based on crossed constellation transformation for M-QAM | |
CN103338175A (zh) | 一种非相干cpm信号解调设备及解调方法 | |
CN111628950A (zh) | 基于差分星座轨迹图的oqpsk信号盲频偏估计方法 | |
CN103152294B (zh) | 基于信号消除进行噪声估计的方法及系统 | |
CN101764773A (zh) | 格雷编码m-qam调制的并行软比特信息计算的实现方法 | |
Şahin et al. | Over-the-air computation over balanced numerals | |
CN106101044B (zh) | 一种用于相干光通信系统的基于多步插值的频偏估计方法 | |
CN103905370B (zh) | 正交振幅调制信号的软解调方法和装置、移动终端 | |
CN102315835B (zh) | 一种成形滤波器滚降系数估计方法 | |
CN111314027A (zh) | 一种双向中继信道物理层网络编码的符号检测方法 | |
Wang et al. | Non-data-aided frequency offset estimation for binary CPM signals | |
Benkeser et al. | Efficient channel shortening for higher order modulation: Algorithm and architecture | |
Pei et al. | Data-driven bandpass filter design for estimating symbol rate of sporadic signal at low SNR | |
CN114531325B (zh) | 一种无线光通信信道估计方法 | |
CN115695124B (zh) | 一种uqpsk相干解调方法和系统 | |
Teronpi et al. | DTW based Modulation Detection-Verification using Software Defined Radio | |
Simoens et al. | Monte Carlo solutions for blind phase noise estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |