搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

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

基于激光器阵列后处理的混沌熵源获取高品质随机数

吴佳辰 宋峥 谢溢锋 周心雨 周沛 穆鹏华 李念强

引用本文:
Citation:

基于激光器阵列后处理的混沌熵源获取高品质随机数

吴佳辰, 宋峥, 谢溢锋, 周心雨, 周沛, 穆鹏华, 李念强

High-quality random number sequences extracted from chaos post-processed by phased-array semiconductor laser

Wu Jia-Chen, Song Zheng, Xie Yi-Feng, Zhou Xin-Yu, Zhou Pei, Mu Peng-Hua, Li Nian-Qiang
PDF
HTML
导出引用
  • 本文提出采用可集成的激光器阵列后处理光反馈半导体激光器的输出, 进而获得无时延特征的优质混沌熵源, 进一步获取高速高品质随机数序列. 方案中采用常规的8位模数转换采样量化和多位最低有效位异或提取处理, 采用国际公认的随机数行业测试标准(NIST SP 800-22)来检验产生的序列. 结果表明, 通过激光器阵列后处理的混沌熵源所获取的随机数序列具有均匀的分布特性, 散点图无明显图案, 可以成功通过NIST SP 800-22的全部测试. 另外, 基于激光器阵列的可扩展性, 本方案可以拓展为可实现同时产生多路并行的高速高品质随机数发生器.
    With the rapid development of the computer technology and communication technology, as well as the popularization of the Internet, information security has received much attention of all fields. To ensure the information security, a large number of random numbers must be generated. It is well accepted that random numbers can be divided into physical random numbers and pseudo random numbers. The pseudo random numbers are mainly generated based on algorithms, which can be reproduced once the seed is decoded. The physical random numbers are extracted from physical entropies. While the bandwidth of the traditional physical entropy source is quite small, the bit rate of generated physical random numbers is limited. In the literature, a lot of methods have been proposed to produce high-quality and high-speed random number sequences with the chaotic entropy source, which exhibits wide bandwidth, large amplitude and random fluctuations. Usually, a semiconductor laser with optical feedback, i.e, an external-cavity semiconductor laser (ECSL), is chosen as a chaotic entropy source to generate a chaotic signal output. However, the chaotic signal output has a high time delay characteristic, which is not conducive to the production of high-quality random numbers. In this paper, to produce high-quality chaos with time-delay signature (TDS) being well suppressed, we propose to employ an integration-oriented phased-array semiconductor laser to post-process the original chaos generated by an ECSL. It is shown that the proposed laser array is effective in TDS suppression, which improves the quality of optical chaos. After certain necessary post-processing, high-speed and high-quality random number sequences can be achieved. In this paper, we employ the conventional post-processing techniques, which include an 8-bit analog-to-digital converter (ADC) for sampling and quantization, and m-bits least significant bit (m-LSB) and exclusive OR (XOR) for removing bias. The simulation results show that the random number sequences obtained from the chaotic entropy source comprised of an ECSL and phased-array semiconductor lasers have uniform distribution characteristic and their scatter diagram contains no obvious pattern. Meanwhile, the obtained random number sequences can pass all tests of the standard randomness benchmark, NIST SP 800-22. Additionally, based on the extensibility of phased-array semiconductor lasers, random number generators that can generate parallel random numbers are achievable.
      通信作者: 周沛, peizhou@suda.edu.cn ; 李念强, nli@suda.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 62004135, 62001317)、江苏省高等学校自然科学研究重大项目(批准号: 20KJA416001)和苏州大学科研启动经费(批准号: Q415900119)资助的课题
      Corresponding author: Zhou Pei, peizhou@suda.edu.cn ; Li Nian-Qiang, nli@suda.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 62004135, 62001317), the Natural Science Research Project of Jiangsu Higher Education Institutions, China (Grant No. 20KJA416001), and the Startup Funding of Soochow University, China (Grant No. Q415900119).
    [1]

    Shannon C E 1949 Bell Syst. Tech. J. 28 656Google Scholar

    [2]

    Durt T, Beimonte C, Lamoureux L P, Panajotov K, van den Berghe F, Thienpont H 2013 Phys. Rev. A. 87 022339Google Scholar

    [3]

    Williams C R S, Salevan J C, Li X W, Roy R, Murphy T E 2010 Opt. Express 18 23584Google Scholar

    [4]

    Guo H, Tang W Z, Liu Y, Wei W 2010 Phys. Rev. E. 81 051137Google Scholar

    [5]

    Ma X F, Xu F H, Xu H, Tan X Q, Qi B, Lo H K 2013 Phys. Rev. A. 87 062327Google Scholar

    [6]

    David P. Rosin, Damien Rontani, Daniel J. Gauthier 2013 Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 87 040902Google Scholar

    [7]

    Li N Q, Kim B, Chizhevsky V N, Locquet A, Bloch M, Citrin D S, Pan W 2014 Opt. Express 22 6634Google Scholar

    [8]

    郭弘, 刘钰, 党安红, 韦韦 2009 科学通报 54 3651Google Scholar

    Guo H, Liu Y, Dang A H, Wei W 2009 Chin. Sci. Bull. 54 3651Google Scholar

    [9]

    Ren M, Wu E, Liang Y, Jian Y, Wu G, Zeng H 2011 Phys. Rev. A 83 023820Google Scholar

    [10]

    周庆, 胡月, 廖晓峰 2008 物理学报 57 5413Google Scholar

    Zhou Q, Hu Y, Liao X F 2008 Acta Phys. Sin. 57 5413Google Scholar

    [11]

    李念强 2016 博士学位论文 (成都: 西南交通大学)

    Li N Q 2016 Ph. D. Dissertation (Chengdu: Southwest Jiaotong University) (in Chinese)

    [12]

    Uchida A, Amano K, Inoue M, Hirano K, Naito S, Someya H, Oowada I, Kurashige T, Shiki M, Yoshimori S, Yoshimura K, Davis P 2008 Nat. Photon. 2 728Google Scholar

    [13]

    Reidier I, Aviad Y, Rosenblush M, Kanter I 2009 Phys. Rev. Lett. 103 024102Google Scholar

    [14]

    Kanter I, Aviad Y, Reidler I, Cohen E, Rosenbluh M 2010 Nat. Photon. 4 58Google Scholar

    [15]

    Harayama T, Sunada S, Yoshimura K, Davis P, Tsuzuki K, Uchida A 2011 Phys. Rev. A. 83 031803Google Scholar

    [16]

    Argyris A, Deligiannidis S, Pikasis E, Bogris A, Syvridis D 2010 Opt. Express 18 18763Google Scholar

    [17]

    Zhang J Z, Wang Y C, Liu M, Xue L G, Li P, Wang A B, Zhang M J 2012 Opt. Express 20 7496Google Scholar

    [18]

    Li P, Wang Y C, Zhang J Z 2010 Opt. Express 18 20360Google Scholar

    [19]

    Wu J G, Tang X, Wu Z M, Xia G Q, Feng G Y 2012 Laser Phys. 22 1476Google Scholar

    [20]

    Li X Z, Chan S C 2013 IEEE J. Quantum Electron. 49 829Google Scholar

    [21]

    Wang A B, Li P, Zhang J G, Zhang J Z, Li L, Wang Y C 2013 Opt. Express 21 20452Google Scholar

    [22]

    Li P, Zhang J G, Sang L X, Liu X L, Guo Y Q, Guo X M, Wang A B, Shore K A, Wang Y C 2017 Opt. Lett. 42 2699Google Scholar

    [23]

    Li P, Sun Y Y, Liu X L, Yi X G, Zhang J G, Guo X M, Guo Y Q, Wang Y C 2016 Opt. Lett. 41 3347Google Scholar

    [24]

    韩韬, 刘香莲, 李璞, 郭晓敏, 郭龑强, 王云才 2017 物理学报 66 124203Google Scholar

    Han T, Liu X L, Li P, Guo X M, Guo Y Q, Wang Y C 2017 Acta Phys. Sin. 66 124203Google Scholar

    [25]

    赵东亮, 李璞, 刘香莲, 郭晓敏, 郭龑强, 张建国, 王云才 2017 物理学报 66 050501Google Scholar

    Zhao D L, Li P, Liu X L, Guo X M, Guo Y Q, Zhang J G, Wang Y C 2017 Acta Phys. Sin. 66 050501Google Scholar

    [26]

    Tang X, Wu Z M, Wu J G, Deng T, Zhong Z Q, Chen J J, Xia G Q 2014 Laser Phys. Lett. 12 015003Google Scholar

    [27]

    Ran C, Tang X, Wu Z M, Xia G Q 2018 Laser Phys. 28 126202Google Scholar

    [28]

    姚晓洁, 唐曦, 吴正茂, 夏光琼 2018 物理学报 67 024204Google Scholar

    Yao X J, Tang X, Wu Z M, Xia G Q 2018 Acta Phys. Sin. 67 024204Google Scholar

    [29]

    Li N Q, Pan W, Xiang S Y, Zhao Q C, Zhang L Y 2014 IEEE Photon. Technol. Lett. 26 1886Google Scholar

    [30]

    Mu P H, Pan W, Xiang S Y, Li N Q, Liu X K, Zou X H 2015 Mod. Phys. Lett. B 29 1550142

    [31]

    Wang Y, Xiang S Y, Wang B, Cao X Y, Wen A J, Hao Y 2019 Opt. Express 27 8446Google Scholar

    [32]

    Xiang S Y, Wang B, Wang Y, Han Y N, Wen A J, Hao Y 2019 J. Light. Technol. 37 3987Google Scholar

    [33]

    Xue C P, Jiang N, Qiu K, Lv Y X 2015 Opt. Express 23 14510Google Scholar

    [34]

    Zhao A K, Jiang N, Wang Y J, Liu S Q, Li B C, Qiu K 2019 Opt. Lett. 44 5957Google Scholar

    [35]

    Li N Q, Pan W, Locquet A, Citrin D S 2015 Opt. Lett. 40 4416Google Scholar

    [36]

    Li S S, Li X Z, Chan S C 2018 Opt. Lett. 43 4751Google Scholar

    [37]

    Xiang S Y, Wen A J, Pan W, Lin L, Zhang H X, Zhang H, Guo X X, Li J F 2016 J. Light. Technol. 34 4221Google Scholar

    [38]

    Jiang N, Wang C, Xue C P, Li G L, Lin S Q, Qiu K 2017 Opt. Express 25 14359Google Scholar

    [39]

    Jiang X X, Liu D M, Cheng M F, Deng L, Fu S N, Zhang M M, Tang M, Shum P 2016 Opt. Lett. 41 1157Google Scholar

    [40]

    Ma Y T, Xiang S Y, Guo X X, Song Z W, Wen A J, Hao Y 2020 Opt. Express 28 1665Google Scholar

    [41]

    Zhou P, Fang Q, Li N Q 2020 Opt. Lett. 45 399Google Scholar

    [42]

    Rukhin A, Soto J, Nechvatal J, Smid M, Barker E, Leigh S, Levenson M, Vangel M, Banks D, Heckert A, Dary J, Vo S 2001 http://csrc.nist.gov/groups/ST/toolkit/rng/documentation_software.html [2020-11-21]

    [43]

    Adams M J, Li N Q, Cemlyn B R, Susanto H, Henning I D 2017 Phys. Rev. A 95 053869Google Scholar

    [44]

    Fang Q, Zhou P, Mu P H, Li N Q 2021 IEEE J. Quantum Electron. 57 1200109.

  • 图 1  基于激光器阵列后处理的混沌熵源获取高品质随机数的示意图($ \lambda /4 $$ 1/4 $波片, PD1、PD2为光电转换器, ADC为模数转换器, LSB为最低有效位, XOR为异或处理)

    Fig. 1.  Schematic diagram of high quality random number generation based on the chaotic entropy source generated by ECSL and post-processed by phased-array semiconductor lasers (λ/4, 1/4 wave plate; PD1 and PD2, photo detector; ADC, analog-to-digital converter; LSB, least significant bit; XOR, exclusive OR).

    图 2  激光器输出混沌信号的时间序列(左列), 自相关函数谱(中列), 功率谱(右列) (a) 光反馈半导体激光器; (b) 注入激光器; (c) 注入激光器阵列

    Fig. 2.  Time series (left column), autocorrelation function (middle column), and power spectra (right column) of the chaotic signal output by laser: (a) ECSL; (b) injection to a single laser A; (c) injection to phased-array lasers.

    图 3  经过激光器阵列后处理混沌熵源的ACF时延处峰值随着注入参数和激光器分离比d/a的演化情况 (a) d/a = 0.2; (b) d/a = 0.4; (c) d/a = 0.6; (d) d/a = 1.0

    Fig. 3.  The evaluation of the ACF peak value located around the feedback delay of the chaotic entropy source that is processed by the phased-array in the plane of injection parameters for several values of laser separation: (a) d/a = 0.2, (b) d/a = 0.4, (c) d/a = 0.6, (d) d/a = 1.0.

    图 4  激光器B输出的混沌信号量化后的统计直方图 (a) 8位ADC输出; (b) 3-LSB输出; (c) XOR输出

    Fig. 4.  Statistical histogram of the quantized chaotic signal of the laser B: (a) The output of 8 bit ADC; (b) the output of 3-LSB; (c) the output of XOR.

    图 5  散点图

    Fig. 5.  Scatter diagram.

    图 6  激光器输出的时间序列与自相关函数 (a) A激光器输出的时间序列; (b) A激光器输出的自相关函数; (c) B激光器输出的时间序列; (d) B激光器输出的自相关函数

    Fig. 6.  Time series and autocorrelation function of the lasers: (a) Time series of laser A; (b) autocorrelation function of laser A; (c) time series of laser B; (d) autocorrelation function of laser. B.

    表 1  NIST统计测试结果

    Table 1.  Result of NIST statistical tests.

    测试名称P-value概率结果
    频数0.5381820.992通过
    块内频数0.2392660.982通过
    累加0.7558190.994通过
    游程0.1404530.988通过
    块内最长游程0.9658600.988通过
    矩阵秩0.2812320.990通过
    离散傅里叶变换0.2066290.982通过
    非重叠模块匹配0.0208310.982通过
    重叠模块匹配0.6993130.984通过
    通用统计0.5101530.994通过
    近似熵0.6993130.994通过
    随机游动0.4436650.986通过
    随机游动变量0.2901580.983通过
    连续性0.0965780.984通过
    线性复杂度0.3408580.986通过
    下载: 导出CSV
  • [1]

    Shannon C E 1949 Bell Syst. Tech. J. 28 656Google Scholar

    [2]

    Durt T, Beimonte C, Lamoureux L P, Panajotov K, van den Berghe F, Thienpont H 2013 Phys. Rev. A. 87 022339Google Scholar

    [3]

    Williams C R S, Salevan J C, Li X W, Roy R, Murphy T E 2010 Opt. Express 18 23584Google Scholar

    [4]

    Guo H, Tang W Z, Liu Y, Wei W 2010 Phys. Rev. E. 81 051137Google Scholar

    [5]

    Ma X F, Xu F H, Xu H, Tan X Q, Qi B, Lo H K 2013 Phys. Rev. A. 87 062327Google Scholar

    [6]

    David P. Rosin, Damien Rontani, Daniel J. Gauthier 2013 Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 87 040902Google Scholar

    [7]

    Li N Q, Kim B, Chizhevsky V N, Locquet A, Bloch M, Citrin D S, Pan W 2014 Opt. Express 22 6634Google Scholar

    [8]

    郭弘, 刘钰, 党安红, 韦韦 2009 科学通报 54 3651Google Scholar

    Guo H, Liu Y, Dang A H, Wei W 2009 Chin. Sci. Bull. 54 3651Google Scholar

    [9]

    Ren M, Wu E, Liang Y, Jian Y, Wu G, Zeng H 2011 Phys. Rev. A 83 023820Google Scholar

    [10]

    周庆, 胡月, 廖晓峰 2008 物理学报 57 5413Google Scholar

    Zhou Q, Hu Y, Liao X F 2008 Acta Phys. Sin. 57 5413Google Scholar

    [11]

    李念强 2016 博士学位论文 (成都: 西南交通大学)

    Li N Q 2016 Ph. D. Dissertation (Chengdu: Southwest Jiaotong University) (in Chinese)

    [12]

    Uchida A, Amano K, Inoue M, Hirano K, Naito S, Someya H, Oowada I, Kurashige T, Shiki M, Yoshimori S, Yoshimura K, Davis P 2008 Nat. Photon. 2 728Google Scholar

    [13]

    Reidier I, Aviad Y, Rosenblush M, Kanter I 2009 Phys. Rev. Lett. 103 024102Google Scholar

    [14]

    Kanter I, Aviad Y, Reidler I, Cohen E, Rosenbluh M 2010 Nat. Photon. 4 58Google Scholar

    [15]

    Harayama T, Sunada S, Yoshimura K, Davis P, Tsuzuki K, Uchida A 2011 Phys. Rev. A. 83 031803Google Scholar

    [16]

    Argyris A, Deligiannidis S, Pikasis E, Bogris A, Syvridis D 2010 Opt. Express 18 18763Google Scholar

    [17]

    Zhang J Z, Wang Y C, Liu M, Xue L G, Li P, Wang A B, Zhang M J 2012 Opt. Express 20 7496Google Scholar

    [18]

    Li P, Wang Y C, Zhang J Z 2010 Opt. Express 18 20360Google Scholar

    [19]

    Wu J G, Tang X, Wu Z M, Xia G Q, Feng G Y 2012 Laser Phys. 22 1476Google Scholar

    [20]

    Li X Z, Chan S C 2013 IEEE J. Quantum Electron. 49 829Google Scholar

    [21]

    Wang A B, Li P, Zhang J G, Zhang J Z, Li L, Wang Y C 2013 Opt. Express 21 20452Google Scholar

    [22]

    Li P, Zhang J G, Sang L X, Liu X L, Guo Y Q, Guo X M, Wang A B, Shore K A, Wang Y C 2017 Opt. Lett. 42 2699Google Scholar

    [23]

    Li P, Sun Y Y, Liu X L, Yi X G, Zhang J G, Guo X M, Guo Y Q, Wang Y C 2016 Opt. Lett. 41 3347Google Scholar

    [24]

    韩韬, 刘香莲, 李璞, 郭晓敏, 郭龑强, 王云才 2017 物理学报 66 124203Google Scholar

    Han T, Liu X L, Li P, Guo X M, Guo Y Q, Wang Y C 2017 Acta Phys. Sin. 66 124203Google Scholar

    [25]

    赵东亮, 李璞, 刘香莲, 郭晓敏, 郭龑强, 张建国, 王云才 2017 物理学报 66 050501Google Scholar

    Zhao D L, Li P, Liu X L, Guo X M, Guo Y Q, Zhang J G, Wang Y C 2017 Acta Phys. Sin. 66 050501Google Scholar

    [26]

    Tang X, Wu Z M, Wu J G, Deng T, Zhong Z Q, Chen J J, Xia G Q 2014 Laser Phys. Lett. 12 015003Google Scholar

    [27]

    Ran C, Tang X, Wu Z M, Xia G Q 2018 Laser Phys. 28 126202Google Scholar

    [28]

    姚晓洁, 唐曦, 吴正茂, 夏光琼 2018 物理学报 67 024204Google Scholar

    Yao X J, Tang X, Wu Z M, Xia G Q 2018 Acta Phys. Sin. 67 024204Google Scholar

    [29]

    Li N Q, Pan W, Xiang S Y, Zhao Q C, Zhang L Y 2014 IEEE Photon. Technol. Lett. 26 1886Google Scholar

    [30]

    Mu P H, Pan W, Xiang S Y, Li N Q, Liu X K, Zou X H 2015 Mod. Phys. Lett. B 29 1550142

    [31]

    Wang Y, Xiang S Y, Wang B, Cao X Y, Wen A J, Hao Y 2019 Opt. Express 27 8446Google Scholar

    [32]

    Xiang S Y, Wang B, Wang Y, Han Y N, Wen A J, Hao Y 2019 J. Light. Technol. 37 3987Google Scholar

    [33]

    Xue C P, Jiang N, Qiu K, Lv Y X 2015 Opt. Express 23 14510Google Scholar

    [34]

    Zhao A K, Jiang N, Wang Y J, Liu S Q, Li B C, Qiu K 2019 Opt. Lett. 44 5957Google Scholar

    [35]

    Li N Q, Pan W, Locquet A, Citrin D S 2015 Opt. Lett. 40 4416Google Scholar

    [36]

    Li S S, Li X Z, Chan S C 2018 Opt. Lett. 43 4751Google Scholar

    [37]

    Xiang S Y, Wen A J, Pan W, Lin L, Zhang H X, Zhang H, Guo X X, Li J F 2016 J. Light. Technol. 34 4221Google Scholar

    [38]

    Jiang N, Wang C, Xue C P, Li G L, Lin S Q, Qiu K 2017 Opt. Express 25 14359Google Scholar

    [39]

    Jiang X X, Liu D M, Cheng M F, Deng L, Fu S N, Zhang M M, Tang M, Shum P 2016 Opt. Lett. 41 1157Google Scholar

    [40]

    Ma Y T, Xiang S Y, Guo X X, Song Z W, Wen A J, Hao Y 2020 Opt. Express 28 1665Google Scholar

    [41]

    Zhou P, Fang Q, Li N Q 2020 Opt. Lett. 45 399Google Scholar

    [42]

    Rukhin A, Soto J, Nechvatal J, Smid M, Barker E, Leigh S, Levenson M, Vangel M, Banks D, Heckert A, Dary J, Vo S 2001 http://csrc.nist.gov/groups/ST/toolkit/rng/documentation_software.html [2020-11-21]

    [43]

    Adams M J, Li N Q, Cemlyn B R, Susanto H, Henning I D 2017 Phys. Rev. A 95 053869Google Scholar

    [44]

    Fang Q, Zhou P, Mu P H, Li N Q 2021 IEEE J. Quantum Electron. 57 1200109.

  • [1] 张依宁, 冯玉玲, 王晓茜, 赵振明, 高超, 姚治海. 半导体激光器混沌输出的延时特征和带宽. 物理学报, 2020, 69(9): 090501. doi: 10.7498/aps.69.20191881
    [2] 赵东亮, 李璞, 刘香莲, 郭晓敏, 郭龑强, 张建国, 王云才. 利用混沌激光脉冲在线实时产生7 Gbit/s物理随机数. 物理学报, 2017, 66(5): 050501. doi: 10.7498/aps.66.050501
    [3] 孙媛媛, 李璞, 郭龑强, 郭晓敏, 刘香莲, 张建国, 桑鲁骁, 王云才. 基于混沌激光的无后处理多位物理随机数高速产生技术研究. 物理学报, 2017, 66(3): 030503. doi: 10.7498/aps.66.030503
    [4] 韩韬, 刘香莲, 李璞, 郭晓敏, 郭龑强, 王云才. 线宽增强因子对光反馈半导体激光器混沌信号生成随机数性能的影响. 物理学报, 2017, 66(12): 124203. doi: 10.7498/aps.66.124203
    [5] 王龙生, 赵彤, 王大铭, 吴旦昱, 周磊, 武锦, 刘新宇, 王安帮. 利用混沌激光多位量化实时产生14 Gb/s的物理随机数. 物理学报, 2017, 66(23): 234205. doi: 10.7498/aps.66.234205
    [6] 杨海波, 吴正茂, 唐曦, 吴加贵, 夏光琼. 反馈强度对外腔反馈半导体激光器混沌熵源生成的随机数序列性能的影响. 物理学报, 2015, 64(8): 084204. doi: 10.7498/aps.64.084204
    [7] 李璞, 江镭, 孙媛媛, 张建国, 王云才. 面向全光物理随机数发生器的混沌实时光采样研究. 物理学报, 2015, 64(23): 230502. doi: 10.7498/aps.64.230502
    [8] 刘莹莹, 潘炜, 江宁, 项水英, 林煜东. 链式互耦合半导体激光器的实时混沌同步. 物理学报, 2013, 62(2): 024208. doi: 10.7498/aps.62.024208
    [9] 黄毅泽, 李毅, 王海方, 俞晓静, 张虎, 张伟, 朱慧群, 孙若曦, 周晟, 张宇明. 双光纤光栅外腔半导体激光器相干失效研究. 物理学报, 2012, 61(1): 014201. doi: 10.7498/aps.61.014201
    [10] 萧宝瑾, 侯佳音, 张建忠, 薛路刚, 王云才. 混沌半导体激光器的弛豫振荡频率对随机序列速率的影响. 物理学报, 2012, 61(15): 150502. doi: 10.7498/aps.61.150502
    [11] 唐曦, 吴加贵, 夏光琼, 吴正茂. 基于互注入半导体激光器的混沌输出产生17.5 Gbit/s随机码. 物理学报, 2011, 60(11): 110509. doi: 10.7498/aps.60.110509
    [12] 张建忠, 王安帮, 张明江, 李晓春, 王云才. 反馈相位随机调制消除混沌半导体激光器的外腔长信息. 物理学报, 2011, 60(9): 094207. doi: 10.7498/aps.60.094207
    [13] 孟丽娜, 张明江, 郑建宇, 张朝霞, 王云才. 外部光注入混沌激光器产生超宽带微波信号的研究. 物理学报, 2011, 60(12): 124212. doi: 10.7498/aps.60.124212
    [14] 陈莎莎, 张建忠, 杨玲珍, 梁君生, 王云才. 基于混沌激光产生1 Gbit/s的随机数. 物理学报, 2011, 60(1): 010501. doi: 10.7498/aps.60.010501
    [15] 操良平, 夏光琼, 邓涛, 林晓东, 吴正茂. 基于非相干光反馈半导体激光器的双向混沌通信研究. 物理学报, 2010, 59(8): 5541-5546. doi: 10.7498/aps.59.5541
    [16] 张继兵, 张建忠, 杨毅彪, 梁君生, 王云才. 外腔半导体激光器随机数熵源的腔长分析. 物理学报, 2010, 59(11): 7679-7685. doi: 10.7498/aps.59.7679
    [17] 牛生晓, 王云才, 贺虎成, 张明江. 光注入半导体激光器产生可调谐高频微波. 物理学报, 2009, 58(10): 7241-7245. doi: 10.7498/aps.58.7241
    [18] 刘四平, 张玉驰, 张鹏飞, 李刚, 王军民, 张天才. 减反膜外腔半导体激光器特性的研究. 物理学报, 2009, 58(1): 285-289. doi: 10.7498/aps.58.285.1
    [19] 赵严峰. 双反馈半导体激光器的混沌特性研究. 物理学报, 2009, 58(9): 6058-6062. doi: 10.7498/aps.58.6058
    [20] 王云才. 增益开关半导体激光器在外光注入下脉冲抖动的实验研究. 物理学报, 2003, 52(9): 2190-2193. doi: 10.7498/aps.52.2190
计量
  • 文章访问数:  1706
  • PDF下载量:  71
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-02
  • 修回日期:  2020-12-19
  • 上网日期:  2021-05-13
  • 刊出日期:  2021-05-20

基于激光器阵列后处理的混沌熵源获取高品质随机数

  • 1. 光电科学与工程学院, 苏州纳米科技协同创新中心, 苏州大学, 苏州 215006
  • 2. 江苏省先进光学制造技术重点实验室, 教育部/江苏省现代光学技术重点实验室, 苏州大学, 苏州 215006
  • 3. 光电信息科学技术学院, 烟台大学, 烟台 264005
  • 通信作者: 周沛, peizhou@suda.edu.cn ; 李念强, nli@suda.edu.cn
    基金项目: 国家自然科学基金(批准号: 62004135, 62001317)、江苏省高等学校自然科学研究重大项目(批准号: 20KJA416001)和苏州大学科研启动经费(批准号: Q415900119)资助的课题

摘要: 本文提出采用可集成的激光器阵列后处理光反馈半导体激光器的输出, 进而获得无时延特征的优质混沌熵源, 进一步获取高速高品质随机数序列. 方案中采用常规的8位模数转换采样量化和多位最低有效位异或提取处理, 采用国际公认的随机数行业测试标准(NIST SP 800-22)来检验产生的序列. 结果表明, 通过激光器阵列后处理的混沌熵源所获取的随机数序列具有均匀的分布特性, 散点图无明显图案, 可以成功通过NIST SP 800-22的全部测试. 另外, 基于激光器阵列的可扩展性, 本方案可以拓展为可实现同时产生多路并行的高速高品质随机数发生器.

English Abstract

    • 随着计算机技术、通信技术的迅速发展, 特别是互联网的普及导致了信息量的爆炸式增长, 信息安全受到了各界的广泛关注. 在信息安全领域, 随机数有着至关重要的地位, 密钥管理、密码学协议、数字签名及身份验证等众多安全技术都需要用到随机数. 而且, 在目前规则下, 绝对安全的保密通讯需要满足Shannon[1]提出的“一次一密”理论, 这就要求大量、高速、安全随机数能实时、快速地产生. 近年来, 随机数的相关研究备受国内外学者的关注, 因此各类随机数发生器被相继提出和验证[2-6].

      随机数可分为真随机数(或物理随机数)和伪随机数[7]. 伪随机数主要是基于算法产生的, 但是因为其固有的周期性使其长度有限, 而且只要获取随机源种子就可复制此类随机数, 不足以保证通信或信息交换、传输的绝对安全. 随着计算机运算能力的提升, 伪随机数用于加密通讯时被破解的可能性大幅度增加, 难以确保通信系统的安全性. 与之不同的是, 真随机数是由物理熵源产生的, 具有高度不可预测性, 使得通信系统的安全性更高. 真随机数的物理熵源主要有电阻热噪声、电子振荡器的频率抖动、电路混沌和激光器相位噪声等[7]. 此外, 利用量子力学基本量的完全随机性及采集生物的无规律行为也可以作为熵源, 通过后处理来提取真随机数[8-10]. 但此类物理熵源的带宽很小, 获取的随机数速率不高, 无法满足当前高速、大容量通信或高速计算模拟的需求. 因此, 寻找新的物理熵源, 通过后处理以实现高速高品质物理随机数发生器的研究成为了当前的研究热点.

      幸运的是, 国内外学者研究发现普通商用半导体激光器在引入一个或多个附加自由度, 如光反馈、光注入或者光电反馈, 可以实现丰富的动力学行为(单周期、多周期、类周期及混沌), 进一步通过优化参数配置即可获取大带宽、高复杂度的混沌熵源[11]. 2008年, 日本Uchida教授的课题组[12]首次通过后处理两路混沌激光器信号, 实现了速率可达1.7 Gb/s的高速随机数生成. 2009年, 以色列著名学者Reidler领导的团队[13]采用8位数模转化器(analog-to-digital converter, ADC)对混沌熵源进行采样量化, 获取了速率为12.5 Gb/s的随机数, 紧接着又采用多级差分处理技术实现了速率可达300 Gb/s的随机数[14]. 上述工作证明了混沌熵源产生高速随机数的可行性, 也把物理随机数的速率提高了多个量级, 掀起了国内外研究采用混沌激光熵源产生高速随机数的热潮[7,15-20]. 特别地, 在国内高校中, 太原理工大学提出了多种产生实时、高速随机数的研究方案, 也提出了基于激光混沌的全光随机数的概念, 最终还实现了随机数发生器样机[18,21-25]; 西南大学在随机数相关研究方面也走在了国际前列, 实现了多种并行随机数发生器[19,26-28]; 西南交通大学提出了采用信息理论区分混沌熵源产生的随机数类型, 并首次实现了基于混沌激光熵源的速率达到2.2 Tb/s量级的随机数[7,29,30]; 成都电子科技大学和西安电子科技大学课题组也在随机数发生器和安全密钥分发方面做了大量优秀的工作[31-34]. 值得注意的是, 在上述混沌熵源的多种产生方案中, 由于光反馈半导体激光器具有系统结构简单、成本低以及动力学丰富的特点, 成为了学者们关注的焦点, 其产生的混沌激光具有大带宽、大幅度和类随机起伏等优点, 因此常用于高速随机数发生器和保密通信领域[11]. 可是光反馈半导体激光器输出的混沌信号具有较高的时间延迟特征(time-delay signature, TDS), 它会阻止生成的序列通过随机数测试标准[7]. 为此, 国内外学者提出了许多可行的方案来削弱或消除这些TDS, 如光注入后处理、光纤传输后处理、互注入或复杂反馈结构等[35-40]. 通过研究发现, 采用上述方案后, 混沌熵源的品质得以提高, 进而进行简单后处理即可产生高速、安全的随机数.

      基于此, 本文采用集成化的激光器阵列作为后处理单元来优化外光反馈激光器产生的原始混沌信号, 消除了混沌信号中的弱周期性[41]. 该方法具有成本低、集成度高及可扩展性高的优点. 需要说明的是, 仿真模型产生的混沌信号与实验条件下的不同, 实验获得的混沌信号统计分布大都是类高斯分布的, 更利于获取随机数, 但在仿真条件下, 由于模型未考虑增益饱和器件或仪器带宽受限带来的滤波效应, 产生的混沌信号统计分布是类指数分布的. 因此, 在仿真条件下, 我们采用了稍微复杂的后处理来实现随机数提取, 其实在实验条件下, 只需采用更加简单的后处理即可. 在本文中, 我们首先将其原始混沌熵源与其自身延迟的信号作差, 然后通过8位ADC量化采样优化后的混沌熵源, 再采用同时提取多位最低有效位(least significant bit, LSB)和异或(exclusive OR, XOR)处理偏差, 获取了随机比特序列. 采用国际公认的随机数测试套件NIST Special Publication 800-22[42]对上述比特序列进行测试, 证明了本方案产生的随机数可通过此测试标准, 进一步论证了采用激光器阵列后处理混沌熵源实现高品质随机数生成的可行性. 另外, 本方案中仅考虑了两节点激光器阵列, 如果引入更多的节点或实现大型激光器阵列, 通过合理设计, 可以保证每个节点输出混沌信号关联度很低, 最后分别通过后处理, 可实现多路并行的高速高品质随机数生成.

    • 基于激光器阵列后处理的混沌熵源获取随机数的装置如图1所示, 它包含了产生混沌熵源和随机数提取两部分. 其中, 混沌熵源产生部分由光反馈半导体激光器和激光器阵列组成. 光反馈半导体激光器输出光作为原始混沌信号单向注入激光器阵列的节点激光器A中, 经激光器阵列中的两个激光器A、B处理后得到高质量混沌激光. 本文主要采用激光器阵列中的节点B输出作为混沌熵源产生高速高品质随机数. 在我们的仿真实验中, 采用了8位ADC采样量化, m-LSB提取和XOR处理. 本文以在每个样本点中提取3-LSB来获取随机数为例说明此方案产生高速高品质随机数的可行性.

      图  1  基于激光器阵列后处理的混沌熵源获取高品质随机数的示意图($ \lambda /4 $$ 1/4 $波片, PD1、PD2为光电转换器, ADC为模数转换器, LSB为最低有效位, XOR为异或处理)

      Figure 1.  Schematic diagram of high quality random number generation based on the chaotic entropy source generated by ECSL and post-processed by phased-array semiconductor lasers (λ/4, 1/4 wave plate; PD1 and PD2, photo detector; ADC, analog-to-digital converter; LSB, least significant bit; XOR, exclusive OR).

      根据图1所示装置, 混沌熵源的速率方程可写为[41]:

      $\begin{split} {{\dot E_{\rm{M}}}} = \;&\varGamma \frac{c}{{2{n_{\rm{g}}}}}{a_{{\rm{diff}}}}\left( {{N_{\rm{M}}} - {N_{{\rm{Mth}}}}} \right)\left( {1 - {\rm{i}}{\alpha _{\rm{H}}}} \right){E_{\rm{M}}}\\ &+ {k_{\rm{f}}}{E_{\rm{M}}}(t - {\tau _{\rm{f}}})\exp ( - {\rm{i}}{\omega _{\rm{M}}}{\tau _{\rm{f}}}), \end{split}$

      $\begin{split} {{\dot E_{{\rm{A,B}}}}} =\;& \varGamma \frac{c}{{2{n_{\rm{g}}}}}{a_{{\rm{diff}}}}\left( {{N_{{\rm{A,B}}}} - {N_{{\rm{Ath,Bth}}}}} \right)\left( {1 - {\rm{i}}{\alpha _{\rm{H}}}} \right){E_{{\rm{A,B}}}} \\ &+ {\rm{i}}\left( {\omega - {\varOmega _{{\rm{A,B}}}}} \right){E_{{\rm{A,B}}}} + {\rm{i}}\eta {E_{{\rm{B,A}}}}\\ &+ {k_{{\rm{inj}}}}{E_{\rm{M}}}(t - {\tau _{{\rm{inj}}}})\exp \left[ { - {\rm{i}}\left( {{\omega _{\rm{M}}} - \omega } \right)t} \right],\\[-8pt] \end{split}$

      $ {{\dot N_j}} = {P_j} - {N_j}{\gamma _{\rm{N}}} - \frac{c}{n}\left[ {{g_{{\rm{th}}}} + {a_{{\rm{diff}}}}\left( {{N_j} - {N_{j{\rm{th}}}}} \right)} \right]{\left| {{E_j}} \right|^2}, $

      式中: 下标M表示光反馈半导体激光器, A和B分别表示激光器阵列中的节点激光器A和B, j = M, A, B; E(t)表示电场; N(t)表示载流子浓度(N0为透明载流子密度); 下标‘th’表征阈值; Г为光场限制因子; c为光速; ng为群折射率; adiff为微分增益; αH为线宽增强因子; ΩA,B为腔谐振频率; PA,B为抽运率(Pth为阈值电流); γN为腔衰减速率; τp为光子寿命; n为折射率; gth为以及益阈值(Гgth = ng/(p)). (1)式表示光反馈半导体激光器中的反馈项, 其中kfτf分别为反馈强度和外腔中的反馈时间. (2)式中右边第三项表征相控阵列中两全同激光器节点之间的横向耦合, 其中η为复耦合强度, 它的表达式可以在文献[43]中的(1)式中找到. (2)式中最后一项表征从外光反馈半导体激光器的注入效应, 其中kinj为注入强度, τinj为光从外光反馈激光器到激光器阵列的传播时间(不失一般性, 我们假定τinj = 0), ωM是外光反馈激光器的角频率, ω是激光器阵列自由运行的角频率. 因此, 频率失谐可以表示为Δf = (ωMω)/2π. 本论文中, 注入项仅存在于激光器阵列中的节点A方程.

    • 利用四阶龙格-库塔算法对该系统速率方程进行数值求解, 得到激光器输出的混沌信号. 在本文数值模拟中, 相关参数取值如下[41]: αH = 5, a = 4 μm, adiff = 2.5 × 10–16 cm2, γN = 1.0 ns–1, τp = 1.53 ps, N0 = 1 × 1018 cm–3, n = 3.4, P = 1.5Pth. 除非特别说明, 我们选择激光器阵列中A、B节点间的分离比(d/a)为0.5, 其中d为A、B节点间距离的1/2, a为激光器节点宽度的1/2, 波导参数选择带增益引导的反折射率引导, 具体定义和参数可参照我们前期的工作[43]. 不失一般性, 我们选取反馈强度kf = 5 ns–1和反馈时延τf = 1 ns, 此时光反馈半导体激光器工作在混沌状态, 其强度时间序列如图2(a1)所示. 通过计算强度时间序列的自相关(autocorrelation function, ACF), 我们发现在反馈时延τf = 1 ns及其倍数处ACF出现峰值, 如图2(a2)所示. 通过观察图2(a3)给出的频谱, 同样可以发现等间隔的峰值, 频率间隔等于反馈时延τf的倒数. 它们表明此原始混沌信号存在周期性, 不利于获取高品质随机数. 研究发现通过注入到单个激光器或激光器阵列可有效改善此时延特性, 特别地, 在同等注入条件下, 激光器阵列可在更大的参数空间内实现时延隐藏. 图2(b)2(c)分别为光反馈半导体激光器产生的混沌信号注入单个激光器和两节点激光器阵列后的混沌熵源及其ACF与频谱特征. 上述图中选择以注入参数kinj = 30 ns–1与Δf = –30 GHz为例. 通过比较可以发现, 同等条件下, 激光器阵列更适合用于后处理混沌熵源, 其ACF和频谱均无明显峰值, 此结果与我们之前报道的结果具有一致性[41].

      图  2  激光器输出混沌信号的时间序列(左列), 自相关函数谱(中列), 功率谱(右列) (a) 光反馈半导体激光器; (b) 注入激光器; (c) 注入激光器阵列

      Figure 2.  Time series (left column), autocorrelation function (middle column), and power spectra (right column) of the chaotic signal output by laser: (a) ECSL; (b) injection to a single laser A; (c) injection to phased-array lasers.

      在我们另外的工作中, 详细研究了不同波导参数、注入参数和耦合参数对于混沌熵源ACF特征的影响[44]. 图3以带增益引导的反折射率引导波导为例, 给出了时延处的ACF峰值随着注入参数和阵列中激光器分离比d/a的演化情况. 红色代表时延处ACF峰值明显的情况, 而蓝色则表示时延被抑制或消除. 从图3可见, 当分离比较小时, 注入强度不宜过大, 在负频率失谐区域更易实现时延隐藏; 随着分离比的增大, 激光器阵列更易实现时延隐藏. 因此, 通过合理设计, 激光器阵列能够有效提升光反馈半导体激光器产生的混沌熵源的性能. 本文重点证明采用激光器阵列后处理的混沌熵源获取高品质随机数的可行性, 对于时延特征的详细分析不再赘述, 可参照我们的其他同步工作.

      图  3  经过激光器阵列后处理混沌熵源的ACF时延处峰值随着注入参数和激光器分离比d/a的演化情况 (a) d/a = 0.2; (b) d/a = 0.4; (c) d/a = 0.6; (d) d/a = 1.0

      Figure 3.  The evaluation of the ACF peak value located around the feedback delay of the chaotic entropy source that is processed by the phased-array in the plane of injection parameters for several values of laser separation: (a) d/a = 0.2, (b) d/a = 0.4, (c) d/a = 0.6, (d) d/a = 1.0.

      接着, 分析利用激光器阵列处理后的混沌熵源输出经过图1所示后处理产生的二进制序列的特性. 正如前面提到的, 在给定仿真参数下, 由于没考虑增益饱和效应和器件或仪器带宽受限, 混沌信号的统计直方图服从近似的指数分布, 远离理想的高斯分布, 并不利于随机数的直接提取, 这里采用混沌熵源与其延迟特定时间后的混沌信号作差. 所得的混沌熵源的统计直方图如图4(a)所示, 其分布的两边存在较长的尾巴. 将激光器输出的混沌信号经过8位ADC后转变为8位二进制序列, 其中ADC的采样速率为20 GHz. 这里以8位二进制量化序列中提取3位LSB拼接为例, 其随机数生产速率可达60 Gb/s (采用高阶差分和拓展激光器阵列节点数目, 可轻松实现Tb/s量级速率的随机数产生), 所得序列的分布如图4(b)所示, 此时分布的均匀性有所改善, 但仍能看到一些量化比特位出现的概率比其他的高一些. 最后采用常规的XOR处理, 得到如图4(c)所示的近似理想的均匀分布, 适合直接用于获取随机数.

      图  4  激光器B输出的混沌信号量化后的统计直方图 (a) 8位ADC输出; (b) 3-LSB输出; (c) XOR输出

      Figure 4.  Statistical histogram of the quantized chaotic signal of the laser B: (a) The output of 8 bit ADC; (b) the output of 3-LSB; (c) the output of XOR.

      为了进一步说明经过以上后处理的序列具有随机数分布特性, 这里采用散点图来表示. 如图5所示, 横纵坐标分别采用了500个“0” “1”比特位进行统计, 可以看出散点图中无明显的特殊图案, 具有随机分布的特点. 虽然这里仅采用了500 × 500比特位来作图, 但是所提取的随机序列都满足本图所给的特征.

      图  5  散点图

      Figure 5.  Scatter diagram.

      进一步, 我们运用美国国家标准技术研究所提供的NIST SP 800-22测试软件对所获得的3位XOR二进制序列进行随机性测试. NIST测试项目共有15项, 每项测试结果用p表示. 若p值大于显著水平值α = 0.01, 则说明该随机数序列通过了相应的测试[39]. 本测试最终结果是多组多次测试的统计, 采用p的分布P-value来表征, 如果P-value大于$q-3\sqrt{ {q\left(1-q\right)}/{m}}$ (q = 1 – α, m表示测试序列的组数字), 则表示通过该项随机数测试. 在本文的工作中, 我们采集了500组数据进行测试, 则要求每项测试通过率大于0.977. 表1给出了测试结果, 从中可以看出激光器阵列作为混沌熵源获取的随机数能够通过NIST的全部随机数测试标准. 值得强调的是, 我们对原混沌熵源(图1(a))和注入单个激光器后的混沌熵源(图1(b))采用了上述一样的后处理, 得到的随机数序列均未完全通过NIST测试, 这可以说明激光器阵列后处理的有效性和实用性. 虽然这里只给出了一个典型结果, 其实如图3所示的蓝色区域基本均表示激光器阵列有效改善了混沌熵源, 通过一系列后处理得到的随机数序列都达到以上测试条件.

      测试名称P-value概率结果
      频数0.5381820.992通过
      块内频数0.2392660.982通过
      累加0.7558190.994通过
      游程0.1404530.988通过
      块内最长游程0.9658600.988通过
      矩阵秩0.2812320.990通过
      离散傅里叶变换0.2066290.982通过
      非重叠模块匹配0.0208310.982通过
      重叠模块匹配0.6993130.984通过
      通用统计0.5101530.994通过
      近似熵0.6993130.994通过
      随机游动0.4436650.986通过
      随机游动变量0.2901580.983通过
      连续性0.0965780.984通过
      线性复杂度0.3408580.986通过

      表 1  NIST统计测试结果

      Table 1.  Result of NIST statistical tests.

      最后, 我们强调激光器阵列后处理光反馈激光器产生混沌信号的另一个优势, 即是其可同时获取多路相关或不相关的随机数序列. 通过选择参数, 激光器阵列中A、B节点的混沌输出均不具有时延特征. 如以kinj = 7 ns–1, Δf = –30 GHz为例, 结果如图6所示, 激光器阵列中节点A、B均实现高维混沌输出, 而且在ACF图中无时延特征峰值. 通过必要后处理, 很容易得到两路高品质的随机数序列. 其实激光器阵列的可扩展性强, 可以实现多节点甚至大型节点的激光器阵列, 因此本文结果可为实时产生多路并行的高速高品质的随机数序列提供思路.

      图  6  激光器输出的时间序列与自相关函数 (a) A激光器输出的时间序列; (b) A激光器输出的自相关函数; (c) B激光器输出的时间序列; (d) B激光器输出的自相关函数

      Figure 6.  Time series and autocorrelation function of the lasers: (a) Time series of laser A; (b) autocorrelation function of laser A; (c) time series of laser B; (d) autocorrelation function of laser. B.

    • 本文采用激光器阵列后处理光反馈半导体激光器产生的混沌熵源, 再经过光电转换、模数转换采样量化及m-LSB提取和XOR处理, 最终产生随机数序列. 研究结果表明, 激光器阵列使得光反馈半导体激光器产生的原始混沌信号的弱周期性得到有效抑制, 可作为混沌熵源来获取随机数. 经过常规后处理, 随机数序列分布均匀, 通过了随机数行业测试标准(NIST SP 800-22)中的全部15项测试, 如果扩展激光器阵列中的节点数, 则有望实现同时获取多路并行的高速高品质随机数.

参考文献 (44)

目录

    /

    返回文章
    返回