【元胞自动机】基于matlab元胞自动机晶体生长【含Matlab源码 232期】 |
您所在的位置:网站首页 › 自动机操作视频 › 【元胞自动机】基于matlab元胞自动机晶体生长【含Matlab源码 232期】 |
✅博主简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,Matlab项目合作可私信。 🍎个人主页:海神之光 🏆代码获取方式: 海神之光Matlab王者学习之路—代码获取方式 ⛳️座右铭:行百里者,半于九十。 更多Matlab仿真内容点击👇 Matlab图像处理(进阶版) 路径规划(Matlab) 神经网络预测与分类(Matlab) 优化求解(Matlab) 语音处理(Matlab) 信号处理(Matlab) 车间调度(Matlab) ⛄一、元胞自动机简介1 元胞自动机发展历程 最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物 细胞的自我复制而提出的. 但是并未受到学术界重视. 1970 年, 剑桥大学的约翰 · 何顿 · 康威设计了一个电脑游戏 “生命游戏” 后, 元胞自动机才吸引了科学家们的注意. 1983 年 S.Wolfram 发表了一系列论文. 对初等元胞机 256 种 规则所产生的模型进行了深入研究, 并用熵来描述其演化行 为, 将细胞自动机分为平稳型, 周期型, 混沌型和复杂型. 2 对元胞自动机的初步认识 元胞自动机(CA)是一种用来仿真局部规则和局部联系的方法。典型的元胞自动机是定义在网格上的,每一个点上的网格代表一个元胞与一种有限的状态。变化规则适用于每一个元胞并且同时进行。典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。 3 元胞的变化规则&元胞状态 典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。 4 元胞自动机的应用 元胞自动机已被应用于物理模拟,生物模拟等领域。 5 元胞自动机的matlab编程 结合以上,我们可以理解元胞自动机仿真需要理解三点。一是元胞,在matlab中可以理解为矩阵中的一点或多点组成的方形块,一般我们用矩阵中的一点代表一个元胞。二是变化规则,元胞的变化规则决定元胞下一刻的状态。三是元胞的状态,元胞的状态是自定义的,通常是对立的状态,比如生物的存活状态或死亡状态,红灯或绿灯,该点有障碍物或者没有障碍物等等。 6 一维元胞自动机——交通规则 定义: 6.1 元胞分布于一维线性网格上. 6.2 元胞仅具有车和空两种状态. 8 什么是元胞自动机 离散的系统: 元胞是定义在有限的时间和空间上的, 并且元 胞的状态是有限. 动力学系统: 元胞自动机的举止行为具有动力学特征. 简单与复杂: 元胞自动机用简单规则控制相互作用的元胞 模拟复杂世界.
(5)规则(状态转移函数) 定义:根据元胞当前状态及其邻居状况确定下一时刻该元胞状态的动力学函数, 简单讲, 就是一个状态转移函数. 分类 : 总和型: 某元胞下时刻的状态取决于且仅取决于它所有邻居 的当前状态以及自身的当前状态. 合法型: 总和型规则属于合法型规则. 但如果把元胞自动机 的规则限制为总和型, 会使元胞自动机具有局限性. (6)森林火灾 宏观连续模型: 1 加速规则:不能超过v m a x ( 2 格 / s ) v_{max}(2格/s)v max(2格/s) 2 防止碰撞:不能超过车距 理论分析: 结果分析: 时空轨迹 % DLA clc;clear;close all; S=zeros(400,500); % 生成状态矩阵 S(end,:)=1; % 设置状态矩阵中最下面一行元素等于1 A=1;B=1;X=0.8; rand(‘state’,0); % 设置随机数的状态数 subplot(121);Ii=imshow(1-S,[]); % 显示状态矩阵 T1=title([‘times = 1’,‘, total particle=’,num2str(sum(S(😃))],… ‘Fontname’,‘times new roman’,‘fontsize’,14); % 显示时间与粒子总数 r=rand(1,500); subplot(122);P1=plot(sum(S,2)/size(S,2),1:size(S,1),‘r’);% 绘制各行的密度值曲线 set(gca,‘Position’,[0.57,0.35,0.33,0.36],‘YDir’,‘reverse’); % 设置坐标轴属性 xlim([0,max(sum(S,2)/size(S,2))]); % 设置x轴的范围 ylabel(‘\ith’,‘fontname’,‘times new roman’,‘fontsize’,14); % y轴标注 xlabel(‘\it\rho’,‘fontname’,‘times new roman’,‘fontsize’,14); % x轴标注 title(‘{\it\rho} ({\ith})’,‘fontname’,‘times new roman’,‘fontsize’,14); % 加注图题 set(gcf,‘DoubleBuffer’,‘on’); % 设置图形窗口的渲染效果 [L1,L2]=size(S); % 返回状态矩阵的行数L1和列数L2 N=500;H=1; % 初始化参数:粒子总数N和时间参数H h=150; % 设置截顶高度 scale=0.5; % 设置剪切系数 while N0.5; % 判断近邻中有固定粒子 S(R1,R2)=1; % 运动粒子被吸附 flag=1; % 标记粒子已经被吸附 else ra=rand; % 粒子进行随机移动的分量 rb=rand; % 粒子进行随机移动的分量 R1=R1+(ra>=0.5)-(ra=0.5)-(rb |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |