多主体仿真平台

多主体仿真平台第1页，课件共44页，创作于2023年2月一、NetLogo简介1.基本情况NetLogo是一个多主体建模仿真集成环境由美国西北大学连接学习与计算机建模中心（CenterforConnectedLearningandComputer-BasedModeling，CCL）开发。2002年发布了1.0版本，最新4.0第2页，课件共44页，创作于2023年2月2.主要功能多主体建模多个移动Agent分布在二维空间中，每个Agent自主行动，所有主体并行异步更新，整个系统随着时间推进而动态变化。运行控制仿真输出提供了多种手段实现仿真运行监视和结果输出实验管理BahaviorSpace，自动管理仿真运行，并记录结果。系统动力学仿真参与式仿真HubNet模型库第3页，课件共44页，创作于2023年2月二、NetLogo仿真框架1.NetLogo软件界面(见软件)：模型部分 三个TabPage

命令行窗口

第4页，课件共44页，创作于2023年2月2.模型的抽象总体：大量的可移动主体在二维空间中交互作用，随着时间推进，微观个体的属性不断发生变化，系统的宏观特征也因此而变化。从三个方面理解：主体空间表达仿真推进第5页，课件共44页，创作于2023年2月（1）主体虚拟世界由主体构成，主体能够接受命令，进行活动，所有主体的行为并行发生。NetLogo中共有三类主体，turtles（小海龟）patches（瓦片）observer（观察者）在4.0中将Links作为一类主体

第6页，课件共44页，创作于2023年2月主体类型turtles指能够在世界中移动的主体。世界是二维的，划分为由patches组成的网格，每个patch占据一个矩形小块。patch不能移动patch也是主体patch和turtle一样可以有自己的属性和行为observer是一个全局主体，它观察着由turtles和patches构成的世界，能够执行指令获取世界全部或部分的状态，或实现对世界的控制。第7页，课件共44页，创作于2023年2月虚拟世界第8页，课件共44页，创作于2023年2月（2）空间表达每个patch有二维坐标（pxcor,pycor）坐标值为整数。默认情况下，二维世界的水平、垂直坐标范围为（-17，17）每个turtle也有坐标（xcor,ycor）turtle坐标不必是整数，因此turtle不一定正好位于某个patch的中心。一个patch上也可以同时有多个turtles。实际上对turtle而言，NetLogo的空间是连续的。第9页，课件共44页，创作于2023年2月拓扑结构根据在水平和垂直方向边界是否进行回绕，形成四种结构环面（torus），默认结构盒子（box）垂直柱面（verticalcylinder）水平柱面（horizontalcylinder）回绕影响主体移动时是否穿越边界软件操作，在View上点击Edit第10页，课件共44页，创作于2023年2月（3）仿真推进没有明确的仿真时钟变量，也没有提供特定的事件处理机制仿真推进是通过不断重复执行某个例程实现的模型中至少要有初始化例程和仿真执行例程初始化例程实现对模型初始状态的设置，生成所需的turtles，设置其状态，以及其它工作。仿真的执行通过例程go实现，在go例程中编写所需执行的各种指令，完成一个仿真步的工作。需要在Interface页中建立一个按钮与go例程相联系，该按钮是一个永久（forever)按钮，点击后将不断重复执行go例程，直到遇到stop指令或用户再次点击该按钮则仿真终止。第11页，课件共44页，创作于2023年2月2.建模基本过程NetLogo模型包括可视化部件和例程两部分，二者具有紧密联系。先在Interface中创建可视化控件，然后在Procedures中实现相应的代码，通过设置控件的属性将二者联系起来。第12页，课件共44页，创作于2023年2月Interface中主要有三类部件运行控制参数控制仿真显示Procedure中的例程分为两类：命令（command）例程报告（reporter）例程第13页，课件共44页，创作于2023年2月基本过程初始化确定活动顺序定义主体属性与行为仿真过程监视图形输出仿真参数控制第14页，课件共44页，创作于2023年2月三、生态系统建模示例假设要模拟一个简单的生态系统，该系统中有一种生物以青草为食，通过吃草获取能量、该类生物经历成长、繁殖、死亡过程。第15页，课件共44页，创作于2023年2月1.初始化创建生物群体，并将它们随机分布在空间中。例程（demo1）tosetup;;定义例程setupclear-all;;设置整个世界为初始状态

;;创建100个turtles，创建后各turtle默认坐标是（0，0）

create-turtles100;;命令所有turtle执行语句setxyrandom-xcorrandom-ycor;;各turtle坐标随机产生，实现turtles在空间中的随机分布

askturtles[setxyrandom-xcorrandom-ycor]end第16页，课件共44页，创作于2023年2月2.仿真执行例程实现turtle随机移动实现go例程，与forever按钮联系go调用move-turtles例程第17页，课件共44页，创作于2023年2月增加代码:(demo2)togo;;定义仿真执行例程gomove-turtles;;调用例程move-turtlesendtomove-turtles;;定义例程move-turtles;;所有turtle执行[]中的命令序列

askturtles[rightrandom360;;右转一个角度，度数随机产生

forward1;;前进距离1]end第18页，课件共44页，创作于2023年2月3.青草的模拟为模拟青草的存在设置patches为绿色改写setup例程，其中调用了两个新的例程setup-turtles和setup-patches，分别设置turtle和patch的初始状态第19页，课件共44页，创作于2023年2月demo3tosetupclear-allsetup-patches;;调用例程setup-patchessetup-turtles;;调用例程setup-turtlesendtosetup-patches;;命令所有patches执行指令setpcolorgreen;;该指令将patch颜色设置为绿色

askpatches[setpcolorgreen]endtosetup-turtlescreate-turtles100askturtles[setxyrandom-xcorrandom-ycor]end第20页，课件共44页，创作于2023年2月4.主体行为现在增加一些行为，假设turtle以青草为食，通过吃草获得能量，移动时要消耗能量。patch代表青草，绿色表示有，黑色表示无。为实现这样的模型，需要重新定义turtles的属性和行为，然后重写go例程。首先给turtle增加变量energy以存储当前能量值，另外添加吃草获取能量的例程，还要修改移动例程以反映能量消耗。第21页，课件共44页，创作于2023年2月自定义turtle变量，修改goturtles-own[energy];;声明turtle变量energytogomove-turtles;;turtle随机移动，消耗能量

eat-grass;;吃草获取能量end第22页，课件共44页，创作于2023年2月定义eatgrassturtle吃草获取能量修改patch颜色表示草的有无toeat-grass;;如果turtle所在patch颜色为绿色，表示有草，则吃草，

;;令该patch颜色变为黑色，表示已无草，然后自身能量增加10askturtles[ifpcolor=green[setpcolorblacksetenergy(energy+10)]]end第23页，课件共44页，创作于2023年2月修改原来的move-turtles例程，添加能量消耗指令。tomove-turtlesaskturtles[rightrandom360forward1setenergyenergy-1;;移动后能量减少1]end第24页，课件共44页，创作于2023年2月假设turtle能量小于等于0就死亡，当能量大于50就繁殖；青草以一定的恢复率再生。由于主体行为发生了改变，需重新定义go例程，togomove-turtles;;移动

eat-grass;;吃草

reproduce;;繁殖

check-death;;死亡

regrow-grass;;青草再生end第25页，课件共44页，创作于2023年2月toreproduceaskturtles[ifenergy>50[;;如果能量大于50则繁殖

setenergyenergy-50;;母体能量减少50hatch1[setenergy50];;产生一个后代，初始能量50]]endtocheck-deathaskturtles[ifenergy<=0[die];;如果能量小于等于0则死亡]endtoregrow-grassaskpatches[;;青草以0.03的概率再生

ifrandom100<3[setpcolorgreen]]end第26页，课件共44页，创作于2023年2月5.监视仿真运行监视器控件(monitor)显示turtle数量，monitor表达式“countturtles”主体标签每个主体有一个变量lable第27页，课件共44页，创作于2023年2月6.图形输出在Interface中创建Plot控件，在Procedures中编制绘图例程。绘图概念每个Plot控件必须指定一个唯一名，在绘图时通过Plot名指定在哪个Plot上绘图绘图时必须使用某个画笔，画笔默认是黑色实线，也可以创建自定义画笔 第28页，课件共44页，创作于2023年2月绘制turtles数量和青草数量曲线添加Plot控件创建两个画笔，分别指定颜色添加do-plots例程修改setup和go例程第29页，课件共44页，创作于2023年2月tosetupclear-allsetup-patchessetup-turtles

do-plots;;绘制初始点endtogomove-turtleseat-grassreproducecheck-deathregrow-grass

do-plots;;随仿真运行绘制图形endtodo-plotsset-current-plot“totals”;;选定plot控件

set-current-plot-pen“turtles”;;选定画笔turtlesplotcountturtles;;绘制turtles数量

set-current-plot-pen“grass”;;选定画笔grassplotcountpatcheswith[pcolor=green];;绘制青草数量end第30页，课件共44页，创作于2023年2月7.仿真参数控制参数控件开关(switch)、滑动条(slider)、选择器(chooser)这些控件都对应一个全局变量,在Procedures中将这些全局变量作为参数使用在程序中，就能实现仿真参数的控制。 第31页，课件共44页，创作于2023年2月例如：选择仿真时长用滑动条设置仿真时长，对应变量名terminate-time添加全局变量ticks修改go例程,随着仿真修改ticks当ticks=terminate-time时停止第32页，课件共44页，创作于2023年2月代码globals[ticks];;声明全局变量ticks…togo

ifticks>=terminate-time[stop];;判断是否应停止

move-turtleseat-grassreproducecheck-deathregrow-grassdo-plots

setticksticks+1;;时钟推进end第33页，课件共44页，创作于2023年2月四、建模技术1.访问邻域原语Neighbors，Moore邻域neighbors4，vonNeumann邻域<agentset>in-radius<number><agentset>at-points<list>演示第34页，课件共44页，创作于2023年2月2.主体交互T-P交互turtle能够直接访问所在之处的patch，对该patch的属性进行读写askturtles[setpcolorblue]turtle还可以利用空间相关操作获取所需的patches，然后对这些patches的属性进行读写askturtles[setpcolor-ofpatch-at10blue] 第35页，课件共44页，创作于2023年2月patch可以通过一些操作获取相关的turtle例如turtles-here就返回当前patch处的turtle集合。也可以通过空间相关操作获取相应patch上的turtles。例如”turtlesatdxdy”返回与当前patch相对距离（dx,dy）处的turtles集合。第36页，课件共44页，创作于2023年2月T-T交互实现T-T交互的第一步是得到目标turtle的句柄，然后进行操作。获取目标turtle句柄的常用方式有三种：随机选取、根据特定条件、空间相关。随机选取是指在特定agent集合中以随机方式选取一个或n个agent。选取一个agent的原语为one-of例如“setcolor-ofone-ofturtlesred”在所有turtles中随机选择一个turtle，将其颜色设为红色。第37页，课件共44页，创作于2023年2月获得目标agent或agentset的第二种方式是根据特定条件。方法是使用with原语，语法为agentsetwith[condition]，返回满足条件condition的agent集合。例如“turtleswith[color=red]”返回红色的turtle集合；第38页，课件共44页，创作于2023年2月获取目标agent或agent集合的第三种方式是运用空间相关操作。比较直接的一种是获取当前patch上的turtle集合，原语有turtles-here和other-turtles-here，二者的区别仅在于是否包含调用者自身。例如“askturtle0[askother-turtles-here[fd10]]”表示ID=0的turtle令处于同一patch上的其他turtle前进10。另一种方式是获取特定patch上的所有turtles。原语有<breeds>-atdxdy，<breeds>-onagentset等。例如“turtles-at10”返回右侧紧邻patch上的所有turtle，“turtles-onpatch-ahead”返回前方patch上的所有turtle。第39页，课件共44页，创作于2023年2月3.持久关系的建立基本方法是在个体中相互保持对方的引用，这样可以随时对对方进行操作。例如在个体之间建立简单的伙伴关系。假设还没有找到伙伴的个体随机移动，直到与另一个尚没有伙伴的个体建立伙伴关系，一旦建立伙伴关系后将长期保持。为实现这一模型，首先为turtle增加变量partner,用于保持对伙伴的引用，然后在移动时从相遇的一些无伙伴的turtle中选择一个，相互建立伙伴关系。第40页，课件共44页，创作于2023年2月4.多类