一种mppm调制的软解调算法及误符号率性能估计方法
【专利摘要】本发明公开了一种MPPM调制的软解调算法及误符号率性能估计方法，软解调算法包括：产生随机的0、1序列，进行(N,M)MPPM信号调制；产生指数韦伯信道衰落和高斯白噪声，得到受高斯白噪声和信道衰落干扰的MPPM信号块，抽样，找到其中M个最大值，解调，仿真足够次数得到误符号率。该算法复杂度低且误符号率性能优于硬解调算法。误符号率估计方法包括：建立大气信道的数学模型，计算信号与非信号时隙抽样值的累积分布函数，分别求出M和N?M个非空时隙随机变量最小值和最大值的概率密度函数，求出其差值ΔY的概率密度函数，计算ΔY的概率密度函数的积分，所得结果即为(N,M)MPPM的误符号率。解决了无法通过蒙特卡洛仿真算法直接对(N,M)MPPM的误符号率进行理论分析的问题。
【专利说明】
一种MPPM调制的软解调算法及误符号率性能估计方法
技术领域
[0001 ]本发明涉及无线光通信系统中的(N，M)MPPM调制方式在指数韦伯信道下的软解调 计算方法，属于无线光通信技术领域。
【背景技术】
[0002] 无线光通信又称自由空间光通信(free-space optical，FS0)，其采用大气信道作 为传输媒介，而大气湍流导致的光学折射率随机起伏会使光信号在传输过程中产生光强闪 烁、光束漂移和光束扩展等效应。这些效应引起的信道衰落会导致传输的光束质量下降，误 码率劣化，进而影响系统的稳定性和可靠性。
[0003] 为了研究大气湍流对FS0通信的影响，国内外科研人员对大气湍流信道模型进行 了深入研究，先后提出了对数正态分布、K分布、gamma-gamma分布等模型，结果表明：对数正 态分布模型适合于弱端流，K分布适合于强端流，而gamma-gamma分布则适合于强、中、弱三 种湍流情况。但上述分布都不能很好地模拟采用有限孔径接收机无线光通信系统的光信号 强度波动，于是2012年首次从可靠性工程领域引入指数韦伯（exponentiated Weibull，EW) 分布以描述孔径平均效应存在情况下的光强衰落。
[0004] 基于上述信道模型，各国研究人员对无线光通信的调制解调方式也进行了深入探 讨和研究。对于调制方式，开关键控(On-Off Keying，00K)作为最常用的调制方案，实现简 单，但其抗干扰性能较差。而新型的多脉冲位置调制（Multipulse Pulse Position Modulation，MPPM)方式，其抗干扰能力相比于00K具有很大的优势，并且带宽利用率较高。
[0005] 对于解调方式，通常分为硬解调和软解调。硬解调通过固定、最佳或动态判决门限 等方法把含噪声和衰落的信号判定为相应的调制电平，技术简单、易实现，但性能往往不够 理想。软解调对于接收到的信号，不急于做判决，通过计算相应比特的软信息，提高判决的 准确性，降低系统误码率，其性能远远优于硬解调。
[0006] 目前的问题在于：由于MPPM的复杂性，关于其软解调特性的报道较少。目前已知软 解调系统算法大多是基于条件信道下的计算，即只考虑存在加性高斯白噪声，不考虑大气 导致的衰落。该算法不能用于指导实际的FS0通信系统接收机设计。
[0007] 因此基于EW分布的无线光通信MPPM调制的软解调算法及其误符号率性能的估计 实属当前十分重要的研究方向。
【发明内容】
[0008] 本发明的目的是提出一种基于指数韦伯(EW)信道的无线光通信(N，M)MPPM调制的 软解调算法并对其误符号率进行性能估计。其中软解调算法的优点在于无需判决门限，且 判决的性能优于一般解调算法。之后提出的一种误符号率性能估计方法的数学模型，对上 述(N，M)MPPM的软解调算法误符号率进行理论分析，进而正确有效的计算出理论误符号率。
[0009] 本发明的目的是通过下述技术方案来实现的。
[0010] -种MPPM调制的软解调算法，包括下述步骤：
[0011] 1)由MATLAB中rand函数产生0、1随机序列作为仿真输入比特流，并根据映射表进 行(N，Μ) MPPM信号调制，其中，N为一个MPPM信号块中的时隙个数，Μ为N个时隙中的脉冲个 数；
[0012] 2)在MATLAB中由逆函数法产生符合指数韦伯分布的衰落信道系数并由randn函数 产生高斯白噪声随机数；
[0013] 3)将步骤2)产生的信道衰落系数和高斯白噪声随机数加入到步骤1)获得的(N，M) MPPM调制信号中，得到受到高斯白噪声和信道衰落干扰的MPPM信号块；
[0014] 4)对步骤3)得到的(N，Μ )MPPM信号块中的每个时隙进行抽样；
[0015] 5)利用MATLAB中的MAX函数找到抽样值中的最大值，将抽样值中Μ个最大值对应的 时隙号判为1，其余均判为〇;
[0016] 6)利用步骤1)中（Ν，Μ)MPPM信号调制时的映射表，将步骤5)经过判决的信号块恢 复为〇、1序列，并将该序列与步骤1)中产生的信号对照，判断是否与其相同，若相同，则为正 确解调;若不相同，则错误符号个数加1;
[0017] 7)重复步骤1)到6)，仿真足够多的次数，得到错误符号个数，计算仿真的误符号 率。
[0018] 相应地，本发明给出了基于上述(N，M)MPPM调制的软解调算法进行误符号率性能 估计的方法，包括下述步骤：
[0019] 1)建立点对点无线光通信大气信道的数学模型；
[0020] 2)计算信号时隙与非信号时隙抽样值的概率密度函数；
[0021] 3)将接收到的(Ν，M)MPPM调制信号的每个时隙抽样值作为互相独立的随机变量， 求出各自的累积分布函数，在此，N为一个MPPM信号块中的时隙个数，Μ为N个时隙中的脉冲 个数；；
[0022] 4)求出Μ个非空时隙随机变量最小值Ysmin的概率密度函数，再求出N-Μ个空时隙随 机变量中最大值Yr?aX的概率密度函数；
[0023] 5)得到非空时隙随机变量的最小值与空时隙随机变量的最大值的差值Δ Y，求出 差值的概率密度函数；
[0024] 6)计算ΔΥ的概率密度函数从负无穷到零处的积分，所得结果即为(N，M)MPPM的误 符号率。
[0025]本发明具有以下优点：
[0026] 1)本发明中首先通过蒙特卡洛仿真实现了一种EW信道下MPPM调制的软解调算法， 该方法主要通过在接收端对(N，M)MPPM块的每个时隙进行抽样得到前Μ个较大值，完成对脉 冲所在时隙的判断。该方法的算法复杂度低，误符号率系统性能上优于硬解调算法；
[0027] 2)提出一种适于理论计算软解调MPPM无线光通信系统误符号率性能的方法，解决 了无法通过上述蒙特卡洛仿真算法直接对(N，M)MPPM的误符号率进行理论分析的问题。
【附图说明】
[0028]图1是本发明总体的流程步骤；
[0029] 图2是弱湍流情况下误符号率拟合值、理论值与仿真值的对比图；
[0030] 图3是中湍流情况下误符号率拟合值、理论值与仿真值的对比图；
[0031] 图4是强湍流情况下误符号率拟合值、理论值与仿真值的对比图；
[0032] 图5是中湍流情况下采用动态门限硬解调和软解调误符号率的对比图。
【具体实施方式】
[0033]为使本发明的目的、技术方案和优点更加清楚，下面结合附图和【具体实施方式】进 一步详细说明。本实施例仅表示对本发明的原理性说明，不代表对本发明的任何限制。 [0034]如图1所示，本发明以（5，2 )MPPM为例，利用MATLAB软件实现蒙特卡洛仿真。详细介 绍指数韦伯(EW)信道下的无线光通信(N，M)MPPM调制的软解调算法。
[0035] 步骤1利用MATLAB中rand函数产生随机的0，1序列，并进行(5，2) MPPM调制
[0036] la)首先需要设定仿真次数T，T的取值需要足够大以确保仿真的一般性，在本实例 中我们取τ = 108。
[0037] 根据MPPM调制原理可知，由于2:' <C\2 = 10<24，所以（5,2)ΜΡΡΜ调制可以对3比特 (bit)为一组的比特流进行调制，其二进制比特流与ΜΡΡΜ符号之间的映射关系如表1所示， 其中（11，12 )表不有脉冲时隙所在位置，此表格表不的映射关系仅为一种实例。
[0038] lb)利用MATLAB中rand函数产生3个0，1随机数，作为一组对应映射表格进行调制， 得到有5个时隙的MPPM信号块[ Χ1，χ2，…，Xi ]，其中i = 5。
[0039] 表1(5,2)MPPM调制映射表
[0040] ±±± UU ± ± U
\0 y J
[0041 ]步骤2产生符合EW的信道衰落和高斯白噪声
[0042] 2a)已知EW型衰落下光强波动I的概率密度函数和累积分布函数分别为：
[0043]
[0044]
[0045]
[0046] 其中，形状参数β>0,其取值与闪烁指数相关。尺度参数n>0,其取值由β决定并与 辐照度均值相关。定义α>0为在给定观测空间内，准直传播的光束与非准直传播的光束被 成功接收的平均量。
[0047] 根据以上公式，利用逆函数法，可以产生符合EW分布的随机数EW_rand。
[0048] 2b)关于加性高斯白噪声，已知平均信噪比为：
[0049] sm = μ2 Ι·2σ]
[0050] 式中，μ为系数，y2=(RNPt/M)2，R表示光电转换效率，P t表示发射端的平均发射功 率。本实例中设定R，Pt的值均为1，N与Μ分别等于5、2。则加性高斯白噪声的方差可以表示 为：
[0051]
[0052]根据上述公式得出的噪声方差利用randn函数即可产生符号高斯白噪声分布的随 机数AWGN。
[0053]步骤3将信号通过信道，得到受噪声和信道衰落干扰的MPPM信号块
[0054] 3a)对于接收到的MPPM信号而言，信号被抽样之后的输出模型可用以下式子表示：
[0055]
[0056]式中，由上到下依次为MPPM信号块中有脉冲的时隙和无脉冲的时隙，η表示的是均 值为〇,方差为h = 乂/2的加性高斯白噪声，No为其单边功率谱密度。h表示信道衰落。
[0057] 3b)将步骤2所产生的符合EW分布的5个随机数依次与一个(5,2)MPPM信号块中5个 时隙的脉冲信号相乘。之后再将产生的5个高斯白噪声随机数依次与脉冲信号相加，最终得 到接收的一个MPPM信号块。即，
[0058] γ? = μ · xi · Eff_rand+AWGN
[0059] 其中μ为系数，EW_rand为步骤2产生的信道随机数，AWGN为步骤2产生的高斯白噪 声。
[0060] 步骤4在接收端对MPPM信号进行抽样
[00611对接收到的信号进行抽样，得到5个时隙的抽样值，即Y = [ Yi，Y2，Y3，Y4，Y5 ]。
[0062]步骤5在抽样值中利用MATLAB函数找到最大值，进行判决 [0063] 在Y3，Y4，Y5]中利用MATLAB中的MAX函数找到前两个最大值(不分次序） Yi，Yj，将其对应的时隙号i、j的信号值判为1，其余均判为0。如果有相同的数值而产生两个 以上的最大值，则直接判断为解调错误，错误符号个数加1，返回步骤1，进行下一次仿真。
[0064] 步骤6根据判决得到的脉冲所在时隙号解调信号
[0065] 由步骤5可知脉冲时隙所在位置(1^12)为（i，j)(或（j，i)，按序号从小到大排列） 对照表1中的映射表，恢复信号为3bit-组的比特流。对照步骤1中产生的信号，判断传输是 否相同，若相同则正确解调;若不同或映射表中没有解调出来的情况则判断为解调错误，错 误符号数加1。
[0066]步骤7得到仿真的误符号率
[0067]重复步骤1到6足够的次数T。完成所有解调后，通过错误个数和解调MPPM信号块的 总个数，计算出误符号率。
[0068]注意，当取不同的平均信噪比snr进行计算时，可以得到不同信噪比下的误符号 率。
[0069]如图1所示，在上述步骤的基础上，指数韦伯(EW)信道下无线光通信(Ν,Μ)ΜΡΡΜ软 解调的误符号率性能估计方法，通过下述具体步骤进行计算：
[0070]步骤1建立点对点无线光通信大气信道的数学模型
[0071] 设X为发送信号序列，y为接收信号序列，则信道的数学模型为：
[0072] y = hx+n
[0073] η表示均值为0、方差为σΗ2的加性高斯白噪声。h为服从EW分布的信道衰落。
[0074]步骤2计算信号时隙与非信号时隙抽样值的概率密度函数
[0075] 根据已知EW型衰落下光强波动的概率密度函数和累积分布函数及步骤1，可以计 算出信号时隙与非信号时隙抽样值的概率密度函数分别为：
[0076]
[0077]
[0078] 其中f(y |yi)为信号时隙抽样值的概率密度函数，f(y|yo)为非信号时隙抽样值的 概率密度函数。
[0079]步骤3将接收端每个时隙抽样值看作互相独立的随机变量，求出各自的累积分布 函数
[0080] 3a)将接收到的(N，Μ) MPPM调制信号划分为Μ个非空时隙以及N-M个空时隙，这些时 隙之间不仅独立而且同分布，设接收到的(Ν，Μ)ΜΡΡΜ信号在Ν个时隙内的抽样变量为：
[0081 ] Y=[Ysl,Ys2,··· ,YsM,Ynl,Yn2,-·· ,Υη(Ν-Μ)]
[0082] 这里为了方便表示，不考虑二进制比特流与MPPM符号之间的映射关系。其中，Ysi表 示的是Μ个信号中第i个时隙的抽样变量，Υ ηι表示的是N-M个非信号时隙中第i个的抽样变 量。则：
[0083] Ysmin=min{Ysi,Ys2, ··· ,Ysm}
[0084] Ynmax = max { Ynl, Yn2 , · · · , Yn(N-M) }；
[0085] 3b)根据步骤2,可分别求出￥51的概率密度函数，进而可计算出其对应的累计分布 函数(.0为：
[0086]
[0087] 式中，F(y I yi)表不信号时隙抽样值的累积分布函数为Q函数。
[0088] 3c)同理，Yni的累积分布函数厂(.V)为：
[0089]
[0090]式中，F(y |yo)表不非信号时隙抽样值的累积分布函数。
[0091] 步骤4求出Μ个非空时隙随机变量最小值(记为￥_^)的概率密度函数，再求出N-M 个空时隙随机变量中最大值(记为Yr?ax)的概率密度函数
[0092] 4a)由步骤3知非空时隙与空时隙抽样值的累积概率密度函数，根据概率论中求解 最大最小分布的相关知识以及Q函数的性质1 -Q ( X ) = Q ( -X )，Ysmin和Ynmax的累积分布函数 和?)可以分别表示为：
[0095] 4b)对YsmiJPY_x的累积分布函数关于y求导，可得到它们的概率密度函数乂_ ( r) 和分别为：
[0093]
[0094]
[0096]
[0097]
[0098] 步骤5求出非空时隙随机变量最小值与空时隙随机变量最大值的差值Δ Y，即Δ Y = Ysmin-Yr?ax，求出ΔΥ的概率密度函数。
[0099] 5a)将Δ Υ看成两个随机变量相WYsmin+(-Ynmax)，则Δ Υ的概率密度函数等于Ysmin 与-Y_x概率密度函数的卷积，其中-Y_x的概率密度函数/iD为：
[0100]
[0101] 5b)再将Λ (.y)与乂r 〇/)进行卷积运算，可求得ΔΥ的概率密度函数fAY(y)为：
[0102]
[0103] 由于Δ Y为非空时隙随机变量最小值与空时隙随机变量最大值的差值，可知当Δ Y 为正，则非空时隙随机变量最小值大于空时隙随机变量最大值，信号可以正确解调；反之， 当Δ Υ为负表示非空时隙随机变量最小值小于空时隙随机变量最大值，信号解调将发生错 误。
[0104] 步骤6计算ΔΥ的概率密度函数从负无穷到零处的积分，所得结果即为(Ν，Μ)ΜΡΡΜ 的误符号率。
[0105] 根据步骤5，对fAY(y)从负无穷到零对y进行积分，即
[0106]
[0107] 式中，Ps为误符号率。
[0108] 将Μ = 2，N=5，代入上式即可获得(5，2 )MPPM在软解调时的误符号率。
[0109] 进一步，为了化简上述误符号率计算过程，可以通过拟合近似代替fAy(y)的表达 式，缩小计算量，进而求得Ps。
[0110] 由于多重积分和高幂次数以及Q函数的存在，f~(y)的表达式非常复杂，导致求解 Ps的解析闭合表达式非常困难。为了方便计算误符号率表达式，引入一种函数拟合的方法 去近似代替hY(y)的表达式，从而避免了多次积分和高幂次计算，大大缩小了计算量。由于 Ps的表达式中是对fM(y)从负无穷到零这一范围进行积分，因此拟合时只需考虑这一部分 的数据，对其进行拟合，再对所获得拟合表达式进行积分，最终得到误符号率的近似值。 [0 111] 借助于MATLAB软件中的曲线拟合工具箱(curve fitting tools)完成拟合。这里， 使用高斯函数作为拟合函数对其进行拟合，高斯函数表达式为：
[0112]
\ ^ J
[0113] 其中a，b，c为需要通过拟合获得的参数。为了提高拟合的准确性与稳定性，通常使 用多个高斯函数相加进行拟合，用叠加后的表达式去表示所需拟合的曲线，其表达式为：
[0114]
[0115] 式中，&1々1，(：1为需要通过拟合获得的参数;1表示序号。
[0116] 对两项叠加的高斯函数进行积分运算，即：
[0117]
[0118] 进行变量替换，令Z = (y-bi)/ci，得：
[0119]
[0120]本发明的正确性和优点可通过以下仿真对比进一步说明：
[0121]通过MATLAB进行解析计算和蒙卡仿真，首先，确定拟合值的准确性;其次，通过理 论值和仿真值的对比证明软解调算法的正确性;最后，再与中湍条件下的硬解调误符号率 特性进行对比，凸显软解调算法的优越性。
[0122] 仿真结果
[0123] 图2、3和4分别给出了弱湍、中湍和强湍条件下(5,2)MPPM采用软解调时，误符号率 与平均信噪比之间的关系曲线。计算用到的指数韦伯信道参数均来自已有文献。图中包含 理论值、拟合值和仿真值。可以看出，三种湍流情况下，拟合值与理论值吻合很好，而且具有 较高的准确性。另外，理论值也与仿真值吻合的很好，证明了理论值推导的正确。
[0124] 图5给出了中湍条件下，（5,2)MPPM采用硬解调和软解调时误符号率与平均信噪比 的关系曲线。硬解调选取的是性能最好的动态判决门限与软解调对比。对比发现，25mm孔径 (a = 2.61，β = 1.37，η = 〇. 68)时，若要达到10-7数量级的误符号率，动态判决门限和软解调 所需的平均信噪比分别约为49.5和46dB;对于60πιπι(α = 1.52β = 4.02，η = 0.99)的孔径而 言，若要达到同等数量级的误符号率，只需31.5和28dB的平均信噪比即可。由此可知，两种 不同孔径时软解调比动态判决门限的硬解调有大约3.5dB的信噪比性能增益。因此，软解调 方法在系统性能上优于硬解调。
[0125] 以上所揭露的仅是本发明的较佳实施例而已，然而不能以此来限定本发明的权利 范围，本领域技术人员利用上述揭露的技术内容做出些许简单修改、等同变化或修饰，仍属 于本发明的保护范围之内。
【主权项】
1. 一种MPPM调制的软解调算法，其特征在于，包括下述步骤： 1) 由MATLAB中rand函数产生0、1随机序列作为仿真输入比特流，并根据映射表进行(N, M)MPPM信号调制，其中，N为一个MPPM信号块中的时隙个数，M为N个时隙中的脉冲个数； 2) 在MATLAB中由逆函数法产生符合指数韦伯分布的衰落信道系数并由randn函数产生 高斯白噪声随机数； 3) 将步骤2)产生的信道衰落系数和高斯白噪声随机数加入到步骤1)获得的(N，M)MPPM 调制信号中，得到受到高斯白噪声和信道衰落干扰的MPPM信号块； 4) 对步骤3)得到的(N，M)MPPM信号块中的每个时隙进行抽样； 5) 利用MATLAB中的MAX函数找到抽样值中的最大值，将抽样值中M个最大值对应的时隙 号判为1，其余均判为〇; 6) 利用步骤1)中（N，M)MPPM信号调制时的映射表，将步骤5)经过判决的信号块恢复为 〇、1序列，并将该序列与步骤1)中产生的信号对照，判断是否与其相同，若相同，则为正确解 调;若不相同，则错误符号个数加1; 7) 重复步骤1)到6)，仿真足够多的次数，得到错误符号个数，计算仿真的误符号率。2. 根据权利要求1所述的方法，其特征在于，所述步骤1)中，进行(N，M)MPPM信号调制， 是通过下述方法来实现的： I a)已知(N，M )MPPM调制的传信能力0为：式中，N为一个MPPM信号块中的时隙个数，M为N个时隙中的脉冲个数;设一个MPPM信号 块以b个比特信号进行调制，则其传信严+等式： 这里，b取满足不等式的最大自然数； Ib)将产生的二进制0，1随机序列以b个比特为一组，依据调制映射表进行(N，M)MPPM调 制。3. 根据权利要求1所述的方法，其特征在于，所述步骤2)中，产生符合指数韦伯分布的 信道衰落系数和高斯白噪声随机数，是通过下述方法来实现的：2a)已知指数韦伯型衰落下光强波动I的概率密度函数fEW(I)和累积分布函数Few(I)分 别为： 和 式中，形状参数β>〇,其取值与闪烁指数相关;尺度参数n>〇,其取值由β决定并与辐照 度均值相关;定义α>0为给定观测空间内，准直传播的光束与非准直传播的光束被成功接 收的平均量； 利用逆函数法，产生符合指数韦伯分布的随机数； 2b)关于加性高斯白噪声，已知平均信噪比为：式中，μ为系数，μ2 = (RNPt/M)2，R表示光电转换效率，Pt表示发射端的平均发射功率； 则加性高斯白噪声的方差σ"2可以表示为：根据得出的噪声方差，利用MATLAB中的randn函数即可产生符合高斯白噪声分布的随 机数。4. 根据权利要求1所述的方法，其特征在于，所述步骤3)中，得到受到高斯白噪声和信 道衰落干扰的MPPM信号块，是通过下述方法来实现的： 3a)根据信号被抽样之后的输出栌W -式中，由上到下依次为MPPM信号块中有脉冲的时隙和没有脉冲的时隙，η表示的是均值 为〇,方差为< =乂/2的加性高斯白噪声，No为其单边功率谱密度，h表示的是信道衰落； 3b)将步骤2)所产生的符合指数韦伯分布的N个随机数依次与一个(N，M)MPPM信号块中 N个时隙的信号相乘，之后再将产生的N个高斯白噪声随机数依次与信号相加，最终得到接 收的一个MPPM信号块。5. -种权利要求1所述的MPPM调制的软解调算法进行误符号率性能估计方法，其特征 在于，包括下述步骤： 1) 建立点对点无线光通信大气信道的数学模型； 2) 计算信号时隙与非信号时隙抽样值的概率密度函数； 3) 将接收到的（N，M)MPPM调制信号的每个时隙抽样值作为互相独立的随机变量，求出 各自的累积分布函数，在此，N为一个MPPM信号块中的时隙个数，M为N个时隙中的脉冲个数； 4) 求出M个非空时隙随机变量最小值Ys-的概率密度函数，再求出N-M个空时隙随机变 量中最大值Y?ax的概率密度函数； 5) 得到非空时隙随机变量的最小值与空时隙随机变量的最大值的差值ΛΥ，求出差值 的概率密度函数； 6) 计算△ Y的概率密度函数从负无穷到零处的积分，所得结果即为(N，M)MPPM的误符号 率。6. 根据权利要求5所述的方法，其特征在于，所述步骤2)中，计算信号时隙与非信号时 隙抽样值的概率密度函数，是通过下述方法来实现的： 根据已知指数韦伯型衰落下光强波动的概率密度函数fEw(I)和累积分布函数Few(I)及 步骤1)中的数学模型，计算出信号时隙与非信号时隙抽样值的概率密度函数分别为：式中，β>〇为形状参数，n>〇为尺度参数，定义α>0为在给定观测空间内，准直传播的 光束与非准直传播的光束被成功接收的平均量，h表示信道衰落，μ为系数，且y = RNPt/M;f (y Iy1)为有信号时隙抽样值的概率密度函数，f(y I y〇)为无信号时隙抽样值的概率密度函 数。7.根据权利要求5所述的方法，其特征在于，所述步骤3)中，求出各自的累积分布函数， 是通过下述方法来实现的： 3a)将接收到的(N，M)MPPM调制信号划分为M个非空时隙以及N-M个空时隙，这些时隙之 间不仅独立而且同分布，设接收到的(N，M)MPPM信号在N个时隙内的抽样变量为： Y= [ Ysl, Ys2-·· YsM, Ynl, Yn2-·· Yn(N-M)] 其中，Ysl表示的是M个信号中第i个时隙的抽样变量，Yni表示的是N-M个非信号时隙中第 i个的抽样变量； 则M个抽样值中的最小值Ysmi4PN-M个抽样值中的最大值Y_x分别为： Ysmin - IIlin{Ysl ,Ys2 , ,Ysm} Ynmax - max { Ynl , Yn2 , , Yn(N-M) }； 3b)根据抽样值的概率密度函数，分别求出Ysl的概率密度函数，则可计算出其对应的累 计分布函数（J)为：其中，of为加性高斯白噪声的方差，ση为加性高斯白噪声的标准差： 3函数； 3c)同理Yni的累积分布函数&()_)为：8.根据权利要求5所述的方法，其特征在于，所述步骤4)中，求解M个非空时隙随机变量 最小值Ysmin的概率密度函数和N-M个空时隙随机变量中最大值YnmaJ^概率密度函数，通过下 述方法实现： 4a)由得到的非空时隙与空时隙抽样值的累积概率密度函数，利用最大最小分布以及Q 函数的性质卜Q(x)=Q(-x)，Ys_和Y_x的累积分布函数尽_>)和及J)，)可以分别表示 为：4b WiYsmir^PYnmax的累积分布函数关于y求导，可得到其概率密度函数/}_(兄)和 分别为：9. 根据权利要求5所述的方法，其特征在于，所述步骤5)中，求出非空时隙随机变量的 最小值与空时隙随机变量的最大值的差值A Y的概率密度函数，通过下述方法实现： 5a)将Δ Y看成两个随机变量相WYsmin+(-Ynmax)，则Δ Y的概率密度函数等于1_与-Y_x概率密度函数的卷积，其中-Ynmax的概率密度函数为：5b)再将.Ziu 与/<_ ()；)进行卷积运算，可求得Δ Y的概率密度函数为：10. 根据权利要求5所述的方法，其特征在于，所述步骤6)中，得到(N，M)MPPM的误符号 率，通过下述方法实现： 根据A Y的概率密度函数求解误符号率Ps，对fM(y)从负无穷到零对y进行积分，BP :
【文档编号】H04L25/03GK105897626SQ201610225127
【公开日】2016年8月24日
【申请日】2016年4月12日
【发明人】王平, 王悦, 曹天, 杨本圣, 杨银堂
【申请人】西安电子科技大学