数学实验 分形.docVIP

用计算机绘制出Koch曲线的图形。你能否自己构造一些生成元，并由此绘出相应的分形图形？任取分形图形的一个局部并将它放大，它同原来的分形图形有什么联系？ 解答：Koch曲线的图形如下(n=5)： 图表 1Koch曲线(n=5) 构造一个和Koch曲线类似的生存元，如下图所示： 图表 2生存元 相应的分形图如下： 图表 3局部放大图形 可知局部放大后和原来的图形一样。 主要的代码如下所示： function output = Koch(n) if n 0 error(n must be positive number!); end x = [0 1 -1 0]; y =[sqrt(3) 0 0 sqrt(3)]; for i =1 : n myLength = numel(x); verX=[]; verY=[];index=1; for j = 1:myLength-1 verX(index) = x(j); verY(index) = y(j); verX(index+1) = 2/3*x(j) + 1/3*x(j+1); verY(index+1) = 2/3*y(j) + 1/3*y(j+1); tmpSol = mySolve(x(j), y(j), x(j+1), y(j+1)); verX(index+2) = tmpSol(1); verY(index+2) = tmpSol(2); verX(index+3) = 1/3*x(j) + 2/3*x(j+1); verY(index+3) = 1/3*y(j) + 2/3*y(j+1); index = index + 4; end verX(index) = verX(1); verY(index) = verY(1); x = verX; y = verY; end fill(x,y,r); title([n= num2str(n)]); axis equal output=[x;y]; function y = MySolve(x1,y1,x2,y2) a1 = 1/2*x1-1/6*y2*3^(1/2)+1/6*y1*3^(1/2)+1/2*x2; b1 = -1/6*3^(1/2)*(-y1*3^(1/2)-y2*3^(1/2)-x2+x1); a2 = 1/2*x1+1/6*y2*3^(1/2)-1/6*y1*3^(1/2)+1/2*x2; b2 = 1/6*3^(1/2)*(y1*3^(1/2)+y2*3^(1/2)-x2+x1); if (a1^2 + (b1 - sqrt(3) * 2/3)^2) (a2^2 + (b2 - sqrt(3) * 2/3)^2) y(1) = a1; y(2) = b1; else y(1) = a2; y(2) = b2; end 2. 定义Weierstrass函数如下： 对不同的s值，画出函数的图像，观察图像的不规则性与ｓ的关系，由此猜测Weierstrass函数图像的维数与ｓ的关系。 主要的程序代码如下： x = linspace(-4,4,100); lambda = 2; i = 1; for s = 1.2: 0.1: 1.5 y = Weierstrass(x, lambda, s); subplot(2,2,i) plot(x,y) title([s= num2str(s)]); i = i + 1; end function y = Weierstrass(x, lambda, s) n = 100; y=zeros(size(x)); for i=1:n y=y + (lambda^((s - 2) * i ) ) * sin(lambda ^ i .* x); end 3.编写Julia集的程序。对不同的参数 观察Julia集的变化。取Julia集的不同局部放大，你能看出某种自相似现象吗？ 产生图像如下所示： 图表 4 c=0 + 1i 图表 5 c = 1 + 0i 图表 6 c = 0.11 + 0.66i 图表 7 c = 0.1028 + 0.95723i 图表 8 c = 1.25 + 0.01i 主要的程序的代码： clear clc p=[0 -1 0.11 -0