1.一种物理层网络编码方法，其特征在于，包括：
步骤1，设置第一源节点、第二源节点和中继节点；
步骤2，将所述中继节点接收到的混叠星座图中，欧氏距离小于预设阈值的点映射到新星座图的同一个星座点，使满足：
[|(h1ai+h2bj)-(h1ap+h2bq)|2-εR]|sij-spq|2≥0
其中，所述第一源节点、所述第二源节点及所述中继节点所用的星座图分别为A＝{a1,...,aM},B＝{b1,...,bN},S＝{sij|i＝1,...,M,j＝1,...,N}，M和N分别表示所述第一源节点的星座图的阶数和所述第二源节点的星座图的阶数，h1和h2表示信道系数，εR为预设阈值，用于调整中继节点处网络编码映射的可靠性；sij＝C(ai,bj)表示当所述第一源节点发射信号ai，所述第二源节点发射信号bj时，所述中继节点发射信号sij，C表示网络编码映射方式：
步骤3，利用两个欧氏距离约束参数εA,εB来分别控制所述第一源节点、所述第二源节点的解码可靠性，使满足：
|sij-sik|2≥εA,sij≠sik∈S
|sji-ski|2≥εB,sji≠ski∈S
步骤4，以 为目标函数，建立优化关系：
[|(h1ai+h2bj)-(h1ap+h2bq)|2-εR]|sij-spq|2≥0
ai,ap∈A,bj,bq∈B,sij,spq∈S,i≠p,j≠q
其中，(.)H表示矩阵的共轭转置；
步骤5，对所述优化关系进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，步骤5具体包括：
步骤51，对所述优化关系进行线性变换和第一变量代换，以获得第一中间优化关系式，其中，n1、n2和n3分别表示相应约束的数量，n1、n2和n3均为阶数M和N的函数，表示相应系数矩阵；z1和z2分别为需要设计的上行星座图和下
行星座图；
第一变量代换具体为：
步骤52，利用第二变量代换将所述第一中间优化关系式转换为二次约束的二次优化问题，所述二次约束的二次优化问题表达式为：
其中， cj表示系数向量；
所述第二变量代换为：
其中，pj,qj为中间变量；
步骤53，对所述二次优化问题进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，步骤53具体包括：
步骤531，利用半定松弛算法对所述二次优化问题进行求解，以获得最优解X*；
步骤532，利用高斯向量法进行随机近似构造法，从所述最优解X*中提取所述二次优化问题的次优解，所述次优解即为近似最优解。
2.根据权利要求1所述的物理层网络编码方法，其特征在于，步骤531具体包括：
根据矩阵乘积的迹的性质，对所述二次优化问题进行如下变量替换：
xHΔx＝tr(xHΔx)＝tr(ΔxxH)
其中，Δ表示任意系数矩阵，xxT与秩为1的半正定矩阵X等效；
将所述二次优化问题经过第三变量代换，转换为第二优化关系式，所述第二优化关系式的表达式为：
s.t.tr(H1iZ1)≥1,i＝1,2,...,n1
tr(H2iZ2)≥1,i＝1,2,...,n2
tr(DjZ1)-cjz3＝1,tr(EjZ2)-djz3＝0
tr(FjZ3)≥0,j＝1,2,...,n3
rank(Z1)＝1,rank(Z2)＝1
其中， rank(.)为求矩阵的秩，tr(.)为矩阵的迹操作；
通过松弛掉rank(Z1)＝1,rank(Z2)＝1，并将 松弛为 将所述第二优化
关系式松弛为第一半定规划问题，其中所述第一半定规划问题为：
s.t.tr(H1iZ1)≥1,i＝1,2,...,n1
tr(H2iZ2)≥1,i＝1,2,...,n2
tr(DjZ1)-cjz3＝1,tr(EjZ2)-djz3＝0
tr(FjZ3)≥0,j＝1,2,...,n3
利用凸优化算法，求解所述第一半定规划问题，以获得最优解X*。
3.一种物理层网络编码方法，其特征在于，包括：
步骤1，设置第一源节点、第二源节点和中继节点；
步骤2，将所述中继节点接收到的混叠星座图中，欧氏距离小于预设阈值的点映射到新星座图的同一个星座点，使满足：
2 2
[|(h1ai+h2bj)-(h1ap+h2bq)|-εR]|sij-spq|≥0
其中，所述第一源节点、所述第二源节点及所述中继节点所用的星座图分别为A＝{a1,...,aM},B＝{b1,...,bN},S＝{sij|i＝1,...,M,j＝1,...,N}，M和N分别表示所述第一源节点的星座图的阶数和所述第二源节点的星座图的阶数，h1和h2表示信道系数，εR为预设阈值，用于调整中继节点处网络编码映射的可靠性；sij＝C(ai,bj)表示当所述第一源节点发射信号ai，所述第二源节点发射信号bj时，所述中继节点发射信号sij，C表示网络编码映射方式：
步骤3，利用两个欧氏距离约束参数εA,εB来分别控制所述第一源节点、所述第二源节点的解码可靠性，使满足：
|sij-sik|2≥εA,sij≠sik∈S
|sji-ski|2≥εB,sji≠ski∈S
步骤4，以 为目标函数，建立优化关系：
[|(h1ai+h2bj)-(h1ap+h2bq)|2-εR]|sij-spq|2≥0
ai,ap∈A,bj,bq∈B,sij,spq∈S,i≠p,j≠q
其中，(.)H表示矩阵的共轭转置；
步骤5，对所述优化关系进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，步骤5具体包括：
步骤51，对所述优化关系进行线性变换和第一变量代换，以获得第一中间优化关系式，其中，n1、n2和n3分别表示相应约束的数量，n1、n2和n3均为阶数M和N的函数，表示相应系数矩阵；z1和z2分别为需要设计的上行星座图和下
行星座图；
第一变量代换具体为：
步骤52，利用第二变量代换将所述第一中间优化关系式转换为二次约束的二次优化问题，所述二次约束的二次优化问题表达式为：
其中， cj表示系数向量；
所述第二变量代换为：
其中，pj,qj为中间变量；
步骤53，对所述二次优化问题进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，步骤53具体包括：
步骤533，利用第四变量代换、第五变量代换和第六变量代换，将所述二次优化问题转换为第三优化关系式，其中所述第三优化关系式为：
min tr(Z1)+tr(Z2)+tr(Z3)
s.t.tr(FiZ1)≥1,i＝1,2,...,m1
tr(GjZ2)≥1,j＝1,2,...,m2
tr(HiZ3)≥1,i＝1,2,...,m3
rank(Z1)＝1,rank(Z2)＝1,rank(Z3)＝1
其中，Fi,Gj,Hi,Dk,fi,gj,hi,dk为相应地系数矩阵和系数向量；
所述第四变量代换为：
所述第五变量代换为：
所述第六变量代换为：
步骤534，松弛掉所述第三优化关系式中的秩约束rank(Z1)＝1,rank(Z2)＝1,rank(Z3)＝1和 以获得松弛问题，所述松弛问题表达式为
min tr(Z1)+tr(Z2)+tr(Z3)
s.t.tr(FiZ1)≥1,i＝1,2,...,m1
tr(GjZ2)≥1,j＝1,2,...,m2
tr(HiZ3)≥1,i＝1,2,...,m3
步骤535，将所述松弛问题的求解转化为第二半定规划问题的求解，所述第二半定规划问题为：
min tr(Z1)+tr(Z2)+tr(Z3)
s.t.tr(FiZ1)≥1,i＝1,2,...,m1
tr(GjZ2)≥1,j＝1,2,...,m2
tr(HiZ3)≥1,i＝1,2,...,m3
步骤536，利用高斯近似法从所述第二半定规划问题的最优解中提取所述二次优化问题的次优解，所述次优解即为近似最优解。
4.一种物理层网络编码装置，其特征在于，包括：
设置模块，用于设置第一源节点、第二源节点和中继节点；
映射模块，用于将所述中继节点接收到的混叠星座图中，欧氏距离小于预设阈值的点映射到新星座图的同一个星座点，使满足：
2 2
[|(h1ai+h2bj)-(h1ap+h2bq)|-εR]|sij-spq|≥0
其中，所述第一源节点、所述第二源节点及所述中继节点所用的星座图分别为A＝{a1,...,aM},B＝{b1,...,bN},S＝{sij|i＝1,...,M,j＝1,...,N}，M和N分别表示所述第一源节点的星座图的阶数和所述第二源节点的星座图的阶数，h1和h2表示信道系数，εR为预设阈值，用于调整中继节点处网络编码映射的可靠性；sij＝C(ai,bj)表示当所述第一源节点发射信号ai，所述第二源节点发射信号bj时，所述中继节点发射信号sij，C表示网络编码映射方式：
控制模块，用于利用两个欧氏距离约束参数εA,εB来分别控制所述第一源节点、所述第二源节点的解码可靠性，使满足：
|sij-sik|2≥εA,sij≠sik∈S
|sji-ski|2≥εB,sji≠ski∈S
优化关系建立模块，用于以 为目标函数，建立如下
优化关系：
[|(h1ai+h2bj)-(h1ap+h2bq)|2-εR]|sij-spq|2≥0
ai,ap∈A,bj,bq∈B,sij,spq∈S,i≠p,j≠q
H
其中，(.) 表示矩阵的共轭转置；
求解模块，用于对所述优化关系进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，求解模块具体包括：
第一中间优化关系式获取子模块，用于对所述优化关系进行线性变换和第一变量代换，以获得第一中间优化关系式，
其中，n1、n2和n3分别表示相应约束的数量，n1、n2和n3均为阶数M和N的函数，表示相应系数矩阵；z1和z2分别为需要设计的上行星座图和下
行星座图；
第一变量代换具体为：
二次优化问题表达式获取子模块，利用第二变量代换将所述第一中间优化关系式转换为二次约束的二次优化问题，所述二次约束的二次优化问题表达式为：
其中， cj表示系数向量；
所述第二变量代换为：
其中，pj,qj为中间变量；
求解子模块，用于对所述二次优化问题进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，求解子模块具体包括：
最优解获取子模块，用于利用半定松弛算法对所述二次优化问题进行求解，以获得最*
优解X；
第一次优解获取子模块，用于利用高斯向量法进行随机近似构造法，从所述最优解X*中提取所述二次优化问题的次优解，所述次优解即为近似最优解。
5.根据权利要求4所述的物理层网络编码装置，其特征在于，最优解获取子模块具体用于：
根据矩阵乘积的迹的性质，对所述二次优化问题进行如下变量替换：
xHΔx＝tr(xHΔx)＝tr(ΔxxH)
其中，Δ表示任意系数矩阵，xxT与秩为1的半正定矩阵X等效；
将所述二次优化问题经过第三变量代换，转换为第二优化关系式，所述第二优化关系式的表达式为：
s.t.tr(H1iZ1)≥1,i＝1,2,...,n1
tr(H2iZ2)≥1,i＝1,2,...,n2
tr(DjZ1)-cjz3＝1,tr(EjZ2)-djz3＝0
tr(FjZ3)≥0,j＝1,2,...,n3
rank(Z1)＝1,rank(Z2)＝1
其中， rank(.)为求矩阵的秩，tr(.)为矩阵的迹操作；
通过松弛掉rank(Z1)＝1,rank(Z2)＝1，并将 松弛为 将所述第二优化
关系式松弛为第一半定规划问题，其中所述第一半定规划问题为：
s.t.tr(H1iZ1)≥1,i＝1,2,...,n1
tr(H2iZ2)≥1,i＝1,2,...,n2
tr(DjZ1)-cjz3＝1,tr(EjZ2)-djz3＝0
tr(FjZ3)≥0,j＝1,2,...,n3
利用凸优化算法，求解所述第一半定规划问题，以获得最优解X*。
6.一种物理层网络编码装置，其特征在于，包括：
设置模块，用于设置第一源节点、第二源节点和中继节点；
映射模块，用于将所述中继节点接收到的混叠星座图中，欧氏距离小于预设阈值的点映射到新星座图的同一个星座点，使满足：
[|(h1ai+h2bj)-(h1ap+h2bq)|2-εR]|sij-spq|2≥0
其中，所述第一源节点、所述第二源节点及所述中继节点所用的星座图分别为A＝{a1,...,aM},B＝{b1,...,bN},S＝{sij|i＝1,...,M,j＝1,...,N}，M和N分别表示所述第一源节点的星座图的阶数和所述第二源节点的星座图的阶数，h1和h2表示信道系数，εR为预设阈值，用于调整中继节点处网络编码映射的可靠性；sij＝C(ai,bj)表示当所述第一源节点发射信号ai，所述第二源节点发射信号bj时，所述中继节点发射信号sij，C表示网络编码映射方式：
控制模块，用于利用两个欧氏距离约束参数εA,εB来分别控制所述第一源节点、所述第二源节点的解码可靠性，使满足：
2
|sij-sik|≥εA,sij≠sik∈S
|sji-ski|2≥εB,sji≠ski∈S
优化关系建立模块，用于以 为目标函数，建立如下
优化关系：
[|(h1ai+h2bj)-(h1ap+h2bq)|2-εR]|sij-spq|2≥0
ai,ap∈A,bj,bq∈B,sij,spq∈S,i≠p,j≠q
其中，(.)H表示矩阵的共轭转置；
求解模块，用于对所述优化关系进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，求解模块具体包括：
第一中间优化关系式获取子模块，用于对所述优化关系进行线性变换和第一变量代换，以获得第一中间优化关系式，
其中，n1、n2和n3分别表示相应约束的数量，n1、n2和n3均为阶数M和N的函数，表示相应系数矩阵；z1和z2分别为需要设计的上行星座图和下
行星座图；
第一变量代换具体为：
二次优化问题表达式获取子模块，利用第二变量代换将所述第一中间优化关系式转换为二次约束的二次优化问题，所述二次约束的二次优化问题表达式为：
其中， cj表示系数向量；
所述第二变量代换为：
其中，pj,qj为中间变量；
求解子模块，用于对所述二次优化问题进行求解，以获得物理层网络编码机制下双向中继网络的新星座图和网络编码映射方式的近似最优解；
其中，求解子模块具体包括：
第三优化关系式获取子模块，用于利用第四变量代换、第五变量代换和第六变量代换，将所述二次优化问题转换为第三优化关系式，其中所述第三优化关系式为：
min tr(Z1)+tr(Z2)+tr(Z3)
s.t.tr(FiZ1)≥1,i＝1,2,...,m1
tr(GjZ2)≥1,j＝1,2,...,m2
tr(HiZ3)≥1,i＝1,2,...,m3
rank(Z1)＝1,rank(Z2)＝1,rank(Z3)＝1
其中，Fi,Gj,Hi,Dk,fi,gj,hi,dk为相应地系数矩阵和系数向量；
所述第四变量代换为：
所述第五变量代换为：
所述第六变量代换为：
松弛问题表达式获取子模块，用于松弛掉所述第三优化关系式中的秩约束rank(Z1)＝
1,rank(Z2)＝1,rank(Z3)＝1和 以获得松弛问题，所述松弛问题表达式为：
min tr(Z1)+tr(Z2)+tr(Z3)
s.t.tr(FiZ1)≥1,i＝1,2,...,m1
tr(GjZ2)≥1,j＝1,2,...,m2
tr(HiZ3)≥1,i＝1,2,...,m3
第二半定规划问题获取子模块，用于将所述松弛问题的求解转化为第二半定规划问题的求解，所述第二半定规划问题为：
min tr(Z1)+tr(Z2)+tr(Z3)
s.t.tr(FiZ1)≥1,i＝1,2,...,m1
tr(GjZ2)≥1,j＝1,2,...,m2
tr(HiZ3)≥1,i＝1,2,...,m3
第二次优解获取子模块，用于利用高斯近似法从所述第二半定规划问题的最优解中提取所述二次优化问题的次优解，所述次优解即为近似最优解。