马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
<FONT face="Courier New">function Mars; <BR>% 标 题: 分岔:迭代函数系统 <BR>% <BR>% <BR>% 发信人: Mars (混沌·分形·周期三), 信区: MATH <BR>% 标 题: 分岔:迭代函数系统 <BR>% 发信站: 一网深情 (Sat Nov 7 09:09:27 1998), 转信 <BR>% <BR>% /*本程序为一个简单的迭代函数系统y=a*x*(1-x) */ <BR>% /*该函数是一个单峰函数,在此条件下可以改为其*/ <BR>% /*他单峰函数。注意:参数a的变化范围可以调整,*/ <BR>% /*但上限小于4. */ <BR>% /*通过不同的a值,你可观察到二、四、八等分岔 */ <BR>% /*也可出现三分岔、五分岔等,中间模糊区域为混沌区*/ <BR>% <BR>% /*Mars写于1996*/ <BR>STARA=2; <BR>ENDA=4; <BR>TA= 0.01; <BR>L=620-((ENDA-STARA)/DA); <BR>STARX=0.2; <BR>K1=1200/ENDA; <BR><BR>K2=200; <BR>J1=STARA*K1; <BR>x0=20; <BR>y0=300; <BR>Z=[]; <BR>for a=STARA:TA:ENDA; <BR> x=STARX; <BR> for k=1:101; <BR> x=a*x*(1-x); <BR> end <BR> for k=1:201; <BR> x=a*x*(1-x); <BR> Z=[Z,a-x*i]; % 简化算式 <BR>% Z=[Z,x0+a*K1-J1+(y0-x*K2)*i]; % 原始式子 <BR> end <BR>end <BR>plot(Z,'.','markersize',2) <BR>axis equal </FONT><BR> |