CN103297989A - 一种高速铁路高架桥场景下时变k因子模型构建方法 - Google Patents
一种高速铁路高架桥场景下时变k因子模型构建方法 Download PDFInfo
- Publication number
- CN103297989A CN103297989A CN2013101810497A CN201310181049A CN103297989A CN 103297989 A CN103297989 A CN 103297989A CN 2013101810497 A CN2013101810497 A CN 2013101810497A CN 201310181049 A CN201310181049 A CN 201310181049A CN 103297989 A CN103297989 A CN 103297989A
- Authority
- CN
- China
- Prior art keywords
- factor
- speed railway
- become
- factor model
- construction method
- 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.)
- Granted
Links
Images
Landscapes
- Monitoring And Testing Of Transmission In General (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种高速铁路高架桥场景下时变K因子模型构建方法,该方法包括将高速铁路高架桥场景分为完全遮挡、部分遮挡和无遮挡区域,并确定各区域相对基站的水平距离和宽度;根据实测信道冲激响应,消除其误差后估计时变K因子;采用串行高斯函数拟合实测K因子确定各区域K因子的峰值;根据各区域的水平距离、宽度和K因子峰值构建K因子模型,并计算与实测结果的标准差。本发明能够准确的描述列车运行过程中无线信道K因子的变化,可再现实际的传播环境,提高通信***测试和仿真时的准确性和真实性。
Description
技术领域:
本申请涉及无线通信领域,特别涉及一种高速铁路高架桥场景下时变K因子模型构建方法。
背景技术:
近几年来,随着高速铁路的飞速发展,高速列车与地面的通信成为学术界和产业界越来越关注的问题。为了满足乘客对高带宽、高质量的通信业务需求,在保证列车运行安全的前提下,国际铁路联盟(UIC)已经提出现有的窄带铁路列控***将向未来的宽带LTE-R平滑演进。准确地认知无线电波传播特性和信道衰落特征,是高速铁路宽带无线通信***设计的基础。一方面,***部署和网络规划需要电波传播模型进行覆盖预测和干扰分析,另一方面,***评估和物理层优化需要信道衰落模型进行性能仿真与算法验证。无线信道测量是获取信道参数最直接、准确的方式,因此高速铁路典型场景的信道测量是研究高铁无线通信的前提。
列车一般运行于空旷地带,具有较明显的直射环境,高速铁路场景下的无线衰落信号通常服从莱斯分布。莱斯K因子是表征信道衰落程度的参数,它表示直射分量相对散射分量的功率之比。K因子对链路预算、发送分集功率分配和自适应传输等方面有重要影响,除此之外,在信道建模中,K因子也是非常重要的参数。因此,在高速铁路无线通信环境中,K因子需要进行准确地分析与建模。
高架桥是高速铁路的一种典型场景,目前研究该场景下的K因子模型较为理想化和简单化。高架桥场景一般采用两种K因子模型,线性模型和统计模型。线性模型是K因子随距离线性变化的模型,它是基于无线信号传播环境处于完全直射情况下的理想假设,即高架桥周围没有任何障碍物遮挡。统计模型是把K因子建模为服从对数正态分布,仅给出列车运行路线的K因子均值和方差,无法体现真实的无线信号传播环境。但是真实的高架桥场景周围存在非连续的障碍物如树木和建筑的遮挡,然而现有模型无法准确描述该场景下的K因子变化。
发明内容:
为了解决上述问题,本申请公开了一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于,包括:
步骤101、将高速铁路高架桥场景分为完全遮挡、部分遮挡和无遮挡区域,并确定各区域相对基站的水平距离和宽度;
步骤102、根据实测信道冲激响应,消除其误差后估计时变K因子;
步骤103、采用串行高斯函数拟合实变K因子确定各区域K因子的峰值;
步骤104、根据各区域的水平距离、宽度和K因子峰值构建K因子模型,并计算与实测结果的标准差。
优选的,所述高速铁路高架桥场景遮挡类型划分依据为,当遮挡物距离高架桥15m之内并高于列车车顶1m以上时,属于完全遮挡情况;当遮挡物高于高架桥并低于列车车顶时,属于部分遮挡情况;当遮挡物不存在或其高度低于高架桥时,属于无遮挡情况。
优选的,所述误差是指测量设备***误差、信道冲激响应中的噪声分量和信道大尺度衰落。
优选的,所述估计K因子的方法包括1,2阶矩估计法和2,4阶矩估计法。
优选的,所述K因子不是一个固定的常数,而是随时间或相对基站距离变化的变量。
优选的,所述采用串行高斯函数拟合实测K因子,得到各区域K因子的峰值,其中完全遮挡区域为0dB,部分遮挡区域为8dB,无遮挡区域为12dB。
优选的,所述K因子模型与实测结果的标准差为4.5dB,该标准差作为模型的波动参数。
与现有技术相比,本申请包括以下优点:
本申请通过将高速铁路高架桥场景分为完全遮挡、部分遮挡和无遮挡区域,基于交替遮挡的时变K因子模型,能够准确的描述列车运行过程中无线信道K因子的变化,可再现实际的传播环境,构建了一种高速铁路高架桥场景下时变K因子模型构建方法,提高了通信***测试和仿真时的准确性和真实性。
附图说明:
图1为本申请一种高速铁路高架桥场景下时变K因子模型构建的方法流程图;
图2是K因子随距离变化的示意图;
图3是本申请K因子模型验证的计算机仿真图。
具体实施方式:
为使本申请的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本申请作进一步详细的说明。
参照图1,示出了本申请一种高速铁路高架桥场景下时变K因子模型构建方法的流程示意图,包括:
步骤101、将高速铁路高架桥场景分为完全遮挡、部分遮挡和无遮挡区域,并确定各区域相对基站的水平距离和宽度。
由于高架桥铁轨沿线存在不同类型的树木和建筑,这些散射体会影响无线电波的传播。例如,高大茂密的树木会完全遮挡基站与列车之间的直射路径,而矮小稀疏的树木会部分遮挡基站与列车之间的直射路径。因此,本申请把遮挡情况进分类:当遮挡物距离高架桥15m之内并高于列车车顶1m以上时,属于完全遮挡情况;当遮挡物高于高架桥并低于列车车顶时,属于部分遮挡情况;当遮挡物不存在或其高度低于高架桥时,属于无遮挡情况。因此,整个高架桥无线传播环境被划分成三种状态,直射状态(LOS)、有遮挡的直射状态(OLOS)和非直射状态(NLOS)。这样,把传播环境看的“更细”了。然后,根据实际地理测量,获得各区域相对基站的水平距离和宽度,作为后面建模的参数。
步骤102、根据实测信道冲激响应,消除其误差后估计时变K因子。
一般宽带信道测量***测量得到的原始信道冲激响应(CIR)存在***误差、噪声和大尺度衰落的影响,无法直接用于提取K因子。因此,首先要对原始CIR进行处理。理想的宽带信道测量***在整个测量带宽内应该具有频率非选择性,然而这是很难实现的,因为发射和接收滤波器具有不理想的脉冲压缩特性,即
PTx(f)·PRx(f)≠1
其中,PTx(f)和PRx(f)分别为发射和接收滤波器的频率响应。为了消除脉冲压缩误差,需要对原始CIR进行校准,即
其中,Hc(t,f)和Hr(t,f)分别为校准的和原始的时变信道频率响应(CTF)。校准后的CIR可以通过对Hc(t,f)进行IFFT变换得到。
通常实测CIR包含有效的多径分量和无效的噪声分量,为了提高K因子的估计精度,我们需要设定合适的门限阈值来区分多径分量和噪声分量,并去除噪声分量。然而,对于噪声门限,如果采用一个常数值,则在信噪比较低或噪声波动较大的情况下,会影响判决的准确性。因此,采用一种动态的噪声门限计算方法。对于低于噪声门限的多径抽头,对它进行迫零。
为了估计K因子,还需要去除信道大尺度衰落的影响包括路径损耗、阴影衰落和天线增益。这里对CIR进行RMS归一化,即
其中,Pa(t)表示平均功率,可通过长度为N的滑动窗对CIR进行平均得到:
一般滑动窗的长度取10个波长。
消除***误差、噪声和大尺度衰落的影响后,通过只存在小尺度衰落的信道冲激响应,可以估计窄带K因子。首先把宽带CIR转换为窄带CIR,即
假设窄带CIR服从莱斯分布,其概率密度函数(PDF)可以表示为
其中,I0(·)为第一类零阶修正贝塞尔函数,K为莱斯K因子,r为窄带信道冲激响应包络,Ω=E(r2)。然后,采用经典的矩估计方法提取窄带K因子:
其中,E(·)和Var(·)分别表示求均值和方差,Ph为
Ph=|hnarr(t)|2
步骤103、采用串行高斯函数拟合实测K因子确定各区域K因子的峰值。
高斯函数的表达式为
f(n)≡f(n;c1,c2,a1,…,a5)=c1+c2g(n;a1,…,a5)
其中
其中,线性参数c1和c2决定了高斯函数的峰值,a1表示高斯函数极值对应的位置,a2和a3决定了右半边函数的坡度,a4和a5决定了左半边函数的坡度。
根据串行级联的高斯函数,采用最小二乘法(LS)拟合K因子变化如图2所示,其中两个高斯函数拟合LOS区域,一个高斯函数拟合OLOS区域,高斯函数的拼接部分拟合NLOS区域。通过拟合结果,得到各区域K因子的峰值,其中NLOS区域为0dB,OLOS区域为8dB,LOS区域为12dB。
步骤104、根据各区域的水平距离、宽度和K因子峰值构建K因子模型,并计算与实测结果的标准差。
最后构建的K因子模型如下:
KdB(d)=g(d;xm,ym,zm)+Kσ,d∈[-600,600]
其中,
其中,d为基站与列车天线的相对距离,M为LOS和OLOS区域的个数,xm表示LOS和OLOS的峰值K因子,ym表示LOS和OLOS区域的中间位置与基站的相对距离,zm表示LOS和OLOS区域的宽度,这里定义zm等于高斯函数的波峰半高度的宽度。Kσ表示与实测结果的波动值,服从均值为0,标准差为σK的正态分布。表1给出了K因子模型参数:
为了验证上述K因子模型,首先仿真服从莱斯分布的窄带CIR,如下式
hsim(d)=μ(d)+σ(d)[XI(d)+jXQ(d)]
其中,XI(d)和XQ(d)为独立的服从标准正态分布的变量,μ(d)和σ(d)分别为仿真的窄带CIR的均值和标准差,如下两式
仿真数据与实测数据进行对比,如图3所示。这里采用PDF作为衡量标准。从图中可以看出,模型数据与实测数据的PDF基本一致,因此验证了本申请K因子模型的正确性。
以上对本申请所提供的一种高速铁路高架桥场景下时变K因子模型构建方法,进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本申请的限制。
Claims (8)
1.一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于,包括:
步骤101、将高速铁路高架桥场景分为完全遮挡、部分遮挡和无遮挡区域,并确定各区域相对基站的水平距离和宽度;
步骤102、根据实测信道冲激响应,消除实测信道冲激响应的误差后估计时变K因子;
步骤103、采用串行高斯函数对时变K因子进行拟合,确定各区域K因子的峰值;
步骤104、根据各区域的水平距离、宽度和K因子峰值构建K因子模型,并计算与实测结果的标准差。
2.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:
所述高速铁路高架桥场景遮挡类型划分依据为,当遮挡物距离高架桥15m之内并高于列车车顶1m以上时,属于完全遮挡情况;当遮挡物高于高架桥并低于列车车顶时,属于部分遮挡情况;当遮挡物不存在或其高度低于高架桥时,属于无遮挡情况。
3.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:
所述误差是指测量设备***误差、信道冲激响应中的噪声分量和信道大尺度衰落。
4.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:
所述估计时变K因子的方法包括1,2阶矩估计法和2,4阶矩估计法。
5.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:
所述K因子不是一个固定的常数,而是随时间或相对基站距离变化的变量。
6.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:
所述采用串行高斯函数拟合实测K因子,得到各区域K因子的峰值,其中完全遮挡区域K因子的峰值为0dB,部分遮挡区域K因子的峰值为8dB,无遮挡区域K因子的峰值为12dB。
7.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:
所述K因子模型与实测结果的标准差为4.5dB,该标准差作为模型的波动参量。
8.根据权利要求1所述的一种高速铁路高架桥场景下时变K因子模型构建方法,其特征在于:所述构建K因子模型为:
KdB(d)=g(d;xm,ym,zm)+Kσ,d∈[-600,600]
其中,
其中,d为基站与列车天线的相对距离,M为LOS和OLOS区域的个数,xm表示LOS和OLOS的峰值K因子,ym表示LOS和OLOS区域的中间位置与基站的相对距离,zm表示LOS和OLOS区域的宽度,这里定义zm等于高斯函数的波峰半高度的宽度。Kσ表示与实测结果的波动值,服从均值为0,标准差为σK的正态分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310181049.7A CN103297989B (zh) | 2013-05-16 | 2013-05-16 | 一种高速铁路高架桥场景下时变k因子模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310181049.7A CN103297989B (zh) | 2013-05-16 | 2013-05-16 | 一种高速铁路高架桥场景下时变k因子模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103297989A true CN103297989A (zh) | 2013-09-11 |
CN103297989B CN103297989B (zh) | 2016-01-20 |
Family
ID=49098166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310181049.7A Expired - Fee Related CN103297989B (zh) | 2013-05-16 | 2013-05-16 | 一种高速铁路高架桥场景下时变k因子模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103297989B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103702336A (zh) * | 2013-12-27 | 2014-04-02 | 中国科学院自动化研究所 | 一种车路交互无线网络部署的测评方法及装置 |
CN106375045A (zh) * | 2016-09-08 | 2017-02-01 | 北京交通大学 | 一种高速铁路场景下多链路信道探测***与方法 |
CN107949883A (zh) * | 2015-08-21 | 2018-04-20 | 美敦力迷你迈德公司 | 个性化参数建模方法及相关设备和*** |
CN112487653A (zh) * | 2020-12-09 | 2021-03-12 | 同济大学 | 一种基于直射路径出现概率的传播图论的信道仿真方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050053125A1 (en) * | 2003-09-10 | 2005-03-10 | Ganesh Balamurugan | Adaptive equalization using a conditional update sign-sign least mean square algorithm |
CN101801073A (zh) * | 2009-02-11 | 2010-08-11 | Ut斯达康通讯有限公司 | 一种莱斯衰落信道k因子的估计方法和设备 |
CN102223192A (zh) * | 2011-06-15 | 2011-10-19 | 北京交通大学 | 一种高速铁路复合小尺度无线信道模型构建的方法和装置 |
CN102710350A (zh) * | 2012-06-13 | 2012-10-03 | 北京交通大学 | 一种马尔科夫链的高铁平原无线信道模型构建方法和装置 |
-
2013
- 2013-05-16 CN CN201310181049.7A patent/CN103297989B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050053125A1 (en) * | 2003-09-10 | 2005-03-10 | Ganesh Balamurugan | Adaptive equalization using a conditional update sign-sign least mean square algorithm |
CN101801073A (zh) * | 2009-02-11 | 2010-08-11 | Ut斯达康通讯有限公司 | 一种莱斯衰落信道k因子的估计方法和设备 |
CN102223192A (zh) * | 2011-06-15 | 2011-10-19 | 北京交通大学 | 一种高速铁路复合小尺度无线信道模型构建的方法和装置 |
CN102710350A (zh) * | 2012-06-13 | 2012-10-03 | 北京交通大学 | 一种马尔科夫链的高铁平原无线信道模型构建方法和装置 |
Non-Patent Citations (3)
Title |
---|
LARRY J.GREENSTEIN: "Ricean K-Factors in Narrow-Band Fixed Wireless Channels:Theory, Experiments, and Statistical Models", 《VEHICULAR TECHNOLOGY, IEEE TRANSACTIONS ON》 * |
刘留: "高速铁路平原场景无线信道小尺度衰落特征的研究", 《铁道学报》 * |
周涛: "基于测量的高铁宽带无线信道莱斯K因子研究", 《铁道学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103702336A (zh) * | 2013-12-27 | 2014-04-02 | 中国科学院自动化研究所 | 一种车路交互无线网络部署的测评方法及装置 |
CN103702336B (zh) * | 2013-12-27 | 2017-03-29 | 中国科学院自动化研究所 | 一种车路交互无线网络部署的测评方法及装置 |
CN107949883A (zh) * | 2015-08-21 | 2018-04-20 | 美敦力迷你迈德公司 | 个性化参数建模方法及相关设备和*** |
US11484651B2 (en) | 2015-08-21 | 2022-11-01 | Medtronic Minimed, Inc. | Personalized parameter modeling methods and related devices and systems |
US11857765B2 (en) | 2015-08-21 | 2024-01-02 | Medtronic Minimed, Inc. | Personalized parameter modeling methods and related devices and systems |
CN106375045A (zh) * | 2016-09-08 | 2017-02-01 | 北京交通大学 | 一种高速铁路场景下多链路信道探测***与方法 |
CN106375045B (zh) * | 2016-09-08 | 2019-01-29 | 北京交通大学 | 一种高速铁路场景下多链路信道探测***与方法 |
CN112487653A (zh) * | 2020-12-09 | 2021-03-12 | 同济大学 | 一种基于直射路径出现概率的传播图论的信道仿真方法 |
CN112487653B (zh) * | 2020-12-09 | 2023-03-17 | 同济大学 | 一种基于直射路径出现概率的传播图论的信道仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103297989B (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Channel measurements and models for high-speed train communication systems: A survey | |
Luan et al. | Fading characteristics of wireless channel on high‐speed railway in hilly terrain scenario | |
Weiler et al. | Measuring the busy urban 60 GHz outdoor access radio channel | |
Zhang et al. | Measurement‐Based Delay and Doppler Characterizations for High‐Speed Railway Hilly Scenario | |
Sood et al. | Integrating physics-based wireless propagation models and network protocol design for train communication systems | |
Müller et al. | Performance analysis of radio propagation models for smart grid applications | |
Yue et al. | Measurements and ray tracing simulations for non-line-of-sight millimeter-wave channels in a confined corridor environment | |
CN103297989A (zh) | 一种高速铁路高架桥场景下时变k因子模型构建方法 | |
CN102223192B (zh) | 一种高速铁路复合小尺度无线信道模型构建的方法和装置 | |
He et al. | Stochastic channel modeling for railway tunnel scenarios at 25 GHz | |
Guan et al. | Excess propagation loss modeling of semiclosed obstacles for intelligent transportation system | |
CN104185189A (zh) | Lte***干扰定位方法及装置 | |
Zelenbaba et al. | Multi-node vehicular wireless channels: Measurements, large vehicle modeling, and hardware-in-the-loop evaluation | |
CN103561412B (zh) | 一种基于平稳随机过程的信道相关阴影衰落构建的方法 | |
Paschalidis et al. | Investigation of MPC correlation and angular characteristics in the vehicular urban intersection channel using channel sounding and ray tracing | |
Zhao et al. | Channel Characteristics of Rail Traffic Tunnel Scenarios Based on Ray‐Tracing Simulator | |
Ding et al. | Broadband Wireless Channel in Composite High‐Speed Railway Scenario: Measurements, Simulation, and Analysis | |
CN115865709A (zh) | 智能反射面辅助传输的变电站混合中继通信性能分析方法 | |
Nisirat et al. | A Hata based model utilizing terrain roughness correction formula | |
Liu et al. | Investigation of Shadowing Effects in Typical Propagation Scenarios for High‐Speed Railway at 2350 MHz | |
Lu et al. | Fading analysis for the high speed railway viaduct and terrain cutting scenarios | |
Abbas et al. | On understanding the nature of slow fading in LOS microcellular channels | |
Perez et al. | Localization approach based on ray-tracing including the effect of human shadowing | |
Sadinov et al. | Approach for MIMO Wireless Channel Modelling and System Characterization for an Indoor Environment | |
Pascual-García et al. | Fast and accurate electric field estimation from a single ray tracing simulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20160516 |
|
CF01 | Termination of patent right due to non-payment of annual fee |