车体垂向振动对高铁FSO通信性能影响的研究 [PDF全文]
(浙江科技学院 a.机械与能源工程学院; b.信息与电子工程学院,杭州 310023)
为了将自由空间光(free space optical,FSO)通信技术应用于高速铁路通信,研究车体垂向振动对FSO通信性能的影响。首先,根据轨道不平顺功率谱密度函数推导其在频域的关系表达式,并利用逆傅里叶变换法对轨道不平顺进行仿真; 然后,根据1/4车体4自由度模型的微分方程,利用拉普拉斯变换推导传递函数搭建基于传递函数的Simulink仿真模型; 最后,通过对车体垂向位移的仿真,推导其概率密度函数表达式,并在高速铁路FSO通信链路中应用动态束腰控制方案,即在发射端实时调整束腰。仿真结果验证了车体垂向位移的概率密度函数表达式,证明动态束腰控制方案在高速铁路FSO通信链路中性能优于传统的固定束腰控制方案。因此,车体垂向位移的概率密度函数和动态束腰控制方案可为高铁FSO通信研究提供理论基础。
Study on impact of carriage vertical vibration on high-speed train FSO communication performance
MIAO Huimina, GUO Jitaoa, MENG Tinga, WU Mingweib, LIU Xi'angb
(a.School of Mechanical and Energy Engineering; b.School of Information and Electronic Engineering, Zhejiang University of Science and Technology, Hangzhou 310023, Zhejiang, China)
To apply FSO(free space optical)communication technology to HST(high-speed train)communication, the impact of carriage vertical vibration on FSO communication performance has been explored. Firstly, the expression of track irregularity in frequency domain was derived in accordance by virtue of the spectral density function of track irregularity, and the inverse Fourier transform method was used to simulate the track irregularity. Secondly, based on the differential equation of the 1/4 carriage model with 4 degrees of freedom, the transfer function was derived using Laplace transform, thereupon establishing the Simulink simulation model. Finally, the PDF(probability density function)was derived by simulatingthe vertical displacement of carriage. And the dynamic beam waist adjustment was applied to the transmitter of HST FSO communication link, where the beam waist was adjusted in real-time at the transmitter. Simulation results verify the PDF expression of the carriage vertical displacement and prove the performance advantage of the dynamic beam waist adjustment scheme over the conventional fixed beam waist adjustment scheme in HST FSO communication link. The PDF of carriage vertical displacement and the dynamic beam waist adjustment scheme can provide a theoretical basis for the study of high-speed train FSO communication.
引言

面对日益增长的交通需求,高速铁路(以下简称高铁)可作为缓解需求的一种重要解决方案,它具有绿色、环保、节能的优点[1]。截至2020年7月底,中国高铁运营里程达到14.14万km[2],根据2016年修订的《中长期铁路网规划》,预计2025年铁路网规模达到17.5万km左右[3]。随着高铁的快速发展,传统的无线射频(radio frequency,RF)技术因存在多普勒频移、带宽低、干扰强等缺点,已经不能满足人们高质量通信的需求。而自由空间光(free space optical,FSO)通信技术可对现有RF技术进行补充,它具有保密性好、抗电磁波干扰[4-5]、部署快捷[6]、不需要频谱许可证[7-8]等优点,不受多普勒频移的影响,并且有足够的带宽来满足数据传输需求。2007年,Kotake等[9]提出了一种在列车和地面之间使用FSO的地对车通信系统,并证明该系统是一种很有发展前景的列车通信系统。Fan等[10]和Mabrouk等[11]分别提出了旋转收发机和多发射机方案,可减少沿轨道部署的基站总数,从而降低了建设成本。Fathi-Kazerooni等[12]和Mohan等[13]分别提出了双波长和分段多波束2种覆盖模型,为列车与基站之间提供无缝连接和切换,因而提高了地对车通信系统性能。高铁地对车FSO通信链路与传统FSO通信链路最大的区别在于列车一直处在前进的运动状态中,并伴随着不可忽略且难以准确估计的机械振动。由于列车的运动方向和速度都是可控且已知的,所以可利用捕获、瞄准和跟踪(acquisition, tracking and pointing,ATP)系统[14-15]修正由于列车的前进运动导致的发射机与探测器在水平方向上的偏移,因此可假设长时间内水平方向上探测器处在发射机发出光束的中心。然而高铁车体垂向振动对通信链路所带来的影响鲜有研究提及,故本研究探索高铁车体移动过程中垂向位移的统计特性,并利用动态束腰控制方案来降低车体垂向振动对通信链路性能的影响。

1 地对车FSO通信链路

图1 高铁地对车FSO通信链路<br/>Fig.1 HST ground-to-train FSO communication link

图1 高铁地对车FSO通信链路
Fig.1 HST ground-to-train FSO communication link

一个典型的高铁地对车FSO通信链路如图1所示,在铁轨的一侧建立基站安装发射机,在列车顶部放置接收机,发射机发出高斯光束被接收机接收。图1中:z为发射机与接收机之间的传输距离; ω0为z=0时高斯光束的光斑半径,即束腰; ωz为高斯光束的光斑半径。

根据高斯光束传输的特性,随着传输距离z的变化,高斯光束的光斑半径ωz表达式[16]11952

图2 高斯光束截面<br/>Fig.2 Gaussian beam cross section

图2 高斯光束截面
Fig.2 Gaussian beam cross section

式(1)中:λ为发射激光波长。图2为高斯光束截面,a为光学探测器C的半径; d为光学探测器C中心与高斯光束中心的偏移距离。在车体静止的情况下,d=0; 在车体运动的情况下,列车前进会导致水平方向上的偏移,该偏移仅与前进时间t和车速v有关,为非随机变量,可利用ATP系统进行校准; 相比车体垂向振动所导致的偏差,车体横向振动所导致的偏差较小[17],可忽略不计,因此d等价于高速列车车体垂向位移。将接收端光学探测器C接收功率和光斑总功率的比例定义为指向性误差hp,可用Marcum Q函数表示[18]

2 车体垂向振动模型2.1 轨道高低不平顺仿真

高铁轨道的轨面并非处于理想的平顺状态,列车行驶过程中带来的轨面的损耗、铁轨的焊接、路基的弹性不均匀等问题都会造成轨道不平顺[19]。通过大量实际数据的测量,已有针对不同轨道拟合出的功率谱密度(power spectral density,PSD)函数来描述轨道不平顺状态。以高速无砟轨道为例,其轨道高低不平顺表达式为S1(F)=A/Fm[20]1,F为空间频率,A、m为拟合系数,参数值及PSD函数曲线参见TB/T 3352—2014《中华人民共和国铁道行业标准》[20]2,4。为了进一步研究PSD函数,我们需要得到以频率f为变量的轨道高低不平顺表达式S2(f)。利用PSD函数在相应单位带宽中均方值相等(S1(F)dF=S2(f)df)的原理[21],再根据空间频率F、频率f及车速v之间的关系f=Fv,可得高速无砟轨道PSD函数的频域表达式为

S2(f)=Af-mvm-1。 (3)

通过估计功率谱的Blackman-Turkey法,有频谱模值与PSD函数之间的关系[22]140如下:

式(4)中:zv(n)为轨道高低不平顺时间序列; Zv(k)为zv(n)的频谱; DFT为离散傅里叶变换; Δf为频率采样间隔; N为时域与频域的采样点数。由于zv(n)为随机过程,其振幅相位皆具有随机性。设n(k)为复高斯白噪声,在式(4)频谱模值基础上乘以n(k),可得轨道高低不平顺频谱

利用离散傅里叶变换的对称性Zv(k)=Z*v(N-K),对Zv(k)进行共轭对称变换,可得Zv(k)(k=0,1,…,N-1)。对Zv(k)进行逆傅里叶变换可得zv(n)。

图3 高速无砟轨道高低不平顺时域样本<br/>Fig.3 Time domain samples of high-speed ballastless track vertical profile irregularity

图3 高速无砟轨道高低不平顺时域样本
Fig.3 Time domain samples of high-speed ballastless track vertical profile irregularity

取车速v=100 m/s,波长范围为2~200 m(即频率范围0.5~50 Hz),时间采样间隔Δt=0.001 s,采样点数N=214进行仿真。图3为高速无砟轨道的高低不平顺时域样本。

2.2 车体模型

图4 1/4车体4自由度模型<br/>Fig.4 1/4 carriage model with 4 degrees of freedom

图4 1/4车体4自由度模型
Fig.4 1/4 carriage model with 4 degrees of freedom

图4为1/4车体4自由度模型。m1为单个转向架质量的1/2; m2为每节满载车体质量的1/4; K1、K2为一系及二系悬架弹簧对应的刚度,用来缓和轨道不平顺所带来的冲击; C1、C2为一系及二系垂向减振器对应的阻尼系数,用来减弱弹簧震荡; Kd1、Kd2为垂向减振器端部连接对应的刚度; zv为轨道高低不平顺输入; zd1、zd2为一系及二系垂向减振器中的活塞在缸筒内的垂向位移; z1为转向架的垂向位移; z2为车体垂向位移。评价人体舒适性的标准为加权加速度均方根值,通过最小化车体垂向振动加权加速度均方根值得最优阻尼比:ξ1=0.33,ξ2=0.22[23]58。最优阻尼比与各悬架系数之间的关系为:

根据D'Alembert's原理,对车体所受动力进行分析,得到垂向振动微分方程组[23]56:

各悬架的详细参数见文献[23]58

2.3 Simulink仿真模型

利用MATLAB软件中的Simulink工具搭建车体仿真模型(图5),其中图5(a)为根据式(6)所搭建的微分方程Simulink仿真模型。由于搭建模型的过程比较繁琐,当车体模型较复杂时,搭建模型的难度将会大大增加,因此我们利用车体垂向位移与轨道不平顺之间的传递函数来搭建车体的Simulink仿真模型,可降低搭建模型的难度。

图5 Simulink仿真模型<br/>Fig.5 Simulink simulation model

图5 Simulink仿真模型
Fig.5 Simulink simulation model

利用拉普拉斯变换,将式(6)从实数域转换到复数域:

式(7)中,s为复变量。求解式(7)中Z2(s)与Zv(s)之比,即为车体垂向位移与轨道不平顺之间的传递函数

式(8)中:

图6 轨道不平顺和车体垂向位移对比<br/>Fig.6 Comparison between track irregularity and carriage vertical displacement

图6 轨道不平顺和车体垂向位移对比
Fig.6 Comparison between track irregularity and carriage vertical displacement

图5(b)为根据式(8)建立的车体传递函数Simulink仿真模型。将2.1节所得的轨道高低不平顺仿真数值zv作为输入,经过车体传递函数Simulink仿真模型,可得输出车体垂向位移z2,轨道不平顺和车体垂向位移对比如图6所示。由图6可知,微分方程Simulink仿真模型与传递函数Simulink仿真模型输出的z2吻合,所以传递函数法也可得到z2,同时大大降低仿真复杂度。另外,z2比zv曲线变得光滑,这是为了乘客的舒适性,以最小化加速度为目的进行的模型参数优化。z2曲线变得光滑表明加速度减少,乘客舒适性提升。但在最小化加速度的情况下,z2幅值变化范围大于zv,因此需要分析z2对通信链路的影响。

2.4 垂向位移仿真数据的统计特性

对于2.3节的车体垂向位移离散序列z2(n),估计其均值[24]90

估计其方差[24]159

通过统计大量仿真数据,μ估计值为0,σ2估计值为1.05×10-5

图7 车体垂向位移概率密度函数<br/>Fig.7 PDF of carriage vertical displacement

图7 车体垂向位移概率密度函数
Fig.7 PDF of carriage vertical displacement

对车体垂向位移z2进行概率密度函数(probability density function,PDF)统计特性分析,车体垂向位移概率密度函数如图7所示。由图7可知,车体垂向位移PDF与高斯分布曲线重叠,由此可发现z2服从高斯分布。因此可定义z2的PDF:

在2.1~2.3节的研究中,为了求得z2需要经过轨道不平顺仿真和车体模型仿真过程,中间的细节繁琐不便。但由z2的PDF可以通过仿真快速得到z2,中间过程大大缩短,减少了仿真时间。以N=225数量级为例,按照传统方法经轨道不平顺仿真、车体模型仿真过程得z2需约95 s,直接由z2的PDF仿真得z2需约0.5 s,仅为前者的0.53%,由此可见z2的PDF仿真法的优越性。另外,在高铁激光通信领域,通过z2的PDF可以推导指向性误差的PDF,这是信道建模的关键步骤,故有重要的研究意义。

3 动态束腰控制方案性能仿真

采用开关键控(on-off keying,OOK)调制方式,系统瞬时误码率[25]

式(12)中:Q(·)为高斯Q函数; R为光学接收器响应度,取0.5; Pt为发射功率,取39.810 7 mW(16分贝毫瓦); N0为接收端噪声单边功率谱密度,取2×10-14 A/Hz; Rdata为数据传输速率,取1 Gb/s[26]

在考虑车体垂向位移z2的所有可能值情况下,利用式(11)~(12)可得系统平均误码率[16]11954

图8 固定及动态束腰控制方案下的系统平均误码率<br/>Fig.8 Average bit error rates under schemes of fixed and dynamic beam waist adjustment systems

图8 固定及动态束腰控制方案下的系统平均误码率
Fig.8 Average bit error rates under schemes of fixed and dynamic beam waist adjustment systems

在传统的固定束腰控制方案中,ω0为固定值,在确定的传输距离下,ωz不变。根据式(13),取仿真参数z=1 km,a=10 cm,λ=1 550 nm[26],固定束腰控制方案下的系统平均误码率如图8所示。

在动态束腰控制方案中,ω0为动态数值,ωz随ω0而变化,文献[16]11956推导其ωz的最优值:

将ωoptz代入式(13),可使得Pe最小化。根据式(13)~(14),取与固定束腰控制方案仿真相同的参数,动态束腰控制方案下的系统平均误码率如图8所示。由图8可知,动态束腰控制方案的系统性能明显优于传统的固定束腰控制方案,并且动态束腰控制方案的性能较稳定。因此动态束腰控制方案可以应用于高铁地对车通信系统来提高系统性能。

4 结 语

本文研究了高铁地对车通信系统,通过对车体进行垂向受力分析,利用拉普拉斯变换推导出车体垂向位移与轨道不平顺之间的传递函数,可以方便快捷地搭建车体传递函数Simulink模型,生成车体垂向位移数据。通过统计,进一步得出了车体垂向位移的概率密度函数表达式,简化了车体垂向振动位移的仿真。采用动态束腰控制方案,将其应用于高铁地对车通信系统,并通过仿真验证了本方案性能的优越性。

参考文献