1.一种用于D2D异构无线通信网络中基于时间反演的抗干扰方法，包括以下步骤：步骤101：利用时间反演镜对蜂窝用户和D2D用户执行信道签名，并分别获得蜂窝用户信干噪比和D2D用户信干噪比；

所述步骤101利用时间反演镜对蜂窝用户和D2D用户执行信道签名包括：基站接收端利用导频信号即直接序列扩频信号进行信道探测，获取信道冲击响应，然后得到信道冲击响应的时间反演共轭形式，最后利用信道冲击响应的时间反演的共轭形式作为时间反演处理器的传输函数；设定蜂窝用户发送信号到基站接收端的时间反演镜为 即定义信道签名： 其中*表示共轭， 表示第i个蜂窝用户到基站归一化时间反演波形即信道签名； 表示第i个蜂窝用户到基站具有长度l的信道延迟分布响应； 表示第i个蜂窝用户发送的符号序列{X(k)}到基站的信道延迟分布响应；L表示信道长度；同时，第j个D2D发送端到基站所获得的时间反演信道签名为 表示第j个D2D发送端到基站归一化时间反演波形即信道签名； 表示第j个D2D用户到基站具有长度l的信道延迟分布响应； 表示第j个D2D蜂窝用户发送得符号序列{X(k)}到基站的信道延迟分布响应，k表示蜂窝用户序号；

步骤102：采用启发式算法为蜂窝用户设置蜂窝用户信干噪比门限和为D2D用户设置D2D用户信干噪比门限，包括：为蜂窝用户设置蜂窝用户信干噪比门限：初始化t＝0，容量增益因子

步骤102A：基站遍历所有用户N，依次增加信干噪比值Δ，并分别计算功率增量以及容量增量

其 表示蜂窝用户此刻信干噪比；并计算容量增益因子 其表示消耗单位功率所获得的容量；Di表示蜂窝用户i的对角矩阵；

2

表示小区其他蜂窝用户对发送用户i的干扰；Pdt表示D2D用户发送功率矩阵；σ表示方差；αbsi为第i个蜂窝用户到基站的路径增益标量系数；I表示单位矩阵； 表示D2D通信对蜂窝用户i造成的共道干扰；

(t) (t)

步骤102B：找出容量增益最大的用户UE ＝argmax{bi }，如果用户容量增益因子差距(t)

甚小 则选择信道增益最大的用户UE ＝argmax{D1/

2 2

αbs1 ,...,Di/αbsi}i＝1,...,M，其中M表示蜂窝用户数；否则选择容量增益最大的用户；

步骤102C：更新信干噪比门限pk (t)

步骤102D：依据香农原理，若系统总容量满足约束目标C ≤C ，则退出循环，返回信干pk (t)

噪比门限，其中C 表示预先设定的容量目标，C 表示t时刻系统总容量；否则重复步骤

102A‑102C；

同理，对于为D2D用户设置D2D用户信干噪比门限与上述一致；

步骤103：根据功率分配的规划目标及所述蜂窝用户信干噪比门限和D2D用户信干噪比门限，采用卡罗需‑库恩‑塔克KKT条件分配用户功率。

2.根据权利要求1所述的基于时间反演的抗干扰方法，其特征在于，步骤101中蜂窝用户信干噪比、D2D用户信干噪比分别表示为：其中， 表示蜂窝用户信干噪比 ， 表示D2D用户信干噪比，为第i个蜂窝用户到基站的路径增益标量系数，其中dbs_i表示蜂窝用户i到基站的距离，αc表示用户i到基站的路径损耗因子，χbs_i表示蜂窝用户i到基站产生的阴影衰落，其服从对数正态分布；相应地，αbsk为第k个蜂窝用户到基站的路径增益标量系数，αdj为第j个D2D用户到基站的路径增益标量系数，αdk为第k对D2D用户到第j对D2D接收端的路径增益标量系数，αdi为第i个蜂窝用户到第j对D2D用户接收端的路径增益标量系数； 表示第i个蜂窝用户的发送功率，如果蜂窝用户k≠i， 表示第k个蜂窝用户发送功率， 表示第j对D2D用户的发送功率，如果D2D用户对k≠j， 表示第k对D2D用户发送功率； 表示第i个蜂窝用户到基站时间反演信道签名矩阵； 表示第k个蜂窝用户到基站时间反演信道签名矩阵； 表示第j对D2D用户时间反演信道签名矩阵； 表示第k对D2D用户时间反演信道签名矩阵；其中 上标H表示上述对应信道签名矩阵共轭转置；

表示第i个蜂窝用户到基站所形成信道拓普利兹矩阵， 表示第i个蜂窝用户到基站所形成信道拓普利兹矩阵的第L行向量；且基站接收到的加性高斯白噪声均值为0，方差2

为σ ，为了简化表达式，第i个蜂窝用户到基站信道协方差矩阵表示为第i个蜂窝用户到基站信道自相关矩阵为表示第i个蜂窝用户到基站产生

码间干扰信道协方差矩阵， 表示D2D干扰用户到基站信道协方差矩阵， 表示除发送用户之外其他干扰用户到基站信道协方差矩阵；其中上标H表示共轭转置， 表示第j个D2D用户到基站所形成信道托普利兹矩阵；同理分别表示第j对D2D用户之间信道自相关矩阵，第j对D2D用户之间码间干扰信道协方差矩阵，除发送D2D用户之外其他干扰用户到第j对接收用户信道协方差矩阵，蜂窝用户到第j对D2D接收用户干扰信道协方差矩阵。

3.根据权利要求2所述基于时间反演的抗干扰方法，其特征在于，所述步骤101蜂窝用户信干噪比和D2D用户信干噪比的规范定义简化式分别如下：其中，上标T表示转置，Φ为串扰矩阵，如果蜂窝用户k≠i，表示小区其他蜂窝用户对发送用户i的干扰，如果k＝i，则 Pct表示蜂窝用户发送功率矩阵，Pdt表示D2D用户发送功率矩阵；A表示D2D用户对蜂窝用户产生的干扰，其取值为 此外，定义蜂窝用户i和D2D用户j的对角矩阵分别为 和 表示蜂

窝用户对D2D接收用户的同层干扰矩阵；相应地，如果D2D用户对k≠j，则串扰矩阵表示小区内其他D2D用户对发送用户的干扰；如果D2D用户对k＝j， 表示D2D发送用户之间的码间干扰； 表示蜂窝用户到第j对D2D接收用户干扰信道协方差矩阵。

4.根据权利要求1所述基于时间反演的抗干扰方法，其特征在于，所述步骤103采用卡罗需‑库恩‑塔克KKT条件分配用户功率，包括以最小化系统总功率为优化目标，最小系统容量和信干噪比门限为约束条件，构建功率分配目标函数：min∑iPi

s.t.C1

C2

sum

其中，i表示蜂窝用户，Pi表示蜂窝用户i发送功率，C 表示当前设定的系统容量，表示信干噪比门限，γi(P)表示当前功率下的信干噪比， 表示在信干噪比门限下所获得的系统容量；约束条件C1，表示在信干噪比门限下的所有用户的容量之和不小于当前设定的系统容量，约束条件C2，表示在当前功率条件下所获得的信干噪比值不小于信干噪比门限。

5.根据权利要求4所述基于时间反演的抗干扰方法，其特征在于，所述采用KKT条件对用户分配功率进一步包括：根据蜂窝用户功率分配的优化目标函数和约束条件，引用KKT条件得到：

其中， α和β是拉格朗日乘子，为常数；

为注水线。

6.根据权利要求4所述基于时间反演的抗干扰方法，其特征在于，所述采用KKT条件对用户分配功率进一步包括：根据D2D用户功率分配的优化目标函数和约束条件，引用KKT条件得到：

其中， θ和ε是拉格朗日乘子，为常数；

为注水线。