《基于迭代译码和判决反馈的MAPSK信噪比估计方法.pdf》由会员分享，可在线阅读，更多相关《基于迭代译码和判决反馈的MAPSK信噪比估计方法.pdf（14页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102882652 A (43)申请公布日 2013.01.16 C N 1 0 2 8 8 2 6 5 2 A *CN102882652A* (21)申请号 201210379846.1 (22)申请日 2012.10.09 H04L 1/00(2006.01) H04L 27/38(2006.01) H04L 27/36(2006.01) (71)申请人南京六九零二科技有限公司 地址 210009 江苏省南京市鼓楼区三步两桥 12号 申请人中国人民解放军理工大学 (72)发明人吴久银 童新海 徐勇 陈伏州 李智深 (74)专利代理机构南京纵横知识产权代理有限 公司。
2、 32224 代理人董建林 (54) 发明名称 基于迭代译码和判决反馈的M-APSK信噪比 估计方法 (57) 摘要 本发明公开了一种基于迭代译码和判决反馈 的M-APSK信噪比估计方法，首先将传送信号利用 变码率译码器得到反馈结果，再调制得到不含噪 声的接收信号，并找出接收信号在星座图上的相 对位置，将其分为若干个圆环，然后通过解旋转获 得各个圆环接收信号，计算出各个圆环接收信号 的幅度均值，并根据幅度均值计算出有用信号总 功率和传送信号总功率，从而计算出信噪比的估 计值。本发明采用迭代估计的思想，利用迭代译码 输出信息比特进行逆调制和星座图分割，然后采 用数据辅助算法完成信噪比的估计，从而。
3、解决星 座图分割错误较大问题，提高信噪比估计精度，同 时避免最大似然估计法带来的实现性和复杂性问 题，具有良好的应用前景。 (51)Int.Cl. 权利要求书2页 说明书6页 附图5页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 6 页 附图 5 页 1/2页 2 1.基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特征在于：包括以下步骤， 步骤（1）将传送信号r(n)利用变码率译码器迭代译码得到反馈结果b(n)，将反馈结果 b(n)根据调制阶数再调制得到接收信号c(k)，并找出接收信号c(k)在星座图上的相对位 置，将其分为m个圆环； 步骤（。
4、2）根据步骤（1）的m个圆环的半径构成恒包络的PSK信号，计算得到各接收信 号c(k)的调制相位 及所处的圆环l，将相同圆环上的接收信号分别作为一个符号集 合，且接收信号分为m个符号集合，统计得到第l个集合的接收信号样本数目为N(l)，满足 N为接收信号的总样本数； 步骤（3）计算出各圆环幅度的均值A l 利用步骤（2）的各接收信号c(k)的调制相位 分别对各接收信号进行解旋转，分别 计算出各圆环上幅度的均值A l ； 步骤（4）计算接受信号c(k)的有用信号总功率PI 有用信号总功率PI为m个圆环上的N(l)个接收信号功率的总和，根据公式（1）得到， 其中m为圆环的总个数，A l 为对应圆环。
5、幅度的均值，l1,2，.m； 步骤（5）计算含噪声传送信号的总功率PT 根据公式（2）得到含噪声传送信号的总功率PT， 其中，N为接收信号的总样本数，|r k | 2 传送信号r(n)的为第k个传送信号的功率； 步骤（6）根据步骤（4）得到的接受信号c(k)的有用信号的总功率PI和步骤（5）得到 的含噪声传送信号的总功率PT，代入信噪比公式（3）， 即得到信噪比估计值SNR， 2.根据权利要求书1所述的基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特 征在于：步骤（1）所述传送信号r(n)为发送信号和噪声的叠加，所述各接收接收信号c(k) 为调制后的信号，不包含噪声，其中n的范围为1到。
6、N，N为接收信号的总样本数。 3.根据权利要求书1所述的基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特 征在于：所述变码率译码器采用Turbo码译码。 4.根据权利要求书1所述的基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特 征在于：所述步骤（3）计算各圆环幅度的均值A l 的步骤如下， 1）分别各圆环上的接收信号c(k)进行解旋转，计算旋转后各接收信号c(k)的幅值； 权 利 要 求 书CN 102882652 A 2/2页 3 2）各圆环幅度的均值A l 为各接收信号c(k)的幅值的平均值，根据公式（5）得到， (5) 其中，r k 为传送信号r(n)中处在第l个圆环上。
7、的传送信号，各接收信号c(k)的调制相 位 为传送信号的噪声相位。 权 利 要 求 书CN 102882652 A 1/6页 4 基于迭代译码和判决反馈的 M-APSK 信噪比估计方法 技术领域 0001 本发明涉及卫星通信技术领域，具体涉及一种用于提高卫星通信的传输速率和频 谱利用率的基于迭代译码和判决反馈的M-APSK信噪比估计方法。 背景技术 0002 随着高速数据业务需求的增加，卫星通信的发展趋势也将采用自适应编码调制技 术，自适应编码调制技术（ACM）是通过改变编码码率和调制的阶数，使系统在功率恒定时提 高传输速率，但自适应调制技术需要要求较高精度的信噪比估计算法，且估计的精度对调 。
8、制编码方式的选择影响重大，如图1所示的自适应编码调制技术的整体框图，由变码率译 码器的译码结果输入编码调制方案选择器后，通过反馈信道，控制发送端的编码器及调制 器。目前，比较可行的编码方案有Turbo码、LDPC码等，调制方式有M-PSK、M-QAM和M-APSK 调制等，在设计适合卫星信道的调制星座时，要求频谱利用率高，尽量减少信号幅度的起 伏，APSK(幅相键控)调制能使星座形状呈圆形、圆周个数少的调制方式，成为极具潜力的 调制方式，APSK调制对信号的幅度和相位同时产生影响，用多个同心圆环上的不同相位的 点表示信息比特。显然，相同进制数的APSK调制的信号可以由不同的圆环和相位组合方式 。
9、获得，较好星座排布方式能够使系统在相同的发射功率下获得更低的误码率，如图2所示， 为经优化为16APSK调制星座图，APSK调制的可分为内外若干个同心圆环，因此，现有技术 中大容量的卫星系统多采用基于APSK调制方式的自适应编码调制技术方案，但是其的关 键技术之一是信噪比精确估计，针对APSK调制信号信噪比估计的解决方案不多，均存在各 自的缺点，具体如下： 0003 1）基于数据辅助方法(Data-Aided,DA)，在DVB-S2标准中采用在有效载荷之前插 入导频信号辅助信噪比估计和信道增益估计，具有较强的实用性、估计性能好，但是插入的 导频信号，不但增加了系统开销，还降低了传输效率； 00。
10、04 2）非数据辅助方法(Non-Data-Aided,NDA)，如最大似然算法和矩方法(Memont ethod,MM)，最大似然算法低信噪比下性能较差，矩方法高信噪比下性能较差； 0005 3）适合多进制振幅移相键控信号信噪比估计的迭代算法，文中定义为MLA算法 (当先验信息为零时定义为ML算法)，在给定先验信息的条件下，由经过平坦衰落的多进制 振幅移相键控信号的概率分布出发，根据迭代最大期望(Expectation Maximum,EM)准则， 由信号最大似然函数的对数期望求导得到信噪比参数和数据的关系，引入数据可信度更新 信噪比参数，从而推导出迭代信噪比估计算法，在高信噪比下性能优于矩。
11、方法，该方法利用 译码得到的发送信息比特的先验信息来辅助实现最大似然信噪比估计，但该算法复杂度很 大，且由于基于先验信息的信噪比估计与先验信息可信度密切相关，对同样的调制星座图， 不同的星座标识和译码器会影响到估计性能，算法实现与接收机其它模块密切相关，需要 预先对验信息建模，实用性较差； 0006 4）基于判决反馈的集分分割信噪比估计算法，即利用判决输出比特将接收信号星 座图进行分割成若干个独立的PSK星座图，再通过M2M4矩方法实现信噪比估计，能够解决 说 明 书CN 102882652 A 2/6页 5 M2M4算法对非恒包络调制信号的信噪比估计问题，但由于基于解调硬判决输出的结果在信 。
12、噪比较低时误码较大，导致星座图分割错误较大，严重影响低信噪下的估计精度。 发明内容 0007 本发明的目的为了克服现有技术中的现有技术中针对APSK调制信号信噪比估计 的解决方法，系统开销大，传输效率低或、实用较差、信噪的估计精度较差等问题。本发明采 用迭代估计的思想，利用迭代译码输出信息比特进行逆调制和星座图分割，然后采用数据 辅助算法完成信噪比的估计，从而解决星座图分割错误较大问题，提高信噪比估计精度，同 时避免最大似然估计法带来的实现性和复杂性问题，具有良好的应用前景。 0008 为了解决上述技术问题，本发明的解决方案是： 0009 一种基于迭代译码和判决反馈的M-APSK信噪比估计方法。
13、，其特征在于：包括以下 步骤， 0010 步骤（1）将传送信号r(n)利用变码率译码器迭代译码得到反馈结果b(n)，将反馈 结果b(n)根据调制阶数再调制得到接收信号c(k)，并找出接收信号c(k)在星座图上的相 对位置，将其分为m个圆环； 0011 步骤（2）根据步骤（1）的m个圆环的半径构成恒包络的PSK信号，计算得到各接 收信号c(k)的调制相位及所处的圆环l，将相同圆环上的接收信号分别作为一个符号集 合，且接收信号分为m个符号集合，统计得到第l个集合的接收信号样本数目为N(l)，满足 N为接收信号的总样本数； 0012 步骤（3）计算出各圆环幅度的均值A l 0013 利用步骤（2）的。
14、各接收信号c(k)的调制相位分别对各接收信号进行解旋转， 分别计算出各圆环上幅度的均值A l ； 0014 步骤（4）计算接受信号c(k)的有用信号总功率PI 0015 有用信号总功率PI为m个圆环上的N(l)个接收信号功率的总和，根据公式（1）得 到， 0016 0017 其中m为圆环的总个数，A l 为对应圆环幅度的均值，l1,2，.m； 0018 步骤（5）计算含噪声传送信号的总功率PT 0019 根据公式（2）得到含噪声传送信号的总功率PT， 0020 0021 其中，N为接收信号的总样本数，|r k | 2 传送信号r(n)的为第k个传送信号的功率； 0022 步骤（6）根据步骤（4。
15、）得到的接受信号c(k)的有用信号的总功率PI和步骤（5） 得到的含噪声传送信号的总功率PT，代入信噪比公式（3）， 0023 0024 即得到信噪比估计值SNR， 说 明 书CN 102882652 A 3/6页 6 0025 0026 前述的基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特征在于：步骤 （1）所述传送信号r(n)为发送信号和噪声的叠加，所述各接收接收信号c(k)为调制后的信 号，不包含噪声，其中n的范围为1到N，N为接收信号的总样本数； 0027 前述的基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特征在于：所述变 码率译码器采用Turbo码译码。 0028。
16、 前述的基于迭代译码和判决反馈的M-APSK信噪比估计方法，其特征在于：所述步 骤（3）计算各圆环幅度的均值A l 的步骤如下， 0029 1）分别各圆环上的接收信号c(k)进行解旋转，计算旋转后各接收信号c(k)的幅 值； 0030 2）各圆环幅度的均值A l 为各接收信号c(k)的幅值的平均值，根据公式（5）得到， 0031 (5) 0032 0033 其中，r k 为传送信号r(n)中处在第l个圆环上的传送信号，各接收信号c(k)的调 制相位为传送信号的噪声相位。 0034 本发明的有益效果是：本发明的基于迭代译码和判决反馈的M-APSK信噪比估计 方法，首先将传送信号利用变码率译码器得。
17、到反馈结果，再调制得到不含噪声的接收信号， 并接收信号在星座图上的找出相对位置，将其分为若干个圆环，然后通过解旋转获得各个 圆环接收信号，计算出各个圆环接收信号的幅度均值，并根据幅度均值计算出有用信号总 功率和传送信号总功率，从而计算出信噪比的估计值。本发明采用迭代估计的思想，利用 迭代译码输出信息比特进行逆调制和星座图分割，然后采用数据辅助算法完成信噪比的估 计，从而解决星座图分割错误较大问题，提高信噪比估计精度，同时避免最大似然估计法带 来的实现性和复杂性问题，具有良好的应用前景。 附图说明 0035 图1是自适应编码调制系统的系统框图。 0036 图2是加噪的16-APSK调制后分割的星。
18、座图。 0037 图3是本发明的基于迭代译码和判决反馈的M-APSK信噪比估计方法的实施图。 0038 图4是本发明的1/2Turbo码率下的APSK调制的性能的仿真图。 0039 图5是本发明的不同进制调制下与DA方式的性能拟合情况的仿真图。 0040 图6是本发明的不同进制调制下的信噪比估计的误差性能的仿真图。 0041 图7是本发明的不同进制调制下的信噪比估计值的标准差的仿真图。 0042 图8是本发明的不同进制调制下的发送一帧数据的信噪比估计性能的仿真图。 说 明 书CN 102882652 A 4/6页 7 具体实施方式 0043 下面将结合说明书附图，对本发明作进一步的说明。 00。
19、44 本发明的基于迭代译码和判决反馈的M-APSK信噪比估计方法，首先将传送信号 利用变码率译码器得到反馈结果，再调制得到不含噪声的接收信号，并接收信号在星座图 上的找出相对位置，将其分为若干个圆环，然后通过解旋转获得各个圆环接收信号，计算出 各个圆环接收信号的幅度均值，并根据幅度均值计算出有用信号总功率和传送信号总功 率，从而计算出信噪比的估计值，采用迭代估计的思想，利用迭代译码输出信息比特进行逆 调制和星座图分割，然后采用数据辅助算法完成信噪比的估计，从而解决星座图分割错误 较大问题，提高信噪比估计精度，同时避免最大似然估计法带来的实现性和复杂性问题，如 图3所示，其具体实现步骤以下步骤：。
20、 0045 第一步，将传送信号r(n)利用变码率译码器迭代译码得到反馈结果b(n)，将反馈 结果b(n)再根据调制阶数M调制得到接收信号c(k)，并找出接收信号c(k)在星座图上的 相对位置，将其分m个圆环，m为圆环的数量，其中传送信号r(n)为发送信号和噪声的叠加， 各接收接收信号c(k)为调制后的信号，不包含噪声，变码率译码器采用Turbo码译码，这里 n的范围为1到N，N为接收信号的总样本数； 0046 第二步，根据第一步的m个圆环半径构成的恒包络的PSK信号，计算得到各接收信 号c(k)的调制相位及所处的圆环l，将相同圆环上的接收信号分别作为一个集合，所有 接收信号分为m个符号集合，统。
21、计得到第l个集合的接收信号数目为N(l)，满足 N为接收信号的总样本数； 0047 第三步，计算出各圆环幅度的均值A l 0048 1)分别各圆环上的接收信号c(k)进行解旋转，计算旋转后各接收信号c(k)的幅 值； 0049 2)各圆环幅度的均值A l 为各接收信号c(k)的幅值的平均值，根据公式(5)得到， 0050 (5) 0051 0052 其中，r k 为传送信号r(n)中处在第l个圆环上的传送信号，各接收信号c(k)的调 制相位为传送信号的噪声相位。 0053 第四步，计算接受信号c(k)的有用信号的总功率PI 0054 根据第三步得到的各圆环幅度的均值A l ，计算有用信号的总功。
22、率PI为m个圆环上 的N(l)个接收信号功率的总和，根据公式(1)得到， 0055 0056 其中m为圆环的总个数，A l 为对应圆环幅度的均值，l1，2，.m； 0057 第五步，计算含噪声传送信号的总功率PT 0058 根据公式(2)得到含噪声传送信号的总功率PT， 0059 说 明 书CN 102882652 A 5/6页 8 0060 其中，N为传送信号r(n)的信号样本数量，|r k | 2 传送信号r(n)的为第k个传送 信号的功率； 0061 第六步，根据第四步得到的接受信号c(k)的有用信号总功率PI和第五步得到的 传送信号总功率PT，代入信噪比估计公式(3)， 0062 00。
23、63 即得到信噪比估计值SNR， 0064 0065 根据公式（4），即可计算出信噪比估计值SNR。 0066 下面结合推理方式的推导出本发明的基于迭代译码和判决反馈的M-APSK信噪比 估计方法的优点，利用计算机仿真的方法给出其估计性能， 0067 仿真采用M进制自适应APSK调制，结合1/2码率的Turbo码，每帧帧长为 1200bits，迭代4次后输出，并将这些比特反馈回信噪比的估计器，一个APSK调制阶数可 以由多种不同的圆环和相位的组合实现，采用经优化后的组合方式，具体为：8APSK：4+4； 16APSK：4+12；32APSK：4+12+16；64APSK：4+12+20+28；。
24、 0068 M-APSK调制结合1/2码率迭代4次Turbo码的误码率性能曲线，如图4所示，若系 统要求传输的译码器输出误比特率小于某一预定值，如10-6，根据图中误码曲线就可以确 定各阶调制方式的门限比特信噪比，在这些门限值附近的信噪比估计值，将影响不同调制 进制数之间的切换，进一步影响系统的传输速率。 0069 卫星通信系统为可靠传输，一般要求误比特率小于10 -3 ，随着信噪比变化，系统在 各种进制APSK信号之间切换，如图5所示，给出了不同进制APSK进制数调制的信噪比与DA 方式的性能拟合情况，如图6所示，给出了各调制进制方式下的信噪比估计的误差性能，从 图5、6中可以看出，当达到各。
25、种调制阶数的切换门限后，AMC系统采用各种进制APSK调制， 此时信噪比估计值SNR的估计误差均小于0.2dB。受计算误差的影响和限制，图中曲线出现 一定程度的抖动，且各误差曲线有交叉，最终误差值小于0.1dB。 0070 如图7所示，给出了算法估计的信噪比偏离均值的程度，即估计误差的标准差曲 线，信噪比较小时，误码得不到纠正，信噪比估计值SNR估计与真实值较大误差，但误差值 变化不大(即标准差较小)，随着信噪比的增大，误码将可能够得到纠正，且能被纠正的概 率不增加，在变码率译码器输出误码率为10-3时，信噪比估计值受误码比特能否被纠正的 影响，波动较大，这是图7中的峰值（即估计值偏离期望值较。
26、大）产生的主要原因。 0071 如图8所示，给出了各调制进制数下发送一帧数据的信噪比估计性能，可以看出 在AMC系统切换到各调制方式后，各调制下的瞬时估计误差小于0.4dB。因此，估计方案可 以用于突发帧传输的SNR估计，图8中的虚线为各调制切换点的信噪比门限值，此时的信道 误码率大于10 -1 。 0072 综上所述，本发明提出了的基于迭代译码和判决反馈的M-APSK信噪比估计方法， 当信噪比大于1.5dB后，在各进制的APSK调制的工作模式下切换，在低信噪比条件下，各工 作模式所对应的估计误差均小于0.2dB，能够大大提高信噪比估计精度，同时避免最大似然 说 明 书CN 102882652。
27、 A 6/6页 9 估计法带来的实现性和复杂性问题，具有良好的应用前景。 0073 以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该 了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原 理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进 都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界 定。 说 明 书CN 102882652 A 1/5页 10 图1 图2 说 明 书 附 图CN 102882652 A 10 2/5页 11 图3 图4 说 明 书 附 图CN 102882652 A 11 3/5页 12 图5 说 明 书 附 图CN 102882652 A 12 4/5页 13 图6 说 明 书 附 图CN 102882652 A 13 5/5页 14 图7 图8 说 明 书 附 图CN 102882652 A 14 。