函数主
清除
L=32； %基础大小设置；
步骤1=500； %成长层数的设定是
niu=1； deltt=0.05； deltx=1； dd=16； aa=(2*DD/Deltx ) ^ )1/2)； %基本参数设置；
cs=100； %重复计算次数
h=Zeros(L2，step1)； %基底高度初始值；
for k=1:step1-1
ww(k )=0；
结束
for cc=1:cs； %重复计算
for t=2:step1
for i=3:L 1
%以下是分数增长方程的离散形式
h(1，t )=h ) L1，t )； %边界条件赋值
h(L2，t )=H ) 2，t )； %边界条件赋值
random=unifrnd(-1/2、1/2、1、1 )； %随机噪声； 随机数；
h(I，t 1 )=h ) I，t )-niu*deltt* ) h(I2，t ) h ) I1，t )5) h ) I，t )-4 ) h ) I-1，t ) h ) I-2，t ) h
结束
结束
% h (1: L1，1: step1)；
for k=1:step1-1；
w(k )=STD ) h ) 2:l1，k 2 )； %surface width； 表面粗糙度； 表面粗糙度
ww(k ) ww(k ) w ) k；
t(k )=) k2 ) *deltt； %生长时间=生长层数step1*基底尺寸L*时间间隔deltt
结束
ccc=cc
结束
save date
plot(log2 ) t )、log2 ) ww/cs )、' square ' ) )
p=polyfit(log2(t )，log2 (ww/cs )，1 ) ) ) ) ) ) ) ) ) ) )。
%plot(h ) 2:L，100 )，'.- ' )