混沌伪随机序列发生器matlab,基于变结构混沌的伪随机序列发生器测试

内容摘要：为产生随机性能良好的伪随机序列，提出了一个新的变结构混沌系统。该混沌系统在一个开关函数控制下其系统结构随时间随机地转换，所产生的混沌信号是两个不同的混沌信号的混合，具有良好的复杂性。基于该变结构混沌系统设计了一种伪随机序列发生器，采用NIST标准和STS-2．0b测试套件对其产生的伪随机序列进行了统计性能测试，测试结果表明该伪随机序列发生器具有良好的随机性，可应用于计算机、通信、信息加密等领域中。

关键词：混沌；变结构混沌；伪随机序列；随机性

引言

伪随机序列在数字通信、密码系统、计算机仿真等领域有着广泛的应用。一个伪随机序列发生器包括随机信号源(种)和一系列的离散、量化及其实现技术，其中良好的随机信号源是伪随机序列设计的关键问题。混沌与传统密码学之间存住着一种自然的联系，混沌动力学特性基本对应着高强度密码系统的某些安全特征，而具有良好混合特性的传统密码又蕴涵着混沌现象。以混沌作为信号源为伪随机序列发生器的设计提供了一种新的途径。

利用连续和离散混沌系统进行伪随机序列发生器的设计已有研究。离散混沌由于算法简单致使其运算速率快，序列码率较高，但缺点是系统参数和初值条件在一般情况下较少，密钥空间小，序列的安全性较低。连续混沌一股情况下是几个非线性微分方程的耦合，其系统参数和初始条件较多，产生伪随机序列的密钥空间较大，缺点是运算复杂，在数字系统实现时运算速率相对较慢。但如果采取合理的量化方法，会较好地弥补这种慢的运算速率。如在抽位量化方法中，如果一次抽取混沌数字迭代值的多位作为0，1序列，可大大提高其码率。因此采用复杂的连续混沌系统作为伪随机序列的源将是混沌序列应用的一个方向。

另一方面，数字系统的编码理论表明，在数字系统中处理非周期的混沌时，由于系统本身的有限位数致使混沌出现周期现象，即短周期或动力学退化问题。为改善这种短周期问题，可通过对混沌系统的变量或参数进行扰动以提高其数字PN序列的统计性能，增大序列的周期。为了提高混沌伪随机序列的复杂性和改善其动力性退化问题，本文设计了一个变结构混沌系统，以期获得性能更好的伪随机序列。所谓变结构混沌系统，是指该系统的代数结构不断地自动变化，而实现这种变化的控制函数是一个开关函数，该函数在自身变量控制下自动地在0，1之间转换。在提出变结构混沌系统之后，对基于该混沌系统的伪随机序列发生器进行了设计，对产生的伪随机序列进行了NIST(National Ins titute of Standards and Technology)测试。测试结果验证了该数字序列具有良好的随机性能。

1 变结构混沌系统构造

首先构造了一个三维连续混沌系统：

式中：a，b，c为可变的系统参数。在Matlab软件平台上计算表明，在较大的a，b，c参数范围内系统(1)都是混沌的，取a=0．8，b=1．5和c=1．5时系统(1)的时域波形和y-z平面上的轨迹(相图)如图1所示。

又构造了另一个三维连续混沌系统：

式中：a，b，c和是为可变的系统参数。在Matlab软件平台上计算表明，在较大的a，b，c和k参数范围内系统(2)都是混沌的，取a=0．8 b=1．5，c=1．5和k=0．32时系统(2)的时域波形和y-z平面上的轨迹(相图)如图2所示。比较图1和图2发现，两者的时域波形和对应的平面轨迹不同。

混沌系统(1)、(2)除第一个方程不同外，其余两个方程完全相同。虽然它们有相似的结构，但其代数结构不同，平衡点也不同(见下一节)，因而它们是非拓扑等价的，即它们在理论上是不同的两个系统。根据两个系统有相似结构但本质不同的特点，采用一个开关控制函数u构造了一个变结构混沌系统：

式中：m为开关控制函数的门限，m∈x取m=0．2，其他参数同前。对变结构混沌系统(3)进行仿真计算，所获得的时域波形x-t和y-z平面上的轨迹如图3所示。

图3中，实线和虚线分别为为系统(1)和(2)的波形或轨迹。

从图3看出，该系统的信号波形或解的轨迹由两个不同的部分构成。当系统的解x≥m=0．2时，u(x-m)=1，混沌系统(3)为混沌系统(2)的结构；当系统的解x

这种由两个不同的混沌信号按时间随机地混杂在一起而形成的一个完整的混沌信号，比之由单一混沌系统产牛的信号要复杂得多，且门限参数本身又是一种密钥参数，它扩展了混沌伪随机序列的密钥空间，使其提高了安全性。

2 伪随机序列发生器设计及性能分析

基于上述的变结构混沌系统可设计一种新的伪随机序列发牛器。主要思路是以变结构混沌系统作为随机信号源，采用一定的方法对其离散、量化，获得一系列的伪随饥序列。

这里研究的变结构混沌系统是一个非线性常微分方程组，在数字系统中对其进行数值解就是一种离散的方法。常微分方程近似求解的数值方法有欧拉算法、改进型的欧拉算法和龙格库塔法等，这都是将连续系统进行近似离散化的方法。其中，欧拉算法速率最快，本文采用欧拉算法将连续混沌离散化。对于一个连续的混沌系统，有：

当τ足够小时，经过欧拉算法离散化后的系统具有与式(3)所示的连续混沌系统相同的动力学特性，此处选择τ=0．004。

在数字系统中迭代求解式(8)所示的离散化系统，迭代过程中的每一个解变量xn，yn和zn都可以通过二进制数据的方式来表示。以xn为例：

式中：b1n，b2n，…，b(k+1+l)n分别为二进制数的所有位(0或1)，混沌系统的解xn随时间不断变化，其二进制表达式中的每一位bm(“0”或“1”)也随时间小断变化。如果抽取随时间变化的一位或多位，可构成一个由“0”或“1”组成的伪随机序列。为了保证提取的序列具有较好的随机性，可以严格地从小数部分中提取其中一位作为随机序列，也可以从{b1n，b2n，…，b(k+1+l)n}中选取随机性能较好的多位作为随机序列，从而增加随机序列的提取速度。这种量化方法可用图4表示。

式(5)～式(9)描述了混沌伪随机序列发生器设计的核心算法。实现一个混沌伪随机序列发生器可借助于软件和硬件平台。如果为计算机或其他软件提供伪随机序列，可借助数字计算机这个性能完善的平台实现式(5)～式(9)的运算，如可用Matlab，C语言等软件实现一个混沌伪随机序列发生器。也可结合实际应用在相关信号处理软硬件平台上实现混沌伪随机序列发生器，如利用DSP芯片对语音或视频信号进行混沌加密，可在DSP内进行上述运算而实现混沌伪随机序列发生器，也可利用FPGA硬件平台实现这种伪随机序列发生器。本文不侧重利用何种平台，如何实现混沌伪随机序列发生器，而是着重基于上述变结构混沌系统的伪随机序列发生器性能的测试。为此，选择Matlab求解变结构混沌系统，通过实现式(5)～式(9)的运算产生一系列伪随机序列，提取序列并进行序列的随机性统计测试。

描述一个序列随机性统计性能的指标有多种，但目前应用最广的是NIST(National Institute of Standardsand Technology，美国国家技术与标准局)标准。NIST推出2．0版本的测试软件包STS是当前最具权威的一种随机性检测工具，它为研究人员提供了一种量化的报告，显式地说明一个伪随机序列性能的好坏。STS-2．0b是当前最新的软件包版本，由十五项核心测试指标组成。

该测试包评价序列性能好坏有两项指标：其一是通过率，另一项是P-value分布的均匀性。测试独立生成的m组随机序列，依据各组每次测试的P-value值是否大于测试水平α=0．01来计算通过率。若各次测试的通过率在可信性区间

内，其中

1-a，则可说明此次测试算法的信任度高。对于P-value分布均匀性，若P-valueT>0．000 1，则说明P-value值是均匀分布的。

在Linux操作系统环境下进行测试。通过编程将变结构混沌系统进行离散迭代运算来产生数字混沌序列，然后将产生的二进制数字序列保存为txt文档，并通过测试指令调用软件包对txt文档中的序列进行测试，测试由STS软件包自动完成，并生成测试报告。基于变结构混沌系统产生的伪随机序列的测试结果如表1所示，序列共有100000000 b，以每组100000 b分为1 000组。

从表1中P-value这一列看出，序列仅在FFT这一项中的P-value值测试不满足P-valueT>0．000 1的条件，这说明序列在该项测试中的P-value值分布不均匀，在其余14项测试中表现为分布均匀。若从通过率来分析，取显著水平α=0．01，那么根据通过率可信区间的计算公式

可得，当PROPORTION的值落在(0．980 560 8，0．999 439 2)区间内时，表明序列通过该测试项，反之则为不通过。表1测试结果显示序列在所有测试中其结果均落在可信区间之内，所有指标均通过该项测试。

3 结论

为产生性能良好的伪随机序列，本文构造了一个新的变结构混沌系统。该系统在一个开关函数控制下自动地在两个混沌子系统之间随时问随机地转换，所产生的混沌信号是两个不同的混沌信号的混合，因而具有较好的复杂性。利用该变结构混沌系统设计了一种伪随机序列发生器，基于NIST标准和STS-2．0b测试套件对其产生的伪随机序列进行了测试，序列通过率全部通过了测试，序列的均匀性只有一项未通过测试。测试结果表明，该伪随机序列发生器具有良好的随机性能，可应用于计算机、通信、信息加密等领域之中。