MATLAB实验报告李犁.docx
MAT1.AB数学试验报告姓名：夏天琦学号：1304120221指导老师；易昆南专业班级：统计1201学号I130412072【ft|91203|姓名Kl|指导老.试目IImK为的丽性那分1、设计实习)目的：1. 了解VAT1.AB在实际问题中的应用2. 通过实践加深对这门语言中M文件的J'解3. 熟识简洁程序结构，如循环结构(for循环、while循环)选择结构(if-else-if),分支语句(switch-case-otherwise).2、试验内容：以f(x)=4x(lx)为例，编程脸证泥沌不是随机的3.具体设计，clearclc;x=0:0.01:1;n=length(x):y=zeros(1.11>:for1=1:101y(i)=4*x(i)(l-x(i);endfigureplot(x,y)：初值为0.IK=ZeroS(1,100);x(l>=.l:forn=l:99x(n÷l)M*x(n)*(l-x(n);enda,b=xcorr(x,unbiased')：plot(bta)初值为0.5X=ZerOS(1.100);x(I)=.5:forn=l:99x(n÷l)=4*x(n)*(l-x(n)：enda,b=xcorr(x,unbiased')：plot(b,a)4,cm*>09x101411210-80-60-40-20204060801005»通过matlab的操作，可以清晰地验证混沌的随机性。姓名：李犁2013年11月“日学I1304120724跳计1203姓名指等老麻目探讨帐Mflt评分1、设计实习目的：I.画出帐篷函数的图形：2.编程演示怅篷函数对初值的坡JS性；3.找出帐逐函数的周期点.2、试内套，帐篷函数的定义为叫M3利用MAT1.AB解决以下问题：(1画H；忤逐函数的图形：(2)编程演示帐篷函数时初值的欷感性：(3找Hl帐篷函数的周期点.3.具体设计，clearclc：K=O:0.01:1;n=lcngth(x);y=zcros(l>n);for1=1:101ifx(i)<=0.5y(i)=8*x(i);elsey(i)-8*(lx(i);endendfircplot(x.y):X=ZerOS(1,100):x(l)=99.0:forn=l:99ifx<=0.5x(n+l)=8*x(n):elsex(nM)=*(l-x(n);endenda,b-xcorr(x,*unbiased*);plot(b.a)x=zeros(l>100);X二100.0;forn=l:99ifx<=0.5x(n+l)=8*x(n);elsex(n÷I)=8*(l-x(n);endenda,b=xcorr(xt*unbiased*);PIot(b,a)4tVIIM*t5»Matlab是个强大的工具，帮助我们探讨帐篷函数的相关性质。姓名Iwl2013年11月11日学号I1304120724三跳计1203姓名指导老府利用方波生成元生成分资评分1、设计（实习）目的：1.了解分形几何的有关学问I2、利用MAT1.AB绘制分形图形.2、试It内容：用图9-3-1中的图形（图U，堤一小段长度相同，且相邻的两段垂直）作为4:成元，端程生成分形，图9-*I方波生成元3.具体设计,u=0.1;fork=l:4m-u/4;uu=«.l4*n*(li).mG4÷li4.l2Hi4*me(-li),1/2.12*b*(-li),ml2li4,m*(li)*3/4-1i4.m÷34;subplot(2,2,k);u=uu;plot(u)end4试5.体妁TrIhtl的过大功能.姓笥w学号1304120724|缜计1203|娃卜亨寺指导老师AM目ItoCh曲线|哪一分1、设计(实习)目的】K了解分形几何的有关学向：2、利用WvnAB绘制分形图形.2、试收内容；修改KoCh曲城算法中板转T的特点由(1)四边形四个初始点动身注制分形，<2>利用下面的舞法籍制Mandelbwt集:3.具体设计，“绘制koch雪花图线P=00;100;新为初始两个点的坐标,第一列为X坐标,其次列为y坐标n=2:%n为结点数A=cospi3-Sin(Pi/3);Sin(Pi/3)cos(pi3):%旋转矩阵fork=l:4d=diff(p)3;%diff计尊相邻两个点的坐标之乱得到相知两点确定的向依%则d就计算出每个向破长度的三分之一,叮遨中将线段三等分对应m=4*n-3:%迭代公式q=p(kn-l.:):%以KI点为起点.前n-1个点的坐标为终点形成向It5:4:。,:)=p(2:n,:):、迭代后处于4k+l位置上的点的坐标为迭代前的相应坐标P(2:4:n.:)=q4d：斓向柒方法计算迭代后处于4k÷2位Bt上的点的坐标P(3:4:m,D=q+d+d*A'：铝H向疥方法计律迭代后处广北+3位置上的点的坐标p(4：4:«.J=q÷2*d:%ffl向阮方法计算迭代后处于4k位置上的点的坐标n=三;%迭代后新的结点数目end%绘制龙曲跳p=010;100fl-10；-100;0IOfAP为四边形四个顶点的坐标,其中第五个点与第一个点Ift合.以便于绘图“第一列为X坐标,灰次列为y坐标11=5:%n为结点数A=cos(-pi3)-sin(-pi3)jsin(-i3)COS(-pi3);%旋转矩阵期时针旋转60度fork15d=diff(p)/3;a=4*n-3;%迭代公式q=p(1：n-1,：)；p(5:4:m,:)=p(2:n.:);p(2:4:m.:)=qj;p(34m=q+2,d+<lA,;p(4:4:nv)=q+2'd;n=m;%绘IM花草树木曲殴p=00；1010:为初始两个点的坐标,第一列为X坐标,其次列为y坐标n=2:M为结点数A=coS(W3)-SwpW)SlMPb3)cos(pi);B=fcos(-p3)-siX-py3)isn(p)cos(-pW);a旋转矩阵A对应于第一次逆时针旋转60度,概转矩阵B时位于其次次顺时针旋转60度fork=1:4d=diff(p)3;dl=d(l:2:n.:);%取抵条线段对应的向lit11=5*n:%迭代公式q1=p(1:2：f>1.:)：p(10:10:tn,；H)(2:2:n,:)：P(W:)=p(l:2:n.:>%迭代后处于IOk与10k+1位皿上的点的坐标为迭代前的相应坐标p(2:10:m.:)=q1+d1;%用向弟方法il算迭代后处于Iok+2.1Ok3.1Ok抬位置上的点的坐标,都相同p(3:10:m.:)=p<2：10:m.:);p(4：10:ii.:)=ql+dl4dl*A，:%用向fit方法计算迭代后处于IOkx-Jtt上的点的坐标p(5:10:m,:)=p(2W:m.:);仅6:IOmpq1+2P1;%用向盘方法计算迭代后处于IOk蛇/Ok+7,1Ok吗位置上的点的坐标,郤相同p(7:10:m.:)=p(6:10:m.:);p<8:10:m.:)=q1*2-dbdrB,;p(9:10:m.:)=p(6:10:m.:);n=m;%迭代后新的结点数目end%绘制变换后的Koch曲线p=oO:IO0：P为初始两个点的坐标.第一列为X坐标,其次列为y坐标n=2:%n为结点数A=0-1:10：%旋转矩阵fork三14d=diff(p)3:%diff计算相邻两个点的坐标之空得到相邻两点确定的向垃%则d就计鸵出班个向量尺度的三分之1与庖中将双段三等分对应m=5*n-4;%迭代公式QP(l：n1,%以原点为起点.前>1个点的坐标为终点形成向嵬p(6:5:n.:)=p(2:n.:)：%迭代后处于5k+1位H上的点的坐标为迭代的的相应坐标p(2:5：B.:)=q4d：%用向量方法计算迭代后处于5k+2位置上的点的电标P(3:5:n,D=q+d+d*A'%用向"方法计算迭代后处于5k+3位置上的点的坐标p(4:5:n,:)q2*dd*A,:%用向量方法计"迭代后处于5k+4位度上的点的坐标p(5:5:m.:)=q+2*d:%用向It方法计算迭代后处于5k位置上的点的坐标n=m:%迭代后新的结点数目4tKm*i输入plot(p(:,l),p(:f2)2出鼻相锦两个点阖axis(O10010)到kod雪花图线输入PbI(NRMZ)ax«s(-1010-1010D得到ONot*Mtoolbarbuttons：databru*hn&Iinkedplots乂息PlaVVideo-10-86-4-2O246810输入plot(p(:,l),p(:,2)多绘出保相钻阚个点的连统axcs(010010)得到输入plot(p(:,l),p(=r2)%绘出将相邻两个点的连线axis(010010>得到5:试Ift总结这次作业的难度相对前几次又有了不同程度的变更，我始终暮在主动的0M老师的步伐，主动学习这笠1»作，争取更加尢自己的作业，有不对的地方还老师报正.姓名：李犁
