第五次实验报告(分形混沌)

南京信息工程大学实验(实习)报告

实验课程数学实验实验名称第五次实验实验日期指导老师

专业数学与应用数学年级姓名学号得分

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 实验目的:

了解有关分形和混沌的基本理论,能够用Mathematica 软件绘制出一些简单的分形和混沌图形。

实验内容:

用Mathematica 软件绘制一个分形的图形,图形类别自选。

Koch 雪花曲线

程序如下:

redokoch[ptlist_List] := Block[{tmp = {}, i, pnum = Length[ptlist]},

For[i =1, i

tmp = Join[tmp, {ptlist[[i]], ptlist[[i]]*2/3 + ptlist[[i + 1]]/3,

(ptlist[[i]] + ptlist[[i + 1]])/2 + {

ptlist[[i]][[2]] - ptlist[[i + 1]][[2]],

ptlist[[i + 1]][[1]] - ptlist[[i]][[1]]}*Sqrt[3]/6,

ptlist[[i]]/3 + ptlist[[i + 1]]*2/3, ptlist[[i + 1]]}]];

tmp] ; lnko01 = {{0, 0}, {1, 0}}; Show[

Graphics[Line[Nest[redokoch, lnko01, 5]], AspectRatio -> Sqrt[3]/6]]

图像如下:

令,其中,绘制出相应的IFS 吸引子图形,并取不同的s ,观察图形的变化。

程序如下:

s=0.5+0.5*I;

p1=0.5;f1[z_]:=s*z+1;

p2=0.5 ;f2[z_]:=s*z-1;

f[z_]:=Block[{tmp},tmp=Random[];

Which[tmp

Array[mu,{150,150}];

showIFS[z0_,shrage_List,divi_List,nmax_]:=

Block[{i,j,z=z0,a=divi[[1]],b=divi[[2]],temp1,temp2,mumax=0},

For[i=a,i>=1,i--,For[j=b,j>=1,j--,mu[i,j]=0]];

For[i=nmax,i>=1,i--,temp1=

Floor[a*(Re[z]-shrage[[1]][[1]])/(shrage[[2]][[1]]-shrage[[1]][[1]])]+1;temp2=Floor[b*(Im[z]-shrage[[1]][[2]])/(shrage[[2]][[2]]-shrage[[1]][[2]])]+1;mu[temp1,temp2]++;

z=f[z];];For[i=a,i>=1,i--,For[j=b,j>=1,j--,mumax=Max[mumax,mu[i,j]]]];

mu1=Table[GrayLevel[1-N[mu[j,i]]/mumax],{i,a},{j,b}];

Show[Graphics[RasterArray[mu1]]]] ;

showIFS[0+I 0,{{-0.1,-0.1},{1.1,1.1}},{150,150},10000]

图像如下:

用Mathematica 软件绘制一个混沌的图形,图形类别自选.

用二次迭代序列迭代函数f(x)=ax(1-x)

程序如下:

IterGeo[u_, x0_] := Module[{p1, p2, i,

pointlist = {},

var = x0, fvar = u*x0*(1 - x0)}, p1 = Plot[{u*x*(1 - x), x}, {x, 0, 1}, DisplayFunction -> Identity]; For[i = 1, i

i++, AppendTo[pointlist, {var, fvar}]; AppendTo[

pointlist, {fvar, fvar}]; var = fvar; fvar = u*var*(1 - var)]; p2 = ListPlot[pointlist,

PlotJoined -> True, DisplayFunction -> Identity];

Show[{p1, p2}, DisplayFunction -> $DisplayFunction]]; IterGeo[3.6, 0.8] 图像:

谈谈你所认识的分形和混沌。

答:分形:具有无限嵌套层次的精细结构,且在不同尺度下保持相似属性。 混沌:对初值敏感,而且不是随机的。

实验要求:

撰写实验报告

写出试验过程中所使用的Mathematica 程序或语句和计算结果

南京信息工程大学实验(实习)报告

实验课程数学实验实验名称第五次实验实验日期指导老师

专业数学与应用数学年级姓名学号得分

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 实验目的:

了解有关分形和混沌的基本理论,能够用Mathematica 软件绘制出一些简单的分形和混沌图形。

实验内容:

用Mathematica 软件绘制一个分形的图形,图形类别自选。

Koch 雪花曲线

程序如下:

redokoch[ptlist_List] := Block[{tmp = {}, i, pnum = Length[ptlist]},

For[i =1, i

tmp = Join[tmp, {ptlist[[i]], ptlist[[i]]*2/3 + ptlist[[i + 1]]/3,

(ptlist[[i]] + ptlist[[i + 1]])/2 + {

ptlist[[i]][[2]] - ptlist[[i + 1]][[2]],

ptlist[[i + 1]][[1]] - ptlist[[i]][[1]]}*Sqrt[3]/6,

ptlist[[i]]/3 + ptlist[[i + 1]]*2/3, ptlist[[i + 1]]}]];

tmp] ; lnko01 = {{0, 0}, {1, 0}}; Show[

Graphics[Line[Nest[redokoch, lnko01, 5]], AspectRatio -> Sqrt[3]/6]]

图像如下:

令,其中,绘制出相应的IFS 吸引子图形,并取不同的s ,观察图形的变化。

程序如下:

s=0.5+0.5*I;

p1=0.5;f1[z_]:=s*z+1;

p2=0.5 ;f2[z_]:=s*z-1;

f[z_]:=Block[{tmp},tmp=Random[];

Which[tmp

Array[mu,{150,150}];

showIFS[z0_,shrage_List,divi_List,nmax_]:=

Block[{i,j,z=z0,a=divi[[1]],b=divi[[2]],temp1,temp2,mumax=0},

For[i=a,i>=1,i--,For[j=b,j>=1,j--,mu[i,j]=0]];

For[i=nmax,i>=1,i--,temp1=

Floor[a*(Re[z]-shrage[[1]][[1]])/(shrage[[2]][[1]]-shrage[[1]][[1]])]+1;temp2=Floor[b*(Im[z]-shrage[[1]][[2]])/(shrage[[2]][[2]]-shrage[[1]][[2]])]+1;mu[temp1,temp2]++;

z=f[z];];For[i=a,i>=1,i--,For[j=b,j>=1,j--,mumax=Max[mumax,mu[i,j]]]];

mu1=Table[GrayLevel[1-N[mu[j,i]]/mumax],{i,a},{j,b}];

Show[Graphics[RasterArray[mu1]]]] ;

showIFS[0+I 0,{{-0.1,-0.1},{1.1,1.1}},{150,150},10000]

图像如下:

用Mathematica 软件绘制一个混沌的图形,图形类别自选.

用二次迭代序列迭代函数f(x)=ax(1-x)

程序如下:

IterGeo[u_, x0_] := Module[{p1, p2, i,

pointlist = {},

var = x0, fvar = u*x0*(1 - x0)}, p1 = Plot[{u*x*(1 - x), x}, {x, 0, 1}, DisplayFunction -> Identity]; For[i = 1, i

i++, AppendTo[pointlist, {var, fvar}]; AppendTo[

pointlist, {fvar, fvar}]; var = fvar; fvar = u*var*(1 - var)]; p2 = ListPlot[pointlist,

PlotJoined -> True, DisplayFunction -> Identity];

Show[{p1, p2}, DisplayFunction -> $DisplayFunction]]; IterGeo[3.6, 0.8] 图像:

谈谈你所认识的分形和混沌。

答:分形:具有无限嵌套层次的精细结构,且在不同尺度下保持相似属性。 混沌:对初值敏感,而且不是随机的。

实验要求:

撰写实验报告

写出试验过程中所使用的Mathematica 程序或语句和计算结果


    相关文章

    实验报告的书写

    动脉血压实验报告的书写 动脉血压测试实验目的: 1. 学习并掌握动脉血压测试的基本原理和基本方法. 2. 通过观察体位改变和极限运动对动脉血压的影响,掌握关于植物神经系统对动脉血压调节相关知识,观察不同个体心血管系统和植物神经系统功能状态和 ...

    关于Duffing方程的报告

    关于方程的报告 Duffing 1 Duffing方程的基本概念 混沌系统对微弱信号具有极强的敏感性同时对噪声具有极大的抑制能力,它的这种性质证明了混沌系统具有可应用于小信号检测的潜力,从检测过程中分析混沌运动发生的间歇性.Duffing方 ...

    混沌通信实验

    混沌通讯实验 实验一: 非线性电阻的伏安特性实验 1.实验目的:测绘非线性电阻的伏安特性曲线 2.实验装置:混沌通信实验仪. 3.实验对象:非线性电阻模块. 4.实验原理框图: 图1 非线性电阻伏安特性原理框图 5.实验方法: 第一步:在混 ...

    大学物理实验报告英文版--非线性与混沌

    Physical Lab Report:Nonlinear devices and chaos Writer: [Purpose] 1.To get familiar with the nonlinear devices and the p ...

    蔡氏二极管实现方法研究

    蔡氏二极管实现方法研究 摘要:该文首先基于一种常用的蔡氏二极管实现方法对蔡氏电路进行了仿真研究,分别使用Matlab 和PSpice 进行了仿真实验.接着,该文提出一种更加简化的蔡氏二极管实现方法,并基于这种新的蔡氏二极管实现方法对蔡氏电路 ...

    本科生科研训练课题指南

    本科生科研训练课题指南 北京师范大学物理系 2008.12 I.思路与方案 科学研究最根本的特征是创新.如何从学生进入大学之初由应试教育过渡到研究性和创新教育,是我国当前大学教育面临的严峻课题.在大学里,科研与教学有机地相结合,以科研促教学 ...

    自适应混沌粒子群优化算法

    计 算 机 工 程 第 37 卷 第15期 Computer Engineering V ol.37 No.15 文章编号:文章编号:1000-3428(2011)15-0128-03 ·人工智能及识别技术·人工智能及识别技术· 2011年 ...

    混沌与分形的奥妙

    数学烟云 自然的奥秘:混沌与分形 丁玖 "自然界的大书是以数学符号写的." -- 伽利略 (一) 引子 2009年的春天,新的一期<美国数学会会刊>(Notices of the American Mathe ...

    混沌神经网络研究进展与展望_董军

    第26卷第5期 1997年10月 信息与控制Infor mation and Contr ol V ol. 26, N o. 5 Oct. , 1997 混沌神经网络研究进展与展望 董 军 胡上序 (浙江大学 杭州 310027) 摘 要 ...