CN104539397A - 交分多址接入中继系统的去噪互信息保持量化转发方法 - Google Patents

Abstract

本发明公开了一种交分多址接入中继系统的去噪互信息保持量化转发方法，包括以下步骤：一、信源对发送信号进行编码、扩频、交织和BPSK调制处理，二、信源将调制后的信号经过信道发送给中继站和目的节点，三、中继站对其接收到的信号进行量化、编码和调制后转发给目的节点，四、目的节点联合译码。本发明设计新颖合理，方法步骤简单，能够很好地处理多址接入中继系统中多个信源量化转发，在译码端进联合译码，并且能够获得更好的性能。

Description

交分多址接入中继系统的去噪互信息保持量化转发方法

技术领域

本发明属于无线通信技术领域，具体涉及一种交分多址接入中继系统的去噪互信息保持量化转发方法。

背景技术

半双工交分多址接入中继系统(ID-MARC)，它包括多用户，一个中继和一个目的端。各用户之间是非正交的，即在相同的时域和频域发送信号，系统通过交织器来区分各个用户发送的信息。中继处常用的处理方法有两种：压缩转发(CF:Compress Forward)和译码转发(DF:DecodingForward)。译码转发需要中继对收到信号进行译码，再转发给目的端，这实际上是一个去噪的过程，在中继能够进行无差错译码的时候，译码转发是最优的策略。但译码转发也有它的缺点：一是要在中继处进行一次译码，它的处理复杂很高；其次如果中继处译码有错会存在错误传播。因此，很多情况下选用半双工的多址接入中继信道的CF策略，CF策略中中继分两个时隙收发送信号，相比DF策略，它不需要进行译码，因而复杂度不会很大，且不会有错误传播，在信源到中继站链路信道条件相对较差的时候，CF策略要比DF策略转发更有效，在相同条件下，CF策略能获得更大的速率。但在CF策略中的量化过程就会涉及到失真度量的问题。Rave提出了最大互信息量化方案(MMIQ)，对接收信号的似然比进行最大互信息量化，得到了比幅度值量化更好的性能，由于均方误差比较直观、简单而且分析起来方便，目前最常用于设计量化器和实随机变量的率失真问题的失真度量是均方误差，但对一些问题它缺乏感性意义。Tishby等提出用量化输出变量与另一个变量的相关性作为失真度量。Zeitler以此为基础提出一种以量化序列和信源发送序列之间互信息作为失真的度量，并提出以最大化两者间的互信息来设计量化器，并给出了二用户正交多址接入中继信道的量化方案。但是，到目前为止，大多数关于多址接入中继信道文章都是关于正交多址接入中继信道，或者只是考虑两个用户的情况，对用户数大于2的非正交的多址接入中继信道模型研究很少。而传统的Lloyd-Max量化方案和最大互信息量化方案(MMIQ)在量化的时候都无法做到尽可能地保持与信源信息的互信息量。

发明内容

本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种交分多址接入中继系统的去噪互信息保持量化转发方法，其设计新颖合理，方法步骤简单，能够很好地处理多址接入中继系统中多个信源量化转发，在译码端进联合译码，并且能够获得更好的性能。

为解决上述技术问题，本发明采用的技术方案是：一种交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于该方法包括以下步骤：

步骤一、信源对发送信号进行编码、扩频、交织和BPSK调制处理：多个信源中的各个信源i∈{1,2,...,K}各自随机产生独立的信息序列在发送端对信息序列以码率k_i/n_i进行Turbo-Hadamard编码得到码比特再经过扩频、交织和BPSK调制得到信号其中，K为信源的个数，k_i为第i个信源的序列长度，n_i为对第i个信源进行Turbo-Hadamard编码后的序列长度；

步骤二、各个信源将信号经过信道发送给中继站，中继站接收到信号其中，为零均值高斯白噪声，s为信源，r为中继站，h_i,r为信源到中继站各信道的乘性系数；

同时，各个信源将信号经过信道发送给目的节点，目的节点接收到信号其中，为零均值高斯白噪声，s为信源，d为目的节点，h_i,d为信源到目的节点各信道的乘性系数；

步骤三、中继站对其接收到的信号进行量化、编码和调制后转发给目的节点：中继站依据尽可能保留接收到的信号中关于信源部分的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息的原则选择出优化的量化区间并根据量化器函数q(y)＝z对其接收到的信号进行量化处理得到量化序列且对量化序列进行信源编码和信道编码得到编码序列c_r，再对编码序列c_r进行调制得到发送序列后发送给目的节点，目的节点接收到信号其中，y为信号的一个实现，z是元素集合Z＝{0，1,2,…,N-1}中的一个元素，N为量化电平数，量化器函数q(y)＝z将N个量化区间[y_min,y₁],[y₁,y₂],…,[y_N-1,y_max]与元素集合Z＝{0，1,2,…,N-1}中的元素0，1,2,…,N-1一一对应，y_min＝y₀<y₁<y₂<…<y_N-1<y_N＝y_max，y_min为量化区间边界的最小取值且第0个量化区间边界的取值y₀＝y_min，y_max为量化区间边界的最大取值且第N个量化区间边界的取值y_N＝y_max；n_r为对编码序列c_r进行调制后的序列长度，为零均值高斯白噪声，h_r,d为中继站到目的节点各信道的乘性系数；

步骤四、目的节点联合译码，具体过程为：

步骤401、目的节点对接收到的信号进行信道译码和信源译码，恢复出量化序列

步骤402、目的节点的多用户检测器联合进行信号检测和译码，恢复出信源信息序列U_i，i∈{1,2,...,K}，具体过程为：

步骤4021、第一信号估计器ESE1和第二信号估计器ESE2分别由其先验信息和接收到的和逐比特联合估计输出外信息，其中，第一信号估计器ESE1的先验信息为译码交织送过来的信号估计器的先验信息加上第二信号估计器ESE2最近一次输出的外信息e_ESE2(x_i)，第二信号估计器ESE2的先验信息为译码交织送过来的信号估计器的先验信息加上第一信号估计器ESE1最近一次输出的外信息e_ESE1(x_i)；

步骤4022、对第一信号估计器ESE1输出的外信息和第二信号估计器ESE2输出的外信息求和得到各信源的联合估计外信息e_ESE(x_i)；

步骤4023、对各信源的联合估计外信息e_ESE(x_i)进行解交织得到译码器的先验信息送给用户译码器译码，用户译码器译码后输出译码外信息e_DEC(x_i)；

步骤4024、用户译码器对译码外信息e_DEC(x_i)进行交织得到信号估计器的先验信息并送给第一信号估计器ESE1和第二信号估计器ESE2；

步骤4025、重复步骤4021～4024迭代给定次数后停止，在用户译码器处进行硬判决输出作为译码结果，恢复出信源1信息序列U_i，i∈{1,2,...,K}。

上述的交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于：步骤三中中继站依据尽可能保留接收到的信号中关于信源部分的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息的原则选择出优化的量化区间的具体过程为：

步骤301、参数输入和计算：给中继站输入信号的概率密度函数p(x)、已知信道条件下信号的概率密度函数p(y|x)、量化电平数N和互信息变化阈值ε，中继站根据公式计算出信号的概率密度函数p(y)；其中，x为信号的一个实现；

步骤302、初始化：随机选择一个满足表达式y_min＝y₀<y₁<y₂<…<y_N-1<y_N＝y_max的量化区间中继站根据公式计算得到初始迭代的输出序列与信源发送序列之间的互信息I^(old)；

步骤303、中继站根据公式 $I^{\left(\mathrm{new}\right)}=\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 计算得到当前迭代的输出序列与信源发送序列之间的互信息I^(new)，其中，j∈{1,2,…,N-1}，p(x,y)为信号和信号的联合概率密度函数且p(x,y)＝p(x)p(y|x)；

步骤304、当|I^(new)-I^(old)|≤ε时，执行步骤308；否则，当|I^(new)-I^(old)|>ε时，执行步骤305；

步骤305、将I^(old)更新为I^(new)；

步骤306、对每个j∈{1,2,…,N-1}，保持y₁,y₂,…,y_j-1,y_j+1…,y_N-1不变，中继站根据公式 $y_j=\underset{y=[y_1,y_2,...,y_{N-1}]}{\arg \max }\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 求解量化区间内的分量y_j，对步骤302中选择的量化区间 $\stackrel{\&RightArrow;}{y}=[y_1,y_2,...,y_{N-1}]$ 进行更新；

步骤307、中继站根据公式 $I^{\left(\mathrm{new}\right)}=\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 计算得到当前迭代的输出序列与信源发送序列之间的互信息I^(new)后，返回步骤304；

步骤308、输出量化区间

上述的交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于：步骤402中所述迭代给定次数为20～40次。

上述的交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于：所述零均值高斯白噪声零均值高斯白噪声和零均值高斯白噪声均为相互独立的零均值高斯白噪声。

本发明与现有技术相比具有以下优点：

1、本发明设计新颖合理，方法步骤简单且使用效果好。

2、本发明提出了一种新的目的节点对信源发送来的信号和中继站发送来的信号进行检测和译码的方法，由第一信号估计器ESE1和第二信号估计器ESE2构成一个逐比特软信息输出检测器，在做到低复杂度的同时能够获得接近最优的检测性能。

3、本发明基于最大化量化序列和信源发送序列之间的互信息的思想而提出，能很好地反应出中继器的量化在保留接收信息的同时去掉接收变量中的噪声，对于非正交的交分多址接入中继信道，本发明的量化方法能够很好地处理多个信源量化转发，在译码端进联合译码，并且能够获得更好的性能。

综上所述，本发明设计新颖合理，方法步骤简单，能够很好地处理多址接入中继系统中多个信源量化转发，在译码端进联合译码，并且能够获得更好的性能。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明交分多址接入中继系统的结构示意图。

图2为本发明交分多址接入中继系统的去噪互信息保持量化转发方法的方法流程图。

附图标记说明:

1—信源；2—中继站；3—目的节点。

具体实施方式

如图1和图2所示，本发明的交分多址接入中继系统的去噪互信息保持量化转发方法，包括以下步骤：

步骤一、信源对发送信号进行编码、扩频、交织和BPSK调制处理：多个信源1中的各个信源1i∈{1,2,...,K}各自随机产生独立的信息序列在发送端对信息序列以码率k_i/n_i进行Turbo-Hadamard编码得到码比特再经过扩频、交织和BPSK调制得到信号其中，K为信源1的个数，k_i为第i个信源1的序列长度，n_i为对第i个信源1进行Turbo-Hadamard编码后的序列长度；

步骤二、各个信源1将信号经过信道发送给中继站2，中继站2接收到信号其中，为零均值高斯白噪声，s为信源，r为中继站，h_i,r为信源1到中继站2各信道的乘性系数；

同时，各个信源1将信号经过信道发送给目的节点3，目的节点3接收到信号其中，为零均值高斯白噪声，s为信源，d为目的节点，h_i,d为信源1到目的节点3各信道的乘性系数；

步骤三、中继站对其接收到的信号进行量化、编码和调制后转发给目的节点：中继站2依据尽可能保留接收到的信号中关于信源部分的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息的原则选择出优化的量化区间并根据量化器函数q(y)＝z对其接收到的信号进行量化处理得到量化序列且对量化序列进行信源编码和信道编码得到编码序列c_r，再对编码序列c_r进行调制得到发送序列后发送给目的节点3，目的节点3接收到信号其中，y为信号的一个实现，z是元素集合Z＝{0，1,2,…,N-1}中的一个元素，N为量化电平数，量化器函数q(y)＝z将N个量化区间[y_min,y₁],[y₁,y₂],…,[y_N-1,y_max]与元素集合Z＝{0，1,2,…,N-1}中的元素0，1,2,…,N-1一一对应，N个量化区间[y_min,y₁],[y₁,y₂],…,[y_N-1,y_max]依据量化区间取得，y_min＝y₀<y₁<y₂<…<y_N-1<y_N＝y_max，y_min为量化区间边界的最小取值且第0个量化区间边界的取值y₀＝y_min，y_max为量化区间边界的最大取值且第N个量化区间边界的取值y_N＝y_max；n_r为对编码序列c_r进行调制后的序列长度，为零均值高斯白噪声，h_r,d为中继站2到目的节点3各信道的乘性系数；

量化器函数q(y)＝z设计的目标是要最大化量化输出序列和信源1发送序列之间的互信息，本发明中的量化方法依据“尽可能保留关于接收到的信号的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息”的原则而提出，这正是人们一直要找的中继处理方案，而现有技术中的量化器只是尽可能地保留接收变量Y的信息，而没有去掉接收变量中关于噪声部分的信息，在保留接收变量信息的同时也保留了大部分的信道噪声。

步骤四、目的节点联合译码，具体过程为：

步骤401、目的节点3对接收到的信号进行信道译码和信源译码，恢复出量化序列

步骤402、目的节点3的多用户检测器联合进行信号检测和译码，恢复出信源1信息序列U_i，i∈{1,2,...,K}，具体过程为：

步骤4021、第一信号估计器ESE1和第二信号估计器ESE2分别由其先验信息和接收到的和逐比特联合估计输出外信息，其中，第一信号估计器ESE1的先验信息为译码交织送过来的信号估计器的先验信息加上第二信号估计器ESE2最近一次输出的外信息e_ESE2(x_i)，第二信号估计器ESE2的先验信息为译码交织送过来的信号估计器的先验信息加上第一信号估计器ESE1最近一次输出的外信息e_ESE1(x_i)；

步骤4022、对第一信号估计器ESE1输出的外信息和第二信号估计器ESE2输出的外信息求和得到各信源1的联合估计外信息e_ESE(x_i)；

步骤4023、对各信源1的联合估计外信息e_ESE(x_i)进行解交织得到译码器的先验信息送给用户译码器译码，用户译码器译码后输出译码外信息e_DEC(x_i)；

步骤4024、用户译码器对译码外信息e_DEC(x_i)进行交织得到信号估计器的先验信息并送给第一信号估计器ESE1和第二信号估计器ESE2；

步骤4025、重复步骤4021～4024迭代给定次数后停止，在用户译码器处进行硬判决输出作为译码结果，恢复出信源1信息序列U_i，i∈{1,2,...,K}。

以上步骤中，步骤一和步骤二是信号收发的第一阶段，步骤三和步骤四是信号收发的第二阶段。

本实施例中，步骤三中中继站2依据尽可能保留接收到的信号中关于信源部分的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息的原则选择出优化的量化区间的具体过程为：

步骤301、参数输入和计算：给中继站2输入信号的概率密度函数p(x)、已知信道条件下信号的概率密度函数p(y|x)、量化电平数N和互信息变化阈值ε，中继站2根据公式计算出信号的概率密度函数p(y)；其中，x为信号的一个实现；ε具体解释为当前迭代的输出序列与信源1发送序列之间的互信息与上一次迭代的输出序列与信源1发送序列之间的互信息的变化阈值，当前迭代为首次迭代时，上一次迭代为初始迭代；

步骤302、初始化：随机选择一个满足表达式y_min＝y₀<y₁<y₂<…<y_N-1<y_N＝y_max的量化区间中继站2根据公式计算得到初始迭代的输出序列与信源1发送序列之间的互信息I^(old)；

步骤303、中继站2根据公式 $I^{\left(\mathrm{new}\right)}=\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 计算得到当前迭代的输出序列与信源1发送序列之间的互信息I^(new)，其中，j∈{1,2,…,N-1}，p(x,y)为信号和信号的联合概率密度函数且p(x,y)＝p(x)p(y|x)；

步骤304、当|I^(new)-I^(old)|≤ε时，执行步骤308；否则，当|I^(new)-I^(old)|>ε时，执行步骤305；

步骤305、将I^(old)更新为I^(new)；

步骤306、对每个j∈{1,2,…,N-1}(即j为元素集合{1,2,…,N-1}中的元素)，保持y₁,y₂,…,y_j-1,y_j+1…,y_N-1不变，中继站2根据公式 $y_j=\underset{y=[y_1,y_2,...,y_{N-1}]}{\arg \max }\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 求解量化区间内的分量y_j，对步骤302中选择的量化区间 $\stackrel{\&RightArrow;}{y}=[y_1,y_2,...,y_{N-1}]$ 进行更新；

步骤307、中继站2根据公式 $I^{\left(\mathrm{new}\right)}=\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 计算得到当前迭代的输出序列与信源1发送序列之间的互信息I^(new)后，返回步骤304；

步骤308、输出量化区间

本实施例中，步骤402中所述迭代给定次数为20～40次。

本实施例中，所述零均值高斯白噪声零均值高斯白噪声和零均值高斯白噪声均为相互独立的零均值高斯白噪声。

综上所述，基于最大化量化序列和信源发送序列之间的互信息的思想，本发明提出了一种应用于交分多址接入中继系统的全新的去噪互信息保持量化迭代算法，能很好地反应出中继器的量化在保留接收信息的同时去掉接收变量中的噪声，对于非正交的交分多址接入中继信道，本发明的量化方法能够很好地处理多个信源量化转发，在译码端进联合译码，并且复杂度也不是很大，能够获得更好的性能。

以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (4)

1.一种交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于该方法包括以下步骤：

步骤一、信源对发送信号进行编码、扩频、交织和BPSK调制处理：多个信源(1)中的各个信源(1)i∈{1,2,...,K}各自随机产生独立的信息序列在发送端对信息序列以码率k_i/n_i进行Turbo-Hadamard编码得到码比特再经过扩频、交织和BPSK调制得到信号其中，K为信源(1)的个数，k_i为第i个信源(1)的序列长度，n_i为对第i个信源(1)进行Turbo-Hadamard编码后的序列长度；

步骤二、各个信源(1)将信号经过信道发送给中继站(2)，中继站(2)接收到信号其中，为零均值高斯白噪声，s为信源，r为中继站，h_i,r为信源(1)到中继站(2)各信道的乘性系数；

同时，各个信源(1)将信号经过信道发送给目的节点(3)，目的节点(3)接收到信号其中，为零均值高斯白噪声，s为信源，d为目的节点，h_i,d为信源(1)到目的节点(3)各信道的乘性系数；

步骤三、中继站对其接收到的信号进行量化、编码和调制后转发给目的节点：中继站(2)依据尽可能保留接收到的信号中关于信源部分的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息的原则选择出优化的量化区间并根据量化器函数q(y)＝z对其接收到的信号进行量化处理得到量化序列且对量化序列进行信源编码和信道编码得到编码序列c_r，再对编码序列c_r进行调制得到发送序列后发送给目的节点(3)，目的节点(3)接收到信号其中，y为信号的一个实现，z是元素集合Z＝{0，1,2,…,N-1}中的一个元素，N为量化电平数，量化器函数q(y)＝z将N个量化区间[y_min,y₁],[y₁,y₂],…,[y_N-1,y_max]与元素集合Z＝{0，1,2,…,N-1}中的元素0，1,2,…,N-1一一对应，y_min＝y₀<y₁<y₂<…<y_N-1<y_N＝y_max，y_min为量化区间边界的最小取值且第0个量化区间边界的取值y₀＝y_min，y_max为量化区间边界的最大取值且第N个量化区间边界的取值y_N＝y_max；n_r为对编码序列c_r进行调制后的序列长度，为零均值高斯白噪声，h_r,d为中继站(2)到目的节点(3)各信道的乘性系数；

步骤四、目的节点联合译码，具体过程为：

步骤401、目的节点(3)对接收到的信号进行信道译码和信源译码，恢复出量化序列

步骤402、目的节点(3)的多用户检测器联合 进行信号检测和译码，恢复出信源信息序列U_i，i∈{1,2,...,K}，具体过程为：

步骤4021、第一信号估计器ESE1和第二信号估计器ESE2分别由其先验信息和接收到的和逐比特联合估计输出外信息，其中，第一信号估计器ESE1的先验信息为译码交织送过来的信号估计器的先验信息加上第二信号估计器ESE2最近一次输出的外信息e_ESE2(x_i)，第二信号估计器ESE2的先验信息为译码交织送过来的信号估计器的先验信息加上第一信号估计器ESE1最近一次输出的外信息e_ESE1(x_i)；

步骤4022、对第一信号估计器ESE1输出的外信息和第二信号估计器ESE2输出的外信息求和得到各信源(1)的联合估计外信息e_ESE(x_i)；

步骤4023、对各信源(1)的联合估计外信息e_ESE(x_i)进行解交织得到译码器的先验信息送给用户译码器译码，用户译码器译码后输出译码外信息e_DEC(x_i)；

步骤4024、用户译码器对译码外信息e_DEC(x_i)进行交织得到信号估计器的先验信息并送给第一信号估计器ESE1和第二信号估计器ESE2；

步骤4025、重复步骤4021～4024迭代给定次数后停止，在用户译码器处进行硬判决输出作为译码结果，恢复出信源1信息序列U_i，i∈{1,2,...,K}。

2.按照权利要求1所述的交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于：步骤三中中继站(2)依据尽可能保留接收到的信号中关于信源部分的信息、同时尽可能去掉接收到的信号中关于噪声部分的信息的原则选择出优化的量化区间的具体过程为：

步骤301、参数输入和计算：给中继站(2)输入信号的概率密度函数p(x)、已知信道条件下信号的概率密度函数p(y|x)、量化电平数N和互信息变化阈值ε，中继站(2)根据公式计算出信号的概率密度函数p(y)；其中，x为信号的一个实现；

步骤302、初始化：随机选择一个满足表达式y_min＝y₀<y₁<y₂<…<y_N-1<y_N＝y_max的量化区间中继站(2)根据公式计算得到初始迭代的输出序列与信源(1)发送序列之间的互信息I^(old)；

步骤303、中继站(2)根据公式 $I^{\left(\mathrm{new}\right)}=\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 计算得到当前迭代的输出序列与信源(1)发送序列之间的互信息I^(new)，其中，j∈{1,2,…,N-1}，p(x,y)为信号和信号的联合概率密度函数且p(x,y)＝p(x)p(y|x)；

步骤304、当|I^(new)-I^(old)|≤ε时，执行步骤308；否则，当|I^(new)-I^(old)|>ε时，执行步骤305；

步骤305、将I^(old)更新为I^(new)；

步骤306、对每个j∈{1,2,…,N-1}，保持y₁,y₂,…,y_j-1,y_j+1…,y_N-1不变，中继站(2)根据公式 $y_j=\underset{y=[y_1,y_2,\&CenterDot;\&CenterDot;\&CenterDot;,y_{N-1}]}{\arg \max }\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 求解量化区间内的分量y_j，对步骤302中选择的量化区间 $\stackrel{\&RightArrow;}{y}=[y_1,y_2,\&CenterDot;\&CenterDot;\&CenterDot;,y_{N-1}]$ 进行更新；

步骤307、中继站(2)根据公式 $I^{\left(\mathrm{new}\right)}=\underset{x}{\mathrm{\&Sigma;}}\underset{j=0}{\overset{N-1}{\mathrm{\&Sigma;}}}[\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right|x)\mathrm{dy}-\log \underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p\left(y\right)\mathrm{dy}]\underset{y_j}{\overset{y_{j+1}}{\&Integral;}}p(x,y)\mathrm{dy}$ 计算得到当前迭代的输出序列与信源(1)发送序列之间的互信息I^(new)后，返回步骤304；

步骤308、输出量化区间

3.按照权利要求1所述的交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于：步骤402中所述迭代给定次数为20～40次。

4.按照权利要求1所述的交分多址接入中继系统的去噪互信息保持量化转发方法，其特征在于：所述零均值高斯白噪声零均值高斯白噪声和零均值高斯白噪声均为相互独立的零均值高斯白噪声。