轻量级分组密码算法ESF的相关密钥不可能差分分析

谢敏; 曾琦雅

doi:10.11999/JEIT180576

轻量级分组密码算法ESF的相关密钥不可能差分分析

doi: 10.11999/JEIT180576

谢敏^,,
曾琦雅

西安电子科技大学综合业务网理论及关键技术国家重点实验室西安 710071

基金项目: 国家重点研发计划(2016YFB0800601)，国家自然科学基金委员会-通用联合基金重点项目(U1636209)，“十三五”国家密码发展基金(MMJJ20180219)

详细信息

作者简介:
谢敏：女，1976年生，副教授，研究方向为编码和密码

曾琦雅：女，1993年生，硕士，研究方向为分组密码算法分析

通讯作者:
谢敏　mxie@xidian.edu.cn

中图分类号: TN918.1
计量
- 文章访问数: 2943
- HTML全文浏览量: 1021
- PDF下载量: 90
- 被引次数: 2
出版历程
- 收稿日期: 2018-06-11
- 修回日期: 2018-12-19
- 网络出版日期: 2018-12-26
- 刊出日期: 2019-05-01

Related-key Impossible Differential Cryptanalysis on Lightweight Block Cipher ESF

Min XIE^,,
Qiya ZENG

State Key Laboratory of Integrated Services Networks, Xidian University, Xi’an 710077, China

Funds: The National Key Research and Development Program of China (2016YFB0800601), The Key Project of the General Joint Fund of the National Natural Science of China (U1636209), National Cryptographic Development Fund of the 13th Five-Year Plan (MMJJ20180219)

摘要

摘要:
八阵图算法(ESF)是一种具有广义Feistel结构的轻量级分组密码算法，可用在物联网环境下保护射频识别(RFID)标签等资源受限的环境中，目前对该算法的安全性研究主要为不可能差分分析。该文通过深入研究S盒的特点并结合ESF密钥扩展算法的性质，研究了ESF抵抗相关密钥不可能差分攻击的能力。通过构造11轮相关密钥不可能差分区分器，在此基础上前后各扩展2轮，成功攻击15轮ESF算法。该攻击的时间复杂度为2^40.5次15轮加密，数据复杂度为2^61.5个选择明文，恢复密钥比特数为40 bit。与现有结果相比，攻击轮数提高的情况下，时间复杂度降低，数据复杂度也较为理想。
- 轻量级分组密码 /
- ESF算法 /
- 相关密钥 /
- 不可能差分分析
Abstract:
Eight-Sided Fortress (ESF) is a lightweight block cipher with a generalized Feistel structure, which can be used in resource-constrained environments such as protecting Radio Frequency IDentification (RFID) tags in the internet of things. At present, the research on the security of ESF mainly adopts the impossible differential cryptanalysis. The ability of ESF to resist the related-key impossible differential cryptanalysis is studied based on the characteristics of its S-boxes and key schedule. By constructing an 11-round related-key impossible differential distinguisher, an attack on 15-round ESF is proposed by adding 2-round at the top and 2-round at the bottom. This attack has a time complexity of 2^40.5 15-round encryptions and a data complexity of 2^61.5 chosen plaintexts with 40 recovered key-bit. Compared with published results, the time complexity is decreased and the data complexity is ideal with the number of attack rounds increased.
- Lightweight block cipher /
- ESF algorithm /
- Related-key /
- Impossible differential attack

HTML全文

1. 引言

近年来，自动驾驶、远程医疗和增强现实等一系列新兴业务和产业蓬勃发展，这对无线网络的信息处理能力提出了更高要求，促使无线网络逐步发展成为通信网络、感知网络和计算网络的融合体^[1–3]。与单纯的无线通信技术相比，通信感知一体化技术不仅具有原生的通信功能，而且还具有对环境的感知功能。此外，系统中的资源被通信和感知功能共享，如何协调分配系统资源成为性能优化的新挑战^[4]。

无人机(Unmanned Aerial Vehicle, UAV)具有按需部署、机动可控等优点^[5]。无人机辅助通信感知一体化系统可以充分利用无人机的机动性，动态调整无人机的位置，提高通信用户或者感知目标的信道质量，更好地实现通信和感知性能的优化^[6–10]。文献[6]联合优化无人机位置和发射功率来最大化通信感知一体化系统的网络效用。文献[7]联合设计无人机的飞行轨迹和发射波束赋形，最大化通信用户的通信速率。文献[8]联合优化无人机轨迹、发射波束赋形和感知任务执行时刻，最大化通信用户的通信速率。文献[9]提出一种联合无人机基站和地面基站预测和跟踪用户位置的通信感知一体化方法，并联合优化无人机轨迹以及无人机和用户的关联来最大化用户通信速率和位置克拉美罗下界的加权和。文献[10]考虑概率视距信道模型，在通信和感知服务质量的约束下，联合优化波束赋形和无人机轨迹来提高系统吞吐量。

然而，在上述工作中，空地信道模型大多采用视距信道模型，并忽略了无人机与回程基站之间的回程链路，具有一定的局限性。具体而言，无人机与地面用户或者感知目标的视距信道模型在城市或山区环境下是不准确的，因为它忽略了随机阴影和小尺度衰落^[11]。为了更加准确地建模实际信道，空地信道模型可以采用莱斯信道模型。莱斯信道模型由一个确定性的视距分量和一个由障碍物的反射、散射和衍射引起的随机多路径分量组成，考虑了小尺度衰落等实际因素，对实际信道建模更加准确。同时，当无人机中继回程链路的信道质量较差或者可用带宽受限时，可能会导致用户通信速率下降^[12]。

解决回程链路问题的一种有效方法是在系统中考虑空地协同网络。将无人车(Unmanned Ground Vehicle, UGV)基站作为无人机中继的回程基站和能源补给站，二者协同构建空地通信感知协同网络，提升通信感知一体化系统的性能和服务范围。文献[13]提出一种无人车基站和无人机基站集群组成空地合作应急网络帮助灾区快速重建通信。文献[14]提出了一种空地协同车载组网体系结构。然而，在上述空地协同工作中，空地协同的方式集中在通信或者感知单一方面，因此其无法直接应用在通信感知一体化系统上。

基于上述工作的局限性，本文研究空地协同通信感知一体化系统，由无人车基站和无人机中继集群组成空地协同网络。采用更加精确且与高度角相关的莱斯信道模型建模空地信道。在莱斯信道下，无人机中继与回程基站、用户或者感知目标的高度角越大，传输过程中遇到的障碍物越少^[11]，所以空地协同网络的通信感知性能不仅与传输距离有关，还受高度角的影响。因此，本文提出面向空地协同通信感知一体化系统的资源分配和轨迹联合优化方法，在保证系统感知性能的同时，提高通信性能。

2. 系统模型

如图1所示，本文考虑一个由无人车基站和无人机中继集群组成的空地协同通信感知一体化系统，其中，1个单天线无人车基站作为K个无人机中继的回程基站，无人机中继集群接收无人车基站发送的数据，然后将其发送给G个单天线用户，并感知R个目标区域。假设无人车、无人机和用户天线各个方向的增益相同^[5,11,15]。在任务开始时，无人车基站搭载并发射无人机集群执行通信感知任务，在任务结束后，无人车回收无人机集群，获取无人机集群对目标区域的回声感知信息。

图 1 空地协同通信感知一体化系统

下载: 全尺寸图片幻灯片

为了减少链路之间的干扰，系统采用频分复用的方式，不同链路的信号带宽相互正交且相等^[6]，即 ${B_{{\mathrm{c}},k}} = {B_k} = {B_{{\text{sum}}}}/(2K),\forall k$ ，其中 ${B_{{\text{sum}}}}$ 表示系统总带宽， ${B_{{\mathrm{c}},k}}$ 表示无人车与无人机k链路所使用的信号带宽， ${B_k}$ 表示无人机k与关联用户以及目标区域链路所使用的信号带宽。

将任务周期T离散化为N个时隙，每个时隙的长度为 ${\delta _t} = T/N$ 。设无人车初始时刻和时隙n的水平位置分别为 ${{\boldsymbol{q}}_{\text{I}}}$ 和 ${{\boldsymbol{q}}_{\text{c}}}[n]$ ，无人机k在时隙n的水平位置和垂直高度分别为 ${{\boldsymbol{q}}_k}[n]$ 和 ${h_k}[n]$ ，则无人车和无人机集群的水平位置约束和速度约束分别为

${{\boldsymbol{q}}_k}[0] = {{\boldsymbol{q}}_{\text{c}}}[0] = {{\boldsymbol{q}}_{\text{I}}},{{\boldsymbol{q}}_k}[N] = {{\boldsymbol{q}}_{\text{c}}}[N],\forall k$

(1)

$\left\| {{{\boldsymbol{q}}_i}[n] - {{\boldsymbol{q}}_i}[n - 1]} \right\| \le {\delta _t}V_{i,\max }^{xy},\forall n,i \in \{ {\text{c}},k\}$

(2)

$\left| {{h_k}[n] - {h_k}[n - 1]} \right| \le {\delta _t}V_{k,\max }^z,\forall n,k$

(3)

其中， $\left\| \cdot \right\|$ 表示欧几里得函数， $V_{i,\max }^{xy},i \in \{ {\text{c}},k\}$ 表示无人车c或者无人机k最大水平行驶速度。 $V_{k,\max }^z$ 表示无人机k垂直方向最大行驶速度。

假设用户g、目标区域r和障碍物j的位置已知^[9]，分别为 ${{\boldsymbol{u}}_g}$ , ${{\boldsymbol{u}}_r}$ 和 ${{\boldsymbol{u}}_j}$ ，则无人车防碰撞约束表示为

${d_i} + {d_{{\text{c}},\min }} \le \left\| {{{\boldsymbol{q}}_{\text{c}}}[n] - {{\boldsymbol{u}}_i}} \right\|,\forall n,i \in \{ g,r,j\}$

(4)

其中， ${d_{{\text{c,min}}}}$ 表示无人车的防碰撞安全距离， ${d_i}[n], \forall i \in \{ g,r,j\}$ 表示用户g，目标区域r或者障碍物j边界与其中心位置的最大距离。

设 ${h_{k,{\text{I}}}}$ 为无人机k在时隙0和时隙N的垂直高度，最大最小飞行高度分别为 ${H_{\max }}$ 和 ${H_{\min }}$ ，则

${h_k}[n] = {h_{k,{\text{I}}}},n \in \{ 0,N\} ,\forall k$

(5)

${H_{\min }} \le {h_k}[n] \le {H_{\max }},\forall n,k$

(6)

设 ${d_{\min }}$ 表示无人机集群的最小安全距离，则无人机集群的防碰撞约束为

$\begin{split} {({d_{\min }})^2}\,& \le {\left\| {{{\boldsymbol{q}}_k}[n] - {{\boldsymbol{q}}_{k'}}[n]} \right\|^2} + {\left| {{h_k}[n] - {h_{k'}}[n]} \right|^2},\\ & \,\forall n,k \ne k' \end{split}$

(7)

定义物体集合 ${{I}} \triangleq \{ {\text{c}},g,r\}$ 。本文采用莱斯信道模型。在时隙n，无人机k到物体i的信道功率增益

${h_{k,i}}[n] = \sqrt {{\beta _{k,i}}[n]{f_{k,i}}[n]} ,\forall n,k,i \in {{I}}$

(8)

其中， ${\beta _{k,i}}[n] = {\beta _0}{({d_{k,i}}[n])^{ - 2}}$ ， ${\beta _0}$ 表示参考距离为1 m时的信道功率增益， ${d_{k,i}}[n]$ 表示无人机k到物体i的距离。 ${f_{k,i}}[n]$ 表示无人机k到物体i的有效衰落功率。

因为 ${f_{k,i}}[n]$ 可以用logistic函数(“S”形状)近似^[5,11,15,16]，即

$\begin{split} {f_{k,i}}[n]& = \left({C_1} + \frac{{{C_2}}}{{1 + \exp ( - ({B_1} + {B_2}{v_{k,i}}[n]))}}\right),\\ &\, \forall n,k,i \in {{I}} \end{split}$

(9)

其中， ${B_1}$ , ${B_2}$ , ${C_1}$ , ${C_2}$ 是最大最小莱斯因子决定的常数。 ${v_{k,i}}[n]$ 表示无人机k与物体i高度角的正弦值¹。

设无人车c和无人机k的峰值发射功率和平均发射功率分别为 ${P_{{\text{c}},\max }}$ , ${\bar P_{\text{c}}}$ , ${P_{k,\max }}$ 和 ${\bar P_k}$ 。在时隙n，设无人车c向无人机k发射的信号功率为 ${p_{{\text{c}},k}}[n]$ ，无人机k的发射功率为 ${p_k}[n]$ ，则

$\sum\limits_{k = 1}^K {{p_{{\text{c}},k}}[n]} \le {P_{{\text{c}},\max }},{p_k}[n] \le {P_{k,\max }},\forall n$

(10)

$\frac{1}{N}\sum\limits_{n = 1}^N {\sum\limits_{k = 1}^K {{p_{{\text{c}},k}}[n]} } \le {\bar P_{\text{c}}},\frac{1}{N}\sum\limits_{n = 1}^N {{p_k}[n]} \le {\bar P_k},\forall k$

(11)

设加性高斯白噪声功率谱密度为 ${N_0}$ 。在时隙n，无人机k回程链路可实现的通信速率(bit/(s·Hz))为

${R_{k,{\text{c}}}}[n] = \frac{1}{{2K}}{\log _2}\left(1 + \frac{{{{\left| {{h_{k,{\text{c}}}}[n]} \right|}^2}{p_{{\text{c}},k}}[n]}}{{{B_{{\text{c}},k}}{N_0}}}\right)$

(12)

为了减少不同用户之间的干扰，在每个时隙内，无人机与地面用户为1对1的关系^[6]，即

$\qquad 0 \le \sum\limits_{g = 1}^G {{a_{k,g}}[n]} \le 1,\forall n,k$

(13)

$\qquad 0 \le \sum\limits_{k = 1}^K {{a_{k,g}}[n]} \le 1,\forall n,g$

(14)

$\qquad {a_{k,g}}[n]({a_{k,g}}[n] - 1) = 0,\forall n,k,g$

(15)

其中， ${a_{k,g}}[n]$ 表示无人机k在时隙n与用户g的关联变量。

在不考虑回程链路的情况下，地面用户g在时隙n的通信速率(bit/(s·Hz))为

$\begin{split} {R_g}[n]& = \frac{1}{{2K}}\sum\limits_{k = 1}^K {{a_{k,g}}[n]{{\log }_2}\left(1 + \frac{{{{\left| {{h_{k,g}}[n]} \right|}^2}{p_k}[n]}}{{{B_k}{N_0}}}\right)} ,\\ & \,\forall n,g\\[-1pt] \end{split}$

(16)

在任意时隙内，无人机k与关联的地面用户的通信速率之和不大于无人机k回程链路通信速率，即

${R_g} \le \sum\limits_{k = 1}^K {{a_{k,g}}[n]{R_{k,{\text{c}}}}[n]} ,\forall n,g$

(17)

为了减少无人机的能耗，在每个时隙内，无人机k最多感知1个目标区域^[17]，即

$\qquad 0 \le \sum\limits_{r = 1}^R {{b_{k,r}}[n]} \le 1,\forall n,k$

(18)

$\qquad {b_{k,r}}[n]({b_{k,r}}[n] - 1) = 0,\forall n,k,r$

(19)

其中， ${b_{k,r}}[n]$ 表示无人机k在时隙n与目标区域r的关联变量。

为了提高感知信息的准确性，在整个任务周期内，目标区域至少需要被感知 ${f_r}$ 次，即

${f_r} \le \sum\limits_{n = 1}^N {\sum\limits_{k = 1}^K {{b_{k,r}}[n]} } ,\forall r$

(20)

在时隙n，无人机k与目标区域r关联，则无人机k对目标区域r的有效感知功率

$\begin{split} {p_{k,r}}[n] = \,& \frac{{{\beta _0}{p_k}[n]}}{{d_{k,r}^2[n]}}\\ & \cdot\left({C_1} + \frac{{{C_2}}}{{1 + \exp ( - ({B_1} + {B_2}{v_{k,r}}[n]))}}\right),\\ \forall & n,k,r\\[-1pt] \end{split}$

(21)

设 ${\varGamma _{\min }}$ 表示区域r的最小有效感知功率阈值。为了保证感知链路的建立同时减少其他非关联区域 $r'$ 的回声信号的干扰，需要满足最小有效感知功率约束和感知水平距离约束^[17]，即

$\quad {b_{k,r}}[n]{\varGamma _{\min }} \le {b_{k,r}}[n]{p_{k,r}}[n],\forall n,k,r$

(22)

$\quad {b_{k,r}}[n]\left\| {{{\boldsymbol{q}}_k}[n] - {{\boldsymbol{u}}_r}} \right\| \le {b_{k,r}}[n]{d_{{\text{s}},\min }},\forall n,k,r$

(23)

定义优化变量 ${\boldsymbol{A}} \triangleq \{ {a_{k,g}}[n],\forall n,k,g\}$ , ${\boldsymbol{B}} \triangleq \{ {b_{k,r}}[n],\forall n,k,r\}$ , ${\boldsymbol{P}} \triangleq \{ {p_{{\text{c}},k}}[n],{p_k}[n],\forall n,k\}$ , ${\boldsymbol{Q}} \triangleq \{ {q_{\text{c}}}[n], {q_k}[n],\forall n,k\}$ 和 ${\boldsymbol{H}} \triangleq \{ {h_k}[n],\forall n,k\}$ 。本文通过联合优化无人机集群对用户的通信关联A、对目标区域的感知关联B、无人车与无人机集群的发射功率P、无人车与无人机集群的水平轨迹Q以及无人机集群的垂直轨迹H，在用户通信、目标感知、回程链路、轨迹等约束下，最大化用户最小平均通信速率，保证通信用户的公平性。相应的问题建模为

${\mathrm{P}}1: \mathop {\max }\limits_{{\boldsymbol{A}},{\boldsymbol{B}},{\boldsymbol{P}},{\boldsymbol{Q}},{\boldsymbol{H}},\eta } \eta$

(24a)

${\mathrm{s.t.}} \;\eta \le \frac{1}{N}\sum\limits_n^N {{R_g}[n]} ,\forall g$

(24b)

$式(1)-式(7)、式(10)-式(20)、式(22)、式(23)$

(24c)

其中， $\eta$ 表示用户最小平均通信速率。上述优化问题是一个包含高度耦合变量且非凸的问题，难以求解。为此，下节将利用块坐标下降法和连续凸优化法对其进行求解。

3. 算法设计

为了方便求解，引入松弛变量 ${\boldsymbol{S}} \triangleq \{ {s_{k,{\text{c}}}}[n],{s_{k,g}}[n], {s_{k,r}}[n],\forall n,k,g,r\}$ 和 ${{\boldsymbol{R}}_{{g}}} \triangleq \{ {R_g}[n],\forall n,g\}$ ，即

${\mathrm{P}}2: \mathop {\max }\limits_{{\boldsymbol{A}},{\boldsymbol{B}},{\boldsymbol{P}},{\boldsymbol{Q}},{\boldsymbol{H}},{\boldsymbol{S}},{{\boldsymbol{R}}_{{g}}},\eta } \eta$

(25a)

$\begin{split} & {\mathrm{s.t}}.\; {R_g}[n] \le \frac{1}{{2K}}\sum\limits_{k = 1}^K {a_{k,g}}[n]\\ & \qquad\qquad\quad \cdot{{\log }_2}\left(1 + \frac{{{{\left| {{h_{k,g}}[n]} \right|}^2}{p_k}[n]}}{{{B_k}{N_0}}}\right) ,\forall n,g \end{split}$

(25b)

$\left. \begin{aligned} & {s_{k,i}}[n] \le {B_1} + {B_2}{v_{k,i}}[n],\forall n,k,i \in {{I}}\\ & 式(1)-式(7)、式(10)-式(15)、式(17)-式(20)、\\ & 式(22)、式(23)、式(24{\mathrm{b}}) \end{aligned}\right\}$

(25c)

3.1 通信感知关联优化

对于式(15)和式(19)约束，引入松弛变量 ${\bar a_{k,g}}[n]$ 和 ${\bar b_{k,r}}[n]$ ，等价转化为

${a_{k,g}}[n](1 - {\bar a_{k,g}}[n]) = 0,{a_{k,g}}[n] = {\bar a_{k,g}}[n],\forall n,k,g$

(26)

${b_{k,r}}[n](1 - {\bar b_{k,r}}[n]) = 0,{b_{k,r}}[n] = {\bar b_{k,r}}[n],\forall n,k,r$

(27)

设 ${{\bar {\boldsymbol{A}}}} \triangleq \{ {\bar a_{k,g}}[n],\forall n,k,g\}$ 和 ${{\bar {\boldsymbol{B}}}} \triangleq \{ {b_{k,r}}[n],\forall n,k,r\}$ ，在P, Q, H和S给定的情况下，P2近似表示为

${\mathrm{P}}3: \mathop {\max }\limits_{{\boldsymbol{A}},{\boldsymbol{B}},{{\boldsymbol{R}}_{{g}}},{{\bar {\boldsymbol{A}}}},{{\bar {\boldsymbol{B}}}},\eta } \eta - {\omega _\rho }$

(28a)

$\begin{split} {\mathrm{s.t.}}\;\;\rho = \,&\sum\limits_n^N \sum\limits_k^K \sum\limits_g^G ({{\left| {{a_{k,g}}[n](1 - {{\bar a}_{k,g}}[n])} \right|}^2}\\ & + {{\left| {{a_{k,g}}[n] - {{\bar a}_{k,g}}[n]} \right|}^2}) \\ & + \sum\limits_n^N \sum\limits_k^K \sum\limits_r^R ({{\left| {{b_{k,r}}[n](1 - {{\bar b}_{k,r}}[n])} \right|}^2} \\ &+ {{\left| {{b_{k,r}}[n] - {{\bar b}_{k,r}}[n]} \right|}^2}) \end{split}$

(28b)

${\omega _\rho } = \frac{1}{{2\mu }}\rho$

(28c)

$\begin{split} & 式(13)、式(14)、式(17)、式(18)、式(20)、式(22)、\\ & 式(23)、式(24{\mathrm{b}})、式(25{\mathrm{b}}) \end{split}$

(28d)

其中， $\mu$ 表示惩罚项系数， $\rho$ 表示关联变量惩罚项。当惩罚项 $\rho$ 趋近于0，式(15)和式(19)的约束满足^[17]。

将P3分解成松弛变量 ${{\bar {\boldsymbol{A}}}}$ , ${{\bar {\boldsymbol{B}}}}$ 和关联变量 ${\boldsymbol{A}}$ , ${\boldsymbol{B}}$ 这两个优化子问题，并通过交替优化求解。在 ${\boldsymbol{A}}$ , ${\boldsymbol{B}}$ , ${{\boldsymbol{R}}_{{g}}}$ 和 $\eta$ 给定的情况下，令目标函数对 ${{\bar {\boldsymbol A}}}$ 和 ${{\bar {\boldsymbol B}}}$ 的1阶偏导为0，得到最优的松弛变量形式^[17]

$\bar a_{k,g}^{{\text{opt}}}[n] = \frac{{{a_{k,g}}[n] + a_{k,g}^2[n]}}{{1 + a_{k,g}^2[n]}},\forall n,k,g$

(29)

$\bar b_{k,r}^{{\text{opt}}}[n] = \frac{{{b_{k,r}}[n] + b_{k,r}^2[n]}}{{1 + b_{k,r}^2[n]}},\forall n,k,r$

(30)

在 ${{\bar {\boldsymbol A}}}$ 和 ${{\bar {\boldsymbol B}}}$ 给定的情况下，P3可以简化为

${\mathrm{P}}4: \mathop {\max }\limits_{{\boldsymbol{A}},{\boldsymbol{B}},{{\boldsymbol{R}}_{{g}}},\eta } \eta - {\omega _\rho }\$

(31a)

$\begin{split} & {\mathrm{s.t}}. \;式(13)、式(14)、式(17)、式(18)、式(20)、\\ & \quad\;\; 式(22)、式(23)、式(24{\mathrm{b}})、式(26) \end{split}$

(31b)

问题P4是一个标准的凸优化问题，可以用内点法^[18]求解。

3.2 功率优化

在A,B,Q,H和S给定的情况下，P2可以简化为

$\qquad {\mathrm{P}}5: \mathop {\max }\limits_{{\boldsymbol{P}},{{\boldsymbol{R}}_{{g}}},\eta } \eta$

(32a)

$\begin{split} & {\mathrm{s.t}}. \;式(10)、式(11)、式(17)、式(22)、\\ & \quad\;\; 式(24{\mathrm{b}})、式(25{\mathrm{b}}) \end{split}$

(32b)

显然，P5是一个标准的凸优化问题，可以借助凸优化工具箱求解。

3.3 无人车与无人机集群水平轨迹优化

定理1　当 $2 \le \alpha \le 4$ , ${C_1},{C_2},\gamma ,x,y \gt 0$ 时， ${\psi _1}(x,y) \triangleq \left({C_1} + \dfrac{{{C_2}}}{x}\right)\dfrac{\gamma }{{{y^{\alpha /2}}}}$ 是关于x和y的凸函数。

证明　定义 ${\psi _1}^\prime (x,y) \triangleq \left({C_3} + \dfrac{{{C_4}}}{x}\right)\dfrac{1}{{{y^{\alpha /2}}}}$ , $(2 \le \alpha \le 4)$ , ${C_3} = {C_1}\gamma \ge 0,{C_4} = {C_2}\gamma \ge 0$ 。

令 ${\psi _1}^\prime (x,y)$ 的海森矩阵为 ${\nabla ^2}{\psi _1}^\prime (x,y)$ ，当 $2 \le \alpha \le 4$ 且 $x,y \gt 0$ 时，对于任意 ${\boldsymbol{t}} = {[{t_1},{t_2}]^{\text{T}}}$ ，有 ${{\boldsymbol{t}}^{\text{T}}}{\nabla ^2}{\psi _1}^\prime (x,y){\boldsymbol{t}} \ge 0$ ，即 ${\psi _1}(x,y)$ 是关于 $x,y \gt 0$ 的凸函数。证毕

为了方便处理1阶泰勒展开式，定义 ${X_{k,i}}[n] = 1 + {{\text{e}}^{ - {s_{k,i}}[n]}},\forall n,k,i \in {{I}}$ , ${Y_{k,{\text{c}}}}[n] = {\left\| {{{\boldsymbol{q}}_k}[n] - {{\boldsymbol{q}}_{\text{c}}}[n]} \right\|^2} + {h_k}{[n]^2},\forall n,k$ , ${Y_{k,i}}[n] = {\left\| {{{\boldsymbol{q}}_k}[n] - {{\boldsymbol{u}}_i}} \right\|^2} + {h_k}{[n]^2}, \forall n, k, i \in \{ g,r\}$ , ${\lambda _{k,{\text{c}}}}[n] = {\beta _0}{p_{{\text{c}},k}}[n]/({B_{{\text{c}},k}}{N_0})$ , ${\lambda _{k,g}}[n] = {\beta _0}\cdot {p_k}[n]/({B_k}{N_0})$ 。设 $X_{k,i}^l[n]$ 和 $Y_{k,i}^l[n]$ 分别为 ${X_{k,i}}[n]$ 和 ${Y_{k,i}}[n]$ 第l次迭代的值。由定理1可知， ${p_{k,r}}[n]$ 是关于 ${X_{k,i}}[n]$ 和 ${Y_{k,i}}[n]$ 的凸函数。将 ${p_{k,r}}[n]$ 在 $X_{k,i}^l[n]$ 和 $Y_{k,i}^l[n]$ 处1阶泰勒展开，得到 ${p_{k,r}}[n]$ 的全局下界 $p_{k,r}^{{\text{lb}}}[n]$ ，即

$\begin{split} {p_{k,r}}[n] \ge\,& \frac{{{\beta _0}{p_k}[n]}}{{Y_{k,r}^l[n]}}({C_1} + \frac{{{C_2}}}{{X_{k,r}^l[n]}}) \\ & + \varphi _{k,r}^l[n]({X_{k,r}}[n] - X_{k,r}^l[n]) + \zeta _{k,r}^l[n]\\ & \cdot({Y_{k,r}}[n] - Y_{k,r}^l[n]) \end{split}\tag{33}$

(33)

$\left. \begin{aligned} & \varphi _{k,r}^l[n] = - \frac{{{C_1}{\beta _0}{p_k}[n]}}{{X_{k,r}^l{{[n]}^2}Y_{k,r}^l[n]}}\\ & \zeta _{k,r}^l[n] = - \frac{{{\beta _0}{p_k}[n]}}{{Y_{k,r}^l{{[n]}^2}}}\left({C_1} + \frac{{{C_2}}}{{X_{k,r}^l[n]}}\right) \end{aligned}\right\}$

(34)

${R_{k,i}}[n]$ 是 ${X_{k,i}}[n]$ 和 ${Y_{k,i}}[n]$ 的凸函数^[11]，将 ${R_{k,i}}[n]$ 在 $X_{k,i}^l[n]$ 和 $Y_{k,i}^l[n]$ 处1阶泰勒展开，得到 $R_{k,i}^{{\text{lb}}}[n]$

$\begin{split} {R_{k,i}}[n] \ge\,& \frac{1}{{2K}}({\log _2}(1 + \frac{{{\lambda _{k,i}}[n]}}{{Y_{k,i}^l[n]}}\left({C_1} + \frac{{{C_2}}}{{X_{k,i}^l[n]}}\right) \\ & + \xi _{k,i}^l[n]({X_{k,i}}[n] - X_{k,i}^l[n]) + t_{k,i}^l[n]\\ & \cdot({Y_{k,i}}[n] - Y_{k,i}^l[n]))\\[-1pt] \end{split}$

(35)

$\begin{split} & \xi _{k,i}^l[n] =\\ & - \frac{{({{\log }_2}{\text{e}}){C_2}{\lambda _{k,i}}[n]}}{{X_{k,i}^l[n](X_{k,i}^l[n]Y_{k,i}^l[n] + {\lambda _{k,i}}[n]({C_1}X_{k,i}^l[n] + {C_2}))}} \end{split}$

(36)

$\begin{split} & t_{k,i}^l[n] =\\ & - \frac{{({{\log }_2}{\text{e}})({C_1}X_{k,i}^l[n] + {C_2}){\lambda _{k,i}}[n]}}{{Y_{k,i}^l[n](X_{k,i}^l[n]Y_{k,i}^l[n] + {\lambda _{k,i}}[n]({C_1}X_{k,i}^l[n] + {C_2}))}} \end{split}$

(37)

设 ${v_{k,i}}[n] \triangleq {h_k}[n]/{({Y_{k,i}}[n])^{1/2}},\forall i \in {{I}}$ 。 ${v_{k,i}}[n]$ 是 ${Y_{k,i}}[n]$ 的凸函数^[11]，将 ${v_{k,i}}[n]$ 在 $Y_{k,i}^l[n]$ 处展开，即

${v_{k,i}}[n] \ge \frac{{{h_k}[n]}}{{{{(Y_{k,i}^l[n])}^{1/2}}}} - \frac{{{h_k}[n]}}{{2{{(Y_{k,i}^l[n])}^{3/2}}}}({Y_{k,i}}[n] - Y_{k,i}^l[n])$

(38)

将式(4)和式(7)分别在 ${\boldsymbol{q}}_{\mathrm{c}}^l[n]$ 和 ${\boldsymbol{q}}_k^l[n] - {\boldsymbol{q}}_{k'}^l[n]$ 处1阶泰勒展开，得

$\begin{split} {d_i} + {d_{{\text{c}},\min }} \le \,& - \left\| {{\boldsymbol{q}}_{\text{c}}^l[n] - {{\boldsymbol{u}}_i}} \right\| + 2{({\boldsymbol{q}}_{\text{c}}^l[n] - {{\boldsymbol{u}}_i})^{\mathrm{T}}}\\ & \cdot({{\boldsymbol{q}}_{\text{c}}}[n] - {{\boldsymbol{u}}_i}),\forall n,i \in \{ g,r,j\} \end{split}$

(39)

$\begin{split} d_{\min }^2 \le \,& {\left| {{h_k}[n] - {h_{k'}}[n]} \right|^2} - {\left\| {{\boldsymbol{q}}_k^l[n] - {\boldsymbol{q}}_{k'}^l[n]} \right\|^2} \\ & + 2{({\boldsymbol{q}}_k^l[n] - {\boldsymbol{q}}_{k'}^l[n])^{\mathrm{T}}}({{\boldsymbol{q}}_k}[n] - {{\boldsymbol{q}}_{k'}}[n]),\forall n,k \ne k' \end{split}$

(40)

在A, B, P和H给定的情况下，P2可以近似表示为

${\mathrm{P}}6: \mathop {\max }\limits_{{\boldsymbol{Q}},{\boldsymbol{S}},{{\boldsymbol{R}}_{{g}}},\eta } \eta$

(41a)

${\mathrm{s.t.}} \;{R_g}[n] \le \frac{1}{{2K}}\sum\limits_{k = 1}^K {{a_{k,g}}[n]R_{k,{\text{c}}}^{{\text{lb}}}[n]} ,\forall n,g$

(41b)

$\quad\;\; {b_{k,r}}[n]{\varGamma _{\min }} \le {b_{k,r}}[n]p_{k,r}^{{\text{lb}}}[n],\forall n,k,r$

(41c)

$\quad\;\; {R_g}[n] \le \frac{1}{{2K}}\sum\limits_{k = 1}^K {{a_{k,g}}[n]} R_{k,g}^{{\text{lb}}}[n],\forall n,g$

(41d)

$\quad\;\; {s_{k,i}}[n] \le {B_1} + {B_2}v_{k,i}^{{\text{lb}}}[n],\forall n,k,i \in {{I}}$

(41e)

$式(1)、式(2)、式(23)、式(24{\mathrm{b}})、式(39)、式(40)$

(41f)

显然，P6是一个标准的凸优化问题，可以用内点法求解。

3.4 无人机集群垂直轨迹优化

式(7)中，在 $h_k^l[n] - h_{k'}^l[n]$ 处1阶泰勒展开，得

$\begin{split} d_{\min }^2 &\le {\left\| {{{\boldsymbol{q}}_k}[n] - {{\boldsymbol{q}}_{k'}}[n]} \right\|^2} - {\left| {h_k^l[n] - h_{k'}^l[n]} \right|^2} \\ & \quad + 2{(h_k^l[n] - h_{k'}^l[n])}({h_k}[n] - {h_{k'}}[n]),\\ & \,\forall n,k \ne k' \end{split}$

(42)

对于式(17)、式(22)和式(26)，与3.3节式(33)和式(34)同理，得到 ${p_{k,r}}[n]$ 和 ${R_{k,i}}[n]$ 的全局下界 $\varGamma _{k,r}^{{\mathrm{lb}}}[n]$ 和 $R_{k,i}^{{\mathrm{lb}}}[n]$ 。则在A,B,P和Q给定的情况下，P10可以近似表示为

${\mathrm{P}}7 \mathop {\max }\limits_{{\boldsymbol{H}},{\boldsymbol{S}},{\boldsymbol{R}}_g,\eta } \eta$

(43a)

${\mathrm{s.t.}}\; {R_g}[n] \le \frac{1}{{2K}}\sum\limits_{k = 1}^K {{a_{k,g}}[n]R_{k,{\text{c}}}^{{\text{lb}}}[n]} \forall n,g,$

(43b)

$\quad\;\; {b_{k,r}}[n]{\varGamma _{\min }} \le {b_{k,r}}[n]p_{k,r}^{{\text{lb}}}[n],\forall n,k,r$

(43c)

$\quad\;\; {R_g}[n] \le \frac{1}{{2K}}\sum\limits_{k = 1}^K {{a_{k,g}}[n]} R_{k,g}^{{\text{lb}}}[n],\forall n,g$

(43d)

$\begin{split} & \quad\;\;式(3)、式(5)、式(6)、式(17)、式(22)、\\ & \quad\;\;式(24{\mathrm{b}})、式(25{\mathrm{b}})、式(25{\mathrm{c}})、式(42) \end{split}$

(43e)

其中，式(25c)不等号右边是关于H的凹函数^[11]，因此P7可以用内点法求解。

通过交替优化迭代上面各个子问题直至收敛，完整的优化算法如算法1所示。

表 1 双层迭代算法

(1) 初始化迭代次数 ${l_1}$ =0和 ${l_2}$ =0，关联变量 ${{\boldsymbol{A}}^{{l_1},{l_2}}}$ 和 ${{\boldsymbol{B}}^{{l_1},{l_2}}}$ ，功率变量 ${{\boldsymbol{P}}^{{l_1}}}$ ，水平轨迹变量 ${{\boldsymbol{Q}}^{{l_1}}}$ ，垂直轨迹变量 ${{\boldsymbol{H}}^{{l_1}}}$ ，惩罚项系数　 ${\mu ^{{l_1},{l_2}}}$ ，目标函数 ${\eta ^{{l_1}}}$ ，收敛阈值 ${\varepsilon _1}$ >0, ${\varepsilon _2}$ >0，内层循环最大迭代次数 ${L_2}$
(2) repeat
(3) 　　令内层迭代次数 ${l_2}$ =0，repeat
(4) 　　　　给定{ ${{\boldsymbol{A}}^{{l_1},{l_2}}}$ , ${{\boldsymbol{B}}^{{l_1},{l_2}}}$ , ${{\boldsymbol{P}}^{{l_1}}}$ , ${{\boldsymbol{Q}}^{{l_1}}}$ , ${{\boldsymbol{H}}^{{l_1}}}$ }，通过式(29)和式(30)计算{ ${{{\bar {\boldsymbol A}}}^{{l_1},{l_2} + 1}}$ , ${{{\bar {\boldsymbol B}}}^{{l_1},{l_2} + 1}}$ }
(5) 　　　　给定{ ${{{\bar {\boldsymbol A}}}^{{l_1},{l_2} + 1}}$ , ${{{\bar {\boldsymbol B}}}^{{l_1},{l_2} + 1}}$ , ${\mu ^{{l_1},{l_2}}}$ }，通过求解P4得到{ ${{\boldsymbol{A}}^{{l_1},{l_2} + 1}}$ , ${{\boldsymbol{B}}^{{l_1},{l_2} + 1}}$ , ${\rho ^{{l_1},{l_2} + 1}}$ }
(6) 　　　　更新惩罚项系数 ${\mu ^{{l_1},{l_2} + 1}}$ , ${l_2} = {l_2} + 1$
(7) 　　until ${\rho ^{{l_1},{l_2}}} \le {\varepsilon _2}$ or ${l_2} \ge {L_2}$
(8) 　　更新 ${{\boldsymbol{A}}^{{l_1} + 1}} = {{\boldsymbol{A}}^{{l_1},{l_2}}}$ , ${{\boldsymbol{B}}^{{l_1} + 1}} = {{\boldsymbol{B}}^{{l_1},{l_2}}}$ ，计算 $\eta _1^{{l_1} + 1}$
(9) 　　给定{ ${{\boldsymbol{A}}^{{l_1} + 1}}$ , ${{\boldsymbol{B}}^{{l_1} + 1}}$ , ${{\boldsymbol{Q}}^{{l_1}}}$ , ${{\boldsymbol{H}}^{{l_1}}}$ }，通过求解P5得到{ ${{\boldsymbol{P}}^{{l_1} + 1}}$ , $\eta _2^{{l_1} + 1}$ }
(10) 　给定{ ${{\boldsymbol{A}}^{{l_1} + 1}}$ , ${{\boldsymbol{B}}^{{l_1} + 1}}$ , ${{\boldsymbol{P}}^{{l_1} + 1}}$ , ${{\boldsymbol{H}}^{{l_1}}}$ }，通过求解P6得到{ ${{\boldsymbol{Q}}^{{l_1} + 1}}$ , $\eta _3^{{l_1} + 1,{\text{lb}}}$ }
(11) 　给定{ ${{\boldsymbol{A}}^{{l_1} + 1}}$ , ${{\boldsymbol{B}}^{{l_1} + 1}}$ , ${{\boldsymbol{P}}^{{l_1} + 1}}$ , ${{\boldsymbol{Q}}^{{l_1} + 1}}$ }，通过求解P7得到{ ${{\boldsymbol{H}}^{{l_1} + 1}}$ , $\eta _4^{{l_1} + 1,{\text{lb}}}$ }
(12) 　更新 ${l_1} = {l_1} + 1$ , ${{\boldsymbol{A}}^{{l_1},0}} = {{\boldsymbol{A}}^{{l_1} - 1}}$ , ${{\boldsymbol{B}}^{{l_1},0}} = {{\boldsymbol{B}}^{{l_1} - 1}}$ , ${\mu ^{{l_1},0}}$ ，计算 ${\eta ^{{l_1}}}$
(13) until $\dfrac{{{\eta ^{{l_1}}} - {\eta ^{{l_1} - 1}}}}{{{\eta ^{{l_1} - 1}}}} \lt {\varepsilon _1}$

下载: 导出CSV

| 显示表格

由的步骤(3)–步骤(12)可得， ${\eta ^{{l_1}}}\mathop \le \limits^{\mathrm{a}} \eta _1^{{l_1} + 1}\mathop \le \limits^{\mathrm{b}} \eta _2^{{l_1} + 1}\mathop = \limits^{\mathrm{c}} \eta _2^{{l_1} + 1,{\text{lb}}}\mathop \le \limits^{\mathrm{d}} \eta _3^{{l_1} + 1,{\text{lb}}}\mathop \le \limits^{\mathrm{e}} \eta _4^{{l_1} + 1,{\text{lb}}}\mathop \le \limits^{\mathrm{f}} {\eta ^{{l_1} + 1}}$ 。其中不等号a表示在的步骤(3)–步骤(7)中，惩罚项趋近于0时，目标值具有非减性；不等号b,d和e分别表示在步骤(9)–步骤(11)中，子问题P5, P6和P7的目标函数max(·)具有非减性。等号c表示在步骤(9)中，1阶泰勒展开式在给定局部点具有相同的值。不等号f表示在步骤(11)中1阶泰勒展开得到的 $\eta _4^{{l_1} + 1,{\text{lb}}}$ 是 ${\eta ^{{l_1} + 1}}$ 的下界。因此每次迭代的目标函数值是非减的，且存在上界，因此是收敛的。假设求解器精度为 $\varepsilon$ >0^[18]，算法1的时间复杂度为O((N(1+G+K+KG+KR))log₂(1/ $\varepsilon$ ))。

4. 仿真结果分析

本节给出仿真结果来验证本文提出方案联合优化算法的性能，并且与以下3种方案作对比：

方案1　固定通信感知关联和功率方案，优化无人车和无人机集群的轨迹。将用户和目标区域分成K簇。无人车和无人机在纯通信时隙分配相等的发射功率，在感知时隙，无人机保持峰值发射功率。

方案2　固定无人机集群高度方案，优化无人机集群的通信感知关联、发射功率和水平轨迹以及无人车的发射功率和行进轨迹，即无人机集群保持相同的高度，即 ${h_k}[n] = {h_{k,{\text{I}}}},\forall n$ 。

方案3　固定无人车基站方案，优化无人机集群的通信感知关联、功率和3维轨迹以及无人车的发射功率。在这个方案中，无人车基站保持静止，即， ${{\boldsymbol{q}}_{\text{c}}}[n] = {{\boldsymbol{q}}_{\text{I}}},\forall n$ 。

对比方案1和方案2的优化算法复杂度均为O((N(1+G+K+KG+KR))log₂(1/ $\varepsilon$ ))，方案3优化算法复杂度为O((N(G+max(K+1,R)+KG+KR))log₂(1/ $\varepsilon$ ))。在仿真中，工作区域的长宽均为1 250 m，无人车初始位置为(625,0,0) m，无人机数量K为3，用户数G为6，用户最大半径0.5 m，目标区域数R为4，目标区域最大半径为1 m，障碍物数J为5，障碍物最大半径为30 m。3架无人机的初始高度分别为80 m, 90 m, 100 m，无人车峰值发射功率0.15 W，平均发射功率0.015 W，无人机峰值发射功率0.12 W，平均发射功率0.024 W。无人车和无人机最大水平速度分别15 m/s和30 m/s。无人机最大垂直速度为20 m/s，最大最小飞行高度分别为300 m和50 m，无人机感知目标区域的最大水平距离为200 m。系统总带宽为1 MHz，任务总周期为160 s，时隙长度为0.5 s，高斯白噪声功率谱密度N₀= –169 dBm/Hz，无人车防碰撞距离为20 m，无人机防碰撞距离为10 m，有效感知功率阈值为–95 dBm，任务周期内每个目标区域的感知频数为20次，参考距离1 m的信噪比为-60 dB，莱斯信道参数^[11]B₁=–4.3221, B₂=6.0750, C₁=0, C₂=1。设中 ${\varepsilon _1}$ =1e–4, ${\varepsilon _2}$ =1e–5, ${L_2}$ =10。

给出了算法1目标函数和惩罚项的迭代情况。在不同的 ${\varGamma _{\min }}$ 和 ${f_r}$ 条件下，所提出算法迭代到30次左右，惩罚项趋近于0，且目标函数值收敛，说明在不同参数下，算法1可以达到良好的收敛效果。

图 2 所提算法中目标函数和惩罚项与迭代次数的关系

下载: 全尺寸图片幻灯片

是在 ${\varGamma _{\min }}$ = –95 dBm和 ${f_r}$ =20 次时，算法1得到的无人机中继的通信感知关联情况。可以看出，感知任务较少的无人机会分配通信时隙给感知任务较多的无人机所关联的用户。

图 3 无人机集群通信感知关联

下载: 全尺寸图片幻灯片

给出了当 ${\varGamma _{\min }}$ = –95 dBm和 ${f_r}$ =20 次时，算法1得到的无人车与无人机集群的水平轨迹和3维轨迹。在图4(a)中，无人车依次为无人机提供信道质量较好的回程链路，无人机飞往关联的目标区域执行感知任务，在感知完后，靠近关联的用户和无人车，获得较好的信道质量。从图4(b)中看出，无人车和无人机动态调整位置，改变回程链路、通信链路和感知链路的信道质量，协调通信感知性能之间的权衡。

图 4 无人车和无人机集群的水平轨迹和3维轨迹

下载: 全尺寸图片幻灯片

给出了当 ${f_r}$ =20 次时几种算法用户最小平均通信速率和有效感知功率阈值的关系。随着有效感知功率阈值的减小，用户最小平均通信速率越来越大，但在–105 dBm以后受限于感知频率等因素，减小有效感知功率阈值，用户最小平均通信速率提升不大。所提优化算法的用户最小平均通信速率均比其他3个对比算法明显提高。因为所提优化算法方案拥有更多的优化自由度，实现通信感知性能之间更好的权衡。

图 5 用户最小平均通信速率随有效感知功率阈值的变化

下载: 全尺寸图片幻灯片

给出了当 ${\varGamma _{\min }}$ =–95 dBm时几种算法的用户最小平均通信速率和目标区域感知次数的关系。所提优化算法的用户最小平均通信速率均比其他对比方案高。从图5和图6可以看出，与其他3种优化方案相比，所提优化算法虽然增加了一定的复杂度(仍是可接受的多项式复杂度)，但显著提高了用户最小平均通信速率，从而验证了所提优化算法的有效性和必要性。

图 6 用户最小平均通信速率随感知频率的变化

下载: 全尺寸图片幻灯片

5. 结论

本文提出一种以空地协同的方式辅助通信感知一体化系统的方法，通过优化无人机中继集群对用户的通信关联和目标区域的感知关联、发射功率和飞行轨迹，以及无人车基站的发射功率和行进轨迹，最大化用户最小平均通信速率。基于块坐标下降和连续凸优化法，提出一种双层迭代的算法求解次优解。仿真结果表明，所提优化算法在不同参数下均有良好的收敛性。此外，相比于无人机集群高度固定、通信感知关联和发射功率固定以及回程基站静止这3种基准方案，所提算法可以实现通信感知一体化系统中通信感知任务关联的均衡，并且动态调整空地协同网络的位置和发射功率，使之与通信感知一体化系统的通信感知任务相适应，协调通信感知性能之间的权衡，使得在相同感知性能下，有效提高了用户最小平均通信速率。

图 1 ESF算法加密流程

下载: 全尺寸图片幻灯片

图 2 ESF算法轮函数

下载: 全尺寸图片幻灯片

图 3 ESF算法的11轮相关密钥不可能差分区分器

下载: 全尺寸图片幻灯片

图 4 ESF算法的15轮相关密钥不可能差分路径

下载: 全尺寸图片幻灯片

表 1 符号约定

符号	意义
$K$	80 bit主密钥
${K_i}$	第 $i\,$ 轮的32 bit轮密钥
${K_{i, j}}$	${K_i}$ 的第 $j$ 个半字节
$K_{i, j}^l$	${K_{i, j}}$ 的第 $l$ 位
${L_i}$	第 $i\,$ 轮输出密文的左边32 bit
${R_i}$	第 $i\,$ 轮输出密文的右边32 bit
$< < < 7$	循环左移7位
$\oplus$	按位异或运算符
$\|\|$	二进制字符联接
${[i]_2}$	常数 $i\,$ 的二进制表示

下载: 导出CSV

表 2 15轮相关密钥差分路径

$\Delta K = (00000200000000000000)$
$\Delta {K_1}$	$0000 0200$	$\Delta {K_9}$	$0000 0000$
$\Delta {K_2}$	$0040 0000$	$\Delta {K_{10}}$	$0000 0000$
$\Delta {K_3}$	$0000 0000$	$\Delta {K_{11}}$	$0000 0000$
$\Delta {K_4}$	$0000 0000$	$\Delta {K_{12}}$	$0000 0000$
$\Delta {K_5}$	$0000 0000$	$\Delta {K_{13}}$	$0000 0020$
$\Delta {K_6}$	$0000 0000$	$\Delta {K_{14}}$	$0004 0000$
$\Delta {K_7}$	$0000 0080$	$\Delta {K_{15}}$	$*000 0000$
$\Delta {K_8}$	$0010 0000$	–	–

下载: 导出CSV

表 3 ESF算法的攻击结果比较

攻击方法	轮数	时间复杂度	数据复杂度	文献
不可能差分	11	${2^{75.5}}$	${2^{59}}$	[4]
不可能差分	11	${2^{32}}$	${2^{53}}$	[5]
不可能差分	12	${2^{60.43}}$	${2^{53}}$	[6]
相关密钥差分	13	${2^{66}}$	${2^{47}}$	[8]
相关密钥不可能差分	15	${2^{40.5}}$	${2^{61.5}}$	本文

下载: 导出CSV

参考文献(16)

WU Wenling and ZHANG Lei. LBlock: A lightweight block cipher[C]. Proceedings of 9th International Conference on Applied Cryptography and Network Security, Nerja, Spain, 2011: 327–344. doi: 10.1007/978-3-642-21554-4_19.

IZADI M, SADEGHIYAN B, SADEGHIAN S, et al. MIBS: A new light-weight block cipher[C]. Proceedings of CANS 2009, Ishikawa, Japan, 2009: 334–348. doi: 10.1007/978-3-642-10433-6_22.

BOGDANOV A, KNUDSEN L, LEANDER G, et al. PRESENT: An ultra-lightweight block cipher[C]. Proceedings of Cryptographic Hardware and Embedded Systems, Vienna, Austria, 2007: 450–466. doi: 10.1007/978-3-540-74735-2_31.

刘宣, 刘枫, 孟帅. 轻量级分组密码算法ESF的不可能差分分析[J]. 计算机工程与科学, 2013, 35(9): 89–95. doi: 10.3969/j.issn.1007-130X.2013.09.014

LIU Xuan, LIU Feng, and MENG Shuai. Impossible differential cryptanalysis of lightweight block ciper ESF[J]. Computer and Engineering Science, 2013, 35(9): 89–95. doi: 10.3969/j.issn.1007-130X.2013.09.014

陈玉磊, 卫宏儒. ESF算法的不可能差分密码分析[J]. 计算机科学, 2016, 43(8): 89–91. doi: 10.11896/j.issn.1002-137X.2016.8.018

CHEN Yulei and WEI Hongru. Impossible differential cryptanalysis of ESF[J]. Computer Science, 2016, 43(8): 89–91. doi: 10.11896/j.issn.1002-137X.2016.8.018

高红杰, 卫宏儒. 用不可能差分法分析12轮ESF算法[J]. 计算机科学, 2017, 44(8): 147–150. doi: 10.11896/j.issn.1002-137X.2017.10.028

GAO Hongjie and WEI Hongru. Impossible differential attack on 12-round block cipher ESF[J]. Computer Science, 2017, 44(8): 147–150. doi: 10.11896/j.issn.1002-137X.2017.10.028

尹军, 马楚炎, 宋健, 等. 轻量级分组密码算法ESF的安全性分析[J]. 计算机研究与发展, 2017, 54(10): 2224–2231. doi: 10.7544/issn1000-1239.2017.20170455

YIN Jun, MA Chuyan, SONG Jian, et al. Security analysis of lightweight block cipher ESF[J]. Journal of Computer Research and Development, 2017, 54(10): 2224–2231. doi: 10.7544/issn1000-1239.2017.20170455

尹军, 宋健, 曾光, 等. 轻量级分组密码算法ESF的相关密钥差分分析[J]. 密码学报, 2017, 4(4): 333–344. doi: 10.13868/j.cnki.jcr.000186

YIN Jun, SONG Jian, ZENG Guang, et al. Related-key differential attack on lightweight block cipher ESF[J]. Journal of Cryptologic Research, 2017, 4(4): 333–344. doi: 10.13868/j.cnki.jcr.000186

KNUDSEN L. Crypatanalysis of LOKI[C] Proceedings of Advances in Cryptology, Gold Coast, Australia, 1991: 22–35.

BIHAM E. New types of cryptanalytic attacks using related keys[J]. Journal of Cryptology, 1994, 7(4): 229–246. doi: 10.1007/BF00203965

BIHAM E, BIRYUKOV A, and SHAMIR A. Cryptanalysis of Skipjack reduced to 31 rounds using impossible differentials[C]. Proceedings of Advances in Cryptolog EUROCRYPT'99. Prague, CZ, 1999: 12–23. doi: 10.1007/3-540-48910-x_2.

JIANG Zilong and JIN Chenhui. Impossible differential cryptanalysis of 8-round Deoxys-BC-256[J]. IEEE Access, 2018, 6: 8890–8895. doi: 10.1109/ACCESS.2018.2808484

徐洪, 苏鹏晖, 戚文峰. 减轮SPECK算法的不可能差分分析[J]. 电子与信息学报, 2017, 39(10): 2479–2486. doi: 10.11999/JEIT170049

XU Hong, SU Penghui, and QI Wenfeng. Impossible differential cryptanalysis of reduced-round SPECK[J]. Journal of Electronics &Information Technology, 2017, 39(10): 2479–2486. doi: 10.11999/JEIT170049

付立仕, 金晨辉. MIBS-80的13轮不可能差分分析[J]. 电子与信息学报, 2016, 38(4): 848–855. doi: 10.11999/JEIT150673

FU Lishi and JIN Chenhui. Impossible differential cryptanalysis on 13-round MIBS-80[J]. Journal of Electronics &Information Technology, 2016, 38(4): 848–855. doi: 10.11999/JEIT150673

XIE Min, LI Jingjing, and ZANG Yuechuan. Related-key impossible differential cryptanalysis of LBlock[J]. Chinese Journal of Electronics, 2017, 26(1): 35–41. doi: 10.1049/cje.2016.06.031

CHENG Lu, XU Peng, and WEI Yuechuan. New related-key impossible differential attack on MIBS-80[C]. Proceedings of 2016 International Conference on Intelligent Networking and Collaborative Systems, Ostrawva, CZ, 2016: 203–206. doi: 10.1109/incos.2016.41.

施引文献

期刊类型引用(2)

1.	潘钰，胡航，金虎，雷迎科，冯辉，姜丽，张孟伯. 非授权频段下无人机辅助通信的轨迹与资源分配优化. 电子与信息学报. 2024(11): 4287-4294 . 本站查看
2.	郭丹. 无线传感器网络簇间通信信息自适应节能优化方法. 长江信息通信. 2024(11): 104-106 . 百度学术

其他类型引用(0)

资源附件(0)

访问统计

图(4) / 表(3)

计量

文章访问数: 2943
HTML全文浏览量: 1021
PDF下载量: 90
被引次数: 2

1. 引言
2. 系统模型
3. 算法设计
3.1 通信感知关联优化
3.2 功率优化
3.3 无人车与无人机集群水平轨迹优化
3.4 无人机集群垂直轨迹优化
4. 仿真结果分析
5. 结论

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

轻量级分组密码算法ESF的相关密钥不可能差分分析

doi: 10.11999/JEIT180576

作者简介:
谢敏：女，1976年生，副教授，研究方向为编码和密码

曾琦雅：女，1993年生，硕士，研究方向为分组密码算法分析

通讯作者:
谢敏　mxie@xidian.edu.cn

计量

Related-key Impossible Differential Cryptanalysis on Lightweight Block Cipher ESF

1. 引言

2. 系统模型

3. 算法设计

3.1 通信感知关联优化

3.2 功率优化

3.3 无人车与无人机集群水平轨迹优化

3.4 无人机集群垂直轨迹优化

4. 仿真结果分析

5. 结论

期刊类型引用(2)

其他类型引用(0)

计量

目录

1. 引言

2. 系统模型

3. 算法设计

3.1 通信感知关联优化

3.2 功率优化

3.3 无人车与无人机集群水平轨迹优化

3.4 无人机集群垂直轨迹优化

4. 仿真结果分析

5. 结论

留言板

轻量级分组密码算法ESF的相关密钥不可能差分分析

doi: 10.11999/JEIT180576

作者简介: 谢敏：女，1976年生，副教授，研究方向为编码和密码 曾琦雅：女，1993年生，硕士，研究方向为分组密码算法分析

通讯作者: 谢敏 mxie@xidian.edu.cn

计量

出版历程

Related-key Impossible Differential Cryptanalysis on Lightweight Block Cipher ESF

1. 引言

2. 系统模型

3. 算法设计

3.1 通信感知关联优化

3.2 功率优化

3.3 无人车与无人机集群水平轨迹优化

3.4 无人机集群垂直轨迹优化

4. 仿真结果分析

5. 结论

期刊类型引用(2)

其他类型引用(0)

计量

出版历程

目录

1. 引言

2. 系统模型

3. 算法设计

3.1 通信感知关联优化

3.2 功率优化

3.3 无人车与无人机集群水平轨迹优化

3.4 无人机集群垂直轨迹优化

4. 仿真结果分析

5. 结论

作者简介:
谢敏：女，1976年生，副教授，研究方向为编码和密码

曾琦雅：女，1993年生，硕士，研究方向为分组密码算法分析

通讯作者:
谢敏　mxie@xidian.edu.cn