一维logistic映射研究.docxVIP

Introduction 混沌是一种貌似无规则的运动，指在确定性非线形系统中，不需要附加 任何随机因素亦可出现类似随机的行为（内在随机性），这种过程即非周期 又不收敛，并且对初始值有极其敏感的依赖性，因此从长期意义上讲，系统 的未来行为是不可预测的。 一维Logistic映射从数学形式上來看就是一个非常简单的混沌映射， 早在20世纪50年代，就被生态学家用来描述种群的变化，此外，基于变形 Logistic映射的混沌流密码算法，在保密通讯领域应用也十分广泛。 一维Logistic映射描述如下： X”i书（心）=凡（1-X』 其中V衣4,该模型看起來简单，并且是确定性的, 但具有极为复杂的动力学行为，产生的序列由入和X的初始值控制，这两个 值任何一个出现微小差别时，所产生的序列将会完全不同。本文即以 F0RTRAN90为编程语言，对一维Logistic映射的性质进行了探索。 explain programming 程序包括module.主程序和若干子程序，下面对各部分一一作简单介绍。 2.1 module Module useful中时装了所有全局变量，需要使用这些变量的函数只要 写入use module usefu 1语句即可。Module usefu 1中，习题要求的最大迭代 步数1000、输出最后10步结果等，均以常量的形式写入；由于希望观察在 不同入下系统特性，因此入使用了一个可变大小数组，在程序运行中根据输 入条件动态取得一系列数值。考虑到非线形迭代对数值非常敏感，Module及 程序其它地方声明的实数均为双精度浮点数。 2.2主程序 Program main中主要调用各个子程序，并将部分结果写入输出文件。由 于程序在最初调试时，只设计了一个入，调试成功后，为简单起见，将多个 入的循环部分放在主程序中。 2.3子程序 子程序依次包括input pretreatx iterative和compare,分别叙述如 下： Input ():输入参数入的最小值kO和最大值kn,以及希望在kO到kn之 间划分的多少等份imk；输入x序列的初始值xO。程序中判断，如果取值范 围不满足要求，提示重新输入。 Pretreat ():对输入进行预处理。希望得到当入在kO到kn之间变化的 多个结果，因此根据输入的imk将区间［kO, kn］进行等分，当i mkHZl时， 参数入取kn值。 Iterative ():主程序中循环调用该子程序，根据一维Logistic映射公 式递归计算Xn序列，程序中用实数数组xval()表示。 Compare ():分析Xn序列特性，是程序的关键部分。算法叙述如下： 从序列点X2开始遍历，直至X1000,执行以下操作； 对点Xn,向前搜索20个点(M20时，搜索到XI停止)，按就近原则 逐一比较，若W|Xn -Xm| error,认为可能有周期出现，进入3) 进一步判断，若没有，则继续1) 比较Xm、Xm+1 Xn-1和Xn、Xn+1……X2n-m-l,若误差均小于常 量error,认为是一个完整的周期，周期数为否则继续1)； 程序中通过逻辑变量blstart、blover的操作，明确当某一个周期数 连续的、不变的出现20次以上时，且结束于最后，认为序列是有周 期性的，否则，认为是混沌的。 output 程序有两种输出情况，一种按习题要求希望得到Xn序列最后10个点 的数值，程序中判断，当imk取1时，输出结果。屏幕输出如下： Msg:input k0^ kii imk... 3.3 3.3 1 Msg: input x0.?. 0.7 Msg:output the last points... 0.479427019824234 0.823603283206069 0.479427019824234 0.823603283206069 0.479427019824235 0.823603283206069 0.479427019824234 0.823603283206069 0.479427019824234 0.823603283206069 0.479427019824235 ki= 3300000000 norbit= 2 另一种情况，当imk较大时，仍在屏幕上输岀则数目较多，则将序列 点写入文本，为便于分析，文本中写入了最后200个点。此时，屏幕中只输 出不同ki时的周期数，若没有周期，输岀chaoticT ,输出如下： Msg: input k0 kn、 imk... 2.5 4.0 10 Msg: input xO... 0.6 ki= 2.500000000 norbit= 1 ki= 2.666666667 norbit= 1 ki= 2.833333333 norbit= 1 ki=