【元胞自动机】基于matlab元胞自动机城市规划【含Matlab源码 125期】

✅博主简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,Matlab项目合作可私信。
🍎个人主页:海神之光
🏆代码获取方式:
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭:行百里者,半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理(进阶版)
路径规划(Matlab)
神经网络预测与分类(Matlab)
优化求解(Matlab)
语音处理(Matlab)
信号处理(Matlab)
车间调度(Matlab)

⛄一、案例简介(附课程报告)

城市规模设计
雄安新区占地总面积约为2000平方公里,涉及河北省雄县、容城、安新3个县及周边部分区域,地处北京、天津、保定腹地,通过ArcGIS地图软件搜索该区域并从中提取出来,区域图如下所示。
在这里插入图片描述
图5 雄安新区区域图
为对雄安新区进行更好的仿真模拟,首先先在地图中截取雄安新区地图,然后进行边缘轮廓提取和白洋淀等不可开发地区的剔除,获得预处理图像。最后用MATlAB进行图像灰度化、二值化处理如下图所示。为后续元胞自动机提供演变地图。
在这里插入图片描述
城市规划模型总步骤:
首先确定其组成的主要元素:元胞、元胞空间、元胞状态、元胞邻域及转变规则
分析模拟城市空间结构;
确定模型的参数:繁殖参数、扩散参数、传播参数及受规划约束参数
确定模型所需元胞转换规则的定义
进行城市发展模拟。
①本文提取的雄安新区地图像素为,元胞空间定义为;元胞状态对应的是该地的四种状态:未城市化(即对应能开发还未开发的区域),城市化,扩展中心城市,不能被开发(如白洋淀等区域)。土地状态用编码表示。

⛄二、元胞自动机简介

1 元胞自动机发展历程
最初的元胞自动机是由冯 · 诺依曼在 1950 年代为模拟生物 细胞的自我复制而提出的. 但是并未受到学术界重视.
1970 年, 剑桥大学的约翰 · 何顿 · 康威设计了一个电脑游戏 “生命游戏” 后, 元胞自动机才吸引了科学家们的注意.

1983 年 S.Wolfram 发表了一系列论文. 对初等元胞机 256 种 规则所产生的模型进行了深入研究, 并用熵来描述其演化行 为, 将细胞自动机分为平稳型, 周期型, 混沌型和复杂型.

2 对元胞自动机的初步认识
元胞自动机(CA)是一种用来仿真局部规则和局部联系的方法。典型的元胞自动机是定义在网格上的,每一个点上的网格代表一个元胞与一种有限的状态。变化规则适用于每一个元胞并且同时进行。典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。

3 元胞的变化规则&元胞状态
典型的变化规则,决定于元胞的状态,以及其( 4 或 8 )邻居的状态。

4 元胞自动机的应用
元胞自动机已被应用于物理模拟,生物模拟等领域。

5 元胞自动机的matlab编程
结合以上,我们可以理解元胞自动机仿真需要理解三点。一是元胞,在matlab中可以理解为矩阵中的一点或多点组成的方形块,一般我们用矩阵中的一点代表一个元胞。二是变化规则,元胞的变化规则决定元胞下一刻的状态。三是元胞的状态,元胞的状态是自定义的,通常是对立的状态,比如生物的存活状态或死亡状态,红灯或绿灯,该点有障碍物或者没有障碍物等等。

6 一维元胞自动机——交通规则
定义:
6.1 元胞分布于一维线性网格上.
6.2 元胞仅具有车和空两种状态.
在这里插入图片描述
7 二维元胞自动机——生命游戏
定义:
7.1 元胞分布于二维方型网格上.
7.2 元胞仅具有生和死两种状态.
在这里插入图片描述
元胞状态由周围八邻居决定.
规则:
在这里插入图片描述
骷髅:死亡;笑脸:生存
周围有三个笑脸,则中间变为笑脸
少于两个笑脸或者多于三个,中间则变死亡。

8 什么是元胞自动机
离散的系统: 元胞是定义在有限的时间和空间上的, 并且元 胞的状态是有限.
动力学系统: 元胞自动机的举止行为具有动力学特征.
简单与复杂: 元胞自动机用简单规则控制相互作用的元胞 模拟复杂世界.

在这里插入图片描述
9 构成要素
在这里插入图片描述
(1)元胞 (Cell)
在这里插入图片描述
元胞是元胞自动机基本单元:
状态: 每一个元胞都有记忆贮存状态的功能.
离散: 简单情况下, 元胞只有两种可能状态; 较复杂情况下, 元胞具有多种状态.
更新: 元胞的状态都安照动力规则不断更新.
(2)网格 (Lattice)
不同维网格
在这里插入图片描述
常用二维网格
在这里插入图片描述
(3)邻居 (Neighborhood)
在这里插入图片描述
(4)边界 (Boundary)
在这里插入图片描述
反射型:以自己作为边界的状态
吸收型:不管边界(车开到边界就消失)

(5)规则(状态转移函数)
定义:根据元胞当前状态及其邻居状况确定下一时刻该元胞状态的动力学函数, 简单讲, 就是一个状态转移函数.
分类 :
总和型: 某元胞下时刻的状态取决于且仅取决于它所有邻居 的当前状态以及自身的当前状态.
合法型: 总和型规则属于合法型规则. 但如果把元胞自动机 的规则限制为总和型, 会使元胞自动机具有局限性.
(6)森林火灾
在这里插入图片描述
绿色:树木;红色:火;黑色:空地。
三种状态循环转化:
树:周围有火或者被闪电击中就变成火。
空地:以概率p变为树木
理性分析:红为火;灰为空地;绿是树
在这里插入图片描述
元胞三种状态的密度和为1
在这里插入图片描述
火转化为空地的密度等于空地转换为树的密度(新长出来的树等于烧没的树)
在这里插入图片描述
f是闪电的概率:远远小于树生成的概率;T s m a x T_{smax}T smax
​是一大群树被火烧的时间尺度
程序实现
周期性边界条件
购进啊
在这里插入图片描述
其中的数字为编号
构建邻居矩阵
在这里插入图片描述
上面矩阵中的数字编号,对应原矩阵相同位置编号的上邻居编号,一 一对应
同样道理:
在这里插入图片描述
(7)交通概念
车距和密度
在这里插入图片描述
流量方程
在这里插入图片描述
守恒方程
在这里插入图片描述
时空轨迹(横轴是空间纵轴为时间)
在这里插入图片描述
红线横线与蓝色交点表示每个时间车的位置。
如果是竖线则表示车子在该位置对应的时间

宏观连续模型:
在这里插入图片描述
最常用的规则:
在这里插入图片描述
红色条表示速度是满的。

1 加速规则:不能超过v m a x ( 2 格 / s ) v_{max}(2格/s)v
max(2格/s)
2 防止碰撞:不能超过车距

理论分析:
在这里插入图片描述
结果分析: 密度与流量
在这里插入图片描述
第一个图:横坐标是归一化后的密度,纵坐标是车流量。第二个图:理论值与CA的结果

结果分析: 时空轨迹
在这里插入图片描述
中间的深色区域是交通堵塞的区域。

⛄三、部分源代码

clear;close all
I=imread(‘city.png’);
level=graythresh(I); % 图像灰度处理
bw=im2bw(I,level); % 图像二值化处理
I=bw;
cells=double(I);

%定义 button
plotbutton=uicontrol(‘style’,‘pushbutton’,…
‘string’,‘Run’,…
‘fontsize’,12,…
‘position’,[100,400,50,20],…
‘callback’,‘run=1;’);
erasebutton=uicontrol(‘style’,‘pushbutton’,…
‘string’,‘Stop’,…
‘fontsize’,12,…
‘position’,[300,400,50,20],…
‘callback’,‘freeze=1;’);
number=uicontrol(‘style’,‘text’,…
‘string’,‘1’,…
‘fontsize’,12,…
‘position’,[20,400,50,20]);

%初始化
cells(33,44)=2;
cells(88,31)=2;
cells(33,80)=2;
%参数定义
kuosan=0.2;fanzhi=0.2;chuanbo=0.2;chengshihua=0.0004;sch=[];skz=[];t=1;
[x,y]=size(cells);
run=0;
freeze=0;
stop=0;
while (stop0)
if(run
1)
for i=1:x
for j=1:y
if(cells(i,j)~=1)
if(cells(i,j)==0) %自然增长
if(rand<chengshihua)
cells(i,j)=2;
end
if(aroundcenter(i,j,cells))
if(rand<chuanbo)
cells(i,j)=2;
end
end
end
if(cells(i,j)==2 && rand<kuosan) %边界增长
if(exist1(i,j,cells))
m=1+3*rand;
switch m
case 1
ii=i-1;jj=j;
case 2

                   otherwise 4
                       ii=i+1;jj=j;
               end
           if(cancity(ii,jj,cells))  
               cells(ii,jj)=2;
           end
           end
    end
   if(cells(i,j)==2 && exist1(i,j,cells))      %新扩展中心型
       if(rand<fanzhi)
           if(cancity(i,j,cells))
               cells(i,j)=3;
           end
       end
   end
    end
end

end
ch=0;kzch=0;
for i=1:x
for j=1:y
if(cells(i,j)==2) ch=ch+1;end
if(cells(i,j)==3) kzch=kzch+1;end
end
end
sch(t)=ch;skz(t)=kzch;
t=t+1;
[A,B]=size(cells);
Area(1:A,1:B,1)=zeros(A,B);
Area(1:A,1:B,2)=zeros(A,B);
Area(1:A,1:B,3)=zeros(A,B);
for i=1:A
for j=1:B
if cells(i,j)==1

end
end
end
pause(0.0005);
Area=uint8(Area);
Area=imagesc(Area);

%计步
stepnumber=1+str2num(get(number,‘string’));
set(number,‘string’,num2str(stepnumber));
end
if freeze==1
run=0;
freeze=0;
end
drawnow
end
function result=cancity(i,j,cells)
s=0
if(cells(i-1,j-1)==1) s=s+1;end
if(cells(i-1,j)==1) s=s+1;end
if(cells(i-1,j+1)==1) s=s+1;end
if(cells(i,j-1)==1) s=s+1;end
if(cells(i,j+1)==1) s=s+1;end
if(cells(i+1,j-1)==1) s=s+1;end
if(cells(i+1,j)==1) s=s+1;end
if(cells(i+1,j+1)==1) s=s+1;end
if(s>=4)
result=0;
else
result=1;
end
end

function result= exist1(i,j,cells)

a=0;
if(cells(i-1,j)==2) a=a+1; end

if(cells(i+1,j)==2) a=a+1;end
if(a>=2)
result=1;
else
result=0;
end
end

⛄四、运行结果

城市规模仿真模拟
初始阶段
在这里插入图片描述

图7 元胞自动机演化初始阶段
中期阶段
在这里插入图片描述

在这里插入图片描述

图8 元胞自动机演化中期阶段
远期阶段
在这里插入图片描述
在这里插入图片描述

图9 元胞自动机演化远期阶段

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]党珊,蒋太刚,巫承军.基于元胞自动机方法的消防疏散仿真研究[J].现代电子技术. 2022,45(21)
[2]帅斌,秦梦瑶,许旻昊.基于元胞自动机的高速铁路列车运行仿真研究[J].计算机仿真. 2022,39(08)
[3]张睿洋.元胞自动机在两类传染病模型中的应用[J].现代信息科技. 2022,6(10)

3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用

生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络(CNN)、LSTM、支持向量机(SVM)、最小二乘支持向量机(LSSVM)、极限学习机(ELM)、核极限学习机(KELM)、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题(TSP)、车辆路径问题(VRP、MVRP、CVRP、VRPTW等)、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

海神之光

有机会获得赠送范围1份代码

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值