船过河,分叉与混沌.doc

实验报告 课程 数学实验 系别 实验日期 2011年 5 月 2 日 班号 组别 交报告日期姓名 学号 同组者 学号 实验一：分叉与混沌 分析：先要大致找出出现分叉的点，然后精确的找。其中有两个变量，用for语句 可以实现。 实验程序： clc clear for a=0:0.3:4.2 x(1)=0.1; for k=1:300; y(k)=a*sin(pi*x(k)); x(k+1)=y(k); if k250 hold on subplot(1,2,1) plot(a,y(k),.) end end end %大致找出分叉点 grid c=input(c=) d=input(d=) %确定分叉点的具体位置 for a=c:0.001:d x(1)=0.1; for k=1:300; y(k)=a*sin(pi*x(k)); x(k+1)=y(k); if k250 hold on subplot(1,2,2) plot(a,y(k)) end end end %作图 3.实验结果 4.改进之处： ①在左图上标出相应的数字 实验二：船过河 1.分析：如果船可以过河， 如右图。无论AB,∠ABF如何变化 DE的斜率始终大于CE。那么通过斜率的比 教即可判断船是否可以过河。通过 pause 命令即可模拟船过河的过程。 实验程序 clear clc for l=0.1:0.01:25 for b=0.01:0.01:pi/2-0.01 c=[5*sin(b),l*sin(b)+5*cos(b)]; if (10-c(2))/(12-c(1))=-tan(b) l break end end if (10-c(2))/(12-c(1))=-tan(b) break end end a=input(a=) for b=0.01:0.01:pi/2-0.01 x1=[10,20,20,10]; y1=[12,12,20,20]; x2=[0,0,20,20]; y2=[0,-2,-2,0]; x3=[-2,-2,0,0]; y3=[20,-2,-2,20]; axis equal axis off hold on fill(x1,y1,r) fill(x2,y2,r) fill(x3,y3,r) x=[a*cos(b),a*cos(b)+5*sin(b),5*sin(b),0]; y=[0,5*cos(b),a*sin(b)+5*cos(b),a*sin(b)]; fill(x,y,y) hold off pause(0.2) end 3.实验结果： L= 4.有待改进： ①图像模拟过程的阴影需要改进