仿真
3 8
合肥学院学报 (然科学版 )自
第 1卷 6
程序如下: fnt n f oii; uco=l sc i gt
i 1;=4;=
%建立一个函数名
f1=+; ( ) O 3 w i i 10 he< 0 l
i= i+1:
%定义初值%进入循环迭代 10次 0
fi 1= fi (一 (); (+ ) () 1 fi)
po() l f t
ed n
12二维离散混沌系统 .
考
虑文献[] 4提出的二维离散电路混沌系统,该系统的动力学方程为
『( n+1 )=a ( )一d n yn y( ) t(, )=一6 ( )+c( ) y/+1/ n y n,
由离散动力系统分支理论及数值计算可知,在分岔点的右侧邻域内,系统出现离散意义下的 hp分 of
支,即在平衡点附近系统出现不变闭合曲线,呈现准周期运动状态;对方程 ( ) 2进行 Lauo指数分析发 ypnv现,继续增大时, a系统将从准周期运动通过倍周期分岔进入混沌状态. 13在 .5
个正的 Lauo指数, ypnv表明系统进入混沌状态; a> .6时,当 14系统同时有两个正的 Lau o指数, ypnv说明系统进入超混沌状态. ei命令建立自定义函数名为 c ci cat .1然后利用 dbg用 d t i u— hoc 1 r t i 1, eu命令直接运
行,混沌图形如图 2所示.
程序如下: f co fg u tn[,]=cci cat; ni i u _hoc r t i
i=1 a=1 4 b=1 1 C=0 1 d=0 2 ; . 5; .; .; .;
%参数赋值
%定义初值
f1= .82 ( ) 1 13; g 1= .96 () O51;
wh l<1 o i i O o e
%进入循环迭代 5 0 0次
fi 1= g i d ( ()^; (+ ) a ()一 g i) !: g i 1= b fi+ g i; (+ ) (一 ) () C ()
i=i+1 :
po(,b ltg ;
h l n; odo
po(,r l f; t
e nd
13离散时滞混沌系统 .
以上给出了一阶差分方程确定的混沌系统,下面讨论离散时滞系统的数值仿真.考虑生物动力学中的