1.一种中继功率分配比例和能量价格的协同优化方法，其特征在于：包括以下步骤：S1、获取第一初始参数、第二初始参数、第三初始参数；

S2、依据第一初始参数、各个电力公司的支付方案和电信运营商分配给中继的传输功率比例，以电力公司总成本最小为目标，以电力公司的支付方案为决策变量，建立第一优化模型；

S3、依据第二初始参数、各个电力公司的支付方案、电信运营商分配给中继的传输功率比例和电信运营商能量需求量，以电信运营商收益最大为目标，以电信运营商分配给中继的传输功率比例和电信运营商能量需求量为决策变量，建立第二优化模型；

S4、依据第三初始参数、能量供应点能量供应量和能量供应点单位能量价格，以能量供应点收益最大为目标，以能量供应点能量供应量和能量供应点单位能量价格为决策变量，建立第三优化模型；

S5、采用迭代算法对所述第一优化模型和所述第二优化模型求解，得到第一组最优解；

S6、将第一组最优解代入到第二优化模型中，采用迭代算法对所述第二优化模型和所述第三优化模型求解，得到第二组最优解；

S7、确定最优调度结果，所述最优调度结果包括第一组最优解和第二组最优解。

2.根据权利要求1所述的一种中继功率分配比例和能量价格的协同优化方法，其特征在于：S1中，所述第一初始参数包括DAU发送速率、DAU目标传输速率、每个DAU分配的带宽、DAU到网关以及中继到网关信道传输链路的信道噪声比、电力公司的预留成本和由风险管理策略决定的成本系数；

所述第二初始参数包括电信运营商的总带宽、电信运营商电信业务基价、基站中DAU数量、每个DAU的带宽；

所述第三初始参数包括能量供应点的成本系数。

3.根据权利要求1所述的一种中继功率分配比例和能量价格的协同优化方法，其特征在于：S2中，在电力公司和电信运营商互动过程中，对于电力公司i来说，决策变量为电力公司i的支付方案mi，建立的所述第一优化模型，包括以下具体步骤：首先利用Taguchi损失函数来描述电力公司i(i∈N＝{1,2,......,N})的平均的电力成本：

预测需求 是Taguchi损失函数的目标值， 是Taguchi损失函数的平均值，实际需求Di的分布为 为实际需求的方差，参数γ由电力市场的风险管理策略决定，为电力公司的预测误差， 是电力公司i的预留电力成本；

假定实际需求的方差与丢失的数据包总数成线性关系：数据包丢失总数表示为：

其中，θ表示由于丢失一个数据包而引起的方差增加， 和 表示DAUi的实际传输速率和目标传输速率，ts为采样周期；

采用DF中继策略时，DAUi的传输速率为：其中，W为DAUi分配的带宽， 为DAUi的传输功率， 为分配给DAUi中继的传输功率，和 分别为DAUi到网关以及中继到网关的信道传输链路的信道噪声比；

其中，mi为电力公司i给电信运营商的支付价格， 为电力公司i的支付方案，N是电力公司集(i∈N＝{1,2,......,i,......,N})，α为电信运营商分配给中继业务的功率比例，p为电信运营商的总功率；

综上，所述第一优化模型为电力公司i的成本模型为：所述第一优化模型约束条件为：

0<α<1

其中， 是分配给电力公司i的中继最大传输功率。

4.根据权利要求1所述的一种中继功率分配比例和能量价格的协同优化方法，其特征在于：S3中，电信运营商的总功率其中有αp分配给中继业务帮助电力公司更高效的收集用户的信息，(1‑α)p分配给电信业务来维持电信运营商自身电信业务的工作；

对于电信运营商来说，决策变量为电信运营商分配给中继的传输功率比例α和电信运营商能量需求量Q，所述第二优化模型为电信运营商的收益模型：max ZT(α,p)＝max(Zu+Zr‑Ue)其中，

Zr＝∑i∈Nmi

Ue＝min[q(E)·Q，q(E)·E]Q＝μp

所述第二优化模型约束条件为：

0<α<1

式中，ZT是电信运营商的收益，第一部分Zu是电信运营商电信业务的利润，第二部分Zr是所有电力公司给电信运营商中继业务的利润，第三部分Ue是电信运营商给能量供应点的能量报酬，ωr是电信业务的基价，Wr是电信运营商的总传输带宽，E为能量供应点产生的能量，q(E)为能量供应点单位能量价格，q(E)＝εE+q0为能量价格函数，ε是价格调节因子，q0是能量基价，Q是电信运营商能量需求量，μ是传输功率与能量之间的转换率。

5.根据权利要求1所述的一种中继功率分配比例和能量价格的协同优化方法，其特征在于：S4中，决策变量为能量供应点能量供应量E和能量供应点单位能量价格q，所述第三优化模型为能量供应点的收益模型为：max ZE(E,ε)＝max(Ze‑Ue)其中，

*

Ze＝min[q(E)·Q，q(E)·E]2

Ue＝aE+bE+c

q(E)＝εE+q0

*

当供大于求时：ZE(E,ε)＝q(E)·Q‑Ue当保守供应时：ZE(E,ε)＝q(E)·E‑Ue所述第三优化模型约束条件为：式中，ZE是能量供应点的收益，Ze是能量供应点收获的电信运营商购买能量的报酬，Ue是生产能量自身消耗的成本，a、b、c是能量成本系数。

6.根据权利要求1所述的一种中继功率分配比例和能量价格的协同优化方法，其特征在于：S5中，所述第一组最优解包括第一优化模型中各个电力公司的支付方案、第二优化模型中电信运营商分配给中继的传输功率比例。

7.根据权利要求1所述的一种中继功率分配比例和能量价格的协同优化方法，其特征在于：S6中，所述第二组最优解包括第二优化模型中电信运营商能量需求量、第三优化模型中能量供应点的能量供应量和能量供应点单位能量价格。

8.一种使用如权利要求1‑7任一项所述的中继功率分配比例和能量价格的协同优化方法的系统，其特征在于：包括参数获取模块、第一模型构建模块、第二模型构建模块、第三模型构建模块、第一求解模块、第二求解模块和调度结果确定模块；

所述参数获取模块，用于获取第一初始参数、第二初始参数和第三初始参数；

所述第一模型构建模块，用于依据第一初始参数、各个电力公司的支付方案和电信运营商分配给中继的传输功率比例，以电力公司总成本最小为目标，以电力公司的支付方案为决策变量，建立第一优化模型；

所述第二模型构建模块，用于依据第二初始参数、各个电力公司的支付方案和电信运营商分配给中继的传输功率比例、电信运营商能量需求量，以电信运营商收益最大为目标，所述电信运营商分配给中继的传输功率比例和电信运营商能量需求量为决策变量，建立第二优化模型；

所述第三模型构建模块，用于依据第三初始参数、能量供应点能量供应量以及单位能量价格，以能量供应点收益最大为目标，所述能量供应点能量供应量和能量供应点单位能量价格为决策变量，建立第三优化模型；

所述第一求解模块，采用迭代算法对所述第一优化模型和第二优化模型求解，得到第一组最优解；

所述第二求解模块，将第一组最优解代入到第二优化模型中，采用迭代算法对所述第二优化模型和第三优化模型求解，得到第二组最优解；

所述调度结果确定模块，用于确定最优调度结果，所述最优调度结果包括第一组最优解、第二组最优解。

9.根据权利要求8所述的一种中继功率分配比例和能量价格的协同优化系统，其特征在于：建立所述第一模型构建模块中的所述第一优化模型，包括以下具体步骤：首先利用Taguchi损失函数来描述电力公司i(i∈N＝{1,2,......,N})的平均的电力成本：

预测需求 是Taguchi损失函数的目标值， 是Taguchi损失函数的平均值，实际需求Di的分布为 表示实际需求的方差，参数γ由电力市场的风险管理策略决定，表示电力公司的预测误差， 是电力公司i的预留电力成本；

假定实际需求的方差与丢失的数据包总数成线性关系：数据包丢失总数表示为：

其中，θ表示由于丢失一个数据包而引起的方差增加， 和 表示DAUi的实际传输速率和目标传输速率，ts表示采样周期；

采用DF中继策略时，DAUi的传输速率为：其中，W表示DAUi分配的带宽， 表示DAUi的传输功率， 表示分配给DAUi中继的传输功率， 和 分别表示DAUi到网关以及中继到网关的信道传输链路的信道噪声比；

其中，mi表示电力公司i给电信运营商的支付价格， 为电力公司i的支付方案，N表示电力公司集(i∈N＝{1,2,......,i,......,N})，α表示电信运营商分配给中继业务的功率比例，p表示电信运营商的总功率；

综上，所述第一优化模型为电力公司i的成本模型为：所述第一优化模型约束条件为：

0<α<1

其中， 是分配给电力公司i的中继最大传输功率。

10.根据权利要求8所述的一种中继功率分配比例和能量价格的协同优化系统，其特征在于：所述第二模型构建模块中的所述第二优化模型为：max ZT(α,p)＝max(Zu+Zr‑Ue)其中，

Zr＝∑i∈Nmi

Ue＝min[q(E)·Q，q(E)·E]Q＝μp

所述第二优化模型约束条件为：

0<α<1

式中，ZT是电信运营商的收益，Zu是电信运营商电信业务的利润，Zr是所有电力公司给电信运营商中继业务的利润，Ue是电信运营商向能量供应点购买能量的成本，ωr是电信业务的基价，Wr是电信运营商的总传输带宽，E为能量供应点产生的能量，q(E)为能量供应点单位能量价格，q(E)＝εE+q0表示能量价格函数，ε是价格调节因子，q0是能量基价，Q是电信运营商能量需求量，μ是传输功率与能量之间的转换率。