- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
几个分形得matlab实现
摘要:给出几个分形得实例,并用matlab编程实现方便更好得理解分形,欣赏其带来得数学美感
关键字:Koch曲线 实验 图像
问题描述:
从一条直线段开始,将线段中间得三分之一部分用一个等边三角形得两边代替,形成山丘形图形如下
?图1
在新得图形中,又将图中每一直线段中间得三分之一部分都用一个等边三角形得两条边代替,再次形成新得图形如此迭代,形成Koch分形曲线。
算法分析:
考虑由直线段(2个点)产生第一个图形(5个点)得过程。图1中,设与分别为原始直线段得两个端点,现需要在直线段得中间依次插入三个点,,。显然位于线段三分之一处,位于线段三分之二处,点得位置可瞧成就是由点以点为轴心,逆时针旋转600而得。旋转由正交矩阵
实现。
算法根据初始数据(与点得坐标),产生图1中5个结点得坐标、结点得坐标数组形成一个矩阵,矩阵得第一行为得坐标,第二行为得坐标……,第五行为得坐标。矩阵得第一列元素分别为5个结点得坐标,第二列元素分别为5个结点得坐标。
进一步考虑Koch曲线形成过程中结点数目得变化规律。设第次迭代产生得结点数为,第次迭代产生得结点数为,则与中间得递推关系为。
三、实验程序及注释:
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
五、注记:
1、参照实验方法,可绘制如下生成元得Koch 分形曲线:
图3
此时,旋转矩阵为:
程序与曲线如下:
p=[0 0;10 0]; %P为初始两个点得坐标,第一列为x坐标,第二列为y坐标
n=2; %n为结点数
A=[0 —1;1 0]; %旋转矩阵
for k=1:4
d=diff(p)/3; %diff计算相邻两个点得坐标之差,得到相邻两点确定得向量
%则d就计算出每个向量长度得三分之一,与题中将线段三等分对应
m=5*n—4; %迭代公式
q=p(1:n—1,:); %以原点为起点,前n-1个点得坐标为终点形成向量
p(6:5:m,:)=p(2:n,:); %迭代后处于5k+1位置上得点得坐标为迭代前得相应坐标
p(2:5:m,:)=q+d; %用向量方法计算迭代后处于5k+2位置上得点得坐标
p(3:5:m,:)=q+d+d*A'; %用向量方法计算迭代后处于5k+3位置上得点得坐标
p(4:5:m,:)=q+2*d+d*A'; %用向量方法计算迭代后处于5k+4位置上得点得坐标
p(5:5:m,:)=q+2*d; %用向量方法计算迭代后处于5k位置上得点得坐标
n=m; %迭代后新得结点数目
end
plot(p(:,1),p(:,2)) %绘出每相邻两个点得连线
axis([0 10 0 10])
?图4
由于中间三分之一部分就是一个正方形时,有很多连接得部分。所以我们将高度压缩到原来得0。7倍,即中间部分为一个长与宽之比为1:0。7得矩形时,得到程序与曲线如下:
p=[0 0;10 0]; %P为初始两个点得坐标,第一列为x坐标,第二列为y坐
文档评论(0)