具有超强可扩放性的Tbps全光并行真随机数发生器的制造方法
【专利摘要】一种具有超强可扩放性的Tbps全光并行真随机数发生器是由保偏光纤依次串连有理数谐波锁模光纤激光器、脉冲光放大器、高非线性色散位移光纤及反常色散光纤，构成一超连续谱熵源；超连续谱熵源输出的脉冲序列经阵列波导光栅被切割产生出 N 路窄带子光脉冲序列，进入光衰减器阵列进行调节，后经 N 路光纤进入全光比较器阵列被量化成 N 路高速真随机脉冲序列，再由光滤波器阵列滤出。本发明单路码率高达Tbps量级，将现有并行真随机码发生器的单路速率提高了3个数量级；可扩放性超强，能同时输出至少10000路的独立、并行真随机码，极大地满足了大规模并行计算及高速保密网络通信的当前需要。
【专利说明】具有超强可扩放性的Tbps全光并行真随机数发生器
【技术领域】
[0001] 本发明与一种并行真随机数发生器有关，尤其是一种具有超强可扩放性的Tbps 全光并行真随机数发生器，适用于并行蒙特卡洛仿真、大规模并行计算及保密通信领域。
【背景技术】
[0002] 并行蒙特卡罗（MonteCarlo)仿真在核物理、计算化学、生物医学、金融工程学、宏 观经济学等领域有着重要应用。随机数是并行蒙特卡罗仿真的基石，其面临的首要问题就 是多路并行随机数的高速产生，称作"并行随机数发生器"。
[0003] -个理想的并行随机数发生器需要满足以下条件：1)不含序列内相关性 (intra-streamcorrelation)。即每个处理器上所用的随机数序列必须具有高质量的随机 特性；2)不含序列间相关性（inter-streamcorrelation)。即多个并行处理器所用的多路 随机数序列之间要相互独立；3)具有可扩放性（scalability)。即根据实际需要，随机数发 生器可同时产生出任意N路独立的随机数序列。
[0004]利用计算机，通过一定的并行化算法可方便地获得并行随机数。该一类并行随机 数发生器具有高速率、低成本、易构建等优点，但存在着众多显著问题。其中最致命的是，该 类发生器是基于种子和确定性算法实现的，不可避免地存在周期性一这一 "阿喀琉斯之 踵"严重限制了其大量产生随机数的能力，局限了其在大规模并行计算领域的应用。因此， 人们将这类并行随机数发生器称作"并行伪随机数发生器"。
[0005]基于自然界随机现象构建并行随机数发生器，可提供非周期、不可预测、无限数量 的真随机数，称作"并行真随机数发生器"。1986年，美国罗切斯特大学的J.Marron等人提 出了首个并行真随机数发生器的雏形，是利用一个二维探测器阵列对激光器散斑分布进行 探测和编码实现的[為叩丄办t25(1): 26-30 (1986)]。1987年，法国巴黎大学F.Devos 进一步发展了这一技术[办t 12(3): 152-154 (1987)]。但遗憾的是，传统并行真随 机数发生器受限于传统物理熵源（如热噪声、振荡器抖动等）带宽过低且可扩放性差，码率 处于Mb/s量级,与实际需求相去甚远。
[0006]近年来，随着光子熵源的出现，串行真随机数发生器取得了跨越式发展，速率可达Gbps量级。典型的实现方法有：1)基于放大自发辐射光噪声（ASE)提取真随机数。如，美 国C.R.S.Williams等人2010年利用铲杂光纤放大器中的ASE获得了 12. 5Gbps的真 随机数[办t^>/7/1^18(23)，23584 - 23597，2010]。2)基于混沌激光提取真随机数。 如， 申请人:所在课题组2013年利用混沌半导体激光器实验实现了 4. 5Gbps真随机数的产 生[办t 21(17): 20452-20462,2013]。3)基于量子真空态获取真随机数。如， 澳大利亚T.Symul等人2011年基于真空态构建了 2Gbps随机数产生系统/--?. Lett.98(23): 231103, 2011] 〇
[0007]但是，这些真随机数发生器却属于"串行"随机数发生器，只能输出一路随机码序 列，不符合高质并行真随机数发生器的要求，无法应用于大规模蒙特卡洛仿真及并行计算 领域。此外，这些随机数发生器提取随机数的过程均利用电子器件在电域中完成，其实时速 率最终将面临"电子瓶颈"的限制，约几十个Gbps。
[0008] 综上所述，并行伪随机数发生器能快速产生随机数(典型速率处于Gbps量级)，但 无法克服算法本身固有周期性的限制，不具备产生大量随机数的能力；传统并行真随机数 发生器拥有产生大量无周期随机数的能力，但受限于随机数信号源带宽，无法实现真随机 数的快速产生，典型速率低于1〇〇Mbps;近年来发展起来的真随机数发生器已具有Gbps的 快速随机数产生能力，但却只能输出一路随机数，可扩放性极差。
[0009] 然而，蒙特卡罗仿真在并行环境中的计算量至少是串行情形下的10?IO5倍，要 求相应随机数的产生速度和数量均需大规模提升。根据2014年"国际TOP500组织"公布 的最新全球超级计算机500强榜单，当前并行计算机已拥有并发执行数千、甚至上万只处 理器的能力，要求具有与之相匹配的可扩性能力的并行真随机数发生器。因此，发展与当前 需求相匹配、兼具超强可扩放性及快速产生大量随机数能力的并行真随机数发生器已迫在 眉睫。
【发明内容】
[0010] 本发明的目的是提供一种具有超强可扩放性的Tbps全光并行真随机数发生器， 以解决现有技术码率不足及可扩放性差的问题，从而公开一种同时具有超高码率和超强可 扩放性的新型并行真随机数发生方案，以满足现代蒙特卡洛仿真、大规模并行计算及网络 安全通信等现实需求。
[0011] 一种具有超强可扩放性的Tbps全光并行真随机数发生器，其特征在于： 一保偏光纤中，依次设置有有理数谐波锁模光纤激光器、脉冲光放大器、高非线性色散 位移光纤及反常色散光纤，构成一超连续谱熵源； 一超连续谱熵源将输出的脉冲序列经阵列波导光栅被切割产生出靡各窄带子光脉冲 序列，进入光衰减器阵列进行调节，后经靡各光纤与连续光激光器输出的靡各连续光信号同 时进入全光比较器阵列被量化成靡各高速真随机脉冲序列，再由光滤波器阵列滤出。
[0012] 进一步的技术特征在于：所述超连续谱熵源输出的脉冲序列的重复频率/= 1THz;所述阵列波导光栅具有#输出波长通道；所述光衰减器阵列是由#光衰减器并列 构成；全光比较器阵列是由A个相同的全光比较器并列构成；所述光滤波器阵列是由外 相同的光滤波器并列构成；所述靡各高速真随机脉冲序列的码率是由超连续谱熵源输出脉 冲的重复频率/7夬定，等于ITbps;所述_取值为10000。
[0013] 实现上述本发明所提供的一种具有超强可扩放性的Tbps全光并行真随机数发生 器，与在先并行随机数产生技术相比，其优点与积极效果在于： 第一，不存在周期性，可提供无限数量的真随机数序列，克服了并行伪随机数发生器固 有周期性的局限； 第二，单路码率可达Tbps量级，将现有并行真随机码发生器的单路速率提高了 3个数 量级； 第三，可同时输出至少10000路的独立、并行真随机码，将现有并行真随机数发生器的 可扩放性提高了 3~4个数量级，极大满足了大规模并行计算及高速保密网络通信的当前需 要； 第四，本发明的并行真随机数发生器的信号处理过程均在光域中进行，不需要任何光 电转换装置及电子模数转换设备，突破了"电子瓶颈"的限制。
【专利附图】
【附图说明】
[0014] 图1是本发明Tbps全光并行真随机数发生器的结构示意图。
[0015] 图2是本发明全光比较器阵列中任一单元模块的结构示意图。
[0016] 图3是本发明超连续谱熵源输出超连续谱脉冲序列的光谱图。
[0017] 图4是本发明靡各窄带窄带子光脉冲序列中任意一路的时序图。
[0018] 图5是本发明全光比较器单元的透射传递特性函数。
[0019] 图6是本发明任意一路全光真随机码时序图。
[0020] 图7是本发明任意一路并行随机码的自相关图。
[0021] 图8是本发明任意两路并行随机码的互相关图。
[0022] 图9是实现本真随机数发生器的真随机数产生方法流程图。
[0023] 图中：1 :有理数谐波锁模光纤激光器；2 :脉冲光放大器；3 :高非线性色散位移光 纤；4 :反常色散光纤；5 :阵列波导光栅；6 :光衰减器阵列；7 :全光比较器阵列；7la:3dB耦 合器I;7lb:光隔离器I;7i。:光隔离器II;7ld:耦合器I;7le:耦合器II;7 1/:波分复用器； 7lg:高非线性光子晶体光纤7lh:3dB耦合器II;8 :连续光激光器；9 :光滤波器阵列。
【具体实施方式】
[0024] 实施本发明上述所提供的一种具有超强可扩放性的Tbps全光并行真随机数发生 器，是由超连续谱熵源、阵列波导光栅5、光衰减器阵列6、全光比较器阵列7、连续光激光器 8和光滤波器阵列9构成；其中，所述超连续谱熵源输出的高重复频率/、峰值功率大幅度随 机起伏的超连续谱脉冲序列经阵列波导光栅5光谱切割产生出靡各峰值功率大幅度随机起 伏的窄带子光脉冲序列，进入光衰减器阵列6中，被各个光衰减器单元（6 :、62、63、…、6#)调 节至平均功率一致，而后经靡各光纤进入全光比较器阵列7,被各自对应的全光比较器单元 (7p72、73、…、7#)及连续光激光器8共同作用量化成靡各独立的高速真随机脉冲序列，分别 由光滤波器阵列9中各自对应的光滤波器单元（9 :、92、93、…、9#)滤出。所述超连续谱熵源 由有理数谐波锁模光纤激光器1、脉冲光放大器2、高非线性色散位移光纤3及反常色散光 纤4构成的，各器件之间通过保偏光纤依次连接；所述超连续谱熵源输出的高重复频率八 峰值功率大幅度随机起伏的超连续谱脉冲序列中作ITHz;所述阵列波导光栅5具有 A个输出波长通道；所述光衰减器阵列6由#光衰减器并列构成；所述全光比较器阵列7 由A个相同的全光比较器并列构成；所述光滤波器阵列9由外相同的光滤波器并列构成； 所述靡各独立的高速并行真随机码序列的码率均由超连续谱熵源输出脉冲的重复频率/7夬 定，等于ITbps;所述施]取值为10000。
[0025] 下面将结合附图对本发明的【具体实施方式】作出详细说明。
[0026]1、本发明超连续谱熵源的实现过程 如附图1中所示，有理数谐波锁模光纤激光器1输出脉宽l±〇. 〇5ps、重复频率ITHz、 波长为1550nm的超短光脉冲序列，经脉冲光放大器2作用后，其峰值功率可以增大到6 kW。以该超短脉冲信号作为泵浦源经保偏光纤进入到一段长5m、非线性系数为25/W/km高 非线性色散位移光纤3,零色散点位于1550nm处，受高非线性色散位移光纤3中自聚集、自 相位调制、交叉相位调制、四波混频和受激拉曼散射等效应的共同作用，泵浦脉冲的光谱中 会产生许多新的频率成分，使得输出脉冲序列的光谱宽度远大于入射脉冲的谱宽，最终产 生谱宽可达2.0 以上的超连续谱光脉冲信号，如附图3所示。
[0027] 由于噪声信号的存在，此时的超连续谱光脉冲信号峰值功率会有微弱起伏，但远 不足以满足后续电子量化系统的要求；为了进一步增强上述超连续谱光脉冲峰值功率的起 伏，本发明引入一段长10 m的反常色散光纤4 ;原始的超连续谱光脉冲信号在反常色散光 纤4传输过程中，噪声驱动的调制不稳定性将发挥主导作用，致使超连续谱光脉冲信号的 稳定性大幅劣化，脉冲峰值强度呈现出强烈的随机起伏。这就实现了本发明的高重频、强度 大幅度起伏的超连续谱熵源。
[0028] 2、本发明靡各并行真随机码的提取过程 如附图1所示，超连续谱熵源产生的重复频率I THz、峰值强度大幅度随机起伏的超连 续谱光脉冲序列经通道间隔为〇. 2 nm、通道数浪10000路的阵列波导光栅5光谱切割后， 可分离出浪10000路的窄带子光脉冲序列；这些序列遗传了超连续谱熵源产生的超连续谱 脉冲序列的高重频及峰值功率大幅度起伏特性，且彼此之间相互独立。窄带子光脉冲序列 相互之间完全独立的本质原因在于超连续谱脉冲序列随机起伏起源于量子独立的激光自 发辐射噪声且阵列波导光栅5各个输出通道光谱上无重叠。附图4是所产生的10000路 窄带子光脉冲序列中任意一路的典型时序图：峰值功率大幅度随机起伏、重复频率固定为1 THz。
[0029] 以上产生的浪10000路的窄带子光脉冲序列被光衰减器阵列6中，光衰减器阵列6 由M v光衰减器6i、62、63、…、6#并行构成，经各自对应的光衰减器单兀6i、6 2、63、…、6#调 节至平均功率一致后，被作为子随机数提取源，进入全光比较器阵列7中，经各自对应的全 光比较器单元7 :、72、73、…、连续光激光器8共同作用量化成靡各独立的高速真随机脉 冲序列。本发明全光比较器阵列7中的各个单元7 :、72、73、…、7N是完全相同的装置，工作 过程亦相同。因此，下面将随机抽取其中全光比较器单元为例，对靡各高速真随机脉冲 序列的产生过程予以说明。
[0030] 附图2是附图1中全光比较器阵列7中全光比较器单元的详细结构示意图。输 入端A、C以及输出端B分别与附图1中全光比较器阵列7中全光比较器单元A、B、C 三端口--对应。
[0031] 由A端口输入的窄带子光脉冲序列，浪10000路窄带子脉冲序列中对应全光比较 器单元7 :的那一路，作为控制光I经过附图2中的波分复用器7 ^进入环路。连续光激光 器8输出的连续光信号作为探测光由C端口进入全光比较器单元，通过3dB耦合器I7la 等分成两路，本发明中称其为上臂、下臂两路。上臂连续光信号经光隔离器I7lb通过耦合器 7ld分成两路：一路连续光信号II通过耦合器17 1(1直通臂向前传输，另一路连续光信号III 则通过耦合器171(1的耦合臂进入高非线性光子晶体光纤7 lg构成的环路中，与控制光信号I 同向传输。传输过程中，连续光信号III相位受到信号线性相移，自相位调制以及与控制光 之间的交叉相位调制的影响。连续光信号III在环路里传输一周后与信号II在耦合器I7 ld 叠加形成新的光场信号IV ;同理，下臂连续信号经光隔离器II7lc;通过耦合器117 16分成两 路：一路连续光信号V通过耦合器II7le的直通臂向前传输，另一路连续光信号VI高非线性 光子晶体光纤7 lg构成的环路中，与控制光信号I逆向传输。此时连续光信号VI与控制光 I之间的交叉相位调制效应可忽略，连续光信号VI相位只受到信号线性相移以及自相位调 制的影响。连续光信号VI在环路里传输一周后与直通臂中的信号V在耦合器II7le叠加形 成新的光场信号VII。最后，新生光场信号IV与信号VII在另一 3dB耦合器II7lh处干涉 由端口B输出。从而，实现对A端口输入的窄带子光脉冲序列的全光量化处理。
[0032] 接下来给出定量分析：上、下臂连续光信号II、III在高非线性光子晶体光纤环形 腔中受到的非线性效应的差异会使得两路光信号之间产生相位差。连续光信号III、V与 II、VI在光耦合器I7ld、光耦合器11716叠加产生新的光场IV、VII的相位分别可以表示为： ⑦ "=⑦。+6 31 /?^/』Jeff+4 31 /--/口An= %+6 31 /?^/』Jeff。这里，⑦。和』分别 是上、下两臂连续光信号的线性相移及波长，7、/?及為ff则是高非线性光子晶体光纤7 18构 成的环形腔的长度、非线性折射率及发生非线性效应的有效横截面积。考虑到光耦合器I7ld 和光耦合器II7le的耦合系数均为r，新生光场IV、VII耦合输出时的有效相位差A0#可 以表达为： AC^eff=arctan[ (1+r2)tan(<?IV/2) / (1-r2)] -arctan[ (1+r2)tan(<?VII/2) / ( 1-r2)] 从而，当新的光场IV、VII在3-dB耦合器II7lh处干涉输出时，透射率f将可以表达为 7Kl-C〇s(A%ff)]/2。合理选择光耦合器71(1、716的耦合系数均为r，可以使得有效相位差 也实现了 n"的跳变，由光路干涉透射率r=[l-cos(A￠#) = /2公式知，输出端口B 处的透射率将实现0和1之间的跳变：当透射率为0时，无脉冲输出，产生"0"码；当透射率 为" 1"时，有脉冲输出且输出脉冲功率恒定，产生" 1"码。
[0033] 附图5是全光比较器单元的透射传递函数。所选用的高非线性光子晶体光纤7 18的 长度人非线性折射率/?及发生非线性效应的有效横截面积Aff分别是〇. 2m、4. 95XKT19 m2/W和10iim2,所选用的光耦合器7ld、7le的耦合系数r均为0.99。由附图可见，全光比较 器阵列7中每个全光比较器单元的投射传递函数是一方波函数，具有陡峭的比较阈值Pth。 当A端口输入的窄带子光脉冲峰值功率大于该阈值时，B输出端有脉冲输出，编码为" 1";反 之，B输出端口无脉冲输出，编码为"0"。最后，这些全光真随机脉冲序列进入光滤波器阵列 9中的各个子单元模块滤出。这里指出，光滤波器阵列9由N个相同的光滤波器并列形成， 中心波长与连续光激光器8完全对应。附图6是相应的经过全光比较器单元作用后产生的 重频ITHz、脉宽约Ips的全光真随机脉冲序列，该重复频率由超连续谱熵源发射的超连续 谱脉冲序列重复频率/7夬定。
[0034] 3、靡各并行真随机码的随机性检验 为了验证本发明并行随机数发生器的质量，本发明任意选取了浪10000路并行真随机 码中的任意两路随机码信号进行了自相关和互相关分析。分析结果如附图7和附图8所示。 附图7是任意一路随机码的自相关特性曲线，它呈现出类S函数形状，这就意味着本发明 产生的每路随机码序列内的码元是统计无关的，具有良好的随机特性；附图8是任意两路 随机码的互相关特性曲线，互相关系数均在0附近，这就表明本发明任意两路并行随机码 序列间是统计无关的。
[0035] 此外，还利用国际随机数行业测试标准，Diehard统计测试包，对本发明所产生的 随机码进行了系统测试，结果证明本发明产生的随机码可通过所有测试项，符合行业标准。 表1是一典型测试结果。
[0036]
【权利要求】
1. 一种具有超强可扩放性的Tbps全光并行真随机数发生器，其特征在于： 一保偏光纤中，依次设置有有理数谐波锁模光纤激光器（1)、脉冲光放大器（2)、高非 线性色散位移光纤（3)及反常色散光纤（4)，构成一超连续谱熵源； 一超连续谱熵源将输出的脉冲序列经阵列波导光栅（5)被切割产生出靡各窄带子光脉 冲序列，进入光衰减器阵列（6)进行调节，后经靡各光纤与连续光激光器（8)输出的靡各连 续光信号同时进入全光比较器阵列（7)被量化成靡各高速真随机脉冲序列，再由光滤波器 阵列（9)滤出。
2. 如权利要求1所述的Tbps全光并行真随机数发生器，其特征在于：所述超连续谱熵 源输出的脉冲序列的重复频率/ = 1THz。
3. 如权利要求1所述的Tbps全光并行真随机数发生器，其特征在于：所述阵列波导光 栅（5)具有#输出波长通道。
4. 如权利要求1所述的Tbps全光并行真随机数发生器，其特征在于：所述光衰减器阵 列（6)是由A个光衰减器并列构成。
5. 如权利要求1所述的Tbps全光并行真随机数发生器，其特征在于：所述全光比较器 阵列（7)是由A个相同的全光比较器并列构成。
6. 如权利要求1所述的Tbps全光并行真随机数发生器，其特征在于：所述光滤波器阵 列（9)是由A个相同的光滤波器并列构成。
7. 如权利要求1所述的Tbps全光并行真随机数发生器，其特征在于：所述靡各高速真 随机脉冲序列的码率是由超连续谱熵源输出脉冲的重复频率夬定，等于1Tbps。
8. 如权利要求1、3、4、5、6或7所述的Tbps全光并行真随机数发生器，其特征在于：所 述施]取值为10000。
【文档编号】G06F7/58GK104516715SQ201410831052
【公开日】2015年4月15日 申请日期:2014年12月29日 优先权日:2014年12月29日
【发明者】李璞, 王云才, 张建忠, 张建国, 王冰洁, 王安帮, 张明江 申请人:太原理工大学