分岔:迭代函数系统
<FONT face="Courier New">function Mars; <BR>% 标题: 分岔:迭代函数系统 <BR>% <BR>% <BR>% 发信人: Mars (混沌·分形·周期三), 信区: MATH <BR>% 标题: 分岔:迭代函数系统 <BR>% 发信站: 一网深情 (Sat Nov7 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=; % 简化算式 <BR>% Z=; % 原始式子 <BR> end <BR>end <BR>plot(Z,'.','markersize',2) <BR>axis equal </FONT><BR> 好东西 以前看过<BR>呵呵 好程序不理解迭代程序,麻烦楼主讲解一下
x=STARX;for k=1:101;
x=a*x*(1-x);
end
for k=1:201;
x=a*x*(1-x);
Z=; % 简化算式
% Z=; % 原始式子
end
程序中k是迭代次数吧,为什么循环迭代了两次? 看了类似论文的一些解释,开始的循环目的是为了使迭代稳定。
回复 6楼 zhly 的帖子
对,个人感觉是为了得到一个比较理想的初始条件。使计算的比较稳定7楼你好
能否帮忙看看这个话题http://forum.vibunion.com/forum/thread-70436-1-1.html
也是用到迭代
页:
[1]