分形曲线及matlable算法

0 koch分形曲线在线演示

从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的两边代替，形成山丘形图形如下

在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边代替，再次形成新的图形如此迭代，形成 koch 分形曲线。

算法分析：由一条线段产生四条线段，故算法中由 n 条线段迭代一次后将产生 4n 条线段。算法针对每一条线段逐步进行，将计算新的三个点。第一个点位于线段三分之一处，第三个点位于线段三分之二处，第二个点以第一个点为轴心，将第一和第三个点形成的向量正向旋转 60 0 而得。正向旋转由正交矩阵

实现。

MATLAB 程序如下

clear

p=[0 0;10 0];n=1;

A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)];

for k=1:5

j=0;

for i=1:n

q1=p(i,:);

q2=p(i+1,:);

d=(q2-q1)/3;

j=j+1;r(j,:)=q1;

j=j+1;r(j,:)=q1+d;

j=j+1;r(j,:)=q1+d+d*A';

j=j+1;r(j,:)=q1+2*d;

end

n=4*n;clear p

p=[r;q2];

end

plot(p(:,1),p(:,2))

koch分形图片flash制作源代码

第一祯

die=4;

bi=1

color=0x000000

alpha=100

stop();

第二祯

_root.createEmptyMovieClip("koch",1); a=new Array(1025);

b=new Array(1025);

c=new Array(1025);

d=new Array(1025);

l=0;

n=1;

a[1]=100;

b[1]=200;

a[2]=500;

b[2]=200;

c[1]=100;

d[1]=200;

c[2]=500;

d[2]=200;

for(m=1;m<=die;m++)

{j=0;

with(_root.koch)

{

for(k=1;k<=n;k++)

{x1=c[k];

y1=d[k];

x2=c[k+1];

y2=d[k+1];

j=j+1;

a[j]=x1;

b[j]=y1;

j=j+1;

a[j]=x1+(x2-x1)/3;

b[j]=y1+(y2-y1)/3;

j=j+1;

a[j]=x1+(x2-x1)/3+((x2-x1)/3)*Math.cos(Math.PI/3)-((y2-y1)/3)*Math.sin(Math.PI/3); b[j]=y1+(y2-y1)/3+((x2-x1)/3)*Math.sin(Math.PI/3)+((y2-y1)/3)*Math.cos(Math.PI/3); j=j+1;

a[j]=x1+2*(x2-x1)/3;

b[j]=y1+2*(y2-y1)/3;

j=j+1;

a[j]=x2;

b[j]=y2;

l=j;

}

for(j=1;j

{if(a[j]==a[j+1]&&b[j]==b[j+1])

{ g=j;

for(;j

{a[j]=a[j+1];

b[j]=b[j+1];

}

j=g+1;

}

}

y=l;

for(f=1;f<=y;f++)

{c[f]=a[f];

d[f]=b[f];

}

}

n=4*n;

}

第三祯

//k1=1;

i=1;

_root.koch.onEnterFrame=function()

{

with(_root.koch)

{ //for(;i<=k1*10&&i<=y;)

{lineStyle(bi,color,alpha);

moveTo(c[i-1],450-d[i-1]);

lineTo(c[i],450-d[i]);

trace(i-1);

trace(d[i-1]);

i++;

}

//k1++;

if(i>=y){

delete _root.koch.onEnterFrame;

}

}

}

第四祯

_root.koch.clear();

stop();

1 矩形分形曲线1在线演示顶部

从一条直线段开始，将线段中间的三分之一部分用一个正方形的三边代替，形成几字形图形如下

在新的图形中，又将图中每一直线段中间的三分之一部分都用一个正方形的三边代替替，再次形成新的图形如此迭代，形成矩形分形曲线 1 。

算法分析：由一条线段产生五条线段，故算法中由 n 条线段迭代一次后将产生 5n 条线段。算法针对每一条线段逐步进行，将计算新的四个点。第一个点位于线段三分之一处，第五个点位于线段三分之二处，第二个点以第一个点为轴心，将第一和第三个点形成的向量正向旋转 90 0 而得。正向旋转由正交矩阵

实现。

MATLAB 程序如下

p=[0 0;10 0];n=2;

A=[0 -1;1 0];

for k=1:4

d=diff(p/3);m=5*n-4;

q=p(1:n-1,:);

p(6:5:m,:)=p(2:n,:);

p(2:5:m,:)=q+d;

p(3:5:m,:)=q+d+1*d*A';

p(4:5:m,:)=q+2*d+1*d*A';

p(5:5:m,:)=q+2*d;

end

plot(p(:,1),p(:,2))

flash制作原代码

1----------------------------------------

_root.createEmptyMovieClip("fenxing",1); a=new Array(1025);

b=new Array(1025);

c=new Array(1025);

d=new Array(1025);

e=new Array(1025);

f=new Array(1025);

a[1]=100;

b[1]=100;

a[2]=500;

b[2]=100;

n=2;

for(k=1;k<=4;k++)

{with(_root.fenxing)

{for(w=1;w

{c[w]=(a[w+1]-a[w])/3;

d[w]=(b[w+1]-b[w])/3;

}//微分d

m=5*n-4;

for(z=1;z<=n;z++)

{ e[z]=a[z];

f[z]=b[z];

}//差距q

h=2;

for(g=6;g<=m;g=g+5)

{if(h<=n)

a[g]=e[h];

b[g]=f[h];

h++

}

z=1;

w=1;

for(g=2;g<=m;g=g+5)

{a[g]=e[z]+c[w];

b[g]=f[z]+d[w];

z++;

}

z=1;

w=1;

for(g=3;g<=m;g=g+5)

{a[g]=e[z]+c[w]-1*d[w];

b[g]= f[z]+d[w]+1*c[w];

z++;

w++;

}

z=1;

w=1;

for(g=4;g<=m;g=g+5)

{a[g]=e[z]+2*c[w]-1*d[w];

b[g]= f[z]+2*d[w]+1*c[w];

z++;

w++;

}

z=1;

w=1;

for(g=5;g<=m;g=g+5)

{a[g]=e[z]+2*c[w];

b[g]=f[z]+2*d[w];

z++;

w++;

}

}

n=m;

}

2-----------------------------------------

i=2;

k=1;

_root.fenxing.onEnterFrame=function() { with(_root.fenxing)

{for(;i<=k*10&&i<=m;)

{lineStyle(1,0x000000,100);

moveTo(a[i-1],400-b[i-1]);

lineTo(a[i],400-b[i]);

trace(i);

trace(a[i]);

trace(b[i]);

i++;

}

k++;

if(i==m){

delete _root.fenxing.onEnterFrame;

}}

}

stop();

3----------------------------------------

_root.fenxing.clear();

stop();

2 矩形分形曲线2在线演示顶部原理同矩形分形曲线1。

matlab程序和图形如下

p=[0 0;10 0];n=2;

A=[0 -1;1 0];

for k=1:4

d=diff(p/3);m=5*n-4;

q=p(1:n-1,:);

p(6:5:m,:)=p(2:n,:);

p(2:5:m,:)=q+d;

p(3:5:m,:)=q+d+.7*d*A';

p(4:5:m,:)=q+2*d+.7*d*A';

p(5:5:m,:)=q+2*d;

n=m;

end

plot(p(:,1),p(:,2))

flash制作原代码

1------------------------------

_root.createEmptyMovieClip("fenxing",1); a=new Array(1025);

b=new Array(1025);

c=new Array(1025);

d=new Array(1025);

e=new Array(1025);

f=new Array(1025);

a[1]=100;

b[1]=100;

a[2]=500;

b[2]=100;

n=2;

for(k=1;k<=4;k++)

{with(_root.fenxing)

{for(w=1;w

{c[w]=(a[w+1]-a[w])/3;

d[w]=(b[w+1]-b[w])/3;

}//微分d

m=5*n-4;

for(z=1;z<=n;z++)

{ e[z]=a[z];

f[z]=b[z];

}//差距q

h=2;

for(g=6;g<=m;g=g+5)

{if(h<=n)

a[g]=e[h];

b[g]=f[h];

h++

}

z=1;

w=1;

for(g=2;g<=m;g=g+5)

{a[g]=e[z]+c[w];

b[g]=f[z]+d[w];

z++;

w++;

}

z=1;

w=1;

for(g=3;g<=m;g=g+5)

{a[g]=e[z]+c[w]-0.7*d[w];

b[g]= f[z]+d[w]+0.7*c[w];

z++;

w++;

}

z=1;

w=1;

for(g=4;g<=m;g=g+5)

{a[g]=e[z]+2*c[w]-0.7*d[w];

b[g]= f[z]+2*d[w]+0.7*c[w];

z++;

w++;

}

z=1;

w=1;

for(g=5;g<=m;g=g+5)

{a[g]=e[z]+2*c[w];

b[g]=f[z]+2*d[w];

z++;

w++;

}

}

n=m;

}

2------------------------------------- i=2;

k=1;

_root.fenxing.onEnterFrame=function() { with(_root.fenxing)

{for(;i<=k*10&&i<=m;)

{ lineStyle(1,0x000000,100);

moveTo(a[i-1],400-b[i-1]);

lineTo(a[i],400-b[i]);

trace(i);

trace(a[i]);

trace(b[i]);

i++;

}k++;

if(i==m+1){

delete _root.fenxing.onEnterFrame;

}}

}

stop();

3---------------------------------------

_root.fenxing.clear();

stop();

3 tree分形曲线1在线演示顶部

树枝形的分形曲线，由一条线段开始，将线段三分之一点处添加一线段作为树枝，再将线段三分之二点处添加线段作为另一树枝。

重复进行下去，可产生分形曲线效果。参考 koch 分开曲线算法，绘制如下图形。

matlab程序如下：

p=[0 0;10 10];n=2;

line(p(:,1),p(:,2))

A=[cos(pi/6) -sin(pi/6);sin(pi/6) cos(pi/6)];

for k=1:4

for j=1:2:n

p1=p(j,:);p2=p(j+1,:); d=(p2-p1)/3;

r(i,:)=p1;i=i+1;

q1=p1+d;

r(i,:)=q1;i=i+1;

r(i,:)=q1;i=i+1;

q2=q1+d*A';

r(i,:)=q2;i=i+1;

r(i,:)=q1;i=i+1;

q3=p1+2*d;

r(i,:)=q3;i=i+1;

r(i,:)=q3;i=i+1;

q4=q3+d*A;

r(i,:)=q4;i=i+1;

r(i,:)=q3;i=i+1;

r(i,:)=p2;i=i+1;

xy=[q1;q2];

line(xy(:,1),xy(:,2)); xy=[q3;q4];

line(xy(:,1),xy(:,2)); end

n=5*n;

end

flash程序如下：

1-----------------------

_root.createEmptyMovieClip("tree",1); a=new Array(1025);

b=new Array(1025);

c=new Array(1025);

d=new Array(1025);

e=new Array(1025);

f=new Array(1025);

g=new Array(1025);

h=new Array(1025);

cc=new Array(1025);

dd=new Array(1025);

a[1]=100;

a[2]=300;

b[1]=100;

b[2]=300;

k=1;

with(_root.tree)

{lineStyle(1,0x000000,100); moveTo(a[1],400-b[1]);

lineTo(a[2],400-b[2]);

}

n=2;

_root.tree.onEnterFrame=function() {with(_root.tree)

{ i=1;

for(j=1;j<=n;j=j+2)

{a[1]=a[j];

b[1]=b[j];

a[2]=a[j+1];

b[2]=b[j+1];

//p1&p2

cc[1]=(a[2]-a[1])/3;

dd[1]=(b[2]-b[1])/3;//d

c[i]=a[1];

d[i]=b[1];//r&p1

i=i+1;

e[1]=a[1]+cc[1];

f[1]=b[1]+dd[1];//q1

c[i]=e[1];

d[i]=f[1];//r&q1

i=i+1;

c[i]=e[1];

d[i]=f[1];//r&q1

i=i+1

e[2]=e[1]+cc[1]*Math.cos(Math.PI/6)-dd[1]*Math.sin(Math.PI/6) f[2]=f[1]+cc[1]*Math.sin(Math.PI/6)+dd[1]*Math.cos(Math.PI/6) //A

c[i]=e[2];

d[i]=f[2];

i=i+1;

c[i]=e[1];

d[i]=f[1];

i=i+1;

e[3]=a[1]+2*cc[1];

f[3]=b[1]+2*dd[1];//q3

c[i]=e[3];

d[i]=f[3];

i=i+1;

c[i]=e[3];

d[i]=f[3];

i=i+1;

e[4]=e[3]+cc[1]*Math.cos(Math.PI/6)+dd[1]*Math.sin(Math.PI/6) f[4]=f[3]-cc[1]*Math.sin(Math.PI/6)+dd[1]*Math.cos(Math.PI/6) //q4

c[i]=e[4];

d[i]=f[4];

i=i+1;

c[i]=e[3];

d[i]=f[3];

i=i+1;

c[i]=a[2];

d[i]=b[2];

i=i+1;// xy down

g[1]=e[1];

h[1]=f[1];

g[2]=e[2];

h[2]=f[2];

trace(e[1]);

trace(e[2]);

lineStyle(1,0x000000,100);

moveTo(g[1],400-h[1]);

lineTo(g[2],400-h[2]);

g[3]=e[3];

h[3]=f[3];

g[4]=e[4];

h[4]=f[4];

lineStyle(1,0xffffff,100);

moveTo(g[3],400-h[3]);

lineTo(g[4],400-h[4]);

}

for(m=1;m<=i;m++)

{a[m]=c[m];

b[m]=d[m]

}

n=5*n;

if(k>=4){

delete _root.tree.onEnterFrame;}

k++;

}

}

stop();

2-------------------------

_root.tree.clear();

stop();

4 tree分形曲线2在线演示顶部原理同矩形分形曲线1。

matlab程序如下：

clear

p=[0 0;0 15];

b=pi/8;e=1.5;

a=[cos(b) -sin(b);sin(b) cos(b)];

n=2;

for k=1:4

i=1;

for j=1:2:n

p1=p(j,:);

p2=p(j+1,:);

d=(p2-p1)/3;

new(i,:)=p1;i=i+1;

q1=p1+d;

new(i,:)=q1;i=i+1;

new(i,:)=q1;i=i+1;

new(i,:)=p1+d+e*d*a';i=i+1;

new(i,:)=q1;i=i+1;

new(i,:)=p1+d+e*d*a;i=i+1;

new(i,:)=p1+d;i=i+1;

q2=p1+2*d;

new(i,:)=q2;i=i+1;

new(i,:)=q2;i=i+1;

new(i,:)=p1+2*d+d*a;i=i+1;

new(i,:)=q2;i=i+1;

new(i,:)=p1+2*d+d*a';i=i+1;

new(i,:)=p1+2*d;i=i+1;

new(i,:)=p2;i=i+1;

end

n=n*7;p=new;

end

for m=1:2:n

line([p(m,1) p(m+1,1)],[p(m,2) p(m+1,2)]) end

flash代码如下：

1---------------------------

_root.createEmptyMovieClip("dtree",1);

a=new Array(1025);

b=new Array(1025);

c=new Array(1025);

d=new Array(1025);

e=new Array(1025);

f=new Array(1025);

g=new Array(1025);

h=new Array(1025);

cc=new Array(1025);

dd=new Array(1025);

a[1]=300;

a[2]=300;

b[1]=300;

b[2]=600;

bb=Math.PI/8;

ee=1.5;

n=2;

for(k=1;k<=4;k++)

{i=1;

for(j=1;j<=n;j=j+2)

{a[1]=a[j];

b[1]=b[j];

a[2]=a[j+1];

b[2]=b[j+1];//p1&p2

cc[1]=(a[2]-a[1])/3;

dd[1]=(b[2]-b[1])/3;//d

c[i]=a[1];

d[i]=b[1];

i=i+1;//new

e[1]=a[1]+cc[1];

f[1]=b[1]+dd[1];//q1

c[i]=e[1];

d[i]=f[1];

i=i+1;

c[i]=e[1];

d[i]=f[1];

//trace(c[i])

i=i+1;

c[i]=a[1]+cc[1]+ee*(cc[1]*Math.cos(bb)-dd[1]*Math.sin(bb)); d[i]=b[1]+dd[1]+ee*(cc[1]*Math.sin(bb)+dd[1]*Math.cos(bb)); //trace(Math.cos(bb));

c[i]=e[1];

d[i]=f[1];

i=i+1;

c[i]=a[1]+cc[1]+ee*(cc[1]*Math.cos(bb)+dd[1]*Math.sin(bb)); d[i]=b[1]+dd[1]+ee*(-cc[1]*Math.sin(bb)+dd[1]*Math.cos(bb)); i=i+1;

c[i]=a[1]+cc[1];

d[i]=b[1]+dd[1];

i=i+1;

e[2]=a[1]+2*cc[1];

f[2]=b[1]+2*dd[1];

c[i]=e[2];

d[i]=f[2];

i=i+1;

c[i]=e[2];

d[i]=f[2];

i=i+1;

c[i]=a[1]+2*cc[1]+cc[1]*Math.cos(bb)+dd[1]*Math.sin(bb);

d[i]=b[1]+2*dd[1]-cc[1]*Math.sin(bb)+dd[1]*Math.cos(bb);

i=i+1;

c[i]=e[2];

d[i]=f[2];

i=i+1;

c[i]=a[1]+2*cc[1]+cc[1]*Math.cos(bb)-dd[1]*Math.sin(bb);

d[i]=b[1]+2*dd[1]+cc[1]*Math.sin(bb)+dd[1]*Math.cos(bb);

i=i+1;

c[i]=a[1]+2*cc[1];

d[i]=b[1]+2*dd[1];

i=i+1;

c[i]=a[2];

d[i]=b[2];

i=i+1;

}

for(m=1;m

{a[m]=c[m];

b[m]=d[m]

//trace(a[m]);

}

n=n*7;

}

2-------------------------------

m=1;

i=1;

_root.dtree.onEnterFrame=function()

{ with(_root.dtree)

{ for(;i<20*k&&i

{lineStyle(1,0x000000,100);

moveTo(a[m],650-b[m]);

m=m+2;

lineTo(a[m],650-b[m]);

i++;

}

k++;

}

if(m>=n-3){

delete _root.dtree.onEnterFrame;

}

}

stop();

3---------------------------------------

_root.dtree.clear();

stop();

顶部

5 洛伦xy曲线在线演示顶部

flash代码如下：

1------------------------------

//初始化并获得到要画轨迹的离散坐标

_root.createEmptyMovieClip("luolun",1);

a=new Array();//记录在x轴上的投影点

b=new Array();//记录在y轴上的投影点

c=new Array();//记录在z轴上的投影点

h=0.01;//时间变换的幅度

d=10;

e=28

f=8/3;

a[0]=0.1;//初始数值

b[0]=0;//初始数值

c[0]=0;//初始数值

for(i=0;i<10000;i++)

{a[i+1]=a[i]+h*d*(b[i]-a[i]);

b[i+1]=b[i]+h*(a[i]*(e-c[i])-b[i]);

c[i+1]=c[i]+h*(a[i]*b[i]-f*c[i]);

stop();

2------------------------------

//画图

i=100;

_root.luolun.onEnterFrame=function()

{with(_root.luolun)

{ i++;

lineStyle(0.1,random(0xffffff),100);

moveTo(300+5*a[i-1],300-5*c[i-1]);

lineTo(300+5*a[i],300-5*c[i]);

if(i>=10000){

delete _root.luolun.onEnterFrame;

}

}

}

stop();

3------------------------------

//清除函数所画的线条

_root.luolun.clear();

stop();

6-----10的原理同5，所以在此就不加以介绍，只提供了在线演示

6 洛伦xz曲线在线演示顶部

7 洛伦yz曲线在线演示顶部

8 洛伦yx曲线在线演示顶部

9 洛伦zx曲线在线演示顶部

10 洛伦zx曲线在线演示顶部

11 抛射曲线包洛在线演示顶部

如果一条曲线上每一点都与曲线簇中某一曲线相切，则称这条曲线为该曲线簇的包络。当曲线族的方程

基于分形几何的分形图绘制与分析

基于分形几何的分形图绘制与分析 摘要：基于分形几何的分形图绘制方法源于l系统、迭代函数系统ifs、复动力系统等。在运用分形原理及算法编程绘制多种分形图的基础上，重点对ifs参数进行实验分析，ifs吸引集实现了对原图形的几何变换。分形图的演变具有渐变性。 关键词：分形几何迭代函数系统分形图绘制渐变 1 分形几何学 现代数学的一个新的分支——，它是由美籍法国数学家曼德勃罗（b.b.mandelbrot）1973年在法兰西学院讲课时，首次提出了分形几何的设想。分形（fractal）一词，是曼德勃罗创造出来的，其原意具有不规则、支离破碎等意义，分形几何学是一门以非规则几何形态为研究对象的几何学。由于不规则现象在自然界是普遍存在的，因此分形几何又称为描述大自然的几何学。分形几何的诞生无论是在理论上还是在实践上都具有重要价值。 2 分形的定义 目前分形还没有最终的科学定义，曼德勃罗曾经为分形下过两个定义： （1）分形是hausdorff-besicovitch维数严格大于拓扑维数的集合。因为它把许多hausdorff维数是整数的分形集合排除在外，例如，经典分形集合peano曲线分形维数 （2）局部与整体以某种方式自相似的形，称为分形。 然而，经过理论和应用的检验，人们发现这两个定义很难包括分形

如此丰富的内容。实际上，对于什么是分形，到目前为止还不能给出一个确切的定义，正如生物学中对“生命”也没有严格明确的定义一样，人们通常是列出生命体的一系列特征来加以说明。对分形的定义也可同样的处理。 (ⅰ) 分形集合在任意小尺度下，它总有复杂的细节，或者说它具有精细的结构。 (ⅱ) 分形集合是非常不规则的，用传统的几何语言无法来描述它的局部和整体，它既不是满足某些条件的点的轨迹，也不是某些简单方程的解集。 (ⅲ) 分形集具有某种自相似形式，可能是近似的自相似或者统计的自相似。 (ⅳ) 以某种方式定义的分形集合的“分形维数”，严格大于它相应的拓扑维数。 (ⅴ) 在大多数令人感兴趣的情形下，分形集合是以非常简单的递归的方法产生的。 3 分形研究的对象 几何学的研究对象是物体的形状，在自然界中，许多物体的形状是极不规则的，例如：弯弯曲曲的海岸线，起伏不平的山脉，变化无偿的浮云，以及令人眼花缭乱的满天繁星，等等。这些物体的形状有着共同的特点，就是极不规则，极不光滑。但是，所有的经典几何学都是以规则而光滑的形状为其研究对象的，例如：初等平面几何的主要研究对象是直线与圆；平面解析几何的主要研究对象是一

分形维数算法

分形维数算法． 分形维数算法 分形包括规则分形和无规则分形两种。规则分形是指可以由简单的迭代或者是按一定规律所生成的分形，如Cantor集，Koch曲线，Sierpinski海绵等。这些分形图形具有严格的自相似性。无规则分形是指不光滑的，随机生成的分形，

如蜿蜒曲折的海岸线，变换无穷的布朗运动轨迹等。这类曲线的自相似性是近 似的或统计意义上的，这种自相似性只存于标度不变区域。 对于规则分形，其自相似性、标度不变性理论上是无限的（观测尺度可以趋于无限小）。不管我们怎样缩小（或放大）尺度（标度）去观察图形，其组成部分和原来的图形没有区别，也就是说它具有无限的膨胀和收缩对称性。因些对于这类分形，其计算方法比较简单，可以用缩小测量尺度的或者不断放大图形而得到。分形维数 D=lnN(λ)/ln(1/λ) （2-20） 如Cantor集，分数维D=ln2/ln3=0.631；Koch曲线分数维 D=ln4/ln3=1.262; Sierpinski海绵分数维D=ln20/ln3=2.777。 对于不规则分形，它只具有统计意义下的自相似性。不规则分形种类繁多，它可以是离散的点集、粗糙曲线、多枝权的二维图形、粗糙曲面、以至三维的[26]。点 集和多枝权的三维图形，下面介绍一些常用的测定方法（1）尺码法 用某个选定尺码沿曲线以分规方式测量，保持尺码分规两端的落点始终在曲线上。不断改变尺码λ，得到一系列长度N（λ），λ越小、N越大。如果作lnN～lnλ图后得到斜率为负的直线，这表明存在如下的幂函数关系

-D（2-21） N～λ上式也就是Mandelbrot在《分形：形状、机遇与维数》专著中引用的Richardson公式。Richardson是根据挪威、澳大利亚、南非、德国、不列颠西部、葡萄牙的海岸线丈量结果得出此公式的，使用的测量长度单位一般在1公里到4公里之间。海岸线绝对长度L被表示为： 1-D（2-22）L=Nλ～λ 他得到挪威东南部海岸线的分维D≈1.52，而不列颠西部海岸线的分维D≈[27]。。这说明挪威的海岸线更曲折一些1.3． ）小岛法（2面积如果粗糙曲线都是封闭的，例如海洋中的许多小岛，就可以利用周长-关系求分维，因此这个方法又被称为小岛法。则与λ的而面积A对于规则图形的周长与测量单位尺寸λ的一次方成正比， 二次方成正比。通常我们可以把它们写成一个简单的比例关系：1/2 (2-23) AP∝对于二维空间内的不规则分形的周长和面积的关系显然更复杂一些，提出，应该用分形周长曲线来代替原来的光滑周长，从而给出了下Mandelbrot 述关系式：21/??D??1/1/D2）（2-24)]?(?)]?[a?AP[(?)][??a(1?D)/DA(?00的P）式），使1（周长光滑时D=1，上式转化成为（2.23这里的分维D大于??的数1变化减缓，a是和岛的形状有关的常数，为小于是测量尺寸，一般取0/D）（1-D??减小而增大。作随测

几个分形的matlab实现

几个分形得matlab实现 摘要:给出几个分形得实例,并用matlab编程实现方便更好得理解分形,欣赏其带来得数学美感 关键字:Koｃh曲线实验图像 一、问题描述: 从一条直线段开始,将线段中间得三分之一部分用一个等边三角形得两边代替,形成山丘形图形如下 ?图1 在新得图形中,又将图中每一直线段中间得三分之一部分都用一个等边三角形得两条边代替,再次形成新得图形如此迭代,形成Kocｈ分形曲线。 二、算法分析: 考虑由直线段(2个点)产生第一个图形(５个点)得过程。图1中,设与分别为原始直线段得两个端点,现需要在直线段得中间依次插入三个点,,。显然位于线段三分之一处,位于线段三分 之二处,点得位置可瞧成就是由点以点为轴心,逆时针旋转600而得。旋转由正交矩阵 实现。 算法根据初始数据(与点得坐标),产生图1中5个结点得坐标、结点得坐标数组形成一个矩阵,矩阵得第一行为得坐标,第二行为得坐标……,第五行为得坐标。矩阵得第一列元素分别为5个结点得坐标,第二列元素分别为5个结点得坐标。 进一步考虑Kｏｃh曲线形成过程中结点数目得变化规律。设第次迭代产生得结点数为,第次迭代产生得结点数为,则与中间得递推关系为。 三、实验程序及注释: p=［0 0;10 0］; %P为初始两个点得坐标,第一列为ｘ坐标,第二列为y坐标 n=2; %n为结点数 A=[cos(pi/３) —ｓin(pi/3);siｎ(ｐi/３) cos(pi/３)]; ％旋转矩阵 ｆor k＝1:4 ｄ=diｆf(p)/3; ％ｄｉff计算相邻两个点得坐标之差,得到相邻两点确定得向量 %则d就计算出每个向量长度得三分之一,与题中将线段三等分对应 m=4＊n-3; ％迭代公式 q=p(1:ｎ—1,:); %以原点为起点,前n—１个点得坐标为终点形成向量 p(5:4:m,:)=ｐ(2:n,:); %迭代后处于４k+１位置上得点得坐标为迭代前得相应坐标 p(2:4:m,:)=q+d; %用向量方法计算迭代后处于４k+２位置上得点得坐标 p(3:4:m,:)=ｑ+d+ｄ＊A'; %用向量方法计算迭代后处于4ｋ＋3位置上得点得坐标 p(4:４:m,:)＝ｑ＋2*d; %用向量方法计算迭代后处于4k位置上得点得坐标 n=m; %迭代后新得结点数目 end pｌot(ｐ(:,1),p(:,２)) ％绘出每相邻两个点得连线 aｘis([0 10 ０ 10]) 四、实验数据记录: 由第三部分得程序,可得到如下得Ｋoch分形曲线:

经典的分形算法 (1)

经典的分形算法 小宇宙2012-08-11 17:46:33 小宇宙 被誉为大自然的几何学的分形（Fractal）理论，是现代数学的一个新分支，但其本质却是一种新的世界观和方法论。它与动力系统的混沌理论交叉结合，相辅相成。它承认世界的局部可能在一定条件下，在某一方面（形态，结构，信息，功能，时间，能量等）表现出与整体的相似性，它承认空间维数的变化既可以是离散的也可以是连续的，因而拓展了视野。 分形几何的概念是美籍法国数学家曼德布罗（B.B.Mandelbrot）1975年首先提出的，但最早的工作可追朔到1875年，德国数学家维尔斯特拉斯（K.Weierestrass）构造了处处连续但处处不可微的函数，集合论创始人康托（G.Cantor，德国数学家）构造了有许多奇异性质的三分康托集。1890年，意大利数学家皮亚诺（G.Peano）构造了填充空间的曲线。1904年，瑞典数学家科赫（H.von Koch）设计出类似雪花和岛屿边缘的一类曲线。1915年，波兰数学家谢尔宾斯基（W.Sierpinski）设计了象地毯和海绵一样的几何图形。这些都是为解决分析与拓朴学中的问题而提出的反例，但它们正是分形几何思想的源泉。1910年，德国数学家豪斯道夫（F.Hausdorff）开始了奇异集合性质与量的研究，提出分数维概念。1928年布利干（G.Bouligand）将闵可夫斯基容度应用于非整数维，由此能将螺线作很好的分类。1932年庞特里亚金（L.S.Pontryagin）等引入盒维数。1934年，贝塞考维奇（A.S.Besicovitch）更深刻地提示了豪斯道夫测度的性质和奇异集的分数维，他在豪斯道夫测度及其几何的研究领域中作出了主要贡献，从而产生了豪斯道夫-贝塞考维奇维数概念。以后，这一领域的研究工作没有引起更多人的注意，先驱们的工作只是作为分析与拓扑学教科书中的反例而流传开来。 真正令大众了解分形是从计算机的普及肇始，而一开始，分形图的计算机绘制也只是停留在二维平面，但这也足以使人们心驰神往。近来，一个分形体爱好者丹尼尔?怀特（英国一钢琴教师）提出一个大胆的方法，创造出令人称奇的3D分形影像，并将它们命名为芒德球（mandelbulb）。

分形维数算法

分形维数算法

分形维数算法 分形包括规则分形和无规则分形两种。规则分形是指可以由简单的迭代或者是按一定规律所生成的分形，如Cantor集，Koch曲线，Sierpinski海绵等。这些分形图形具有严格的自相似性。无规则分形是指不光滑的，随机生成的分形，如蜿蜒曲折的海岸线，变换无穷的布朗运动轨迹等。这类曲线的自相似性是近似的或统计意义上的，这种自相似性只存于标度不变区域。 对于规则分形，其自相似性、标度不变性理论上是无限的（观测尺度可以趋于无限小）。不管我们怎样缩小（或放大）尺度（标度）去观察图形，其组成部分和原来的图形没有区别，也就是说它具有无限的膨胀和收缩对称性。因些对于这类分形，其计算方法比较简单，可以用缩小测量尺度的或者不断放大图形而得到。分形维数 D=lnN(λ)/ln(1/λ) （2-20） 如Cantor集，分数维D=ln2/ln3=0.631；Koch曲线分数维D=ln4/ln3=1.262; Sierpinski海绵分数维D=ln20/ln3=2.777。 对于不规则分形，它只具有统计意义下的自相似性。不规则分形种类繁多，它可以是离散的点集、粗糙曲线、多枝权的二维图形、粗糙曲面、以至三维的点集和多枝权的三维图形，下面介绍一些常用的测定方法[26]。 （1）尺码法 用某个选定尺码沿曲线以分规方式测量，保持尺码分规两端的落点始终在曲线上。不断改变尺码λ，得到一系列长度N（λ），λ越小、N越大。如果作lnN～lnλ图后得到斜率为负的直线，这表明存在如下的幂函数关系 N～λ-D（2-21） 上式也就是Mandelbrot在《分形：形状、机遇与维数》专著中引用的Richardson公式。Richardson是根据挪威、澳大利亚、南非、德国、不列颠西部、葡萄牙的海岸线丈量结果得出此公式的，使用的测量长度单位一般在1公里到4公里之间。海岸线绝对长度L被表示为： L=Nλ～λ1-D（2-22） 他得到挪威东南部海岸线的分维D≈1.52，而不列颠西部海岸线的分维D≈1.3。这说明挪威的海岸线更曲折一些[27]。

Koch分形雪花图的面积计算

K o c h分形雪花图的面 积计算 内部编号：（YUUT-TBBY-MMUT-URRUY-UOOY-DBUYI-0128）

Koch分形雪花图的面积计算 一、问题叙述 分形几何图形最基本的特征是自相似性，这种自相似性是指局部与整体在形态、功能、信息、时间、空间等方面具有统计意义上的相似。在具有自相似性的图形中，图形局部只是整体的缩影，而整体图形则是局部的放大。而本文我们要分析的是Koch分形雪花图，包含以下三个问题： 1.描述Koch分形雪花 2.证明Koch分形雪花图K n的边数为 3.求Koch分形雪花图的面积（数据），求 二、问题分析 在分析Koch分形雪花图之前，我们首先介绍Koch分形曲线。Koch分形曲线的绘制原理是：从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的两边代替，形成四条线段的折线，如图所示： 图对一条线段进行第一次Koch分形 然后，对形成的四条直线段的每一条的中间的三分之一部分用等边三角形的两边代替，形成十六条线段的折线。这种迭代继续进行下去可以形成Koch分形曲线。在迭代过程中，图形中的点数将越来越多，而曲线的最终显示细节的多少将取决于迭代次数和显示系统的分辨率。 设P1和P2分别是原始的两个端点，现在需要在直线段的中间依次插入点Q1， Q2,Q3以产生第一次迭代图形。显然，Q1位于P1右端直线段的三分之一处，Q3位于P1点右端直线段的三分之二处，而Q2点的位置可以看作由Q3绕Q1逆时针旋转60度而得到的，故可以处理Q Q13经过正交变换而得到Q Q12。算法如下： （1）Q1P1+P P Q P1+P P/3; （2-1）/3;32（2-1） ←←

《高频电子线路》课程设计指导书.doc

《高频电子线路》课程设计指导书 一、课程设计基本信息 核心课程名称（中文）高频电子线路核心课程名称（英文）High-frequency Electronic Circuits 课程设计名称高频电子线路课程设计 课程设计编号课程设计类型实物制作 相关辅助课程电路分析、电子线路（线性部分） 教材及实验指导书教材《电子线路（非线性部分）》,谢嘉奎，高等教育出版 课程设计时间：第五学期18 周 面向专业电子信息科学与技术 二、课程设计的目的 《高频电子线路》课程是电子信息专业继《电路理论》、《电子线路（线性部分）》之后必修的主要技术基础课，同时也是一门工程性和实践性都很强的课程。课程设计是在课程内容学习结束，学生基本掌握了该课程的基本理论和方法后，通过完成特定电子电路的设计、安装和调试，培养学生灵活运用所学理论知识分析、解决实际问题的能力，具有一定的独立进行资料查阅、电路方案设计及组织实验的能力。通过设计，加深对调幅的理解，学会电路的调整；进一步培养学生的动手能力 三、主要仪器设备 序号实验项目名称仪器设备名称仪器设备编号 1调幅收音机设计高频信号发生器、数字示波器、稳压电源 四、课程设计的内容与要求 1、内容：根据所学知识，设计一超外差调幅收音机电路，选择合适的元器件，进行安装和调试电路；应能接收正常广播，且接收的广播节目不少于3套° 序 号 名称目的方式场所要求

1调幅收音机设计加深对调幅的理解，学会 电路的调整；进一步培养 学生的动手能力 实物制作 通信学 院 2、要求 1设计电路图； 2供电电压：直流3V 3 接收频段：535kHz ~ 1605kHz; 4输出功率：P o> 1W。 5为满足偷出功率要求，采用两级放大电路； 6采用互补推挽功率放大器作为输出级。 五、考核与报告 考核内容：1实际操作：包括电路设计、安装、焊接及调试 2设计报告：包括原理、电路图、元器件的选择 成绩评定：实际操作和设计报告各占50%o 六、主要参考文献 1、《电子线路（非线性部分）》，谢嘉奎，高等教育出版社 2、《实用电子电路手册》，孙肖子，高等教育出版社 3、《电子技术技能训练》，张大彪，电子工业出版社七、课程设计报告 1、报告内容 目的、原理、电路图、安装注意事项、调试过程及结果。 2、版面格式 （1）A4纸打印，上、下、左、右边距为2. 5cm,段落间距0,行间距1. 5倍; （2）标题使用四号黑体、居中，正文使用小四号宋体; 一级标题：小四号黑体（如：1、2、3……）；

曲线型组合图形的面积计算方法

曲线型组合图形的面积计算方法姓名对于不规则图形面积的计算问题一般将它转化为若干基本规则图形的组合，分析整体与部分的和、差关系，问题便得到解决.常用的基本方法有： 一、相加法：这种方法是将不规则图形分解转化成几个基本规则图形，分别计 算它们的面积，然后相加求出整个图形的面积。例如下图中，要求整个图形的面积，只要先求出上面半圆的面积，再求出下面正方形的面积，然后把它们相加就可以了。 30厘米 二、相减法：这种方法是将所求的不规则图形的面积看成是若干个基本规则图 形的面积之差。例如下图中，若求阴影部分的面积，只需先求出正方形面积再减去里面圆的面积即可。 三、

四、 重新组合法：这种方法是将不规则图形拆开，根据具体情况和计算上的需要，重新组合成一个新的图形，设法求出这个新图形面积即可.例如，欲求下图中阴影部分面积，可以把它拆开使阴影部分分布在正方形的4个角处，这时采用相减法就可求出其面积了。 五、 辅助线法：这种方法是根据具体情况在图形中添一条或若干条辅助线，使不规则图形转化成若干个基本规则图形，然后再采用相加、相减法解决即可.如下图，求两个正方形中阴影部分的面积.此题虽然可以用相减法解决，但不如添加一条辅助线后用直接法作更简便. 六、 割补法：这种方法是把原图形的一部分切割下来补在图形中的另一部分使之成为基本规则图形，从而使问题得到解决.例如，如下图，欲求阴影部分的面积，只需把右边弓形切割下来补在左边，这样整个阴影部分面积恰是正方形面积的一半。 七、 平移法：这种方法是将图形中某一部分切割下来平行移动到一恰当位置，使之组合成一个新的基本规则图形，便于求出面积。例如下图，欲求阴影部分面积，可先沿中间切开把左边正方形内的阴影部分平行移到右边长方形内，这样整个阴影部分恰是一个长方形。 旋转法：这种方法是将图形中某一部分切割下来之后，使之沿某一点或某一轴旋转一定角度贴补在另一图形的一侧，从而组合成一个新的基本规则的图形，便于求出面积.例如，欲求下左图中阴影部分的面积，可将左半图形绕B 点逆时针方向旋转180°，使A 与C 重合，从而构成如下右图（2）的样子，此时阴影部分的面积可以看成半圆面积减去中间等腰直角三角形的面积。 九、 对称添补法：这种方法是作出原图形的对称图形，从而得到一个新的基本规则图形.原来图形面积就是这个新图形面积的一半。例如，欲求下图中阴影部分的面积，沿AB 在原图下方作关于AB 为对称轴的对称扇形ABD.弓形CBD 的面积的一半就是所求阴影部分的面积。 十、 重叠法：这种方法是将所求的图形看成是两个或两个以上图形的重叠部分，然后运用“容斥原理”（SA ∪B ＝SA ＋SB-SA ∩B ）解决。例如欲求下图中阴影部分的面积，可先求两个扇形面积的和，减去正方形面积，因为阴影部 分的面积恰好是两个扇形重叠的部分。 10厘米 6厘米 4厘米 20厘米 8厘米 10厘米 20厘米 30厘米 10厘米

分形插值算法和MATLAB实验

一，分形插值算法 ——分形图的递归算法1，分形的定义 分形（Fractal）一词，是法国人B.B.Mandelbrot 创造出来的，其原意包含了不规则、支离破碎等意思。Mandelbrot 基于对不规则的几何对象长期地、系统地研究，于1973 年提出了分维数和分形几何的设想。分形几何是一门以非规则几何形状为研究对象的几何学，用以描述自然界中普遍存在着的不规则对象。分形几何有其显明的特征，一是自相似性；分形作为一个数学集合, 其内部具有精细结构, 即在所有比例尺度上其组成部分应包含整体, 而且彼此是相似的。其定义有如下两种描述： 定义 1如果一个集合在欧式空间中的 Hausdorff 维数H D 恒大于其拓扑维数 r D ，则称该集合为分形集，简称分形。 定义 2组成部分以某种方式与整体相似的形体叫分形。 对于定义 1 的理解需要一定的数学基础，不仅要知道什么是Hausdorff 维数，而且要知道什么是拓扑维数，看起来很抽象，也不容易推广。定义 2 比较笼统的说明了自然界中的物质只要局部和局部或者局部和整体之间存在自相似性，那么这个物质就是分形。正是这一比较“模糊”的概念被人们普遍接受，同时也促进了分形的发展。 根据自相似性的程度，分形可分为有规分形和无规分形。有规分形是指具有严格的自相似的分形，比如，三分康托集，Koch 曲线。无规分形是指具有统计意义上的自相似性的分形，比如，曲折的海岸线，漂浮的云等。本文主要研究有规分形。

2. 分形图的递归算法 2.1 三分康托集 1883 年，德国数学家康托(G.Cantor)提出了如今广为人知的三分康托集。三分康托集是很容易构造的，然而，它却显示出许多最典型的分形特征。它是从单位区间出发，再由这个区间不断地去掉部分子区间的过程构造出来的(如图2.1)。 其详细构造过程是：第一步，把闭区间[0，1]平均分为三段，去掉中间的 1/3 部分段，则只剩下两个闭区间[0，1/3]和[2/3，1]。第二步，再将剩下的两个闭区间各自平均分为三段，同样去掉中间的区间段，这时剩下四段闭区间：[0，1/9]，[2/9，1/3]，[2/3，7/9]和[8/9，1]。第三步，重复删除每个小区间中间的 1/3 段。如此不断的分割下去，最后剩下的各个小区间段就构成了三分康托集。三分康托集的 Hausdorff 维数是0.6309。 图2.2 三分康托集的构造过程

分形之Julia集及其算法实现

成绩：课程名称：智能信息处理概论 分形之Julia集及其算法实现 摘要：本文从自然界的几何现象引出分形的概念，再从其定义、几何特征和分形维的计算这三个方面来加以介绍。以Julia集和Mandelbort集为例来具体描述分形。本文主要从Julia集的特点和算法实现来描述分形以及其实现的方法。 关键词：分形、分数维、Julia集、Mandelbort集、算法实现 引言 大自然是个很伟大的造物者，它留给我们一大笔美丽景观：蜿蜒曲折的海岸线、起伏不定的山脉，变幻无常的浮云，粗糙不堪的断面，袅袅上升的烟柱，九曲回肠的河流，纵横交错的血管，令人眼花缭乱的满天繁星……那么，我们又能从这些美妙的自然现象中得到什么有趣的结论呢？ 正文 分形概述 分形的英文单词为fractal，是由美籍法国数学家曼德勃罗（Benoit Mandelbrot）创造出来的。其取自拉丁文词frangere（破碎、产生无规则碎片）之头，撷英文之尾所合成，本意是不规则的、破碎的、分数的。他曾说：分形就是通过将光滑的形状弄成多个小块，反复的碎弄。1975年，曼德勃罗出版了他的法文专著《分形对象：形、机遇与维数》，标志着分形理论正式诞生。【1】 两种定义 其一：具有自相似性结构的叫做分形； 其二：数学定义：豪斯道夫维Df>=拓扑维Dt。 若一有界集合，包含N个不相重叠的子集，当其放大或缩小r倍后，仍与原集合叠合，则称为自相似集合。自相似集合是分形集。具有相似性的系统叫做分形。 当放大或缩小的倍数r不是一个常数，而必须是r（r1，r2，….）的各种不同放大倍数去放大或缩小各子集，才能与原集合重合时，称为自仿射集合。具有自仿射性的系统叫做分形。【2】 特征 1.自相似性：局部与整体的相似，是局部到整体在各个方向上的等比例变换的结果； 2.自仿射性：是自相似性的一种拓展，是局部到整体在不同方向上的不等比例变换的结果； 3.精细结构：即使对该分形图放大无穷多倍，还是能看到与整体相似的结构，表现出无休止的重复； 4.分形集无法用传统几何语言来描述，它不是某些简单方程的解集，也不是满足某些条件的点的轨 迹； 5.分形集一般可以用简单的方法定义和产生，如递归、迭代；分形其实是由一些简单的图形，经过 递归或者迭代产生的复杂、精细的结构； 6.无确定的标度且具有分数维数。【3】

分形图形与分形的产生

分形图形 分形理论是非线性科学的主要分支之一，它在计算机科学、化学、生物学、天文学、地理学等众多自然科学和经济学等社会科学中都有广泛的应用。分形的基本特征是具有标度不变性。其研究的图形是非常不规则和不光滑的已失去了通常的几何对称性；但是，在不同的尺度下进行观测时，分形几何学却具有尺度上的对称性，或称标度不变性。研究图形在标度变换群作用下不变性质和不变量对计算机图形技术的发展有重大的意义。 说到分形（fractal），先来看看分形的定义。分形这个词最早是分形的创始人曼德尔布诺特提来的，他给分形下的定义就是：一个集合形状，可以细分为若干部分，而每一部分都是整体的精确或不精确的相似形。分形这个词也是他创造的，含有“不规则”和“支离破碎”的意思。分形的概念出现很早，从十九世纪末维尔斯特拉斯构造的处处连续但处处不可微的函数，到上个世纪初的康托三分集，科赫曲线和谢尔宾斯基海绵。但是分形作为一个独立的学科被人开始研究，是一直到七十年代曼德尔布诺特提出分形的概念开始。而一直到八十年代，对于分形的研究才真正被大家所关注。 分形通常跟分数维，自相似，自组织，非线性系统，混沌等联系起来出现。它是数学的一个分支。我之前说过很多次，数学就是美。而分形的美，更能够被大众所接受，因为它可以通过图形化的方式表达出来。而更由于它美的直观性，被很多艺术家索青睐。分形在自然界里面也经常可以看到，最多被举出来当作分形的例子，就是海岸线，源自于曼德尔布诺特的著名论文《英国的海岸线有多长》。而在生物界，分形的例子也比比皆是。 近20年来，分形的研究受到非常广泛的重视，其原因在于分形既有深刻的理论意义，又有巨大的实用价值。分形向人们展示了一类具有标度不变对称性的新世界，吸引着人们寻求其中可能存在着的新规律和新特征；分形提供了描述自然形态的几何学方法，使得在计算机上可以从少量数据出发，对复杂的自然景物进行逼真的模拟，并启发人们利用分形技术对信息作大幅度的数据压缩。它以其独特的手段来解决整体与部分的关系问题，利用空间结构的对称性和自相似性，采用各种模拟真实图形的模型，使整个生成的景物呈现出细节的无穷回归的性质，丰富多彩，具有奇妙的艺术魅力。分形对像没有放大极限，无论如何放大，总会看到更详细的结构。借助于分形的计算机生成，从少量的数据生成复杂的自然景物图形，使我们在仿真模拟方面前进了一大步。在分形的诸多研究课题中，分形的计算机生成问题具有明显的挑战性，它使传统数学中无法表达的形态(如山脉、花草等)得以表达，还能生成一个根本“不存在”的图形世界。分形在制造以假乱真的景物方面的进展和潜在的前途，使得无论怎样估计它的影响也不过分。可以肯定，分形图案在自然界真实物体模拟、仿真形体生成、计算机动画、艺术装饰纹理、图案设计和创意制作等具有广泛的应用价值。 分形图形简介一、关于分形与混沌 关于分形的起源，要非常准确的找出来是非常困难的。研究动态系统、非线形数学、函数分析的科学家，已数不胜数。尽管分形的早期线索已非常古老，但这一学科却还很年轻。比如关于动态系统和细胞自动机的大部分工作可以追溯到冯-诺依曼；但是，直到Mandelbrot 才如此清楚地将自然现象和人工现象中的混沌及分形同自相似性联系在一起。大家如果对此感兴趣，可进一步查阅有关资料。下面我们看一看分形的概念。 什么是分形呢？考虑到此文的意图，我们无意给出它严格的定义，就我们的目的而言，一个分形就是一个图象，但这个图象有一个特性，就是无穷自相似性。什么又是自相似呢？在自然和人工现象中，自相似性指的是整体的结构被反映在其中的每一部分中。比如海岸线，常举的例子，你看它10公里的图象（曲线），和一寸的景象（曲线）是相似的，这就是自相似性。 与分形有着千差万屡的关系的，就是混沌。混沌一词来源与希腊词汇，原意即“张开咀”，但是在社会意义上，它又老爱和无序联系在一起。解释分形和混沌的联系，要注意到分形是

各种有趣的分形

各种有趣的分形 我们看到正方形，圆，球等物体时，不仅头脑里会迅速反映出它是什么，同时，只要我们有足够的数学知识，我们头脑中也反映出它的数学概念，如正方形是每边长度相等的四边形，圆是平面上与某一点距离相等的点的集合，等等。 但是，当我们看到一个山的形状时，我们会想到什么?"这是山"，没错，山是如此的不同于其他景象，以至于你如果绘画水平不高，根本画不出象山的东西。可是，山到底是什么?它既不是三角形，也不是球，我们甚至不能说明山具有怎样的几何轮廓，但为什么我们却有如此直观而又强烈的山的印象?分形的创始人是曼德布洛特思考了这个问题。让 图中的风景图片又是说明分形的另一 很好的例子。这张美丽的图片是利用分 形技术生成的。在生成自然真实的景物 中，分形具有独特的优势，因为分形可 以很好地构建自然景物的模型。 这是一棵厥类植物，仔细观察，你会发 现，它的每个枝杈都在外形上和整体相 同，仅仅在尺寸上小了一些。而枝杈的 枝杈也和整体相同，只是变得更加小 了。 Sierpinski三角形具有严格的自相似特 性

Kohn雪花具有严格的自相似特性 分维及分形的定义 分维概念的提出 对于欧几里得几何所描述的整形来说，可以由长度、面积、体积来测度。但用这种办法对分形的层层细节做出测定是不可能的。曼德尔布罗特放弃了这些测定而转向了维数概念。分形的主要几何特征是关于它的结构的不规则性和复杂性，主要特征量应该是关于它的不规则性和复杂性程度的度量，这可用“维数”来表征。维数是几何形体的一种重要性质，有其丰富的内涵。整形几何学描述的都是有整数维的对象：点是零维的，线是一维的，面是二维的，体是三维的。这种几何对象即使做拉伸、压缩、折叠、扭曲等变换，它们的维数也是不变的；这种维数称为“拓扑维”，记为d。例如当把一张地图卷成筒，它仍然是一个二维信息载体；一根绳子团成团，仍然是一维结构。但曼德尔布罗特认为，在分形世界里，维数却不一定是整数的。特别是由于分形几何对象更为不规则，更为粗糙，更为破碎，所以它的分数维（简称“分维”，记为D）不小于它的拓扑维，即D≥d。 维数和测量有密切关系。如为了测一平面图形的面积，就要用一个边长为l、面积为l2的标准面元去覆盖它，所得的数目就是所测的面积。

分形算法与应用

《分形算法与应用》教学大纲 1 课程的基本描述 课程名称：分形算法与应用Algorithm and Application of Fractal 课程编号：5301A36 课程性质：专业课适用专业：计算机专业 教材选用：孙博文编著，《分形算法与程序设计》，科学出版社，2004.11 总学时：32学时理论学时：32学时 实验学时：0学时课程设计：无 学分：2学分开课学期：第七学期 前导课程：算法分析 后续课程：毕业设计 2 教学定位 2.1 能力培养目标 通过本课程的学习，培养学生的认知和理解能力、逻辑思维能力，以及算法设计与分析能力，程序设计和实现能力。一方面使学生掌握非规则图形的计算机绘制的基本方法，以便实现对不规则对象的算法设计。另一方面，学习本课程的过程也是进行复杂程序设计的训练过程。 2.2 课程的主要特点 本课程是一门重要的专业课，有理论性、设计性与实践性的特点。介绍分形的基本概念及算法设计的基本方法。它是介于计算机软件、程序设计和数学三门课程之间的核心课程。不仅为后续专业课提供了必要的知识基础，也为计算机、软件工程的专业人员提供了必要的技能训练。

2.3 教学定位 通过本课程的学习，使学生达到知识和技能两方面的目标： 1．知识方面：从算法设计及其实现这两个层次的相互关系的角度，系统地学习和掌握非规则图形的算法设计方法，了解并掌握分析、比较和选择不同非规则结构的设计方案，不同运算实现的原则和方法。 2．技能方面：系统地学习和掌握在不同非规则对象实现的不同算法及其设计思想，从中体会并掌握结构选择和算法设计的思维方式及技巧，使分析问题和解决问题的能力得到提高。 3 知识点与学时分配 3.1掌握分形的基本概念 分形简介 分形 分维 分形的测量 共2学时 3.2分形图生成算法之一 分形图的递归算法 Cantor三分集、Koch曲线、Sierpinski垫片、 Peano曲线、分形树等的递归算法。 共2学时 3.3分形图生成算法之二 文法构图算法 LS文法、单一规则的LS文法生成、多规则的LS文法生成、 随机LS文法生成。 共2学时 3.4分形图生成算法之三 迭代函数系统

分形图程序

（1）Koch曲线程序koch.m function koch(a1,b1,a2,b2,n) %koch(0,0,9,0,3) %a1,b1,a2,b2为初始线段两端点坐标,n为迭代次数 a1=0;b1=0;a2=9;b2=0;n=3; %第i-1次迭代时由各条线段产生的新四条线段的五点横、纵坐标存储在数组A、B中 [A,B]=sub_koch1(a1,b1,a2,b2); for i=1:n for j=1:length(A)/5; w=sub_koch2(A(1+5*(j-1):5*j),B(1+5*(j-1):5*j)); for k=1:4 [AA(5*4*(j-1)+5*(k-1)+1:5*4*(j-1)+5*(k-1)+5),BB(5*4*(j-1)+5*(k-1)+1:5*4*(j-1)+5*(k-1)+5)]=sub_koch1(w(k,1),w(k,2),w(k,3),w(k,4)); end end A=AA; B=BB; end plot(A,B) hold on axis equal %由以(ax,ay),(bx,by)为端点的线段生成新的中间三点坐标并把这五点横、纵坐标依次分别存%储在数组A,B中 function [A,B]=sub_koch1(ax,ay,bx,by) cx=ax+(bx-ax)/3; cy=ay+(by-ay)/3; ex=bx-(bx-ax)/3;

ey=by-(by-ay)/3; L=sqrt((ex-cx).^2+(ey-cy).^2); alpha=atan((ey-cy)./(ex-cx)); if (ex-cx)<0 alpha=alpha+pi; end dx=cx+cos(alpha+pi/3)*L; dy=cy+sin(alpha+pi/3)*L; A=[ax,cx,dx,ex,bx]; B=[ay,cy,dy,ey,by]; %把由函数sub_koch1生成的五点横、纵坐标A,B顺次划分为四组,分别对应四条折线段中 %每条线段两端点的坐标,并依次分别存储在4*4阶矩阵k中,k中第i(i=1,2,3,4)行数字代表第%i条线段两端点的坐标 function w=sub_koch2(A,B) a11=A(1);b11=B(1); a12=A(2);b12=B(2); a21=A(2);b21=B(2); a22=A(3);b22=B(3); a31=A(3);b31=B(3); a32=A(4);b32=B(4); a41=A(4);b41=B(4); a42=A(5);b42=B(5); w=[a11,b11,a12,b12;a21,b21,a22,b22;a31,b31,a32,b32;a41,b41,a42,b42];

利用定积分求曲线围成的面积

12.9 利用定积分求曲线围成的面积 武汉外国语学校 汪家硕 一．复习回顾： 1.定积分的几何意义：当()0f x ≥时，积分()b a f x dx ?在几何上表示由()y f x =、x a =、x b =与x 轴所围成的曲边梯形的面积。 当()0f x ≤时，由()y f x =、x a =、x b =与x 轴所围成的曲边梯形位于x 轴的下方。 2.牛顿—莱布尼茨公式 定理（微积分基本定理）如果()f x 是区间[,]a b 上的连续函数，并且' ()()F x f x =，则 ()()()b a f x dx F b F a =-? 二．曲线围成的面积 1.设f 和g 是区间[,]a b 上的连续函数且对任意的[,]x a b ∈有()()f x g x ≥，则直线x a =和直线x b =以及曲线间围成的面积可以表示为：()()()()b b b a a a f x dx g x dx f x g x dx -=-? ?? 例1.求抛物线2y x =和直线2y x =所围成的区域面积。 解：先求出P 点坐标。 解方程组22y x y x ?=?=? ? 02x x =??=? ∴ P 点的坐标是(2,4)。 所求的面积= 2 23220 08424333x x x dx x ??-=-=-=????? 例1 例2．计算曲线21y x =+和24y x =-，以及直线1x =和1x =-所围成的区域面积。 解：所求面积= 11132221112144(1)32333x x x dx x dx x ---??--+=-=-=?????? 例2 2.前面的例题都是一个曲线总在另外一个曲线的上方，如果它们交叉会是什么结果？ 考虑区间112233[,],[,],[,],[,]a c c c c c c b ，阴影部分面积可以表示为： 1 23123()()()()()()()()c c c b a c c c f x g x dx g x f x dx f x g x dx g x f x dx -+-+-+-???? ?b a f (x )dx =c a f (x )dx +b c f (x )dx 。

几个分形的matlab实现资料

几个分形的matlab 实现 摘要：给出几个分形的实例，并用matlab 编程实现方便更好的理解分形，欣赏其带来的 数学美感 关键字：Koch 曲线 实验 图像 一、问题描述： 从一条直线段开始，将线段中间的三分之一部分用一个等边三角形的两边代替，形成山丘形图形如下 图1 在新的图形中，又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边代替，再次形成新的图形如此迭代，形成Koch 分形曲线。 二、算法分析： 考虑由直线段（2个点）产生第一个图形（5个点）的过程。图1中，设1P 和5P 分别为原始直线段的两个端点，现需要在直线段的中间依次插入三个点2P ，3P ，4P 。显然2P 位于线段三分之一处，4P 位于线段三分之二处，3P 点的位置可看成是由4P 点以2P 点为轴心，逆时针旋转600 而得。旋转由正交矩阵 ?????? ? ?-=)3cos()3sin()3sin()3cos(ππππA 实现。 算法根据初始数据（1P 和5P 点的坐标），产生图1中5个结点的坐标。结点的坐标数组形成一个25?矩阵，矩阵的第一行为1P 的坐标，第二行为2P 的坐标……，第五行为5P 的坐标。矩阵的第一列元素分别为5个结点的x 坐标，第二列元素分别为5个结点的y 坐标。 进一步考虑Koch 曲线形成过程中结点数目的变化规律。设第k 次迭代产生的结点数为k n ，第1+k 次迭代产生的结点数为1+k n ，则k n 和1+k n 中间的递推关系为341-=+k k n n 。 三、实验程序及注释：

p=[0 0;10 0]; %P为初始两个点的坐标,第一列为x坐标,第二列为y坐标 n=2; %n为结点数 A=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)]; %旋转矩阵 for k=1:4 d=diff(p)/3; %diff计算相邻两个点的坐标之差,得到相邻两点确定的向量 %则d就计算出每个向量长度的三分之一,与题中将线段三等分对应 m=4*n-3; %迭代公式 q=p(1:n-1,:); %以原点为起点,前n-1个点的坐标为终点形成向量 p(5:4:m,:)=p(2:n,:); %迭代后处于4k+1位置上的点的坐标为迭代前的相应坐标 p(2:4:m,:)=q+d; %用向量方法计算迭代后处于4k+2位置上的点的坐标 p(3:4:m,:)=q+d+d*A'; %用向量方法计算迭代后处于4k+3位置上的点的坐标 p(4:4:m,:)=q+2*d; %用向量方法计算迭代后处于4k位置上的点的坐标 n=m; %迭代后新的结点数目 end plot(p(:,1),p(:,2)) %绘出每相邻两个点的连线 axis([0 10 0 10]) 四、实验数据记录： 由第三部分的程序，可得到如下的Koch分形曲线： 图2 五、注记： １．参照实验方法，可绘制如下生成元的Koch 分形曲线： 图3

数学实验之分形图的绘制

钦州学院数学与计算机科学学院 数学实验报告 专业 : 数学与应用数学班级姓名：学号： 实验完成日期 :2010 年 11 月 1 日，第 10 周，星期一 成绩等级（五级分制）评阅教师评阅日期年月日数学实验报告填写要求：思路清晰，中间结果和最终结果真实；字迹工整，报告完整。[实验题目及内容] 实验题目：分形图形的绘制 实验内容：利用二叉树的画法对生成元带参数进行迭代绘制分形图。 [问题描述]（用自己组织的相关数学语言重述现实问题；注意对约定的条件作说明） 分形图是由一个简单的枝杈不断向周围延伸增加枝干而成，由简单元素生成整体，其中包含有旋转、带参数深层迭代等步骤，对生成元的张开角度和线段长度也有所控制才能绘制出多彩的图形，所以就要设计几个能控制生成图的角度的圆，随时改动分形图的伸张。 [模型建立或思路分析]（建立合理，可解释的数学模型，通过公式、表格或图形直观明确地描述模型的结构；无法通过建立模型解决的，给出解题的思路及办法。） 整个分形图就由几个简单的枝杈进行带参数深层迭代而成，所以先做一个作为整棵树的树主干，做线段AB，以一个B端点作为旋转中心，做两个能控制角度旋转的圆，以圆上所选的角度做适当旋转将线段AB及端点A向上旋转得到两条线段，将得到的线段进行缩放到原来的三分之二，三条线段就组成一个树杈，再继续做另外两个圆选好角度将由线段AB旋转得到的两条线段再向上旋转得到另两条线段，将得到的线段进行缩放到原来的一半。新建参数n=1,对AB两点和参数n进行深度迭代，使得旋转得到的线段的起始点对应

线段AB的起始点，改变n值，即可得到一棵参天大树，即分形图完成。 [实验结果]（通过数学表达式、列表或图形图像的方式显示实验结果。） [结果分析及结论]（对实验结果进行定量分析、合理性分析或误差分析；对所讨论的问题重新认识或提出相关类似问题的拓延；给出自己的意见和合理建议。） 得出的分形图伸张程度和倾斜程度都可以由原先做出的角度控制，改变圆上的角度的大小就可以改变树的弯曲倾斜程度，改变三层基层线段的粗细和颜色可以让分形图更形象，分形图的迭代情况有参数n控制，改变n值增加迭代次数，让树的枝丫伸展使得分形图更多彩。也可以改变生成元的构成，可以在基层增加枝干，进行深层迭代后得出不同形象的分形图。 [求解方法或解题步骤]（针对所建模型或解题思路，给出具体的求解方法或解题步骤。对通过编程解决的问题，画出流程图，给出细节部分的算法，给出相关软件的代码；其他方法解决的，给出详细的解题步骤。）