function main
clear
L=32;%基底尺寸的设定;
step1=500;%生长层数的设定,
niu=1;deltt=0.05;deltx=1;dd=16;aa=(2*dd/deltx)^(1/2);%基本参数的设定;
cs=100; %重复计算次数
h=zeros(L+2,step1);%基底高度的初值;
for k=1:step1-1
ww(k)=0;
end
for cc=1:cs;%重复计算
for t=2:step1
for i=3:L+1
%以下是分数生长方程的离散形式
h(1,t)=h(L+1,t);%边界条件赋值
h(L+2,t)=h(2,t);%边界条件赋值
random=unifrnd(-1/2,1/2,1,1);%随机噪声;随机数;
h(i,t+1)=h(i,t)-niu*deltt*(h(i+2,t)-4*h(i+1,t)+5*h(i,t)-4*h(i-1,t)+h(i-2,t))+aa*(12*deltt)^(1/2)*random;
end
end
% h(1:L+1,1:step1);
for k=1:step1-1;
w(k)=std(h(2:L+1,k+2));%surface width;surface roughnenss;表面粗糙度
ww(k)=ww(k)+w(k);
t(k)=(k+2)*deltt; %生长时间=生长层数step1*基底尺寸L*时间间隔deltt
end
ccc=cc
end
save date
plot(log2(t),log2(ww/cs),'square')
p=polyfit(log2(t),log2(ww/cs),1)
% plot(h(2:L,100),'.-')