声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2214|回复: 7

[共享资源] 分岔:迭代函数系统

[复制链接]
发表于 2005-7-12 21:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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>
回复
分享到:

使用道具 举报

发表于 2005-7-14 15:39 | 显示全部楼层
好东西
发表于 2005-9-21 16:34 | 显示全部楼层
以前看过<BR>呵呵
发表于 2005-11-27 02:22 | 显示全部楼层
好程序
发表于 2008-9-26 10:03 | 显示全部楼层

不理解迭代程序,麻烦楼主讲解一下

x=STARX;
for k=1:101;
x=a*x*(1-x);
end
for k=1:201;
x=a*x*(1-x);
Z=[Z,a-x*i]; % 简化算式
% Z=[Z,x0+a*K1-J1+(y0-x*K2)*i]; % 原始式子
end
程序中k是迭代次数吧,为什么循环迭代了两次?
发表于 2008-9-28 21:00 | 显示全部楼层
看了类似论文的一些解释,开始的循环目的是为了使迭代稳定。

评分

1

查看全部评分

发表于 2008-9-29 19:34 | 显示全部楼层

回复 6楼 zhly 的帖子

对,个人感觉是为了得到一个比较理想的初始条件。使计算的比较稳定

评分

1

查看全部评分

发表于 2008-10-15 16:42 | 显示全部楼层

7楼你好

能否帮忙看看这个话题
http://forum.vibunion.com/forum/thread-70436-1-1.html
也是用到迭代
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-23 02:19 , Processed in 0.060779 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表