用于16位音频加密的三维连续时间混沌系统设计与OGY稳定控制

《IEEE Access》：Design and OGY Stabilization of a 3-D Continuous-Time Chaotic System for 16-bit Audio Encryption

　　

摘要翻译：连续时间混沌系统的探索仍是非线性动力学中的一个关键前沿领域，对安全通信、密码学和复杂系统控制具有深远意义。本研究对一种新提出的三维（3D）连续时间混沌系统进行了全面研究，该系统以独特的双曲正切和绝对值非线性为特征。研究人员分析了该系统的数学结构，以确定其基本动力学性质，包括耗散性、相空间拓扑、平衡点结构、李雅普诺夫指数谱和分岔行为。分析揭示了两个满足必要Shilnikov型条件的非平凡鞍-焦点，以及一个自激的、有效单稳态混沌吸引子。为展示可控性，研究人员采用了一种完全数值化的、映射级线性化方法（适用于非光滑流）实现了Ott-Grebogi-Yorke（OGY）方法；混沌轨迹被稳定到目标不稳定周期轨道（UPO）上，控制能量低至8.7952单位。此外，研究人员利用混沌动力学开发了一种基于明文相关的置换-扩散架构的安全数字音频加密算法。对语音、音乐和环境数据集的安全分析显示：相关系数接近零，16位空间中的香农熵为13.14，NPCR和UACI满足相应统计准则，表现出对已知/选择明文攻击的抵抗能力，密钥空间达2¹⁸⁶。该框架还被证明可扩展至8位图像加密。

论文解读文章

一、研究背景与问题

非线性动力学与混沌理论为理解复杂系统提供了数学框架，混沌系统固有的初值敏感性和不可预测性使其在安全通信与密码学中具有重要应用价值。近年来，三维混沌系统的研究已从光滑多项式非线性拓展至包含非光滑与超越非线性的结构，以期增强系统不可预测性。然而，现有用于加密的三维混沌系统均为光滑多项式非线性，其庞加莱映射可解析求解雅可比矩阵，经典OGY控制方法可直接应用。当非线性项不可微且切换流形与庞加莱截面重合时，经典OGY框架的适用性成为悬而未决的方法论问题。此外，混沌控制在加密应用中的双重角色——既需保证混沌发生器稳定运行，又不能在加密时施加控制以免破坏密钥流的随机性——尚未被系统阐述。同时，音频数据体量大、冗余度高，传统分组密码存在延迟和效率问题，亟需基于混沌的轻量级加密方案。本研究正是为解决上述非光滑系统中OGY控制的理论缺口，并构建实用的音频加密框架而开展。

二、研究内容与结论

研究人员提出了一种新的三维连续时间混沌系统，其状态方程包含双曲正切tanh(z)和绝对值|y|非线性。在名义参数p=2下，该系统呈现耗散性（散度为-1.4），具有三个平衡点：原点E0为指数1鞍-焦点，两个非平凡平衡点E1和E2为指数2鞍-焦点，且满足Shilnikov型必要条件（|λ|>ρ）。数值模拟表明系统呈现自激的、有效单稳态混沌吸引子，最大李雅普诺夫指数LE1=0.0523，Kaplan-Yorke维数为2.036。分岔分析显示p=2位于宽混沌带内，远离周期窗口边界，具有参数鲁棒性。研究人员通过完全数值化的映射级OGY线性化方法，克服了非光滑切换流形导致的解析雅可比不可得问题，成功将混沌轨迹稳定到目标UPO上，控制能量仅8.7952单位。基于该系统构建的明文相关置换-扩散音频加密算法，在语音、音乐和环境声音数据集上实现了近零相关、香农熵13.14、NPCR和UACI满足统计准则、抵抗已知/选择明文攻击等安全指标，密钥空间达2¹⁸⁶，并可扩展至8位图像加密。论文发表于《IEEE Access》。

三、关键技术方法

本研究采用的主要技术方法包括：（1）非线性动力学分析，通过雅可比矩阵特征值、散度、Lyapunov指数谱（Benettin/QR算法）和分岔图表征系统动力学特性；（2）数值映射级OGY混沌控制，在y=0庞加莱截面上通过递归搜索定位UPO，用中心差分（ε=10^-5）数值估计线性化矩阵A和B，进而计算控制增益Cx、Cz，施加有界参数脉冲δp使轨迹落入稳定流形；（3）明文相关密钥调度，利用HMAC-SHA-256密钥化摘要扰动初始条件和分岔参数（扰动幅度≤5×10^-4），使每个明文对应不同混沌轨迹；（4）置换-扩散加密架构，使用y序列排序索引置换明文采样，XOR混合x和z投影生成16位密钥流，再经链式XOR（含初始化向量IV）完成扩散。样本来源包括Free Spoken Digit Dataset（语音）、ESC-50（环境声音）和MATLAB内置音乐录音，图像测试使用标准cameraman图像。

四、主要研究结果

A. 平衡点与Shilnikov型分析
通过分支解析（y>0和y<0）和数值多起始点Newton搜索，确定系统在p=2时存在三个平衡点：原点E0和两个非平凡点E1=(0.6511,0.4653,0.3733)、E2=(-0.7023,0.3417,-0.6508)。E1和E2的特征值构型满足Shilnikov型必要条件。研究还解析证明了y<0分支无平衡点：令y=-u（u>0），推导出x(u)和z(u)的显式表达式，并利用tanh(z)≤z（z≥0）的性质证明g(u)>0恒成立（当p>1.6时）。对61×61网格初始条件的积分显示所有轨迹均收敛到同一混沌吸引子，排除了多稳定性，证明系统为有效单稳态自激吸引子。

B. 参数选择与分岔机制
通过计算最大Lyapunov指数LE1(p)（p∈[0.1,2.55]）并与分岔图对齐，揭示了通向混沌的路径：约p≈0.12进入混沌带，随后周期窗口与混沌带交替；最强混沌出现在p≈1.33（LE1=0.0874）。名义值p=2被置于宽混沌带内，距最近混沌边界（p=1.9158）有0.0842的余量，是密钥调度最大扰动幅度（5×10^-4）的168倍，从而保证了密钥调度过程中系统始终处于混沌状态。结构系数鲁棒性测试显示，在±5%扰动下98.7%样本保持混沌。

C. 与经典三维混沌系统的定量比较
在相同数值设置下（同一Benettin/QR实现、相同积分步长和平均窗口），将本系统与Lorenz、Chen、Lü、Sprott B/E系统比较。虽然本系统原始最大Lyapunov指数较小，但归一化到收缩率后，分形维数与经典系统同量级。本系统的独特之处在于同时具备：非光滑绝对值与超越双曲正切非线性、切换流形与庞加莱截面重合、极强收缩（稳定乘子量级10^-10，回映条件数~10⁶），这些性质的组合使得解析雅可比不可得，必须采用数值映射级线性化。

D. OGY混沌控制的实现与性能
在y=0截面定义离散回映Xn+1=P(Xn,p)，通过回归搜索定位周期-1固定点，用ε=10^-5的中心差分估计矩阵A和B。控制律为δpn=Cx(xn-x*)+Cz(zn-z*)，激活窗口为||δXn||-9，故λs=-7.0×10^-10，这解释了极低控制能量——一次参数脉冲将状态投影到稳定流形后，偏差在一个回映内坍缩至UPO。与延迟反馈控制的定量比较显示，延迟反馈在允许增益区间(-1,-0.9748)内最小闭环乘子为0.9874，百倍衰减需363次回映（约5200 s），而OGY控制单次回映内收敛。自适应辨识方法因回映矩阵条件数达10⁶而不可行。滑模控制因作用通道不同未被纳入比较。

E. 音频加密架构与安全分析
加密算法采用明文相关密钥调度：对音频数据计算HMAC-SHA-256密钥化摘要，扰动初始条件(x0,y0,z0)和参数p（|δi|≤5×10^-4）。蒙特卡洛实验（500次）中所有样本保持LE1>0，零失败时95%单侧置信上限为0.60%的脱混沌概率。丢弃瞬态设为3×10⁵积分步（1500 s），满足最坏情况下LE1≥0.024时放大10^-14扰动所需时间。音频信号归一化至[-1,1]后量化为16位无符号整数；用y序列排序索引置换，密钥流为x和z投影取模2¹⁶后异或；最后进行链式XOR扩散（C0由密钥和摘要派生的IV）。在语音、音乐、环境声音三个10000样本数据集上，加密后相邻样本相关系数接近零，香农熵13.14（16位空间理论值16），NPCR和UACI采用精确泊松尾部准则（N=10⁴时允许1个重合）均通过；N=10⁶时同样通过。位级NBCR为49.98±0.14%，密钥灵敏度（扰动10^-14）样品级100%、位级50%。已知/选择明文攻击测试中，从一对明密文恢复的密钥流用于解密其他消息时仅恢复0.00%样本，与随机猜测率1/2¹⁶一致。

F. 密钥空间与攻击抵抗
密钥包含四个实值分量(x0,y0,z0,p)，灵敏度10^-14，密钥空间(10¹⁴)⁴=10⁵⁶≈2¹⁸⁶，远超2¹²⁸基准。明文相关密钥调度使得单对明密文无法揭示其他消息的密钥流，实现了对已知/选择明文攻击的抵抗。该方案提供机密性但不提供认证，如需完整性保护可采用encrypt-then-MAC结构（额外成本4.2%）。

G. 计算性能与特征对比
密钥流生成复杂度：积分O(N)、排序置换O(NlogN)、扩散O(N)，固定瞬态成本与消息长度无关。N=10³时有效速率4412样本/s，低于测试信号8192 Hz采样率，不满足实时性；N≥10⁴时达到数倍实时。与AES-256-CTR相比慢约429倍，但本方案的价值在于控制与生成的耦合而非吞吐量。

H. 图像加密扩展
将密钥流模数从2¹⁶降至2⁸，其余结构不变，加密256×256 cameraman图像：熵7.9973，相邻像素相关系数（水平-0.0051、垂直-0.0054、对角-0.0051），NPCR=99.6399%（允许临界值99.5693%），UACI=33.6027%（接受区间[33.2824,33.6447]%），直方图卡方=248.20<293.25。证明框架可有效迁移至图像数据。

五、讨论与结论

讨论部分指出，OGY控制分析认证了伪随机发生器安全运行的参数区域；加密时故意不施加控制，因为稳定周期-1轨道会生成周期性密钥流。控制实验展示的是同一架构可按需驱动至运行稳定，而在自由运行时保持高熵输出。与替代控制策略比较，OGY因极强体积收缩实现单次回映收敛，优于延迟反馈控制（需数千秒）；自适应辨识因回映矩阵条件数过大而失败，这恰好证明了离线数值线性化的必要性。研究局限性在于：非光滑绝对值和双曲正切非线性的硬件可行性未经模拟电路实验验证，但数字硬件实现是合适的；置换阶段需要缓冲整个消息，是硬件实现中最不友好的部分。未来工作将聚焦FPGA或DSP硬件实现、实时加密性能评估，以及扩展至视频段加密和更广泛的图像加密场景。

研究结论：本研究分析了一种具有非光滑和超越非线性的新型3D混沌系统，验证了其满足Shilnikov条件的非平凡平衡结构、有效单稳态自激吸引子、正最大Lyapunov指数和分形维数。针对非光滑回映开发的数值映射级OGY实现，以高能效成功稳定了混沌状态。最后，研究开发了鲁棒的明文相关置换-扩散音频加密应用，在语音、音乐和环境数据集上实现了卓越的安全指标，展示了对已知/选择明文攻击的抵抗能力和2¹⁸⁶密钥空间，表明该系统在下一代安全通信中的潜力。该框架还被证明可自然扩展至8位图像加密。