技术领域
[0001] 本
发明涉及无线信道建模领域，更具体的说是涉及一种瑞利圆拱复衰落信道仿真方法。
背景技术
[0002] 经典瑞利信道模型描述了小尺度范围内无线信道多径效应和
多普勒效应，是无线通信信道最重要、最
基础的仿真模型。该信道模型能很好地描述市区无直射路径环境下的移动无线通信信道状态，并且在该模型下，信道相
角统计特性服从瑞利分布，相角统计服从均匀分布，且其多普勒
功率谱为经典Jakes功率
密度函数，表现为U型谱。然而，在2.5GHz固定无线通信信道中，由于发射器和接收器的相对运动，造成了非常小的多普勒扩展，在这种情况下，多普勒功率谱的功率主要集中在一定
频率范围的中心，此时，多普勒功率谱表现为圆拱谱。该信道模型也是经典瑞利信道模型的一个重要分支，但鉴于相关文献和
专利并未对瑞利圆拱复衰落信道仿真方法进行相关论述，因此，有必要提出一种瑞利圆拱复衰落信道建模和仿真方法。
发明内容
[0003] 本发明的目的是在于提供一种瑞利圆拱复衰落信道仿真方法，旨在解决现有仿真方法中无法准确描述圆拱衰落信道的问题。
[0004] 本发明解决其技术问题采取的技术方案是：
[0005] 一种对瑞利圆拱复衰落信道仿真方法，包括以下步骤：
[0006] 步骤一，程序开始、进行相关参数初始化；
[0007] 步骤二，产生瑞利圆拱复衰落信道；
[0008] 步骤三，对瑞利圆拱复衰落信道进行统计分析、程序结束；
[0009] 步骤四，设计瑞利圆拱复衰落信道
硬件仿真系统。
[0010] 进一步，所述步骤一包括以下步骤：
[0011] 程序开始后，即可进行参数初始化，设定瑞利圆拱复衰落信道的最大多普勒频移fd，系统
采样频率fs，系统采样点数N，以及多普勒圆拱
滤波器相关参数α0、α2、α4。
[0012] 进一步，所述步骤二包括以下步骤：
[0013] 参数设置完成后，即可产生瑞利圆拱复衰落信道。将两路高斯白噪声输入具有圆拱形功率谱的成形滤波器即多普勒圆拱滤波器，从而产生两路具有圆拱形多普勒功率谱的有色高斯白噪声序列，分别组成瑞利圆拱复衰落信道的
实部和
虚部，该瑞利圆拱复衰落信道的包络服从瑞利分布，相角服从均匀分布。两路
信号即实部和虚部的产生原理一致，以其中一路为例加以说明，通过高斯白噪声生成瑞利圆拱复衰落信道的实部如下：
[0014] 第2-1步，时域复高斯白噪声经过FFT变换后，转换至频域。
[0015] 第2-2步，多普勒圆拱滤波器对频域复高斯白噪声进行滤波，并且该滤波器的数学表达式均可由频率f的多项式给出：
[0016]
[0017] 其中，瑞利圆拱复衰落信道的最大多普勒频移fd、参数α0、α2、α4的数值由参数初始化过程定义，且
[0018] 第2-3步，频域复高斯白噪声经滤波处理后，由IFFT变换将频域信号转换至时域信号，即可得到瑞利圆拱复衰落信道的实部 hI(t)。
[0019] 同理，另一路的频域复高斯白噪声由IFFT变换，并通过换向器换向处理后，即可得到瑞利圆拱复衰落信道的虚部hQ(t)，且实部 hI(t)和虚部hQ(t)相互
正交，实部hI(t)和虚部hQ(t)相加得到瑞利圆拱复衰落信道，即h(t)＝hI(t)+jhQ(t)。
[0020] 进一步，所述步骤三包括以下步骤：
[0021] 产生瑞利圆拱复衰落信道后，即可分别对其进行一阶和二阶统计分析。其中一阶统计分析包括对瑞利圆拱复衰落信道幅值分布和相角分布的分析，二阶统计分析包括瑞利圆拱复衰落信道自相关函数和多普勒功率谱密度函数分析，通过与理论值相比较，验证该方法的准确性与科学性。
[0022] 一阶、二阶统计分析为业内现有通用技术，属于信号分析的范畴，如专利号为201410288810.1公开的《一种Nakagami复衰落信道建模方法》。
[0023] 进一步，所述步骤四中，瑞利圆拱复衰落信道硬件仿真系统，主要包括发射器模
块、
模数转换模块、FPGA信道仿真器、计算机、
数模转换模块、接收器模块。发射器模块通过模数转换模块与FPGA信道仿真器相连，计算机亦与FPGA信道仿真器相连，FPGA信道仿真器通过数模转换模块与接收器模块相连。
[0024] 其中，发射器模块产生基带发射信号s(t)，发送至模数转换模块，并将数字化的基带发射信号s(t)传输至FPGA信道仿真器作为
输入信号。同时，通过计算机的MATLAB和ISE
软件创建并生成瑞利圆拱复衰落信道h(t)的IP核，将IP核导入至FPGA信道仿真器作为模拟信道。经FPGA信道仿真器后，得到离散化的模拟接收器接收信号 r(t)，并由数模转换模块转换为连续信号，最终传输至接收器模块。
[0025] 离散化的模拟接收器接收信号r(t)表示为：
[0026]
[0027] 式中，n为无线信道多径数目；hk(t)＝hI(t)+jhQ(t)为第k条路径下的瑞利圆拱复衰落信道，且hk(t)均由计算机产生的瑞利圆拱复衰落信道h(t)给出，具有时变性；s(t)为基带发射信号，由发射器模块产生；τk为时延，由计算机进行定义。
[0028] 通过FPGA实现多径时变
频率选择性衰落信道时，离散化的模拟接收器接收信号r(t)可表示为多径时变频率选择性衰落信道。
[0029] 其中，所述步骤一至步骤三可在计算机MATLAB完成。所述步骤四可联合计算机和FPGA板卡等实现巴特沃兹复衰落瑞利信道仿真的功能。
[0030] 本发明与现有仿真方法相比较具有以下优点：
[0031] (1)采用成形滤波器即多普勒圆拱滤波器，便于实现瑞利信道圆拱形多普勒功率谱，从而实现瑞利圆拱复衰落信道。
[0032] (2)成形滤波器法是用随机噪声产生信道的，因而比
正弦波叠加法更具随机性，其一阶和二阶仿真结果更为精确，具有明显优势。
[0033] (3)该仿真方法可以通过优
化成形滤波器的相关参数，从而能便捷地改善仿真模型的性能。
[0034] (4)该信道硬件仿真系统装置简洁、易于实现，能应用于城区的移动通信无线信道的仿真和测试中。
附图说明
[0035] 图1为本发明瑞利圆拱复衰落信道仿真方法
流程图。
[0036] 图2为本发明瑞利圆拱复衰落信道的生成流程图。
[0037] 图3为本发明仿真结果之瑞利圆拱复衰落信道时域幅值
波动图。
[0038] 图4为本发明仿真结果之瑞利圆拱复衰落信道一阶统计特性图。
[0039] 图5为本发明仿真结果之瑞利圆拱复衰落信道二阶统计特性图。
[0040] 图6为本发明瑞利圆拱复衰落信道硬件仿真系统结构图。
[0041] 图7为本发明模拟接收器接收信号r(t)生成原理图。
具体实施方式
[0042] 以下结合附图和实施方法对本发明作进一步说明。
[0043] 参照图1，本发明瑞利圆拱复衰落信道仿真方法，包括以下步骤：程序开始、参数初始化U1、产生瑞利圆拱复衰落信道U2、对瑞利圆拱复衰落信道进行统计分析U3、程序结束。
[0044] 程序开始后，即可进行参数初始化U1，不失一般的，我们设定瑞利圆拱复衰落信道的最大多普勒频移fd＝500Hz，系统
采样频率 fs＝15000Hz，系统采样点数N＝106，α0＝1，α2＝-1.72，α4＝0.785。
[0045] 参数设置完成后，即可产生瑞利圆拱复衰落信道U2。将两路高斯白噪声输入具有圆拱形功率谱的成形滤波器，从而分别产生具有圆拱形多普勒功率谱的有色高斯白噪声序列，并分别组成了瑞利圆拱复衰落信道的实部和虚部，该瑞利圆拱复衰落信道的包络服从瑞利分布，相角服从均匀分布。
[0046] 产生瑞利圆拱复衰落信道后，即可分别对其进行一阶和二阶统计分析。其中一阶统计分析包括对瑞利圆拱复衰落信道幅值分布和相角分布分析，二阶统计分析包括对瑞利圆拱复衰落信道自相关函数和多普勒功率谱密度函数分析，通过与理论值相比较，验证该方法的准确性与科学性。
[0047] 参照图2，本发明瑞利圆拱复衰落信道，主要包括复高斯白噪声 U21和U25、FFT变换U22和U26、多普勒圆拱滤波器U23和U27、IFFT 变换U24和U28、换向器U29、瑞利圆拱复衰落信道U30。其中，复高斯白噪声U21、FFT变换U22、多普勒圆拱滤波器U23、IFFT变换 U24产生信道瑞利圆拱复衰落信道U30的实部hI(t)；复高斯白噪声 U25、FFT变换U26、多普勒圆拱滤波器U27、IFFT变换U28、换向器 U29产生瑞利圆拱复衰落信道U30的虚部hQ(t)，并且实部hI(t)和虚部hQ(t)相互正交。两路信号的产生原理一致，以其中一路为例加以说明。
[0048] 时域复高斯白噪声U21经过FFT变换U22后，转换至频域。
[0049] 多普勒圆拱滤波器U23和U27对频域复高斯白噪声进行滤波，并且该滤波器的数学表达式均可由频率f的多项式给出：
[0050]
[0051] 其中，瑞利圆拱复衰落信道的最大多普勒频移fd、α0、α2、α4数值由参数初始化过程定义，且
[0052] 频域复高斯白噪声经滤波处理后，由IFFT变换U24将频域信号转换至时域信号，即可得到瑞利圆拱复衰落信道U30的实部hI(t)。同理，另一路的频域复高斯白噪声由IFFT变换U28，并通过换向器 U29换向处理后，即可得到瑞利圆拱复衰落信道U30的虚部hQ(t)，且实部hI(t)和虚部hQ(t)相互正交，实部hI(t)和虚部hQ(t)相加得到瑞利圆拱复衰落信道U30，即h(t)＝hI(t)+jhQ(t)。
[0053] 本发明对瑞利圆拱复衰落信道进行了仿真，其中，相关参数分别为：最大多普勒频移fd＝500Hz，系统采样频率fs＝15000Hz，系统采样点数N＝106，α0＝1，α2＝-1.72，α4＝0.785，相应的仿真结果如图 3、图4和图5所示。其中，图3为时域瑞利圆拱复衰落信道U30的幅值
波形图，描述了信道在不同时刻的衰落情况。图4描述了瑞利圆拱复衰落信道U30的一阶统计特性，可以看出，其幅值的概率密度分布服从瑞利分布，而相角的概率密度分布服从均匀分布，符合瑞利信道的一阶统计特性。图5描述了瑞利圆拱复衰落信道U30的二阶统计特性，可以看出本发明输出信道的二阶统计特性自相关函数和多普勒功率谱密度函数与理论值均能很好的吻合。
[0054] 完成参数初始化U1、产生瑞利圆拱复衰落信道U2、对瑞利圆拱复衰落信道进行统计分析U3后，还包括设计瑞利圆拱复衰落信道硬件仿真系统。
[0055] 参照图6，瑞利圆拱复衰落信道硬件仿真系统，主要包括发射器模块6-1、模数转换模块6-2、FPGA信道仿真器6-3、计算机6-4、数模转换模块6-5、接收器模块6-6。发射器模块6-1通过模数转换模块6-2与FPGA信道仿真器6-3相连，计算机6-4亦与FPGA信道仿真器6-3相连，FPGA信道仿真器6-3通过数模转换模块6-5与接收器模块6-6相连。
[0056] 其中，发射器模块6-1产生基带发射信号s(t)，发送至模数转换模块6-2，并将数字化的基带发射信号s(t)传输至FPGA信道仿真器6-3作为输入信号。同时，计算机6-4的MATLAB和ISE软件创建并生成瑞利圆拱复衰落信道h(t)的IP核并导入至FPGA信道仿真器 6-3作为模拟信道。经FPGA信道仿真器进行相关数乘运算后，得到离散化的模拟接收器接收信号r(t)，并由数模转换模块6-5转换为连续信号，最终传输至接收器模块6-6进行相关分析。不失一般地， FPGA信道仿真器的板卡采用Xilinx公司Kintex-7系列的KC705型号。
[0057] 参照图7，通过FPGA实现多径时变频率选择性衰落信道时，离散化的模拟接收器接收信号r(t)可表示为多径时变频率选择性衰落信道。离散化的模拟接收器接收信号r(t)表示为：
[0058]
[0059] 式中，n为无线信道多径数目；hk(t)＝hI(t)+jhQ(t)为第k条路径下的瑞利圆拱复衰落信道，且hk(t)均由计算机产生的瑞利圆拱复衰落信道h(t)给出，具有时变性；s(t)为基带发射信号，由发射器模块产生；τk为时延，由计算机进行定义。
[0060] 本发明之一种瑞利圆拱复衰落信道建模方法及仿真器设计能实现对瑞利圆拱复衰落信道进行建模和仿真的功能，具有良好
稳定性和准确性的特点，适用于需要对无线信道建模、仿真和特征参数提取的场合。
[0061] 本领域的技术人员可以对本发明进行各种
修改和变型，倘若这些修改和变型在本发明
权利要求及其等同技术的范围之内，则这些修改和变型也在本发明的保护范围之内。
[0062]
说明书中未详细描述的内容为本领域技术人员公知的
现有技术。
