数学模型7.docVIP

数学模型7.3蛛网图 function f=zwt(a,x0,m,n) x=0:0.01:1; y=a*x.*(1-x); plot(x,x,r,x,y,r); hold on; clear x,y; x(1)=x0; y(1)=a*x(1)*(1-x(1)); x(2)=y(1); if m2 plot([x(1),x(1),x(2)],[0,y(1),y(1)]); end for i=2:n y(i)=a*x(i)*(1-x(i)); x(i+1)=y(i); if im plot([x(i),x(i),x(i+1)],[y(i-1),y(i),y(i)]); end end hold off; 输入zwt(1.7,0.2,1,100) 输入zwt(2.6,0.2,1,100) 输入zwt(3.3,0.2,1,100) 输入zwt(3.45,0.2,1,100) 输入zwt(3.55,0.2,1,100) 输入zwt(3.57,0.2,1,100) 本课其它代码 clear b=3.2; x=0.2; n=100; for i=1:n x=b*x*(1-x);x=b*x*(1-x);y(i)=x;x1(i)=i; end subplot(1,2,1) plot(x1,y,-o) x1(1)=0.2; y1(1)=0; for i=1:n-1 x1(2*i)=x1(2*i-1); x1(2*i+1)=y(i); y1(2*i)=y(i); y1(2*i+1)=y1(2*i); end; x2=0:0.01:1; y2=x2; y3=b.*x2.*(1-x2);y3=b.*y3.*(1-y3); subplot(1,2,2) plot(x1,y1,k-,x2,y2,k-,x2,y3,k-) clear b=3.3; x=0.2; n=100; for i=1:n x=b*x*(1-x);x=b*x*(1-x);y(i)=x;x1(i)=i; end subplot(1,2,1) plot(x1,y,-o) x1(1)=0.2; y1(1)=0; for i=1:n-1 x1(2*i)=x1(2*i-1); x1(2*i+1)=y(i); y1(2*i)=y(i); y1(2*i+1)=y1(2*i); end; x2=0:0.01:1; y2=x2; y3=b.*x2.*(1-x2);y3=b.*y3.*(1-y3); subplot(1,2,2) plot(x1,y1,k-,x2,y2,k-,x2,y3,k-) clear r=3.2; % syms z s=solve(r*r*z*(1-z)*(1-r*z*(1-z))-z,z) sr=subs(s,r,r) ssort=sort(sr) x10(1)=ssort(1)+(ssort(2)-ssort(1))*rand(1); x10(2)=ssort(2)+(ssort(3)-ssort(2))*rand(1); x10(3)=ssort(3)+(ssort(4)-ssort(3))*rand(1); x10(4)=ssort(4)+(1-ssort(4))*rand(1); n=100; for k=1:4 x=x10(k); for i=1:n x=r*x*(1-x); x=r*x*(1-x);y(i,k)=x;x11(i)=i; end subplot(4,2,1+2*(k-1)) plot(x11,y(:,k),-o) x1(1)=x10(k); y1(1)=0; for i=1:n-1 x1(2*i+1)=y(i,k); x1(2*i)=x1(2*i-1); y1(2*i,k)=y(i,k); y1(2*i+1,k)=y1(2*i,k); end; x2=0:0.01:1; y2=x2; y3=r.*x2.*(1-x2);y3=r.*y3.*(1-y3); subplot(4,2,2+2*(k-1)) plot(x1,y1(:,k),k-,x2,y2,k-,x2,y3,k-) end %为了观测r对迭代格式x(k+1)=rx(k)(1-x(k))的影响,将区间(0,4]离散化. %对每个离散的r值进行迭代,忽略前50个迭代值,把点(r,x51),(r,x52),...,(r,x100) %显示在坐标平面上.这样形成的图称为Feigenbaum图,它反映了分叉与混沌的基本特性. clear;