clc;
clear;
close all;
warning off;
addpath(genpath(pwd));
numIt=1000;
x=50;
y=50;
z=50;
grid = zeros(x,y,z);
for a=20:30
for b=20:30
for c=20:30
grid(a,b,c)=1;
end
end
end
for n=1:numIt
for i=2:x-1
for j=2:y-2
for k=2:z-2
if ( grid(i,j,k)== 1 )
grid(i,j,k)=0;
r=rand(1);
%Direction possibility #1
if (r < 0.1666)
if ( grid(i+1,j,k)==1 )
grid(i,j,k)=1;
else
grid(i+1,j,k)=1;
end
end
%Direction possibility #2
if (r >= 0.1666 && r < 0.3332)
if ( grid(i-1,j,k)==1 )
grid(i,j,k)=1;
else
grid(i-1,j,k)=1;
end
end
%Direction possibility #3
if (r >= 0.3332 && r < 0.5000)
if ( grid(i,j+1,k)==1 )
grid(i,j,k)=1;
else
grid(i,j+1,k)=1;
end
end
%Direction possibility #4
if (r >= 0.5000 && r < 0.6666)
if ( grid(i,j-1,k)==1 )
grid(i,j,k)=1;
else
grid(i,j-1,k)=1;
end
end
%Direction possibility #5
if (r >= 0.6666 && r < 0.8333)
if ( grid(i,j,k+1)==1 )
grid(i,j,k)=1;
else
grid(i,j,k+1)=1;
end
end
%Direction possibility #6
if (r >= 0.8333)
if ( grid(i,j,k-1)==1 )
grid(i,j,k)=1;
else
grid(i,j,k-1)=1;
end
end
else
end
end
end
end
%Plot
figure(1)
xslice = [25, 50]; % location of y-z planes
yslice = [25, 50]; % location of x-z plane
zslice = [1, 25]; % location of x-y planes
slice(grid,xslice,yslice,zslice)
drawnow
end
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
基于随机游走算法的气体扩散matlab仿真模拟.rar (2个子文件)
基于随机游走算法的气体扩散matlab仿真模拟
操作录像0031.avi 6.14MB
Runme.m 2KB
共 2 条
- 1
资源评论
fpga和matlab
- 粉丝: 16w+
- 资源: 2561
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功