matlab画分形minkowski,分形系列之matlab绘制koch曲线
从一条直线段开始,将线段中间的三分之一部分用一个等边三角形的两边代替,形成山丘形图形如下,在新的图形中,又将图中每一直线段中间的三分之一部分都用一个等边三角形的两条边代替,再次形成新的图形如此迭代,形成Koch分形曲线。
koch
在matlab中绘制koch曲线,来欣赏下数学之美,开始分形之旅。
1、算法分析:
考虑由直线段(2个点)产生第一个图形(5个点)的过程。图中,设和分别为原始直线段的两个端点,现需要在直线段的中间依次插入三个点P2,P3,P4。显然P2,P4分别位于线段三分之一处,线段三分之二处,P3点的位置可看成是由点以点为轴心,逆时针旋转600而得。旋转由正交矩阵实现。
旋转矩阵
根据初始数据(P1和P5点的坐标),产生图中5个结点的坐标。结点的坐标数组形成一个矩阵,矩阵的第一行为P1的坐标,第二行为P2的坐标……,第五行为P5的坐标。矩阵的第一列元素分别为5个结点的x坐标,第二列元素分别为5个结点的y坐标。
进一步考虑Koch曲线形成过程中结点数目的变化规律。设第k次迭代产生的结点数为n(k),第k+1次迭代产生的结点数为n(k+1),则n(k)和n(k+1)中间的递推关系为n(k+1)=4n(k)-3。
2、代码实现
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:5
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])
%设置坐标显示
3、运行结果
koch曲线
附:
可更改代码中的k值来调节koch曲线的密度,最好不要调的超过10。
参考:<>
matlab画分形minkowski,分形系列之matlab绘制koch曲线相关推荐
- 前端画圆弧html弧线的像素,[js高手之路] html5 canvas系列教程 - arc绘制曲线图形(曲线,弧线,圆形)...
arc:画弧度 cxt.arc( x, y, 半径, 开始角度,结束角度,是否逆时针 ); x, y: 为弧度的中心横坐标和纵坐标,如果这是画一个圆.那么x,y就是圆的圆心. 开始角度与结束角度都是以 ...
- matlab画平面心型线,如何用matlab画出心形线
心形线,是一个圆上的固定一点在它绕着与其相切且半径相同的另外一个圆周滚动时所形成的轨迹线.下面就简单讲解一下如何用matlab画出心形线. 1.心形线的数学定义 2.编制的绘制心形线的matlab程序 ...
- matlab 画电机效率图,如何用matlab绘制电机效率map图或发动机万有特性曲线
如何用matlab绘制电机效率map图或发动机万有特性曲线 前段时间写论文,需要绘制电机效率map图,其实和发动机万有特性曲线一样.. 看了好多资料都不会,问问师兄也没具体画过..困惑中查到貌似有几个 ...
- matlab画y x 2,怎么用MATLAB画x^2+y^2=4*x图像
第一问:不能用通常的三维命令来绘制椭圆球体,应用专用的绘制椭圆球体命令ellipsoid().更改后 由热心网友提供的答案1: 给你图片吧,纯代码常常会被百度认为不符合要求,你自己打这个方程应该是单叶 ...
- matlab 画狗头彩蛋,我用matlab画制动力分配曲线,求程序代码
matlab画函数曲线 c = 3.00e8;h = 6.63e-34;k = 1.38e-23;lambda = (1:10:5000).*1e-9; figure(1)subplot(2,1,1) ...
- 如何用matlab画三维函数相交,如何用matlab绘制三元函数f=x^2 y^3 z^4在
www.mh456.com防采集. 1234567syms x y z;f = (x-y)^2+(x-z)^2+(y-z)^2-100;z = solve(f,'z');for p = 1:lengt ...
- Matlab画圆方向图,如何在matlab里画天线的立体方向图?
如何在matlab里画天线的立体方向图?有没有具体的函数? 有的,你看下一般的matlab教程里都有着方面的介绍,查一下就可以了 有数据,Matlab可以画 归根到底,画极坐标方向图是坐标变换的问题. ...
- matlab画交线曲面平面交线,用matlab如何画以下两面的交线。圆柱面方程x^2+y^2=125^2,平面方程x+z*tan(pi/9)-pi*5/4=0。拜托各位大侠...
答:用AutoCAD可以画出来,我用的2006版 答:clear all [x,y,z]=cylinder(1,100);%生成标准的100个面的圆柱数据,半径为1,高为1,圆心0,0:z1=[z(1 ...
- MATLAB画三维曲面时去掉复数,matlab画三维图怎么除去复数
Q3:matlab 怎么画复数函数? t = -10:0.05:10; fun = @(x) cos(x)+1i*sin(x); figure,plot(fun(t.*pi)), axis squar ...
最新文章
- maven的配置-2019-4-13
- VS添加服务引用和 Web引用的区别
- spring-mvc里的 mvc:resources 及静态资源访问
- 调试实战 | 通过转储文件分析程序无响应之使用 windbg + IDA 逆向篇
- P4593-[TJOI2018]教科书般的亵渎【拉格朗日差值】
- 开源合同管理系统_「物联网架构」最适合物联网的开源数据库
- vscode修改c 项目_windows 10上使用vscode编译运行和调试C/C++
- C++类的定义要注意最后的分号不要丢了
- 卡巴斯基KEY官方生成器
- 税控盘版开票助手开发
- 尺缩钟慢之动钟变慢——思想实验推导狭义相对论(七)
- C/C++手机通信录
- linux美元符号进入井号,MyBatis中井号与美元符号的区别
- Erlang 游戏开发经验总结
- python的数据类型分为复数_Python的基本数据类型
- 报头压缩-ROHC压缩系统
- opencv图像旋转和翻转,cv2.flip,cv2.rotate
- coursera 吴恩达 -- 第一课 神经网络和深度学习 :第三周课后习题 Shallow Neural Networks Quiz, 10 questions
- java例外站点添加不上,java例外站点添加不了
- 【启动盘制作】单/多系统启动盘制作