CN106850109B - 一种莱斯因子估计方法及系统 - Google Patents
一种莱斯因子估计方法及系统 Download PDFInfo
- Publication number
- CN106850109B CN106850109B CN201710040940.7A CN201710040940A CN106850109B CN 106850109 B CN106850109 B CN 106850109B CN 201710040940 A CN201710040940 A CN 201710040940A CN 106850109 B CN106850109 B CN 106850109B
- Authority
- CN
- China
- Prior art keywords
- parameter
- estimation
- received sample
- received
- moment
- 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 31
- 241000209094 Oryza Species 0.000 claims abstract description 47
- 235000007164 Oryza sativa Nutrition 0.000 claims abstract description 47
- 235000009566 rice Nutrition 0.000 claims abstract description 47
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000007480 spreading Effects 0.000 claims description 15
- 238000005562 fading Methods 0.000 claims description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- 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/3911—Fading models or fading generators
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种莱斯因子估计方法及系统。该方法包括步骤:采用Nagakami分布对接收样本的主路径分量进行建模,获得接收样本模型;根据接收样本模型计算接收样本的二阶矩、四阶矩和六阶矩,构建并求解估计参数的方程组,以获取估计参数的数值;根据估计参数的数值对莱斯因子进行估计。该系统包括第一计算模块、第二计算模块和第三计算模块。本发明通过采用Nagakami分布对接收样本的主路径分量进行建模,计算接收样本的二、四和六阶矩以构建估计参数方程组和建立接收样本的似然函数,能够有效提高莱斯因子的估计精度。
Description
技术领域
本发明涉及通信技术领域,尤其涉及一种莱斯因子估计方法及系统。
背景技术
在车载移动通信中,多径严重影响着无线信号的传播。实际上,路边的建筑、树木以及其它障碍会造成信号的随机衰落,当信号强度变小时,车辆之间的通信将变得非常困难。人们采用主路径(信号路径)和辅助路径(所有非信号路径)的能量之比来衡量无线多径信道的衰落程度,称为莱斯(Rician)因子。该名称源于众所周知的莱斯模型,即当接收信号是一个复数常值(对应于主路径)加上一个复高斯随机变量(对应于所有的辅助路径)时,其幅度服从莱斯分布。莱斯因子的引入目的是用来区分主路径和辅助路径的,一旦多径特性被确定,人们就可以采用合适的信号检测策略来完成无线通信或定位。
莱斯因子估计的准确性对无线通信性能的影响至关重要,常用的估计方法是借助接收样本的二阶矩和四阶矩。从数学角度来看,它是一个方程组的解,该方程组具有两个未知数(主路径能量和噪声能量)和两个方程(所测量的二、四矩)。但是,这个模型不适用于典型城市环境的车载移动通信,这是因为这时发射机和接收机的相对运动,将导致随机衰落从而影响接收信号。
发明内容
本发明的目的在于,解决现有技术中存在的上述问题,提供了一种莱斯因子估计方法及系统,通过采用Nagakami分布对接收样本的主路径分量进行建模,计算接收样本的二、四和六阶矩以构建估计参数方程组和建立接收样本的似然函数,能够有效提高莱斯因子的估计精度。
为了实现上述目的,一方面,本发明提供了一种莱斯因子估计方法。该方法包括步骤:采用Nagakami分布对接收样本的主路径分量进行建模,获得接收样本模型;根据接收样本模型计算接收样本的二阶矩、四阶矩和六阶矩,构建并求解估计参数的方程组,以获取估计参数的数值;根据估计参数的数值对莱斯因子进行估计。
优选地,该方法还包括:当估计参数的数值有多组取值时,根据接收样本模型,建立接收样本似然函数;根据估计参数的多组取值确定估计参数的唯一取值,并根据估计参数的唯一取值对莱斯因子进行估计。
优选地,采用Nagakami分布对接收样本的主路径分量进行建模,获得接收样本模型步骤具体包括:
其中,g(sl;ml,Ωl)为sl的概率密度函数;sl为第l个接收样本的主路径分量的幅度,服从Nakagami分布;ml为的第l个接收样本的未知的形状参数,Ωl为第l个接收样本的未知的扩展参数,Γ(ml)为伽马函数。
优选地,根据接收样本模型计算接收样本的二阶矩、四阶矩和六阶矩,构建并求解估计参数的方程组,以获取估计参数的数值步骤具体包括:
优选地,根据接收样本模型,建立接收样本似然函数步骤具体包括:
其中,m表示未知的形状参数,σ2表示方差,r表示接收信号,Ω表示未知的扩展参数,θ表示主路径分量的相位。
另一方面,本发明提供了一种莱斯因子估计系统。该系统包括第一计算模块用于采用Nagakami分布对接收样本的主路径分量进行建模,获得接收样本模型;第二计算模块,用于根据接收样本模型计算接收样本的二阶矩、四阶矩和六阶矩,构建并求解估计参数的方程组,以获取估计参数的数值;第三计算模块,用于根据估计参数的数值对莱斯因子进行估计。
优选地于,第二计算模块还用于:当估计参数的数值有多组取值时,根据接收样本模型,建立接收样本似然函数;根据估计参数的多组取值确定估计参数的唯一取值。
优选地,第一计算模块具体用于:
其中,g(sl;ml,Ωl)为sl的概率密度函数;sl为第l个接收样本的主路径分量的幅度,服从Nakagami分布;ml为的第l个接收样本的未知的形状参数,Ωl为第l个接收样本的未知的扩展参数,Γ(ml)为伽马函数。
优选地,第二计算模块具体用于:
其中,m表示未知的形状参数,σ2表示噪声方差,r表示接收信号,Ω表示未知的扩展参数,θ表示主路径分量的相位。
本发明提供的一种莱斯因子估计方法及系统,通过采用Nagakami分布对接收样本的主路径分量进行建模,计算接收样本的二、四和六阶矩以构建估计参数方程组和建立接收样本的似然函数,能够有效提高莱斯因子的估计精度。
附图说明
图1为本发明实施例提供的一种莱斯因子估计方法的流程示意图;
图4为本发明实施例提供的一种莱斯因子估计系统的结构示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细、清楚、完整的说明。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
图1为本发明实施例提供的一种莱斯因子估计方法的流程示意图。如图1所示,该方法包括步骤101-104:
步骤101:采用Nagakami分布对接收样本的主路径分量进行建模,获得接收样本模型。
将经预处理和数字采样后获得的第l个接收样本记为rl,其中,l=1,…,L,L为接收样本的总数。接收样本中包括主路径分量和辅路径分量:
rl表示第l个接收样本;θl为第l个接收样本的主路径分量的相位,服从在区间[0,2π]上的均匀分布;nl是第l个接收样本的噪声分量,服从方差为σ2的复高斯分布,即L为接收样本总数;sl表示第l个接收样本的主路径分量的幅度,服从Nakagami分布,其概率密度函数为:
其中,ml为第l个接收样本的未知的形状参数,Ωl为第l个接收样本的未知的扩展参数,Γ(ml)为伽马函数。
此外,r1,…,rL是独立同分布的随机变量。莱斯因子为主路径和辅助路径的能量之比,根据式(1)和(2)可以推导出其表达式为:
步骤102,根据接收样本模型计算接收样本的二阶矩、四阶矩和六阶矩,构建并求解估计参数的方程组,以获取估计参数的数值。
其中,|r|=|re-jθ|,结合式(1)进行推导可得:
根据式(4)-(7)可得出估计参数方程组:
为了更加清楚地展现方程组(8)的求解过程,记x1=Ω,x2=σ2,x3=1/m,则相应的方程组(8)可以记为:
根据方程组(9)中的第1个和第2个方程,可得:
根据式(10),可以将式(12)重写为:
根据式(11)可得:
将式(14)代入式(12)中可得:
进一步可以得到:
进一步可以得到:
因式(19)是x1的一元二次方程,求解它我们能够得到x1的两个解。显然,x1的估计值必须是正实数,即其中代表正实数集合。在计算时,如果有充足的接收样本,那么一定存在。求解出后,通过方程组(9)可以得到x2,x3的估计和同样地,和也必须是正实数,即
步骤103,根据估计参数的数值对莱斯因子进行估计。
将估计参数的数值代入式(3)中,以对莱斯因子进行估计:
如果方程组(9)的解均不符合式(20),则需要更多的接收样本重新进行计算以对方程组(9)进行求解。如果根据方程组(9)求解的估计参数的数值有多组取值时,根据接收样本模型,建立接收样本似然函数。根据估计参数的多组取值确定估计参数的唯一取值,并根据估计参数的唯一取值对莱斯因子进行估计。
其中,根据接收样本模型,建立接收样本似然函数的过程如下,且为了简化推导过程,将rl的下角标省略,统一用r表示:
利用全概率定理可以得到
其中hrsθ(r,s,θ;σ2,m,Ω)是r,s和θ的联合概率密度函数,hr(r;σ2,m,Ω)是r的概率密度函数,hr|sθ(r;s,θ,σ2)是在给定s和θ时r的条件概率密度函数。当给定s和θ时,其概率密度函数为
根据上式,式(22)可以被改写为
其中,C=(mmexp{-|r|2/σ2})/((πσ)2Γ(m)Ωm),
需要说明的是,式(25)可以通过数值计算方法来求解。将估计参数的数值的多组取值代入式(25)中,选取使该似然函数的取值最大的那组估计参数的数值为估计参数的唯一取值,并将估计参数的唯一取值代入式(3)中,以对莱斯因子进行估计。
通过以下性能仿真实验进一步证明本发明提出的一种莱斯因子估计方法能达到的有益效果。
基于莱斯信号模型的传统估计方法采用二、四阶矩进行计算可以得到莱斯因子的估计统计量为:
采用蒙特卡洛方法对采用本发明实施例提供的莱斯因子估计法计算出的莱斯因子和采用传统估计方法计算出的莱斯因子进行性能分析,独立仿真的次数是1000。假设随机信号具有Nagakami模型,噪声服从高斯分布,σ2=1。仿真结果如图2和3所示。
图2为和的均值和标准差随K值的变化曲线,其中和的标准差作为置信区间,考虑了两种典型的m值,即m=5和参数m=20。从图中可以看出的均值偏差远小于与真值相差无几。图3给出的是和的均值和置信区间随参数m的变化曲线,考虑了两种K值,即K=5和K=15。图3的结果表明的均值偏差仍小于尤其是在K值比较大的情况。注意到当K=5时,随着m的增大,的均值和真值之间的偏差随之增大,说明所提出的估计方法在K=5时更适用于小m值情况。
本发明实施例提供的一种莱斯因子估计方法,通过采用Nagakami分布对接收样本的主路径分量进行建模,计算接收样本的二、四和六阶矩以构建估计参数方程组和建立接收样本的似然函数,能够有效提高莱斯因子的估计精度。
图4为本发明实施例提供的一种莱斯因子估计系统的结构示意图。该系统200包括第一计算模块201,第二计算模块202和第三计算模块203。
第一计算模块201用于采用Nagakami分布对接收样本的主路径分量进行建模,获得接收样本模型。
第一计算模块具体用于:
其中,g(sl;ml,Ωl)为sl的概率密度函数;sl为第l个接收样本的主路径分量的幅度,服从Nakagami分布;ml为的第l个接收样本的未知的形状参数,Ωl为第l个接收样本的未知的扩展参数,Γ(ml)为伽马函数。
第二计算模块202用于根据接收样本模型计算接收样本的二阶矩、四阶矩和六阶矩,构建并求解估计参数的方程组,以获取估计参数的数值。
第二计算模块具体用于:
第二计算模块还用于当估计参数的数值有多组取值时,根据接收样本模型,建立接收样本似然函数;根据估计参数的多组取值确定估计参数的唯一取值。
第二模块用于根据接收样本模型,建立接收样本似然函数,具体包括:
其中,m表示未知的形状参数,σ2表示噪声方差,r表示接收信号,Ω表示未知的扩展参数,θ表示主路径分量的相位。
第三计算模块203用于根据估计参数的数值对莱斯因子进行估计。
本发明提供的系统实施例与方法实施例相对应,在此不作赘述。
本发明实施例提供的一种莱斯因子估计系统,通过采用Nagakami分布对接收样本的主路径分量进行建模,计算接收样本的二、四和六阶矩以构建估计参数方程组和建立接收样本的似然函数,能够有效提高莱斯因子的估计精度。
以上的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
2.根据权利要求1所述的方法,其特征在于,还包括:
当所述估计参数的数值有多组取值时,根据所述接收样本模型,建立接收样本似然函数;根据所述估计参数的多组取值确定所述估计参数的唯一取值,并根据所述估计参数的唯一取值对莱斯因子进行估计。
6.根据权利要求5所述的系统,其特征在于,所述第二计算模块还用于:
当所述估计参数的数值有多组取值时,根据所述接收样本模型,建立接收样本似然函数;根据所述估计参数的多组取值确定所述估计参数的唯一取值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710040940.7A CN106850109B (zh) | 2017-01-17 | 2017-01-17 | 一种莱斯因子估计方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710040940.7A CN106850109B (zh) | 2017-01-17 | 2017-01-17 | 一种莱斯因子估计方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106850109A CN106850109A (zh) | 2017-06-13 |
CN106850109B true CN106850109B (zh) | 2020-04-28 |
Family
ID=59119171
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710040940.7A Active CN106850109B (zh) | 2017-01-17 | 2017-01-17 | 一种莱斯因子估计方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106850109B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024088919A1 (fr) | 2022-10-27 | 2024-05-02 | Thales | Estimation de paramètres statistiques d'un canal de rice en présence de masquage basée sur du deep learning |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1603360A1 (en) * | 2004-06-04 | 2005-12-07 | Samsung Electronics Co., Ltd. | Velocity estimation apparatus in mobile communication environments |
CN101801073A (zh) * | 2009-02-11 | 2010-08-11 | Ut斯达康通讯有限公司 | 一种莱斯衰落信道k因子的估计方法和设备 |
CN102025680A (zh) * | 2010-12-15 | 2011-04-20 | 华为技术有限公司 | 速度估计方法和装置 |
CN102201876A (zh) * | 2011-05-11 | 2011-09-28 | 南京航空航天大学 | 一种快速的Nakagami信道模型仿真方法 |
-
2017
- 2017-01-17 CN CN201710040940.7A patent/CN106850109B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1603360A1 (en) * | 2004-06-04 | 2005-12-07 | Samsung Electronics Co., Ltd. | Velocity estimation apparatus in mobile communication environments |
CN101801073A (zh) * | 2009-02-11 | 2010-08-11 | Ut斯达康通讯有限公司 | 一种莱斯衰落信道k因子的估计方法和设备 |
CN102025680A (zh) * | 2010-12-15 | 2011-04-20 | 华为技术有限公司 | 速度估计方法和装置 |
CN102201876A (zh) * | 2011-05-11 | 2011-09-28 | 南京航空航天大学 | 一种快速的Nakagami信道模型仿真方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024088919A1 (fr) | 2022-10-27 | 2024-05-02 | Thales | Estimation de paramètres statistiques d'un canal de rice en présence de masquage basée sur du deep learning |
Also Published As
Publication number | Publication date |
---|---|
CN106850109A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi | A new weighted centroid localization algorithm based on RSSI | |
CN105334495B (zh) | 一种无线网络中基于信号到达时间的非视距稳健定位方法 | |
CN108668358B (zh) | 一种应用于无线传感网络的基于到达时间的协作定位方法 | |
CN106851821B (zh) | 一种基于无线通信基站的室内三维定位方法 | |
JP6032462B2 (ja) | 発信源推定方法およびそれを利用した発信源推定装置 | |
CN104066179B (zh) | 一种改进的自适应迭代ukf的wsn节点定位方法 | |
CN106714301B (zh) | 一种无线定位网络中的载波优化方法 | |
CN105898865A (zh) | 非线性非高斯条件下基于ekf和pf的协同定位方法 | |
CN103781095B (zh) | 一种tdoa测量误差的校正方法、传输点和系统 | |
CN101982953B (zh) | 宽带无线通信信道频域多维参数化模型及建模方法 | |
CN104902567A (zh) | 基于最大似然估计的质心定位方法 | |
Meissner et al. | Analysis of position-related information in measured UWB indoor channels | |
CN109521412B (zh) | 基于局部统计量融合的雷达组网空域目标检测方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
CN104661232B (zh) | 基于费歇尔信息矩阵指纹定位精度基本限的ap布置方法 | |
CN112887901A (zh) | 一种基于量化toa量测的凸优化目标定位方法 | |
CN104660356A (zh) | 一种具有可靠虚警性能的半盲协作式频谱感知方法 | |
CN112162244A (zh) | 一种相关噪声和随机丢包环境下的事件触发目标跟踪方法 | |
CN104053234A (zh) | 一种基于rssi的坐标误差补偿定位系统和方法 | |
CN105425206A (zh) | 一种非同步无线网络中的稳健最小二乘定位方法 | |
CN101545964B (zh) | 一种基于无线传感器网络的一维定位方法及装置 | |
CN104076324A (zh) | 一种未知信源数高精度波达方向估计方法 | |
CN106125059B (zh) | 非参数联合估计信号及位置的被动定位方法 | |
CN106850109B (zh) | 一种莱斯因子估计方法及系统 | |
CN104105049A (zh) | 一种减少传声器使用数量的房间冲激响应函数测量方法 |
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 |