一种阵列天线卫星通信系统的多用户检测译码方法
【专利摘要】本发明提供一种阵列天线卫星通信系统的多用户检测译码方法，该方法用于针对阵列天线卫星多用户采用SC?FDMA多址方式接入的系统的上行链路的干扰消除和译码，所述方法为：通过将消息在“译码节点”、“映射/反映射节点”和“多用户检测节点”之间反复迭代，在达到一定迭代次数以后，进而实现阵列天线卫星通信系统多用户之间的联合干扰消除和译码，输出消除干扰的译码结果。本发明涉及阵列天线卫星通信系统多用户SC?FDMA反向链路中基于期望传播算法的联合干扰消除与译码实现方法。本发明结合因子图的结构联合处理多用户系统的干扰消除与译码。并且通过采用期望传播算法进一步降低计算复杂度。本发明具有线性计算复杂度且具有良好的干扰消除性能。
【专利说明】
-种阵列天线卫星通信系统的多用户检测谭码方法
技术领域
[0001] 本发明设及一种基于期望传播算法的阵列天线卫星通信系统SC-抑MA链路实现方 法属于通信技术领域，特别设及阵列天线卫星通信系统多用户SC-抑MA(单载波频分多址） 反向链路中基于期望传播算法化xpec化tion Propagation)的多用户联合干扰消除与译码 实现方法。
【背景技术】
[0002] 随着人们对卫星通信的通信质量需求越来越高和卫星通信技术的不断发展，卫星 通信阵列天线技术W其增益高、宽带宽、低旁瓣等优点得到越来越广泛的应用。为了能够进 一步提高系统的频谱效率和系统容量，可W在阵列天线卫星通信系统中采用SC-抑MA多址 方式，利用下行控制信号，将地面用户逻辑上划分为多个小区，同小区的用户采用不同的子 载波集合复用整个频带;不同小区采用单色频率复用或者分数频率复用，即构成了采用SC- FDMA多址方式的阵列天线卫星多用户通信系统。但当系统中小区用户密集时，小区间的干 扰将严重恶化系统性能。
[0003] 针对上述问题，应用干扰消除技术可W有效改善系统性能约束。其常见实现方法 如利用匪沈(最小均方误差)滤波检测器、最大似然检测器等模块进行。运些方法能够显著 提升系统容量，但存在计算复杂度高的缺点，如基于MMSE滤波的迭代方法，其计算复杂度W 阵列天线单元数的=次方增长。
【发明内容】
[0004] 本发明实施例提供一种基于期望传播算法的多用户联合干扰消除和译码方法，将 其应用在阵列天线卫星通信系统SC-抑MA反向链路实现中，进而实现基于期望传播的算法 实现多用户联合干扰消除与译码，相比于传统的干扰消除方法，系统复杂度得到了极大地 降低。
[0005] 第一方面，提供一种阵列天线卫星通信系统的干扰消除和译码方法，该方法用于 针对阵列天线卫星多用户采用SC-抑MA多址方式接入的系统的上行链路的多用户联合干扰 消除和译码，所述方法为:通过将消息在"译码节点"、"映射/反映射节点"和"多用户检测节 点"之间反复迭代，在达到一定迭代次数W后，进而实现阵列天线卫星通信系统多用户之间 的联合干扰消除和译码，输出消除干扰的译码结果。
[0006] 结合上述第一方面，在第一种可能的实现方式中，所述方法具体包含：
[0007] 步骤101)利用因子图中消息传递的规则，计算并更新各变量节点与函数节点间传 递的消息，并将更新的各变量节点与函数节点间传递的消息传递至多用户检测节点，再由 多用户检测节点将更新的各变量节点与函数节点间传递的消息传递至"映射/反映射"节 占 . '?、、，
[000引在"映射/反映射"节点，利用当前消息计算用户端发送的变量符号对应的各发送 比特的最大似然比LLR并将计算得到的最大似然比LLR发送至"编码/译码"节点；
[0009] "编码/译码"节点进行译码操作得到译码结果并输出更新后的编码比特的最大似 然比LLR;
[0010] 其中，所述变量节点代表各用户在各子载波上的发送符号；函数节点代表信道对 应于各天线阵元的各接收子载波上的分量；
[0011] 步骤102)判断是否满足了迭代要求，如果满足了迭代条件则结束迭代并将译码器 在最后一次迭代中的译码结果作为用户的译码结果进行输出；
[0012] 否则返回步骤101)对接收的所有用户的发送信息开始再一次迭代。
[0013] 结合上述第一方面，和第一种实现方式，在第二种可能的实现方式中所述步骤 101)之前还包含：
[0014] 步骤100)初始化与多用户检测节点相关的消息；
[001引其中，当节点间传递的消息均为高斯近似时，所述与多用户检测节点相关的消息 包含:均值和方差。
[0016] 结合上述第一方面，和第一种及第二种实现方式，在第=种可能的实现方式中，当 节点间传递的消息均为高斯近似时，所述步骤100)进一步包含：
[0017] 步骤100-1)假设诚为t时刻第n个用户发送的第k个子载波上的符号，i为迭代计算 次数，初始化过程包含对W下参数赋初值：
[001 引 i = l
[0019]
[0020]
[0021]
[0022] 其中，明_;、哉二、;,为从函数节点记传递至变量节点苗的消息在第i和第"i-1" 次迭代时高斯分布的均值；为从函数节点传递至变量节点4的消息在第i 和第"i-r次迭代时高斯分布的方差;、.?二为从变量节点4传递至信道转移函数 节点洁的消息在第巧日第"i-i"次迭代时高斯分布的均值;唯三。；、为从变量节点冶 传递至函数节点g文的消息在第i和第"i-r次迭代时的方差；
[0023] 步骤100-2)发送时用户将编码比特信息端映射为发送符号為，巧,为t时刻第n个 用户发送符号姑的第q位编码比特，所述编码比特娩取值为0或1，在第i = l次迭代中编码比 特的概率分布初始化为：
[0024]
[00巧]其中，q = l，2,3...Q，Q为系统调制阶数。
[0026] 结合上述第一方面，和第一种、第二种及第=种实现方式，在第四种可能的实现方 式中所述步骤101)进一步包含：
[0027] 步骤101-1)计算变量节点始与信道转移函数节点篇之间传递的消息片若(4)
[002引 和片吝K)，通过高斯近似的方法采用如下公式计算传递的消息的均值和方差：
[0029]
[0030]
[0031]
[0032] 其中，0<n<N+l、0<m<M+l、0<k<K+l，N为系统用户数量、M为卫星天线阵元数 量、K为子载波数目；社、自为消息自(治；)的均值和方差，为消息 的均值和方差，2技,;,、r器为第"i-i"次迭代时函数节点传递至变量节点 4消息峰芯;,(.《)的均值和方差，端一,'。、嘴一、为第i次迭代时函数节点挺传递至变量节 点碱消息的均值和方差，.?;兰;.、为第"i-r次迭代变量节点4时传递至 函数节点蘇消息/<;九;,(4)的均值和方差，旅为t时刻卫星第m个天线阵元接收的第k个子 载波上的信号，O2为信道加性高斯白噪声的方差，Wmn系统频域信道矩阵的第m行第n列的分 量；
[0033] 步骤101-2)基于更新后的均值和方差从变量节点记传递至函数节点 复文的消息表示为/<;,U;,(苗)砸过高斯近似计算更新消息/4U,;,(刮的均值和方差：
[0034]
[0035]
[0036] 其中，^_^&,^、分别为消息/的均值和方差，2；；^4,、巾为消息 片2一(苗）的均值和方差，W皿系统频域信道矩阵的第m行第n列的分量；
[0037] 步骤101-3)基于更新后的均值自和方差从函数节点結传递至变量节 点4的消息表示为片直_,;. (乂,)，通过高斯近似计算更新消息A',,)的均值和方差：
[00；3 引
[0039]
[0040] 其中，分别为消息自(4；)的均值和方差，?^自、分别为 第i次迭代变量节点诚传递至函数节点g;,消息成-二;,（乂,)的均值和方差，、墙一,',,分别 为消息乂的均值和方差，F为傅里叶变换矩阵，Fki为F中第k行第1列的元素；
[0041] 步骤101-4)变量节点《流向映射节点姑的消息表示为冷1_^4心；：，），计算消息 UJ诚对应的高斯近似在第i迭代中的均值^和方差:
[0042]
[0043]
[0044] 其中，4:一、(。、堪一,,,,,分别为消息地:一、,>;")在第i次迭代时的均值和方差，F为傅里 叶变换矩阵，Fki为F中第k行、第1列的元素；
[0045] 步骤101-5)在步骤101-4)计算结果的基础上，计算第i次迭代符号4对应编码比 特的最大似然比LLR;对于符号兹对应的第q位编码比特瑞，计算出避前最大似然比LLR:
[0046]
[0047] 其中，0<q<Q+l，Q为系统调制阶数，A为2。阶调制符号集合，Aj、分别为第q位 比特为0或1的符号子集，瑞巧第q位比特的概率分布；
[0048] 译码器W计算得到的悼(瑞)}作为输入并输出外信息批>(端亦
[0049] 对于编号q，更新编码比特墙的概率分布，其中：0<q<Q+l:
[(K)加 ]
[0化1 ]
[0化2] 其中，有〇<n<化1和0<1<K+1;
[0053]步骤101斗}利用軒>(瑞)}，计算更新映射节点拖传递向变量节点《的消息片;;;一4片：!；
[0化4]
[OO对其中，0 < n <化1、0 < 1 < K+1;诚为符号4对应的第q位比特；
[0化6]步骤101-7)基于更新后的乂 :)，对于编号n、l计算变量节点這的归一化置信 度嫂(引，并讲一步计貸锭中。的挽值巧)巧力差巧，：
[0化7]
[0化引 [0化9]
[0060] 其中，〇<n<N+l、0<1<K+1、巧)为第i次迭代中消息峰L成片。)的均值和方 差；
[0061] 步骤1〇1_7)对于编号n、k、l，更新从变量节点4传递至函数节点装:的消息片2_;,片,)， 计算消息高斯近似的均值^^^>和方差^^*。：
[0062]
[0063]
[0064] 其中，0<n<N+l、0<k<K+l、0<l<K+l ;巧>、巧i为步骤101-6计算的归一化置信 度片;:'K)的均值和方差，項一、[2-1,分别为消息(乂，)的均值和方差，切傅里叶变 换矩阵，Fki为F中第k行、第1列的元素；
[0065] 步骤101-8)对于编号〇心更新函数节点技传递至变量节点.<:的消息//^；^"。* (4)， 计算消息//占(乂,）高斯近似的均和方差;
[0066；
[0067；
[006引其中，0<n<N+l、0<k<K+l，巧,，一X;,、巧为消息杉,)高斯近似的均值和 方差，F为傅里叶变换矩阵，Fki为F中第k行、第1列的元素。
[0069] 结合第一方面，和/或第一种实现方式，和/或第二种实现方式，和/或第=种的实 现方式，和/或第四种实现方式，在第五种可能的实现方式中，根据迭代输出结果是否满足 系统需求或者迭代次数是否达到规定次数两个原则之一判断是否满足了迭代条件。
[0070] 第二方面，提供一种阵列天线卫星通信系统的多用户联合干扰消除和译码系统， 所述系统包含：
[0071] 多用户联合检测模块，用于通过迭代计算实现卫星接收信号的多用户检测W及干 扰消除;所述多用户联合检测模块的输入包括各天线阵元的接收信号，由映射/反映射处理 模块传递至该模块的消息，且该用户联合检测模块的计算输出为传递至映射/反映射处理 模块的消息；
[0072] 映射/反映射处理模块，用于实现用户编码比特与时域发送符号之间的映射/反映 射功能;其中，所述映射/反映射处理模块的输入信号为所述多用户联合检测模块传递至该 模块的消息、由译码模块计算输出的更新后的各个编码比特的最大似然比化R，所述映射/ 反映射处理模块的输出包括传递至多用户联合检测模块的消息、传递至译码模块的各个编 码比特的最大似然比LLR;
[0073] 译码模块，用于实现卫星对接收信号的译码功能;其中，所述译码模块的输入为所 述映射/反映射处理模块处理得到的时域符号各编码比特的最大似然比化R，所述译码模块 的输出为系统的译码输出、传递向所述映射/反映射处理模块的更新后的各编码比特的最 大似然比LLR。
[0074] 上述多用户联合检测模块的工作原理和过程为:信道转移函数节点乂自W及变量节 点4,其中信道转移函数节点代表了信道对应于各天线阵元的各接收子载波上的分量，而 变量节点代表了各用户在系统各子载波上的发送符号;系统迭代计算过程中消息在信道转 移函数节点W及变量节点之间相互传递。经过一定的迭代计算之后，可W求得各变量节点 的边缘概率分布，从而实现了对于接收混合信号的多用户检测W及干扰消除。所述映射/反 映射处理模块的工作原理和过程为:映射/反映射函数节点.《和变量节点4,其中函数节点 代表了各用户通过编码比特得到各时域符号的映射和反映射变换，变量节点则代表了各用 户发送的各个时域符号;借助于模块的输入消息等，模块内通过对于变量节点和函数节点 之间传递消息的迭代计算可W得到各编码比特的最大似然比化R。所述译码模块根据各用 户编码比特的最大似然比LLR，计算输出各用户译码结果即各用户发送信息。
[0075] 与现有技术相比，本发明中基于期望传播算法的阵列天线卫星通信系统SC-FDMA 链路实现方法相比于传统方法具有W下显著特征：
[0076] 多用户联合干扰消除与译码，提高了系统频谱利用率;本发明中阵列天线卫星通 信系统使用SC-抑MA多址技术，通过联合干扰消除与译码，有效消除了各小区之间的同步的 干扰，极大提高了系统对频谱资源的利用率。本发明中使用基于期望传播的算法实现联合 干扰消除与译码，相比于传统的干扰消除方法，系统复杂度得到了极大地降低。
【附图说明】
[0077] 图1基于期望传播算法的多用户联合干扰消除与译码实现流程图；
[0078] 图2基于期望传播算法的多用户联合干扰消除与译码消息传递因子图；
[0079] 图3基于期望传播算法的阵列天线卫星通信系统SC-FDMA链路实现系统结构图。
【具体实施方式】
[0080] 下面结合附图和实施例对本发明所述方法进行详细说明。
[0081 ] 实施例1
[0082] 本发明提出的基于期望传播算法的阵列天线卫星通信系统SC-抑MA链路实现方法 的流程图如附图1所示，本发明包括W下步骤：
[0083] 算法从"多用户检测节点"开始，首先初始化与其相关的消息（本方法中节点间传 递的消息都采用高斯近似来简化运算，因此只要初始化各消息相应的均值和方差即可）。之 后利用因子图中消息传递的规则，计算更新各变量节点与函数节点间传递的消息。在"映 射/反映射"节点，利用当前消息计算发送变量符号对应的各发送比特的最大似然比化R并 发送至"编码/译码"节点。"编码/译码"节点进行译码操作得到译码结果并输出更新后的编 码比特的LLR。利用新的LLR基于因子图消息传递的规则，计算更新各变量节点与函数节点 间传递的消息至"多用户检测节点"，本次迭代完成。
[0084] 实施例2
[0085] W下实施例假设用户在t时刻发送信息，接收端在t时刻接收用户发送的信息并估 计发送端在t时刻的发送信息，即W下实施例未考虑发送端与接收端传输信息的延时，而当 具体应用考虑延时信息时下述技术方案仅需要结合本领域的公知常识简单变形即可，并不 需要付出创造性的劳动。此外，下述技术方案将发送信号近似为高斯分布，当发送信号不符 合高斯分布时也可W借鉴如下技术方案实现消除干扰的联合译码的效果。
[0086] 本发明使用基于期望传播的算法，通过将消息在"译码节点"、"映射/反映射节点" 和"多用户检测节点"之间反复迭代，在达到一定迭代次数W后，本发明可W实现阵列天线 卫星通信系统多用户之间的联合干扰消除和译码。
[0087] 为了降低复杂度，规定了消息传递顺序：多用户检测节点根据输入接收信号信息 更新消息，并传递至映射节点，映射节点更新消息后传递至译码节点，一旦译码器更新了 化R(最大似然比)并将其传递至多用户检测节点，新一轮迭代开始，消息从因子图底部向上 传递并立即返回。
[0088] 本发明针对阵列天线卫星多用户系统的上行链路，其应用SC-抑MA多址方式。系统 中有N个独立用户，第n个用户在t时刻在时域发送符号记为为Sa=y"，苗,，...，请了乂 EA A 为2。阶调制符号集合，K为子载波的数量。对应的频域发送符号记为X。,=[為，4，...，诚了治 为用户n在t时刻第k个子载波上发送的信号。卫星配备了由M个天线阵元组成的天线阵列， 所有天线阵元在第k个子载波上接收到的频域信号为yk，j4为第m个天线阵元t时刻在第k个 子载波上的接收信号。
[0089] 本发明在上述系统中利用期望传播算法进行联合干扰消除和译码，第i次迭代中 从变量节点據传递至信道转移函数节点记的消息记为(為) > 而在反方向传递的 消息记为本方法中为降低计算复杂度将其近似为高斯分布
Wmn多波束多用 户频域信道矩阵的第m行第n列的分量，哉;_^和为高斯分布的均值，:和嘴为 高斯分布的方差。
[0090] 同理，第i次迭代中从变量节点冶传递至函数节点茲的消息高斯近似记为 ^
................ 3均值，*^<^>^3自为方差)，而在反方向传递的消息高 斯近似记为.j
如匀值，为方差);第i次迭代中 从变量节点這传递至函数节点茲的消息高斯近似记为鐵(《）S 乂;:(妃.繁，於蝶W;、)（ 为均值,?,?为方差），而在反方向传递的消息高斯近似记为 '
;为均值，为方差），瑪为傅里叶变换矩 阵中的元素;从变量节点4传递至映射节点短的消息近似为
((61为均值，嫂为方差）。
[0091] 一种基于期望传播算法的阵列天线卫星通信系统SC-FDMA链路实现方法特征在 于，所述的方法具体过程是依次按照W下步骤实现的：
[0092] 步骤(1)初始化：
[0093] 碱为t时刻第n个用户发送的第k个子载波上的符号，i为迭代计算次数，算法初始 化阶段i = l，设置
[0094]
[0095]
[0096]
[0097] 发送时用户将编码比特信息<L(q=l，2,3…切映射为发送符号:《，瑞,为t时刻第n 个用户发送符号4的第q个编码比特。為/取值为0、1，第i = l次迭代中其概率分布初始化为 [009引
[0099] 步骤(2)传播和更新信息：
[0100] 在规定的迭代次数内，对于时刻t的系统在第i次迭代过程中，执行W下步骤：
[0101] 步骤(2.1)对于编号n、m、k，有 0<n<N+l、0<m<M+l、0<k<K+l，计算变量节点4 传递与信道转移函数节点.谋之间传递的消息A心(.攻)/岩^、?。(為），通过高斯近似的方 法只需计算消息/<;;^；^冶)和片<^.,;>；：,)均值与方差，通过^下步骤更新：
[0102]
[0103]
[0104]
[0105]
[0106] 步骤(2.2)基于更新后的均值评方差巧>对于编号n、k、I，有0 <n<N+1、 0<1^<1(+1、0<1<1(+1，从函数节点滅传递至变量节点4的消息雌:~?;,,(乂，)，通过高斯近 似，计算更新消息堪一,Ja'u )均值方差：
[0107]
[010 引
[0109] 步骤(2.3)对于编号n、1，有0 <n <N+1、0 < 1 <K+1，变量节点攻流向映射节点始的 消息为成Ui (>,'?)，计算消息為UJ 乂I)高斯近似在迭代中的均值巧;和方差巧'，
[0110]
[0111]
[0112] 利用计算得到的均值巧>日方差Cf，计算产生符号社对应编码比特的化R。对于编 号q，有0<q<Q+l，计算出编码比特瑞，的LLR:
[0113]
[0114] 译码器化掉'(端)}作为输入并输出外信息护'(域/)}，对于编号q，有。<q<Q+l， 更新编码比特切，/的概率分布：
[0115]
[0116]
[0117] 步骤(2.4)利用批'(瑞)}对于编号n、1，有0<11<壯1、0<1<1(+1，计算更新映射 节点據传递向变景节点si的消息) _
[011 引
[0119] 步骤(2.4)基于更新后的片^、;,(乂,)>对于编号11、1，有0<11<壯1、0<1<1(+1，计算 变量节点记的扫一化晉信度货'片,）并讲一巧计貸巧挽估巧^日方差巧;i。
[0120]
[0121]
[0122]
[0123] 步骤(2.5)对于编号〇心1，有0<11<町1、0<1^<1(+1、0<1<1(+1，更新从变量节点 4传递至函数节点各:,的消息峰:-4,(乂:)，计算其高斯近似的均值巧,^',和方差巧一苗，
[0124]
[0125]
[0126] 步骤(2.6)对于编号11、4，有0<]1<壯1、0<1^<1(+1，更新函数节点备,^传递至变量节 点的消息冷U (诚)，计算其高斯近似的均值端_^和方差嘴一,t :
[0127]
[012 引
[0129] 少骤U;里发W J：少骤U;直到满足迭代次数，输出译码结果。
[0130] 此外，本发明还提供一种阵列天线卫星通信系统的多用户联合干扰消除和译码系 统，其特征在于，所述系统包含：
[0131] 多用户联合检测模块，用于通过迭代计算实现卫星接收信号的多用户检测W及干 扰消除;所述多用户联合检测模块的输入包括各天线阵元的接收信号，由映射/反映射处理 模块传递至该模块的消息，且该用户联合检测模块的计算输出为传递至映射/反映射处理 模块的消息；
[0132] 映射/反映射处理模块，用于实现用户编码比特与时域发送符号之间的映射/反映 射功能;其中，所述映射/反映射处理模块的输入信号为所述多用户联合检测模块传递至该 模块的消息、由译码模块计算输出的更新后的各个编码比特的最大似然比化R，所述映射/ 反映射处理模块的输出包括传递至多用户联合检测模块的消息、传递至译码模块的各个编 码比特的最大似然比LLR;
[0133] 译码模块，用于实现卫星对接收信号的译码功能;其中，所述译码模块的输入为所 述映射/反映射处理模块处理得到的时域符号各编码比特的最大似然比化R，所述译码模块 的输出为系统的译码输出、传递向所述映射/反映射处理模块的更新后的各编码比特的最 大似然比LLR。
[0134] 实施例3
[0135] 本发明提出的一种基于期望传播算法的阵列天线卫星通信系统SC-抑M链路实现 方法，下面W-个多用户SC-抑MA系统上行链路为例进行说明。
[0136] 考虑一个多用户的SC-抑MA卫星通信系统，系统结构如附图3所示。地面端拥有N = 32个不同用户，每个用户来自不同的小区，卫星配置了由M = 32个天线阵元组成的天线阵 列。系统采用单色频率复用，共有N=32个来自不同小区的独立用户共享某块相同的时频资 源。每一个用户使用长度为Ne = 4096,1/2码率的LDPC编码与QPSK调制，每个用户子载波数K = 16。
[0137] 多用户的SC-抑MA卫星通信系统实现关键在于多用户检测与干扰消除。本方法中 提出在接收端利用期望传播算法进行联合干扰消除与译码。传统方法则利用MMSE检测器迭 代进行干扰消除W及译码。分别利用本发明中提出的方法W及传统方法对于W上方法进行 实现，通过仿真验证了本发明方法可W在线性计算复杂度的情况下，达到良好的干扰消除 效果。
[0138] 总之本发明提出一种基于期望传播算法的阵列天线卫星通信系统多用户联合干 扰消除和译码方法。本发明结合因子图的结构联合处理多用户系统的干扰消除与译码。并 且通过采用期望传播算法进一步降低计算复杂度。本发明具有线性计算复杂度且具有良好 的干扰消除性能。
[0139] 最后所应说明的是，W上实施例仅用W说明本发明的技术方案而非限制。尽管参 照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方 案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明 的权利要求范围当中。
【主权项】
1. 一种阵列天线卫星通信系统的多用户检测译码方法，该方法用于针对阵列天线卫星 多用户采用SC-抑MA多址方式接入的系统的上行链路的多用户联合干扰消除和译码，所述 方法为:通过将消息在"译码节点"、"映射/反映射节点"和"多用户检测节点"之间反复迭 代，在达到一定迭代次数W后，进而实现阵列天线卫星通信系统多用户之间的联合干扰消 除和译码，输出消除干扰的译码结果。2. 根据权利要求1所述的阵列天线卫星通信系统的多用户联合干扰消除和译码方法， 所述方法具体包含： 步骤101)利用因子图中消息传递的规则，计算并更新各变量节点与函数节点间传递的 消息，并将更新的各变量节点与函数节点间传递的消息传递至多用户检测节点，再由多用 户检测节点将更新的各变量节点与函数节点间传递的消息传递至"映射/反映射"节点； 在"映射/反映射"节点，利用当前消息计算用户端发送的变量符号对应的各发送比特 的最大似然比LLR并将计算得到的最大似然比LLR发送至"编码/译码"节点； "编码/译码"节点进行译码操作得到译码结果并输出更新后的编码比特的最大似然比 LLR； 其中，所述变量节点代表各用户在各子载波上的发送符号；函数节点代表信道对应于 各天线阵元的各接收子载波上的分量； 步骤102)判断是否满足了迭代要求，如果满足了迭代条件则结束迭代并将译码器在最 后一次迭代中的译码结果作为用户的译码结果进行输出； 否则返回步骤101)对接收的所有用户的发送信息开始再一次迭代。3. 根据权利要求2所述的阵列天线卫星通信系统的多用户联合干扰消除和译码方法， 其特征在于，所述步骤101)之前还包含： 步骤100)初始化与多用户检测节点相关的消息； 其中，当节点间传递的消息均为高斯近似时，所述与多用户检测节点相关的消息包含： 均值和方差。4. 根据权利要求3所述的阵列天线卫星通信系统的多用户联合干扰消除和译码方法， 其特征在于，当节点间传递的消息均为高斯近似时，所述步骤100)进一步包含： 步骤100-1)假设4为t时刻第η个用户发送的第k个子载波上的符号，i为迭代计算次 数，初始化过程包含对W下参数赋初值：其中，二分别为从函数节点/nt传递至变量节点端的消息在第"i-1"次迭代 时高斯分布的均值和方差;·为从函数节点传递至变量节点诚的消息在第i 次迭代时高斯分布的均值和方差；自分别为从变量节点讀传递至函数节点gf, 的消息在第"i-r次迭代时的均值和方差； 步骤100-2)发送时用户将编码比特信息站f映射为发送符号延，鸿，为t时刻第η个用户 发送符号4的第q位编码比特，所述编码比特項,取值为0或1，在第i = l次迭代中编码比特的 概率分布初始化为：其中，q=l，2,3. . .Q，Q为系统调制阶数。5.根据权利要求4所述的阵列天线卫星通信系统的多用户联合干扰消除和译码方法， 其特征在于，所述步骤101)进一步包含： 步骤lOl-U计算变量节点4与信道转移函数节点总之间传递的消息(诚)和 请)，通过高斯近似的方法采用如下公式计算传递的消息的均值和方差：其中，0<n<N+l、0<m<M+l、0<k<K+l，N为系统用户数量、M为卫星天线阵元数量、K 为子载波数目；.巧^;,、哨为消息/堪(4 )的均值和方差，鸣Ad、为消息 嫂-r,V(诚)的均值和方差，4-1;,、rjli,为第"i-1"次迭代时函数节点记，传递至变量节点 讀消息/货(讀)的均值和方差，.?一、、嘴一4为第i次迭代时函数节点詔传递至变量节 点诚消息端一,,,(4)的均值和方差，.?:过、为第"i-1"次迭代变量节点苗时传递至 函数节点症t消息/<瓦4(冶）的均值和方差，诚为t时刻卫星第m个天线阵元接收的第k个子 载波上的信号，σ2为信道加性高斯白噪声的方差，wmn系统频域信道矩阵的第m行第η列的分 量； 步骤101-2)基于更新后的和方差从变量节点诚传递至函数节点的 消息表示为，通过高斯近似计算更新消息的均值和方差：g中，鸣^自、結分别为消，島均值芽η方^，蹲^'禮为消，島 片<,^；^、^_^4!)的均值和方差，*^^系统频域信道矩阵的第111行第11列的分量； 步骤101-3)基于更新后的均值i?^q;,和方差^54g自，从函数节点g^传递至变量节点妓的 消息表示为，(乂,），通过高斯近似计算更新消息片24,;>;")的均值和方差：其中，·分别为撤島的均值和方差，巧^^' 分别为第i次 迭代变量节点4,传递至函数节点茲消息αΓ；。。自(攻)的均值和方差，<；U.、墙分别为消 息八^_>,^"^；：,)的均值和方差^为傅里叶变换矩阵，。1^为。中第4行第1列的元素； 步骤1〇1_4)变量节点這流向映射节点玻的消息表示为成K)，计算消息成Ui (与，） 对应的高斯近似在第i迭代中的均值和方差;其中，堪一,4、墙一^"分别为消息冷U',, (乂,)在第i次迭代时的均值和方差，切傅里叶变 换矩阵，Fki为F中第k行、第1列的元素； 步骤101-5)在步骤101-4)计算结果的基础上，计算第i次迭代符号4对应编码比特的最 大似然比LLR;对于符号么对应的第q位编码比特瑞，计算出瑞的最大似然比LLR:其中，0<q<Q+l，Q为系统调制阶数，A为29阶调制符号集合，A^e、Agi分别为第q位比特 为0或1的符号子集，碱)为第q位比特的概率分布； 译码器W计算得到的·(巧，(端作为输入并输出外信息·(皆i (瑞)}; 对于编号q，更新编码比特項/的概率分布，其中:〇<q<Q+l:其中，有〇<n<N+巧Ρ0<1<Κ+1; 步骤loi-e)利用·(王…(瑞)}，计算更新映射节点《传递向变量节点4的消息诚'一心片,）；其中，0<n<N+l、0< 1 <Κ+1;瑞为符号达对应的第q位比特； 步骤101-7)基于更新后的(乂,)，对于编号n、l计算变量节点冻的归一化置信度 巧'K)，并进一步计算嫂倒的均值巧和方差嘴:其中，0<n<N+l、0<l<K+l，嘴、蝶为第i次迭代中撤島卷《(乂，)的均值和方差； 步骤101-7)对于编号n、k、l，更新从变量节点攻传递至函数节点公,的消息乂1)， 计算消高斯近似的均值.和方差其中，0<n<N+l、0<k<K+l、0<1<Κ+1;^、为步骤101-6计算的归一化置信度 巧1片）的均值和方差，墙_";,,、墙^、,，分别为消息的均值和方差，F为傅里叶变换 矩阵，Fki为F中第k行、第1列的元素； 步骤101-8)对于编号n、k，更新函数节点传递至变量节点^的消息(.诗)，计算 消息高斯近似的均和方差=其中，0<n<N+l、0<k<K+l，?^_巧5、?U?；^为消息片??一s4K)高斯近似的均值和方 差，F为傅里叶变换矩阵，Fki为F中第k行、第1列的元素。6.根据权利要求1-5中任一项所述的阵列天线卫星通信系统的多用户联合干扰消除和 译码方法，其特征在于，根据迭代输出结果是否满足系统需求或者迭代次数是否达到规定 次数两个原则之一判断是否满足了迭代条件。
【文档编号】H04L25/03GK105978662SQ201610302674
【公开日】2016年9月28日
【申请日】2016年5月9日
【发明人】吴胜, 张剑, 匡麟玲, 倪祖耀
【申请人】清华大学