正在加载图片...
画图3.1的程序为 8x(1)=0 x(+1)=a*x()*(1-x(1), plo(0:20,x) 画图3.2的程序为 a=[32,3.5,3.56443.8284]y=0.2*[1,1,1,1] end x(1,)=y x(计+1,)=a.*x(1:)*(1-x(1:) end belote(2,2,1)plot(020,(x(;1) subplot(2,2,2)plot(020,(x(:,2) subplot(2, 2, 3 ),plot(0: 20, (x( 3))) plot(2,2,4.plot(020,(X(:4) 从线性联结图上,看不出当 44 时, logistic迭代到底是几一周期的?用下面蛛网图就 可以看清楚。 2.蛛网图 先计算出 再做下列工作: (1)画曲线y=ax(1-x),x∈[0,和直线y=x; (2)出发点A(x10x 过A做竖直线交曲线于点B(x100x0o),过B做水平 线交直线于新的点A(x010100):过新A做竖直线 交曲线于新的点B(x0o1xo0),过新B做水平线交直 线于新的点A(x100,x00);重复多次。 画图34中第3个图的程序为 hold =0:0.05:1y=3.5644*X*(1-x) plot(x, y),plot([o, 1][0,1D a=35644x=0.2; =a*x*(1-x) end 3.5644*X.*(1-x) plot(Ix, x].x,y D),plot(lx, yD画 图 3.1 的程序为: a=3.8;x(1)=0.2; for i=1:20 x(i+1)=a*x(i)*(1-x(i)); end plot(0:20,x) 画 图 3.2 的程序为: a=[3.2,3.5,3.5644,3.8284];y =0.2*[1,1,1,1]; for i=1:10000 y=a.*y.*(1-y); end x(1,:)=y; for i=1:20 x(i+1,:)=a.*x(i,:).*(1-x(i,:)); end subplot(2,2,1),plot(0:20,(x(:,1))') subplot(2,2,2),plot(0:20,(x(:,2))') subplot(2,2,3),plot(0:20,(x(:,3))') subplot(2,2,4),plot(0:20,(x(:,4))') 从线性联结图上,看不出当 a=3.5644, a=3.8284 时,Logistic 迭代到底是几---周期的? 用下面蛛网图就 可以看清楚。 2.蛛网图 先计算出 10000 x ,再做下列工作: (1)画曲线 y = ax(1− x), x[0,1] 和直线 y = x ; (2) 出发点 ( , ) 10000 10000 A x x , 过 A 做竖直线交曲线于点 ( , ) 10000 10001 B x x ,过 B 做水平 线交直线于新的点 ( , ) 10001 10001 A x x ;再过新 A 做竖直线 交曲线于新的点 ( , ) 10001 10002 B x x ,过新 B 做水平线交直 线于新的点 ( , ) 10002 10002 A x x ;重复多次。 画 图 3.4 中第 3 个图 的程序为: hold on x=0:0.05:1;y=3.5644*x.*(1-x); plot(x,y),plot([0,1],[0,1]) a=3.5644;x=0.2; for i=1:10000 x=a*x*(1-x); end for i=1:20 y=3.5644*x.*(1-x); plot([x,x],[x,y]),plot([x,y],[y,y]) x=y;
<<向上翻页向下翻页>>
©2008-现在 cucdc.com 高等教育资讯网 版权所有