MATLAB实验报告李犁.docVIP

15 MATLAB数学实验报告 姓名：夏天琦 学号:1304120221 指导老师：易昆南 专业班级：统计1201 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 验证混沌的随机性 I评 分 1、 设计(实习》目的： 了解MATLAB在实际问题中的应用 通过实践加深对这门语言中M文件的了解 熟悉简单程序结构，如循环结构(for循环、while循环)选择结构 (if-else-if)、分支语句(switch-case-otherwise)。 2、 实验内容： 以f(x)=4x(l-x)为例，编程验证混沌不是随机的 3、 详细设计： clear clc; x=0:0. 01:1; n=length(x); 尸zeros (1，n); for i=l:101 y(i)=4*x(i)*(l-x(i)); end figure plot (x, y); 初值为0. 1 x=zeros (1, 100); X ⑴=.1; for n=l:99 x(n+l)=4*x(n)*(l-x(n)); end [a, b]=xcorr (x/ unbiased’)； plot (b, a) 初值为0. 5 x=zeros (1,100); x(l)=. 5; for n=l:99 x (n+1) =4*x (n) * (1-x (n)): end [a, b]=xcorr (x/ unbiased’)； plot (b, a) 4：实验结果： x 10 5：实验总结 通过mat lab的操作，可以清楚地验证混沌的随机性。 姓名：李犁 2013年11月11日 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 研究帐篷函数 评 分 1、设计（实习）目的： 画出帐篷函数的图形； 编程演示帐篷函数对初值的敏感性； 找出帐篷函数的周期点。 2、实验内容： 帐篷函数的定义为 T （x）=错误!未找到引用源。 利用MATLAB解决以下问题： （1） 画山帐篷函数的图形； （2） 编程演示帐篷函数对初值的敏感性； （3） 找出帐篷函数的周期点。 3.详细设计： clear clc; x=0:0.01:1; n=length(x); y=zeros(l, n); for i=l:101 if x(i)=0. 5 y⑴=8*x(i); else y(i)=8*(l-x(i)); end end figure plot (x, y); x二zeros (1, 100); x(l)=99. 0; for n二1:99 if x〈=0. 5 x(n+l)=8*x(n); else x(n+l)=8*(l-x(n)); end end [a, b]=xcorr(x, ’ unbiased’)； plot (b, a) x=zeros(l, 100); x(l) = 100.0; for n=l:99 if x〈=0.5 x (n+l)=8*x (n); else x(n+l)=8*(l-x(n)); end end [a, b]=xcorr (x/ unbiased’)； plot (b, a) 4:实验结果: x 10 100 -80 -60 -40 -20 0 20 40 60 80 100 5：实验总结 Mat lab是个强大的工具， 帮助我们研究帐篷函数的相关性质。 姓名：李犁 2013年11月11日  学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 利用方波生成元生成分形 评 分 1、设计（实习〉目的： 1、 了解分形几何的有关知识； 2、 利用MATLAB绘制分形图形。 2、实验内容： 用图9-3-1中的图形（图中 1每一小段 长度相同，且相邻的两段垂直）作为召 ￡成元，编程生成分形。 图9-3-1方波 生成元 3.详细设计： u=[0，l]; for k=l:4 m=u/4; uu=[m, l/4+m*(li)，m+l/4+li/4, 1/2+1i/4+m*(-1i), 1/2, 1/2+m*(_li)，m+l/2_li/4, m*(li)+3/4-1i/4, m+3/4]: subplot (2, 2, k): u=uu; plot (u) end 4:实验结果： 0.4 0.2 0 -0.2 -0.4 0 0.5 1 5：实验总结 体会到了 Mat lab的强大功能。 姓名：李犁 2013年11月11日 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 Koch曲线 评 分 1、设计(实习》目的： 1、了解分形几何的有关知识； 2、利用MATLAB绘制分形阁形。 2、实验内容： 修改Koch曲线算法中