分形实验报告.docVIP

分形实验报告 实验：分形 1. 用计算机绘制出Koch曲线的图形。你能否自己构造一些生成元，并由此绘出相应的分形图形？任取分形图形的一个局部并将它放大，它同原来的分形图形有什么联系？ 代码： function result=KochSnow if -a)+-a) c=[a+-a)/3 a+-a)/3]; d=[a+2*-a)/3 a+2*-a)/3]; alpha=atan-c)/-c)); l=sqrt-c)+-c)); if e=c+cos*l; e=c+sin*l; result=0; KochSnow;KochSnow; KochSnow;KochSnow; endend clear clc n=input; a=[0 0]; b=[130 0]; c=[b/2 b*sqrt/2]; KochSnow; KochSnow; KochSnow; 2016秋 数学实验实验：分形 页 第 1 页 共 10 实验：分形 运行结果： 请输入次数:6 120160 806040200-20-40-20 020406080100120160 2. 定义Weierstrass函数如下： W???ksin,??1, k?1 ? 1?s?2 对不同的s值，画出函数的图像，观察图像的不规则性与ｓ的关系，由此猜测Weierstrass函数图像的维数与ｓ的关系。 代码： function y = Weierstrass n=1000; y=zeros); for i=1:n 2016秋 数学实验实验：分形第 2 页 共 10 页 实验：分形 y=y+*i))* sin; end end clear clc x = linspace; lambda = 2; i = 1; for s = : : y = Weierstrass; subplot plot title]); i = i + 1; end 运行结果： s= 21 0-1-2-4 -2 02 4 s= 21 0-1-2-4 -2 2 4 2016秋 数学实验实验：分形 页 s= 21 0-1-2-4 -2 02 4 s= 42 0-2-4-4 -2 2 4 第 3 页 共 10 实验：分形 3. 编 写 Julia 集 的 程 序 。 对 不 同 的 参 数 :,,,, 观察Julia集的变化。取Julia集的不同局部放大，你能看出某种自相似现象吗？ 代码： function Julia if nargin r=max,2); l=linspace; A=ones*l+i**l); B=zeros; for s = 1:m B = B+ imagesc; colormap; hold off; axis equal; axis off; clear clc A=[0 1;-1 0; ;- ;- -]; a=;b=; p=0;q=1; n=500; m=14; 2016秋 数学实验实验：分形 页 第 4 页 共 10 实验：分形 for i=1:5 subplot; p=A; q=A; Julia; title) + num2str) i]); end 运行结果： c=0+ 1 i c=-+ i c=-+ i c=-1+ 0 i c=+ i c=-+ - i c=-+ i 2016秋 数学实验实验：分形第 5 页 共 10 页 实验六：分形图的生成