基于Netlogo建模基础知识讲解

Netlogo建模基础知识讲解14302010008秦奕深NetLogo简介1.基本情况NetLogo是一个用来对自然和社会现象进行仿真的可编程建模环境建模仿真集成环境由美国西北大学连接学习与计算机建模中心（CenterforConnectedLearningandComputer-BasedModeling，CCL）开发。2002年发布了1.0版本，用户手册为4.0.2版本，最新为4.1.1版本2.主要功能多主体建模（重点）多个移动Agent分布在二维空间中，每个Agent自主行动，所有主体并行异步更新，整个系统随着时间推进而动态变化。运行控制仿真输出提供了多种手段实现仿真运行监视和结果输出实验管理BahaviorSpace，自动管理仿真运行，并记录结果。系统动力学仿真参与式仿真HubNet模型库3.软件特色完全可编程简单语言结构是LOGO语言的扩展，支持智能体和网络结构可以定义无限个智能体和变量多种内置命令帮助使用支持整型和双精度型浮点数计算跨平台、可复用的应用。二、NetLogo仿真框架1.NetLogo软件2.模型的抽象总体：大量的可移动主体在二维空间中交互作用，随着时间推进，微观个体的属性不断发生变化，系统的宏观特征也因此而变化。从三个方面理解：主体空间表达仿真推进（1）主体（Agents）虚拟世界由主体构成，主体能够接受命令，进行活动，所有主体的行为并行发生。NetLogo中共有三类主体，turtles（海龟）patches（瓦片）observer（观察者）主体类型turtles指能够在世界中移动的主体。世界是二维的，划分为由patches组成的网格，每个patch占据一个矩形小块。patch不能移动patch也是主体patch和turtle一样可以有自己的属性和行为observer是一个全局主体，它观察着由turtles和patches构成的世界，能够执行指令获取世界全部或部分的状态，或实现对世界的控制。虚拟世界ObserverPatchTurtle（2）空间表达每个patch有二维坐标（pxcor,pycor）坐标值为整数。默认情况下，二维世界的水平、垂直坐标范围为（-17，17）每个turtle也有坐标（xcor,ycor）turtle坐标不必是整数，因此turtle不一定正好位于某个patch的中心。一个patch上也可以同时有多个turtles。实际上对turtle而言，NetLogo的空间是连续的。（3）仿真推进没有明确的仿真时钟变量，也没有提供特定的事件处理机制仿真推进是通过不断重复执行某个例程实现的模型中至少要有初始化例程和仿真执行例程初始化例程实现对模型初始状态的设置，生成所需的turtles，设置其状态，以及其它工作。仿真的执行通过例程go实现，在go例程中编写所需执行的各种指令，完成一个仿真步的工作。需要在Interface页中建立一个按钮与go例程相联系，该按钮是一个永久（forever)按钮，点击后将不断重复执行go例程，直到遇到stop指令或用户再次点击该按钮则仿真终止。2.建模基本过程NetLogo模型包括可视化部件和例程两部分，二者具有紧密联系。先在Interface中创建可视化控件，然后在Procedures中实现相应的代码，通过设置控件的属性将二者联系起来。Interface中主要有三类部件运行控制参数控制仿真显示Procedure中的例程分为两类：命令（command）例程报告（reporter）例程仿真的基本框架初始化tosetup每一个仿真周期togo是否结束？结束对所有Turtle循环askturtles每个Turtle做出决策forward1是否三、编程指南breed[]turtles-own[]globals[]tosetup[]endtogo[]endaskturtles[setpcolorwhite]set-default-shapeturtlespersonsetxyrandom-xcorrandom-ycorsetturtles-ownrandomsetage(22+random53)setage(age+1)sethouse0clear-all(ca)forward(fd)1create-turles100sprout-breedsnumber[]sproutnumber[]hatchnumber[]hatch-breedsnumber[]askone-ofmen[setbreed]ififelse四、生态系统建模示例假设要模拟一个简单的生态系统，该系统中有一种生物以青草为食，通过吃草获取能量、该类生物经历成长、繁殖、死亡过程。1.初始化创建生物群体，并将它们随机分布在空间中。例程（demo1）tosetup;;定义例程setupclear-all;;设置整个世界为初始状态;;创建100个turtles，创建后各turtle默认坐标是（0，0）create-turtles100;;命令所有turtle执行语句setxyrandom-xcorrandom-ycor;;各turtle坐标随机产生，实现turtles在空间中的随机分布askturtles[setxyrandom-xcorrandom-ycor]end2.仿真执行例程实现turtle随机移动实现go例程，与forever按钮联系go调用move-turtles例程增加代码:(demo2)togo;;定义仿真执行例程gomove-turtles;;调用例程move-turtlesendtomove-turtles;;定义例程move-turtles;;所有turtle执行[]中的命令序列askturtles[rightrandom360;;右转一个角度，度数随机产生forward1;;前进距离1]end3.青草的模拟为模拟青草的存在设置patches为绿色改写setup例程，其中调用了两个新的例程setup-turtles和setup-patches，分别设置turtle和patch的初始状态demo3tosetupclear-allsetup-patches;;调用例程setup-patchessetup-turtles;;调用例程setup-turtlesendtosetup-patches;;命令所有patches执行指令setpcolorgreen;;该指令将patch颜色设置为绿色askpatches[setpcolorgreen]endtosetup-turtlescreate-turtles100askturtles[setxyrandom-xcorrandom-ycor]end4.主体行为现在增加一些行为，假设turtle以青草为食，通过吃草获得能量，移动时要消耗能量。patch代表青草，绿色表示有，黑色表示无。为实现这样的模型，需要重新定义turtles的属性和行为，然后重写go例程。首先给turtle增加变量energy以存储当前能量值，另外添加吃草获取能量的例程，还要修改移动例程以反映能量消耗。自定义turtle变量，修改goturtles-own[energy];;声明turtle变量energytogomove-turtles;;turtle随机移动，消耗能量eat-grass;;吃草获取能量end定义eat-grassturtle吃草获取能量修改patch颜色表示草的有无toeat-grass;;如果turtle所在patch颜色为绿色，表示有草，则吃草，;;令该patch颜色变为黑色，表示已无草，然后自身能量增加10askturtles[ifpcolor=green[setpcolorblacksetenergy(energy+10)]]end修改原来的move-turtles例程，添加能量消耗指令。tomove-turtlesaskturtles[rightrandom360forward1setenergyenergy-1;;移动后能量减少1]end假设turtle能量小于等于0就死亡，当能量大于50就繁殖；青草以一定的恢复率再生。由于主体行为发生了改变，需重新定义go例程，togomove-turtles;;移动eat-grass;;吃草reproduce;;繁殖check-death;;死亡regrow-grass;;青草再生endtoreproduceaskturtles[ifenergy50[;;如果能量大于50则繁殖setenergyenergy-50;;母体能量减少50hatch1[setenergy50];;产生一个后代，初始能量50]]endtocheck-deathaskturtles[ifenergy=0[die];;如果能量小于等于0则死亡]endtoregrow-grassaskpatches[;;青草以0.03的概率再生ifrandom1003[setpcolorgreen]]end5.监视仿真运行监视器控件(monitor)显示turtle数量，monitor表达式“countturtles”主体标签每个主体有一个变量lable6.图形输出在Interface中创建Plot控件，在Procedures中编制绘图例程。绘图概念每个Plot控件必须指定一个唯一名，在绘图时通过Plot名指定在哪个Plot上绘图绘图时必须使用某个画笔，画笔默认是黑色实线，也可以创建自定义画笔绘制turtles数量和青草数量曲线添加Plot控件创建两个画笔，分别指定颜色添加do-plots例程修改setup和go例程tosetupclear-allsetup-patchessetup-turtlesdo-plots;;绘制初始点endtogomove-turtleseat-grassreproducecheck-deathregrow-grassdo-plots;;随仿真运行绘制图形endtodo-plotsset-current-plot“totals”;;选定plot控件set-current-plot-pen“turtles”;;选定画笔turtlesplotcountturtles;;绘制turtles数量set-current-plot-pen“grass”;;选定画笔grassplotcountpatcheswith[pcolor=green];;绘制青草数量end7.仿真参数控制参数控件开关(switch)、滑动条(slider)、选择器(chooser)这些控件都对应一个全局变量,在Procedures中将这些全局变量作为参数使用在程序中，就能实现仿真参数的控制。例如：选择仿真时长用滑动条设置仿真时长，对应变量名terminate-time添加全局变量ticks修改go例程,随着仿真修改ticks当ticks=terminate-time时停止代码globals[ticks];;声明全局变量ticks…togoifticks=terminate-time[stop];;判断是否应停止move-turtleseat-grassreproducecheck-deathregrow-grassdo-plotssetticksticks+1;;时钟推进end