youbrave 发表于 2009-4-2 10:56

大家帮忙看看周纪卿书后面的 点映射程序?

点映射的基本公式为 x(n+1) = u*x(n) (1-x(n)/4),代码如下:
ax=326.667;
ay = 32;
h = 0.006;
y1=0.4;

title('一维点映射')
% axis off
for i = 1:200
    xx = h.*i + 2.8;
    xi = (xx - 2.8).*ax;
    for j = 1:160
      y2 = xx.*y1.*(1.0 - y1/4.0);
      if j<=120
            yj = y2.*ay;
            hold on
            plot(xi,yj)
            y1 = y2;
      end
    end
end
这个代码里面ax,ay又什么作用,为什么要用xj,yj乘这些数值,这些初始值都好奇怪,是怎么设置的,还有j为什么只算到160,而在120时判断。

[ 本帖最后由 youbrave 于 2009-4-2 11:10 编辑 ]
页: [1]
查看完整版本: 大家帮忙看看周纪卿书后面的 点映射程序?