简介一份围绕光通信物理层安全展开的专业文档面向通信工程、光网络与信息安全领域的研究者、工程师及高年级学生也可作为相关课程设计与毕业设计的参考。资源从5G时代光纤传输面临的安全挑战入手指出上层加密在保护范围与密钥管理上的短板系统阐述基于高非线性光纤简并四波混频效应与超结构光纤布拉格光栅的混合光加密方案涵盖QPSK信号加解密流程、密钥隐匿传输机制、SSFBG编解码原理及系统结构设计并结合相位矢配条件、弱耦合近似等理论分析说明方案可行性。文档还给出了VPI仿真平台的参数配置、误码率测试曲线以及加密前后信号对比便于读者掌握从理论推导、系统搭建到性能验证的完整链路。资源为1个docx文档大小668KB内容精炼且公式图表齐全已有89人学习可作为课题研究、毕业设计或光通信安全技术预研的参考资料。1. 混合光加密正在接管光链路的物理层安全传统密码学保护光网络时加密和解密都发生在数据链路层光脉冲进入光纤后以“明文”物理波形传输一根分光器复制光信号不需要破解任何上层协议。混合光加密的思路是把混沌密钥生成与相位掩膜处理放到光域在发送端同时改变信号的幅度与相位分布合法接收端凭密钥恢复原始符号窃听端只能看到接近噪声的基带图样。这个方案适合想补齐物理层安全短板的传输工程师也适合对光链路做安全评估的网络规划人员。下面的内容会从理论边界起步逐步推进到最小可复现仿真、参数评估和落地前的验证方法用一个连续链路把原理和指标串起来。2. 混合光加密的物理层实现框架与光域安全边界2.1 为什么在电域完成加密后仍要关心物理层安全在真实光网络里窃听者并不需要破译管理面或数据面的密码系统。最常见的手段是在光缆护套开一个口子引入分光器或在某一段制造微弯窃取一小部分光功率这种截取方式几乎不会引发上层告警。此时上层加密保护的是载荷内容但光层信号的时间结构、波长位置、调制格式都以明文暴露给攻击者。更关键的是骨干网节点之间的光信号是连续波逐跳解密再加密的成本极高因此“光纤上走的是密文光纤本身却是透明的”这个矛盾长期存在。电域加密还有一个结构性问题OFDM 或 QAM 这类高阶调制在光纤中传输时星座图、峰均比、循环前缀这些物理特征是明文暴露的。攻击者可以使用机器学习识别调制类型与帧结构辅助流量分析和定向干扰。物理层安全的目标不是替代 AES 这类算法而是抬高窃听信道的代价让合法信道与窃听信道之间形成确定的容量差。混合光加密的价值正在于它把破解难度从“已知密文且已知算法”的标准密码分析问题转换成了“未知信道模型加未知随机映射”的估计问题对想获取明文的一方来说这两类问题所需的工具完全不在一个数量级。2.2 混合光加密的三个组成面混沌密钥、相位掩膜与符号映射“混合”二字的含义是几种手段在同一个发送端叠加使用。我一般把发送端拆成三个面来设计。密钥面由混沌映射产生伪随机序列作为相位掩膜的种子和动态扰码的钥流掩膜面把相位掩膜乘到频域复信号上使星座点的相位分布整体打散映射面将原始比特通过密钥控制的映射关系从固定 QAM 星座切换到随帧变化的星座子集。三者的结合在光域上的表现就是信号的模拟波形随密钥变化而变化静态星座图不复存在。攻击者在有限截获样本上做的聚类、调制识别和符号估计都会因为映射关系逐帧漂移而失效。下表方便对照三个组成面的作用边界组成面典型实现作用边界不设置时的影响混沌密钥Logistic 映射、Chen 混沌决定密钥空间与帧同步要求掩膜与映射均可被识别相位掩膜随机相位板、Hadamard 矩阵对抗幅度统计与频谱估计星座相位可直接观测符号映射动态 QAM 子集选择对抗聚类与模式识别窃听端可直接硬判决在系统实现顺序上混沌密钥先于其他两面生成它既驱动掩膜也决定映射表。映射面的工程实现一般用一张密钥索引的置换表把输入比特组映射到 M 个星座点的非自然顺序。这样即使攻击者截获到一段没有相位掩膜的导频也无法利用导频建立可信的星座骨架。2.3 物理层安全收益怎么算安全容量与窃听信噪比评估“安全”不能只看加密算法强度工程上更多以信道不对称性来度量。Wyner 安全容量的核心表达式为合法信道容量减窃听信道容量安全容量 C_s max{0, C_legit - C_eve}每个容量都按香农公式 C B * log2(1 SNR) 计算。光链路中合法接收机与窃听者从同一根光纤分支出去物理层的唯一差异是接收光功率OSNR和相位噪声累积。混合光加密真正改善的是窃听者解调后的等效信噪比掩膜和动态映射使其恢复星座时产生一个叠加惩罚而这个惩罚转化为安全收益用来抵消加密给合法用户带来的消光比损失。下面这段代码可以直接估算给定合法与窃听 OSNR 时的安全容量% 安全容量估算B 为信号带宽10 GHz snrLegitDb 18; snrEveDb 8; B 10e9; Clegit B * log2(1 10^(snrLegitDb/10)); Ceve B * log2(1 10^(snrEveDb/10)); fprintf(安全容量约 %.2f Gbit/s\n, max(0, Clegit - Ceve));这段代码的参数说明snrLegitDb 和 snrEveDb 分别是合法与窃听接收机在解调前的信噪比B 是子载波聚合后的有效带宽。注意这里把窃听信道容量估算为普通高斯信道容量而混合光加密的实际作用是让窃听者等效到比 C_eve 更差的位置所以该公式给出的是偏保守的下界。我通常以“可实现的保密速率不低于业务速率且链路保持稳定”为验收口径而不仅仅看这个理论容量值。3. 用 MATLAB 本地跑通混合光加密的最小仿真3.1 以 CO-OFDM 为载体的混合光加密链路结构选择 CO-OFDM相干光正交频分复用作为载体是因为它天然是复数信号方便在频域乘相位掩膜做仿真时也容易对照星座图观察加密效果。常见做法是构建如下链路发送端生成随机比特按选定的 QAM 阶数做星座映射在频域上给每个子载波乘上混沌生成的相位掩膜IFFT 回到时域加循环前缀和训练序列并入 ASE 噪声与色散损伤模拟光纤信道接收端做同步、信道估计后用密钥重放掩膜并共轭相乘解映射并统计误码率。相位掩膜放在 IFFT 之前是为了让掩膜的扰动作用均匀展开到时域全部采样点上。如果放在 IFFT 之后掩膜会直接改变时域波形包络对峰均比产生不可控影响放在频域则能保持每个子载波相对独立的相位偏移。3.2 最小仿真代码混沌序列生成与相位掩膜加密下面这段代码把混沌密钥和频域掩膜这两个核心动作放在一个可运行的最小脚本里省略 IFFT 和循环前缀以便单独观察加密对误码率的影响。% 混合光加密最小仿真发送端 接收端 clear; rng(1); N 256; % 子载波数 M 16; % 16QAM bitsPerSym 4; % 1. 混沌密钥logistic 映射迭代 N 次 r 3.99; % 混沌常数接近 4 时序列遍历性更好 x0 0.314159; % 密钥初值 idx zeros(N,1); x x0; for k 1:N x r * x * (1 - x); % logistic 迭代方程 idx(k) x; end % 2. 相位掩膜把混沌值映射到 [0, 2*pi) mask exp(1j * 2 * pi * idx); % 3. 发送端16QAM 符号与频域掩膜相乘 data randi([0 1], N*bitsPerSym, 1); sym qammod(data, M, gray, InputType, bit); symEnc sym .* mask; % 频域加密 % 4. ASE 噪声用加性高斯白噪声近似 snr 15; % 光信噪比在电域用 SNR 近似 rx awgn(symEnc, snr, measured); % 5. 解密持有密钥时用掩膜共轭恢复 symDec rx .* conj(mask); dataRx qamdemod(symDec, M, gray, OutputType, bit); ber sum(dataRx ~ data) / numel(data); fprintf(16QAM 加掩膜后的误码率%.4f\n, ber);代码背后的逻辑logistic 映射在 r 接近 4 时对初值异常敏感两个不同初值在有限精度下迭代约 20 次后就会完全分叉因此 x0 和 r 的组合就是密钥利用该序列生成的掩膜乘以符号后不会改变信号功率只会旋转相位。接收端把掩膜的共轭乘回去就能恢复原星座。qammod 中的 gray 参数表示格雷映射相邻星座点只有一位翻转可以减少解调误码awgn 的 measured 选项表示噪声功率按信号实际功率计算。把 x0 改成相差 1e-9 的值解密后的误码率会接近 0.5这是混沌初值敏感性的直观体现。3.3 解密与同步掩膜共轭恢复与帧同步上面的代码在理想同步前提下工作实际链路还需要解决掩膜同步问题。合法接收端需要知道当前帧用的掩膜一种常用做法是每帧在导频子载波上发送 BPSK 编码的帧序号接收端根据帧序号和本地保存的初值重新生成掩膜。这样只需要在建立连接时同步一次密钥初值之后每帧的掩膜都能靠帧号重放。做法上我会在每帧设置 4 个导频子载波承载相同的帧号编码接收端对导频做相关判决后取多数投票。该方案的代价是每帧多出约 1% 的带宽但在 OSNR 不低于 12dB 时帧号检测成功率可以做到 1-1e-4 以上。需要注意解密本身不能简单在时域做共轭乘法正确顺序是先做信道均衡再把掩膜共轭与均衡后的残余相位误差合并处理。若顺序颠倒信道相位噪声会叠加进掩膜恢复过程造成额外误码。% 帧同步验证的最小循环 syncErr zeros(100,1); for k 1:100 maskFrame genMask(x0, frameNo(k)); % 本地按帧号重放 est recoverMask(rxFrame(k,:), pilot); % 从导频恢复掩膜 syncErr(k) sum(abs(est - maskFrame) 0.1) / N; end fprintf(掩膜对齐错误比例%.4f\n, mean(syncErr));逻辑说明genMask 与 recoverMask 分别表示本地重放和导频提取恢复两者的差值大于 0.1 的样本点视为失配。对齐错误比例若超过某个阈值说明信道相位变化速度超过帧周期这时候不需要换加密算法而是把掩膜更新粒度从帧级别降到块级别让每个资源块使用独立掩膜。3.4 最小可运行版本与扩展方向上述代码在 MATLAB R2016b 及以上版本可直接运行。若手边没有 MATLAB把 qammod 和 qamdemod 替换成 Python 的 numpy 逐元素运算也可以核心的乘法与共轭乘法不依赖通信工具箱。要完整跑通包含 IFFT、循环前缀和色散补偿的 CO-OFDM 链路则需要在此基础上再补约八百行收发代码。扩展方向通常有三类把 logistic 换成更高维的 Chen 或超混沌系统增大密钥空间把单个相位掩膜换成“频域时域”双掩膜抵抗针对单掩膜的频谱分析把固定的 QAM 映射改成密钥控制的动态映射表这一步对抗机器学习解调特别有效。4. 混合光加密的关键参数与双用户 BER 安全评估4.1 密钥初值、混沌常数与掩膜颗粒度的设定这几个参数直接影响安全强度也存在相互制约。下表是我在仿真与实验里常用的参考区间参数推荐范围对安全性的作用设置过小的影响混沌初值 x0随机位宽 52 位以上决定密钥空间大小可被暴力搜索混沌常数 r3.903.999决定序列遍历范围序列退化并进入周期掩膜颗粒度子载波级或 4 子载波块级加密单元越细统计越难窃听者能做星座聚类掩膜更新周期每帧或每 10 帧抗重放和样本积累攻击者积累同分布样本掩膜颗粒度指的是一个相位值覆盖的子载波个数子载波级掩膜给每个子载波独立随机相位安全收益最高但对相位噪声更敏感块级掩膜让相邻 4 个子载波共用同一个相位偏移抗激光器线宽能力强一些。实际网络里我会把 x0 与帧号绑定每一帧的掩膜由“初值 帧号”共同驱动生成。好处是连接建立时只同步一次初值后续掩膜从帧头读出帧号就能在本地重放不需要在数据面额外传输掩膜本身。4.2 合法用户与窃听者双 BER 曲线的仿真对比安全评估不能只看合法用户的误码率要看同一信道上合法用户和窃听者的差异。下面代码在一个循环里同时统计双方误码率窃听者假设没有掩膜信息直接对接收符号做硬判决。% 双用户BER合法用户带掩膜窃听者不做掩膜处理 M 16; N 256; data randi([0 1], N*4, 1); % 混沌掩膜生成与第 3 节相同 r 3.99; x0 0.314159; mask zeros(N,1); x x0; for k 1:N x r * x * (1 - x); mask(k) exp(1j * 2 * pi * x); end sym qammod(data, M, gray, InputType, bit); symEnc sym .* mask; snrVec 5:2:21; berLegit zeros(size(snrVec)); berEve zeros(size(snrVec)); for i 1:numel(snrVec) rx awgn(symEnc, snrVec(i), measured); % 合法用户持有掩膜共轭恢复星座 sym1 rx .* conj(mask); bit1 qamdemod(sym1, M, gray, OutputType, bit); berLegit(i) sum(bit1 ~ data) / numel(data); % 窃听者无掩膜直接硬判决 bit2 qamdemod(rx, M, gray, OutputType, bit); berEve(i) sum(bit2 ~ data) / numel(data); end table(snrVec, berLegit, berEve, ... VariableNames, {snr_dB,ber_legit,ber_eve})逻辑说明合法用户的误码率随 SNR 升高快速下降窃听者的误码率则始终维持在高位因为掩膜把星座旋转成近似圆对称分布对没有掩膜的一方来说有效星座点落在决策边界附近。自适应判决阈值在这里也帮不上忙16QAM 加频域掩膜后任何线性均衡只能乘一个固定复数无法还原逐子载波变化的相位。只有正确的掩膜序列才能逐点对齐。从仿真结果看当合法用户误码率降到 1e-4 时窃听者仍停留在 0.3 左右。这个差距是混合光加密安全性的直观证据。建议在报告中把两条 BER 曲线画在一起横轴统一为 OSNR纵轴为误码率。4.3 从仿真读数到安全边际的换算仿真曲线出来后我习惯先换算成两个可操作的量。第一个是合法用户与窃听者的接收灵敏度差在目标误码率 1e-3 时合法用户需要 OSNR 12dB窃听者可能需要 28dB差值 16dB 就是安全余量含义是攻击者必须把窃取光功率提升 16dB 才能达到与合法用户同等的解码效果。第二个是 BER 曲线间隔正常设计下两条曲线至少应相距 6dB。注意若合法用户与窃听者曲线间隔小于 3dB先检查相位掩膜是否只有有限分辨率或者混沌序列是否已经退化到周期状态不要急着增加加密层数掩盖参数失配的根源。5. 混合光加密装进现有光网络前要验证的三个边界5.1 密钥同步与帧头对准的最小验证混合光加密最容易在同步环节出问题。实际链路里帧号恢复错误会导致整个数据块误码因此在改动信道编码之前先用一个小脚本反复验证帧同步的余量。上面的 syncErr 循环可以扩展为不同 OSNR 下的失配统计把 OSNR 从 8dB 以 1dB 步进扫到 20dB记录失配比例开始升高的拐点。拐点如果高于系统工作 OSNR 2dB 以上说明同步设计有足够余量如果拐点与工作点接近需要增加导频密度或缩短掩膜更新周期。5.2 激光器线宽对加密误码率的压制相干光接收机的相位噪声会直接叠加在掩膜相位上。激光器线宽在 100kHz 以下时子载波级掩膜可行使用 DFB 激光器线宽 1~3MHz时相邻子载波的相位跳变会被线宽噪声抹平掩膜共轭恢复后残留随机旋转需要把掩膜颗粒度从子载波级调整为 4 个子载波一块。这个调整会把安全强度降低约 2dB但能让误码率回到工程可用范围。实验阶段可以用高线宽激光器加色散补偿模块观察固定 OSNR 下误码率变化快速判断当前掩膜粒度是否超出相位噪声容限。5.3 用 OSNR 余量反推安全运行区间最后把安全余量和传输余量合并成一张运行区间表分别记录系统无误码的最低 OSNR、FEC 极限 OSNR 和窃听者要达到相同误码率所需的 OSNR。三者之间的关系决定了加密能否在现有光缆上安全部署。举例来说若合法用户无误码 OSNR 是 12dBFEC 极限是 9dB而窃听者需要 28dB那么系统允许链路劣化 3dB同时还有 16dB 的安全深度。如果实测发现合法用户无误码 OSNR 比理论值劣化了 2dB 以上优先排查相位噪声与掩膜同步而不是怀疑混沌算法本身。验证方法和技巧强调一句。本文还有配套的精品资源点击获取
混合光加密：光网络物理层安全的关键技术与仿真实践
简介一份围绕光通信物理层安全展开的专业文档面向通信工程、光网络与信息安全领域的研究者、工程师及高年级学生也可作为相关课程设计与毕业设计的参考。资源从5G时代光纤传输面临的安全挑战入手指出上层加密在保护范围与密钥管理上的短板系统阐述基于高非线性光纤简并四波混频效应与超结构光纤布拉格光栅的混合光加密方案涵盖QPSK信号加解密流程、密钥隐匿传输机制、SSFBG编解码原理及系统结构设计并结合相位矢配条件、弱耦合近似等理论分析说明方案可行性。文档还给出了VPI仿真平台的参数配置、误码率测试曲线以及加密前后信号对比便于读者掌握从理论推导、系统搭建到性能验证的完整链路。资源为1个docx文档大小668KB内容精炼且公式图表齐全已有89人学习可作为课题研究、毕业设计或光通信安全技术预研的参考资料。1. 混合光加密正在接管光链路的物理层安全传统密码学保护光网络时加密和解密都发生在数据链路层光脉冲进入光纤后以“明文”物理波形传输一根分光器复制光信号不需要破解任何上层协议。混合光加密的思路是把混沌密钥生成与相位掩膜处理放到光域在发送端同时改变信号的幅度与相位分布合法接收端凭密钥恢复原始符号窃听端只能看到接近噪声的基带图样。这个方案适合想补齐物理层安全短板的传输工程师也适合对光链路做安全评估的网络规划人员。下面的内容会从理论边界起步逐步推进到最小可复现仿真、参数评估和落地前的验证方法用一个连续链路把原理和指标串起来。2. 混合光加密的物理层实现框架与光域安全边界2.1 为什么在电域完成加密后仍要关心物理层安全在真实光网络里窃听者并不需要破译管理面或数据面的密码系统。最常见的手段是在光缆护套开一个口子引入分光器或在某一段制造微弯窃取一小部分光功率这种截取方式几乎不会引发上层告警。此时上层加密保护的是载荷内容但光层信号的时间结构、波长位置、调制格式都以明文暴露给攻击者。更关键的是骨干网节点之间的光信号是连续波逐跳解密再加密的成本极高因此“光纤上走的是密文光纤本身却是透明的”这个矛盾长期存在。电域加密还有一个结构性问题OFDM 或 QAM 这类高阶调制在光纤中传输时星座图、峰均比、循环前缀这些物理特征是明文暴露的。攻击者可以使用机器学习识别调制类型与帧结构辅助流量分析和定向干扰。物理层安全的目标不是替代 AES 这类算法而是抬高窃听信道的代价让合法信道与窃听信道之间形成确定的容量差。混合光加密的价值正在于它把破解难度从“已知密文且已知算法”的标准密码分析问题转换成了“未知信道模型加未知随机映射”的估计问题对想获取明文的一方来说这两类问题所需的工具完全不在一个数量级。2.2 混合光加密的三个组成面混沌密钥、相位掩膜与符号映射“混合”二字的含义是几种手段在同一个发送端叠加使用。我一般把发送端拆成三个面来设计。密钥面由混沌映射产生伪随机序列作为相位掩膜的种子和动态扰码的钥流掩膜面把相位掩膜乘到频域复信号上使星座点的相位分布整体打散映射面将原始比特通过密钥控制的映射关系从固定 QAM 星座切换到随帧变化的星座子集。三者的结合在光域上的表现就是信号的模拟波形随密钥变化而变化静态星座图不复存在。攻击者在有限截获样本上做的聚类、调制识别和符号估计都会因为映射关系逐帧漂移而失效。下表方便对照三个组成面的作用边界组成面典型实现作用边界不设置时的影响混沌密钥Logistic 映射、Chen 混沌决定密钥空间与帧同步要求掩膜与映射均可被识别相位掩膜随机相位板、Hadamard 矩阵对抗幅度统计与频谱估计星座相位可直接观测符号映射动态 QAM 子集选择对抗聚类与模式识别窃听端可直接硬判决在系统实现顺序上混沌密钥先于其他两面生成它既驱动掩膜也决定映射表。映射面的工程实现一般用一张密钥索引的置换表把输入比特组映射到 M 个星座点的非自然顺序。这样即使攻击者截获到一段没有相位掩膜的导频也无法利用导频建立可信的星座骨架。2.3 物理层安全收益怎么算安全容量与窃听信噪比评估“安全”不能只看加密算法强度工程上更多以信道不对称性来度量。Wyner 安全容量的核心表达式为合法信道容量减窃听信道容量安全容量 C_s max{0, C_legit - C_eve}每个容量都按香农公式 C B * log2(1 SNR) 计算。光链路中合法接收机与窃听者从同一根光纤分支出去物理层的唯一差异是接收光功率OSNR和相位噪声累积。混合光加密真正改善的是窃听者解调后的等效信噪比掩膜和动态映射使其恢复星座时产生一个叠加惩罚而这个惩罚转化为安全收益用来抵消加密给合法用户带来的消光比损失。下面这段代码可以直接估算给定合法与窃听 OSNR 时的安全容量% 安全容量估算B 为信号带宽10 GHz snrLegitDb 18; snrEveDb 8; B 10e9; Clegit B * log2(1 10^(snrLegitDb/10)); Ceve B * log2(1 10^(snrEveDb/10)); fprintf(安全容量约 %.2f Gbit/s\n, max(0, Clegit - Ceve));这段代码的参数说明snrLegitDb 和 snrEveDb 分别是合法与窃听接收机在解调前的信噪比B 是子载波聚合后的有效带宽。注意这里把窃听信道容量估算为普通高斯信道容量而混合光加密的实际作用是让窃听者等效到比 C_eve 更差的位置所以该公式给出的是偏保守的下界。我通常以“可实现的保密速率不低于业务速率且链路保持稳定”为验收口径而不仅仅看这个理论容量值。3. 用 MATLAB 本地跑通混合光加密的最小仿真3.1 以 CO-OFDM 为载体的混合光加密链路结构选择 CO-OFDM相干光正交频分复用作为载体是因为它天然是复数信号方便在频域乘相位掩膜做仿真时也容易对照星座图观察加密效果。常见做法是构建如下链路发送端生成随机比特按选定的 QAM 阶数做星座映射在频域上给每个子载波乘上混沌生成的相位掩膜IFFT 回到时域加循环前缀和训练序列并入 ASE 噪声与色散损伤模拟光纤信道接收端做同步、信道估计后用密钥重放掩膜并共轭相乘解映射并统计误码率。相位掩膜放在 IFFT 之前是为了让掩膜的扰动作用均匀展开到时域全部采样点上。如果放在 IFFT 之后掩膜会直接改变时域波形包络对峰均比产生不可控影响放在频域则能保持每个子载波相对独立的相位偏移。3.2 最小仿真代码混沌序列生成与相位掩膜加密下面这段代码把混沌密钥和频域掩膜这两个核心动作放在一个可运行的最小脚本里省略 IFFT 和循环前缀以便单独观察加密对误码率的影响。% 混合光加密最小仿真发送端 接收端 clear; rng(1); N 256; % 子载波数 M 16; % 16QAM bitsPerSym 4; % 1. 混沌密钥logistic 映射迭代 N 次 r 3.99; % 混沌常数接近 4 时序列遍历性更好 x0 0.314159; % 密钥初值 idx zeros(N,1); x x0; for k 1:N x r * x * (1 - x); % logistic 迭代方程 idx(k) x; end % 2. 相位掩膜把混沌值映射到 [0, 2*pi) mask exp(1j * 2 * pi * idx); % 3. 发送端16QAM 符号与频域掩膜相乘 data randi([0 1], N*bitsPerSym, 1); sym qammod(data, M, gray, InputType, bit); symEnc sym .* mask; % 频域加密 % 4. ASE 噪声用加性高斯白噪声近似 snr 15; % 光信噪比在电域用 SNR 近似 rx awgn(symEnc, snr, measured); % 5. 解密持有密钥时用掩膜共轭恢复 symDec rx .* conj(mask); dataRx qamdemod(symDec, M, gray, OutputType, bit); ber sum(dataRx ~ data) / numel(data); fprintf(16QAM 加掩膜后的误码率%.4f\n, ber);代码背后的逻辑logistic 映射在 r 接近 4 时对初值异常敏感两个不同初值在有限精度下迭代约 20 次后就会完全分叉因此 x0 和 r 的组合就是密钥利用该序列生成的掩膜乘以符号后不会改变信号功率只会旋转相位。接收端把掩膜的共轭乘回去就能恢复原星座。qammod 中的 gray 参数表示格雷映射相邻星座点只有一位翻转可以减少解调误码awgn 的 measured 选项表示噪声功率按信号实际功率计算。把 x0 改成相差 1e-9 的值解密后的误码率会接近 0.5这是混沌初值敏感性的直观体现。3.3 解密与同步掩膜共轭恢复与帧同步上面的代码在理想同步前提下工作实际链路还需要解决掩膜同步问题。合法接收端需要知道当前帧用的掩膜一种常用做法是每帧在导频子载波上发送 BPSK 编码的帧序号接收端根据帧序号和本地保存的初值重新生成掩膜。这样只需要在建立连接时同步一次密钥初值之后每帧的掩膜都能靠帧号重放。做法上我会在每帧设置 4 个导频子载波承载相同的帧号编码接收端对导频做相关判决后取多数投票。该方案的代价是每帧多出约 1% 的带宽但在 OSNR 不低于 12dB 时帧号检测成功率可以做到 1-1e-4 以上。需要注意解密本身不能简单在时域做共轭乘法正确顺序是先做信道均衡再把掩膜共轭与均衡后的残余相位误差合并处理。若顺序颠倒信道相位噪声会叠加进掩膜恢复过程造成额外误码。% 帧同步验证的最小循环 syncErr zeros(100,1); for k 1:100 maskFrame genMask(x0, frameNo(k)); % 本地按帧号重放 est recoverMask(rxFrame(k,:), pilot); % 从导频恢复掩膜 syncErr(k) sum(abs(est - maskFrame) 0.1) / N; end fprintf(掩膜对齐错误比例%.4f\n, mean(syncErr));逻辑说明genMask 与 recoverMask 分别表示本地重放和导频提取恢复两者的差值大于 0.1 的样本点视为失配。对齐错误比例若超过某个阈值说明信道相位变化速度超过帧周期这时候不需要换加密算法而是把掩膜更新粒度从帧级别降到块级别让每个资源块使用独立掩膜。3.4 最小可运行版本与扩展方向上述代码在 MATLAB R2016b 及以上版本可直接运行。若手边没有 MATLAB把 qammod 和 qamdemod 替换成 Python 的 numpy 逐元素运算也可以核心的乘法与共轭乘法不依赖通信工具箱。要完整跑通包含 IFFT、循环前缀和色散补偿的 CO-OFDM 链路则需要在此基础上再补约八百行收发代码。扩展方向通常有三类把 logistic 换成更高维的 Chen 或超混沌系统增大密钥空间把单个相位掩膜换成“频域时域”双掩膜抵抗针对单掩膜的频谱分析把固定的 QAM 映射改成密钥控制的动态映射表这一步对抗机器学习解调特别有效。4. 混合光加密的关键参数与双用户 BER 安全评估4.1 密钥初值、混沌常数与掩膜颗粒度的设定这几个参数直接影响安全强度也存在相互制约。下表是我在仿真与实验里常用的参考区间参数推荐范围对安全性的作用设置过小的影响混沌初值 x0随机位宽 52 位以上决定密钥空间大小可被暴力搜索混沌常数 r3.903.999决定序列遍历范围序列退化并进入周期掩膜颗粒度子载波级或 4 子载波块级加密单元越细统计越难窃听者能做星座聚类掩膜更新周期每帧或每 10 帧抗重放和样本积累攻击者积累同分布样本掩膜颗粒度指的是一个相位值覆盖的子载波个数子载波级掩膜给每个子载波独立随机相位安全收益最高但对相位噪声更敏感块级掩膜让相邻 4 个子载波共用同一个相位偏移抗激光器线宽能力强一些。实际网络里我会把 x0 与帧号绑定每一帧的掩膜由“初值 帧号”共同驱动生成。好处是连接建立时只同步一次初值后续掩膜从帧头读出帧号就能在本地重放不需要在数据面额外传输掩膜本身。4.2 合法用户与窃听者双 BER 曲线的仿真对比安全评估不能只看合法用户的误码率要看同一信道上合法用户和窃听者的差异。下面代码在一个循环里同时统计双方误码率窃听者假设没有掩膜信息直接对接收符号做硬判决。% 双用户BER合法用户带掩膜窃听者不做掩膜处理 M 16; N 256; data randi([0 1], N*4, 1); % 混沌掩膜生成与第 3 节相同 r 3.99; x0 0.314159; mask zeros(N,1); x x0; for k 1:N x r * x * (1 - x); mask(k) exp(1j * 2 * pi * x); end sym qammod(data, M, gray, InputType, bit); symEnc sym .* mask; snrVec 5:2:21; berLegit zeros(size(snrVec)); berEve zeros(size(snrVec)); for i 1:numel(snrVec) rx awgn(symEnc, snrVec(i), measured); % 合法用户持有掩膜共轭恢复星座 sym1 rx .* conj(mask); bit1 qamdemod(sym1, M, gray, OutputType, bit); berLegit(i) sum(bit1 ~ data) / numel(data); % 窃听者无掩膜直接硬判决 bit2 qamdemod(rx, M, gray, OutputType, bit); berEve(i) sum(bit2 ~ data) / numel(data); end table(snrVec, berLegit, berEve, ... VariableNames, {snr_dB,ber_legit,ber_eve})逻辑说明合法用户的误码率随 SNR 升高快速下降窃听者的误码率则始终维持在高位因为掩膜把星座旋转成近似圆对称分布对没有掩膜的一方来说有效星座点落在决策边界附近。自适应判决阈值在这里也帮不上忙16QAM 加频域掩膜后任何线性均衡只能乘一个固定复数无法还原逐子载波变化的相位。只有正确的掩膜序列才能逐点对齐。从仿真结果看当合法用户误码率降到 1e-4 时窃听者仍停留在 0.3 左右。这个差距是混合光加密安全性的直观证据。建议在报告中把两条 BER 曲线画在一起横轴统一为 OSNR纵轴为误码率。4.3 从仿真读数到安全边际的换算仿真曲线出来后我习惯先换算成两个可操作的量。第一个是合法用户与窃听者的接收灵敏度差在目标误码率 1e-3 时合法用户需要 OSNR 12dB窃听者可能需要 28dB差值 16dB 就是安全余量含义是攻击者必须把窃取光功率提升 16dB 才能达到与合法用户同等的解码效果。第二个是 BER 曲线间隔正常设计下两条曲线至少应相距 6dB。注意若合法用户与窃听者曲线间隔小于 3dB先检查相位掩膜是否只有有限分辨率或者混沌序列是否已经退化到周期状态不要急着增加加密层数掩盖参数失配的根源。5. 混合光加密装进现有光网络前要验证的三个边界5.1 密钥同步与帧头对准的最小验证混合光加密最容易在同步环节出问题。实际链路里帧号恢复错误会导致整个数据块误码因此在改动信道编码之前先用一个小脚本反复验证帧同步的余量。上面的 syncErr 循环可以扩展为不同 OSNR 下的失配统计把 OSNR 从 8dB 以 1dB 步进扫到 20dB记录失配比例开始升高的拐点。拐点如果高于系统工作 OSNR 2dB 以上说明同步设计有足够余量如果拐点与工作点接近需要增加导频密度或缩短掩膜更新周期。5.2 激光器线宽对加密误码率的压制相干光接收机的相位噪声会直接叠加在掩膜相位上。激光器线宽在 100kHz 以下时子载波级掩膜可行使用 DFB 激光器线宽 1~3MHz时相邻子载波的相位跳变会被线宽噪声抹平掩膜共轭恢复后残留随机旋转需要把掩膜颗粒度从子载波级调整为 4 个子载波一块。这个调整会把安全强度降低约 2dB但能让误码率回到工程可用范围。实验阶段可以用高线宽激光器加色散补偿模块观察固定 OSNR 下误码率变化快速判断当前掩膜粒度是否超出相位噪声容限。5.3 用 OSNR 余量反推安全运行区间最后把安全余量和传输余量合并成一张运行区间表分别记录系统无误码的最低 OSNR、FEC 极限 OSNR 和窃听者要达到相同误码率所需的 OSNR。三者之间的关系决定了加密能否在现有光缆上安全部署。举例来说若合法用户无误码 OSNR 是 12dBFEC 极限是 9dB而窃听者需要 28dB那么系统允许链路劣化 3dB同时还有 16dB 的安全深度。如果实测发现合法用户无误码 OSNR 比理论值劣化了 2dB 以上优先排查相位噪声与掩膜同步而不是怀疑混沌算法本身。验证方法和技巧强调一句。本文还有配套的精品资源点击获取
