suiji.m

来自「1． 对薄膜形成与生长中相关物理过程及现象进行分析和建模。 2． 利用蒙特卡罗」· M 代码 · 共 33 行

%suiji.m
d=round(rand(20)*M);%(rand(20)*m)20*20的随机矩阵，每元素在0~1*100间,round对每元素四舍五入取整；
e=round(rand*M);%0~100间的随机数；
f=round(rand*M);%0~100间的随机数；
if(e~=0&f~=0)%若产生的随机数都不为零；
m=d(f);%d矩阵中从竖行数第f个元素值；
n=d(e);
elseif(f==0&e~=0)
    m=d(f+1);
    n=d(e);
elseif(e==0&f~=0)
    m=d(f);
    n=d(e+1);
else
    f=f;e=e;
end
while  m==0|n==0%若在上面(rand(20)*m)20*20的随机矩阵中随机取到元素有一为零，则重新取数,防止脚标取到零；
d=round(rand(20)*M);%重新产生(rand(20)*m)20*20的随机矩阵，每元素在0~1*100间,round对每元素四舍五入取整；
e=round(rand*M);%重新产生e并覆盖前面的e值；
f=round(rand*M);
if(e~=0&f~=0)
m=d(f);
n=d(e);
elseif(f==0&e~=0)
    m=d(f+1);
    n=d(e);
elseif(e==0&f~=0)
    m=d(f);
    n=d(e+1);
else
    f=f;e=e;
end
end