迭代函数系统matlab实现,线性随机IFS迭代Matlab示例程序.doc

线性随机IFS迭代Matlab示例程序

线性随机IFS迭代 Matlab示例程序 Tree Leaf Mountain

%% 线性随机IFS迭代[iterated function systems迭代函数系统]

clear;clc;close

a1 [0 0 0 0.16 0 0 0.01;%树叶 0.85 0.01 -0.04 0.85 0 1.6 0.85 0.2 -0.26 0.23 0.22 0.01 1.6 0.07 -0.15 0.28 0.26 0.24 0 0.44 0.07];

a2 [0.08 -0.031 0.084 0.0306 5.17 7.97 0.03%山脉 0.0801 0.0212 -0.08 0.0212 6.62 9.4 0.025 0.75 0 0 0.53 -0.357 1.106 0.22 0.943 0 0 0.474 -1.98 -0.65 0.245 -0.402 0 0 0.402 15.513 4.588 0.21 0.217 -0.052 0.075 0.15 3 5.74 0.07 0.262 -0.105 0.114 0.241 -0.473 3.045 0.1 0.22 0 0 0.43 14.6 4.286 0.1];

a3 [0.05000.6000.1%松树 0.0500-0.5010.1 0.460.32-0.3860.38300.60.2 0.47-0.1540.1710.423010.2 0.430.275-0.260.476010.2 0.421-0.3570.3540.30700.70.2];

a a3;

N 50000;

x0 1;y0 1;

[m,n] size a ;

b [0 a :,n '];

for i 1:N r rand; for j 1:m if r sum b 1:j &r sum b 1:j+1 ;k j;end end X [a k,1 ,a k,2 ;a k,3 ,a k,4 ]*[x0;y0]+[a k,5 ;a k,6 ]; x1 X 1 ;y1 X 2 ;x0 x1;y0 y1; plot x1,y1,'k' ;hold on; if ~mod i,1000 ;pause eps ;title num2str i ;end

end