move.m

来自「利用matla模拟薄膜生长」· M 代码 · 共 29 行

%原子团迁移
DD=zeros(222);
EE=zeros(222);
EE(:)=100;
r=ceil((100-E)/10); %搜索半径
for i=12:211
    for j=12:211
        a=r(i,j);
        if a<2    %若半径过小则不执行
            continue
        end
        for k=2:a
            for ii=i-k:i+k
                for jj=j-k:j+k
                    if ii==i&jj==j   %判断自身点的时候跳过执行程序
                        continue
                    else if E(ii,jj)<100
                            m=fix(i+((ii-i)*N(ii,jj))/(N(i,j)+N(ii,jj)));
                            n=fix(j+((jj-j)*N(ii,jj))/(N(i,j)+N(ii,jj)));%m,n为迁移后的坐标
                            DD(m,n)=D(i,j)+D(ii,jj);%迁移后的粒子个数
                            EE(m,n)=E(i,j)-D(ii,jj);%迁移后的能量
                        end
                    end
                end
            end
            break
        end
    end
end