简介这是一份基于MATLAB编写的演化博弈仿真代码包面向博弈论初学者、生物与社会经济模型研究者和MATLAB仿真爱好者。资源围绕X与Y两种策略在群体中的动态演化过程展开通过复制动态或Fermi规则等机制模拟策略更新并利用绘图函数输出频率变化与稳定状态。包内共4个m文件均为可直接运行的MATLAB脚本体积仅1KB代码简短精炼适合用于教学演示或二次修改。当前已有2846人学习下载。借助该代码读者可以直观观察不同博弈矩阵和初始策略分布下合作的涌现或竞争的消亡加深对演化稳定策略(ESS)等核心概念的理解同时也能锻炼MATLAB编程与结果分析能力。 演化博弈代码2.zip这个压缩包是我前后反复改了两轮才整理出来的。里面装的不只是几十个.py和.ipynb文件而是一整套做演化博弈仿真时常用的基础库、示例脚本、参数配置还有结果可视化模板。演化博弈evolutionary game theory很多人一开始觉得门槛高其实自己动手跑一遍代码比看十篇综述都管用。这个 zip 包适合刚接触演化博弈的学生、要复现论文的实验者以及想快速在合作、背叛、惩罚这些经典场景里搭一套可调参数的仿真框架的人。先说一个我踩过的经历这个包第一次发给别人时对方反馈“解压就报错failed to copy spatial iop zip”我一开始还以为是代码问题后才发现是压缩包在传输过程中损坏了。类似这种 zip 相关的问题其实比博弈模型本身还劝退新手。所以这次我不光整理代码也会把从解压、加密、分卷到映射 Git 仓库这些最容易翻车的环节一起讲清楚。1. 项目整体思路:这套代码到底在解决什么问题1.1 演化博弈到底跑的是什么演化博弈不是传统博弈论那种“理性人算出最优策略”的静态分析而是把策略看成生物种群里的个体让它们在多轮交互中优胜劣汰。核心就是三个概念种群、策略频率、适应度。每一轮随机配对博弈收益高的策略在群体里占比上升最终可能稳定在一个演化稳定策略ESS上也可能出现周期振荡甚至混沌。我的代码里默认搭好了两类基础模型一类是复制者动态replicator dynamics适合分析无限大群体用一组常微分方程描述频率变化另一类是有限群体的随机过程比如 Moran 过程群体大小固定每次随机选一个个体复制再用另一个个体被淘汰。这两类模型是演化博弈的两大支柱论文里 80% 的场景都能套进去。1.2 为什么是“2”而不是“1”第一版代码功能上没毛病但结构非常散所有脚本堆在一个目录里参数全部写死在源码里想换个收益矩阵得把代码翻个底朝天。这次重写有几个明确目标把博弈矩阵、群体规模、迭代次数、突变率全部抽到配置文件里改参数不用再动源码。把核心算法和具体场景解耦新增一个博弈场景只需要在games/目录加一个类。增加可视化模块实时输出策略频率变化图和空间博弈的网格快照。顺带补了requirements.txt和完整的README新机器上从零跑通不超过十分钟。1.3 压缩包里的内容结构拿到 zip 后建议别急着双击运行先建个文件夹把它完整解压出来完整的目录结构如下演化博弈代码2/ ├── README.md ├── requirements.txt ├── config/ │ ├── pd_config.yaml │ └── moran_config.yaml ├── src/ │ ├── __init__.py │ ├── games/ │ │ ├── base_game.py │ │ ├── prisoner_dilemma.py │ │ └── public_goods.py │ ├── dynamics/ │ │ ├── replicator.py │ │ ├── moran.py │ │ └── spatial.py │ └── utils/ │ ├── payoff.py │ └── plotter.py ├── examples/ │ ├── run_replicator.py │ ├── run_moran.py │ └── run_spatial_pd.py └── data/ └── record_output/第一版的问题就是没有config和utils层所有东西混在一起改一个展示逻辑就可能把仿真主循环带到沟里去。现在每一层只干一件事games负责定义策略和收益dynamics负责演化规则utils负责画图和数据处理。边界清晰之后调试效率提升非常明显这也是我强烈建议所有仿真项目都这么拆的原因。2. 解压前的功课别把时间耗在 zip 的坑里我见过太多人卡在第一步github 上明明写了“Download ZIP”下载回来却怎么都解压不了。这里包含很多与 zip 本身相关的常见问题大家一定要留意。zip 不是代码它只是分发形态你把它的属性搞清楚后面的实验才会一路顺畅。2.1 解压报错“could not find EOCD”是什么鬼EOCD 是 zip 压缩包的结尾目录记录用来告诉解压工具这个压缩包到底有哪些文件、压缩到哪个位置。如果解压时提示invalid zip archive: could not find eocd基本可以断定这个 zip 没有被完整下载或者文件被人为截断过。最简单粗暴的办法是用ls -l检查文件大小再和源网站上的文件大小比对。如果对不上删掉重下不要花时间修复。遇到确实已经下载完成但还是报错的情况可以尝试用压缩工具的修复功能。Linux/macOS 下执行zip -FF damaged.zip --out repaired.zip它会尽可能读取完整文件生成一个新压缩包Windows 上用 7-Zip 可以打开损坏 zip 后选择“修复”。但请记住这种修复只能救回部分文件源码里缺了任何一个.py文件都可能让整个项目跑不起来。2.2 zip 密码加密与“移除密码”的边界代码项目很少有人加密但如果你拿到别人转发的 zip 时恰好有密码这事就有点烦了。zip 压缩包里常见的加密方法有两种传统 ZipCrypto 和 AES 加密。ZipCrypto 强度比较弱网上能搜到很多 zip 密码恢复工具比如用 hashcat 或 john 工具跑字典策略是先把 zip 头信息导出成 hash再拿到爆破工具里跑。不过作为分享代码的人我真心建议不要给代码包加密码尤其是当你希望别人顺利复现实验时。你永远无法预知对方会不会因为忘记密码或不会解密工具干脆放弃你的代码。如果你自己不小心忘记密码也可以试试“百事牛zip密码恢复工具”这类图形化软件但密码复杂度高的话就别抱太大希望。更稳妥的做法是直接找原始分享者要一份未加密的副本。2.3 分卷压缩包 .z01 没有 zip 怎么办分卷压缩常出现在网盘传大文件时比如xxx.zip、xxx.z01、xxx.z02。很多人只下载了.z01发现后缀不对咋都打不开。注意.z01只是第一个分卷它依赖最终的.zip分卷来提供索引信息。完整下载所有分卷之后用 7-Zip 打开.zip那个分卷它会自动识别同目录下的.z01/.z02直接解压即可。如果是 Linux 环境可以用命令行方式合并分卷zip -s 0 演化博弈代码2.zip --out 全量代码.zip这里的-s 0表示把分卷合并成一个完整 zip 再输出。处理完后再用unzip 全量代码.zip正常解压。2.4 GitHub 下载的 zip 怎么和 Git 仓库关联下载 zip 和git clone最大的区别是zip 不包含.git目录所以它只是一个普通文件夹。有人想基于 GitHub 下载的 zip 继续开发于是先解压、再运行git init然后git remote add origin ...结果执行git pull或git rebase origin/main时频繁失败。原因很简单你还没有建立本地分支和远程分支的追踪关系或者远程仓库里已经存在大量完全不同的提交历史。我的推荐流程是这样# 1. 解压后进入目录 cd 演化博弈代码2 # 2. 初始化并关联远程 git init git remote add origin gitgithub.com:yourname/repo.git # 3. 先拉取远程但用 --allow-unrelated-histories git pull origin main --allow-unrelated-histories如果远程仓库本来就有代码合并且产生冲突是很正常的。我一般会先查看git status把冲突文件打开后手动保留需要的版本再提交。注意千万不要在没有任何远程代码的情况下用git push --force它会直接覆盖远端。3. 运行环境准备与依赖安装代码归类 zip 没问题但解压只是前菜真正能否跑起来要看运行环境。3.1 用虚拟环境隔离依赖我强烈建议先建一个独立环境不要让系统级 Python 环境被各种版本的库污染。conda 和 venv 二选一即可我用的是 venvpython -m venv evogame_env source evogame_env/bin/activate pip install --upgrade pip pip install -r requirements.txtrequirements.txt里锁定了numpy1.24.3、scipy1.10.1、matplotlib3.7.1。这几个版本在 Python 3.9~3.11 下都稳定。如果你机器上 Python 版本偏新或偏旧建议先看 README 的兼容说明再决定是否微调版本号。3.2 解压后导入模块报错怎么定位有人解压后直接双击运行run_replicator.py结果出现ModuleNotFoundError: No module named src。这是因为脚本的执行路径不在项目根目录Python 找不到src包。正确做法是确保当前工作目录在项目根目录再运行python examples/run_replicator.py如果还是报错就把执行路径手动加进去或者在脚本开头写一段路径处理逻辑import os, sys ROOT os.path.dirname(os.path.dirname(os.path.abspath(__file__))) if ROOT not in sys.path: sys.path.insert(0, ROOT)这个写法在小项目里很实用不至于每次都要配置 PYTHONPATH。4. 核心模块拆解与代码实现解析4.1 复制者动态无限大群体的主力算法复制者动态的核心公式是策略i的频率增长率等于当前收益和种群平均收益的差。我代码里的replicate函数这样实现def replicator_step(freqs, payoff_matrix): fitness payoff_matrix freqs avg_fitness np.dot(freqs, fitness) return freqs * (fitness - avg_fitness)这段代码的精髓就是那个freqs * (fitness - avg_fitness)。如果某个策略的收益超过平均值它的占比就上升低于平均值就下降。默认的囚徒困境收益矩阵设定为Temptation T 1.6 Reward R 1.0 Punishment P 0.0 Sucker S 0.0这个取法保证 T R P ≥ S所以纯理性下背叛是优势策略。但你在复现时会发现如果初始合作者比例足够高合作策略可以在一段时间内维持甚至在某些迭代次数内占优。这就是演化博弈最有意思的地方演化路径依赖初始频率而不一定收敛到 Nash 均衡。4.2 Moran 过程有限群体里的随机性复制者动态里没有种群规模这个概念它是无限大且确定性的。真实实验和真实社会里群体一定是有限的此时随机漂变不能忽略。Moran 过程每步只更新一个个体频率变化天生带噪声因此结果必须重复多次取统计学上的平均才能判断某个策略是否真的更容易被选中。我的run_moran.py示例里会模拟一个 100 个个体的种群突变率设为 0.01每轮随机挑选一个个体按其适应度比例复制再随机淘汰一个个体。注意这里的适应度通常写成1 - w w * payoffw是选择强度。当w 0时演化完全是中性漂变w越大收益对演化结果的影响越强。经验值是w 0.1~0.3最适合观察策略竞争的中期动态太大会让系统迅速锁定太小则噪声完全淹没信号。4.3 空间演化博弈加一张网格就完全不一样把随机配对博弈放到二维网格上策略只和邻居交互整个系统行为会发生质的改变。比如在方格网络里进行空间囚徒困境合作者可以通过形成簇群来自我保护即使背叛者收益很高也无法轻易侵入合作者抱团形成的“安全区”。空间博弈的每步循环逻辑是for i in range(rows): for j in range(cols): opponent choose_random_neighbor(i, j) if effective_payoff[i, j] effective_payoff[opponent]: strategy_grid[i, j] strategy_grid[opponent]这个简单的“收益低就模仿邻居”规则其实就是一类局部更新规则。实现时注意边界处理我代码里默认使用周期性边界条件也就是网格上下左右循环连接。如果你改成固定边界结果会有明显差异这属于正常现象不是 bug。4.4 可视化把仿真结果变成能进论文的图仿真代码不配一张好图价值直接打对折。我的utils/plotter.py里做了两个主要功能一是画策略频率随时间变化的曲线二是把空间博弈每一步的网格状态输出为 PNG。频率曲线用 Matplotlib 的折线图即可网格快照用imshow加自定义 colormap。在保存图片时有一个很容易踩的坑如果你在循环里连续保存 200 张网格图目录名和文件名一定要格式化补零否则排序会乱。建议用snapshot_%03d.png这种格式这样可以用 FFmpeg 直接把图片序列转成视频做动态演示非常方便。5. 参数设置与实验设计经验5.1 不同模型参数速查表为了方便大家抄作业我整理了一个参数经验表。这些数值不是拍脑袋定的而是权衡了计算量和结果稳定性之后比较合适的默认值。模型群体规模/网格迭代次数关键参数建议取值复制者动态无限大数值积分1000初始频率0.1~0.9 扫参Moran 过程10050000选择强度 w0.1~0.3空间囚徒困境100x1003000邻居半径4 邻或 8 邻重复实验次数上Moran 过程建议至少 50 次独立运行然后取均值空间博弈则依赖网格大小100x100 网格跑 3000 代在普通笔记本上大约需要一到两分钟可以接受。5.2 我的三个实操心得第一肉眼观察策略频率曲线很容易误判收敛。最好写一个滑动窗口统计比如最后 500 代里频率波动小于 0.01才判定为收敛而不是简单看到曲线“好像不涨了”就停止。第二参数扫描一定要避免单点调参。如果你要对某个参数画相变图保证其它参数不变的前提下一口气扫 20 个点以上而不是只试两三个值就开始下结论否则很容易得出局部偶然结果。第三随机数种子务必要固定。演化博弈是随机过程同一个参数跑两次结果不可能完全一样。代码里我默认在main函数开头加一句np.random.seed(42)这既是实验可复现性的底线也是调试时的救星至少报错后你能排除随机性带来的干扰。5.3 性能优化无从下手时怎么办如果网格规模大迭代次数多Python 的纯 for 循环很容易慢到怀疑人生。我的建议是先用numba给更新函数加jit核心循环就能提速几倍到十几倍。如果不想引入新依赖也可以用numpy的向量化操作替代嵌套循环比如计算邻居收益时用切片直接提取周围区域而不是每个位置单独访问数组。我自己实际测试过100x100 网格跑 3000 代用纯 Python 大约需要 90 秒用 numba 的 nopython 模式优化后只要 5 秒。值不值得加库不需要纠结直接说结论值得。6. 常见问题与排查技巧实录我把自己做这套代码以来遇到的高频问题整理成一个速查表避免大家在同样的坑里反复折腾。现象原因解决办法解压提示 could not find eocdzip 文件不完整或传输中断重新下载或用 zip -FF 修复zip 有密码解不开密码或加密方式不匹配联系原始提供者或尝试 hashcat 字典爆破有 .z01 但没 .zip 分卷分卷下载不全补齐所有分卷用 7-Zip 打开 .zipGitHub 下载 zip 后 rebase 失败本地无远程提交历史用 git pull --allow-unrelated-historiesModuleNotFoundError: src运行路径不在项目根目录cd 到项目根目录再执行曲线锯齿严重无法判断趋势单次随机过程噪声大多次独立运行求平均加滑动窗口内存被网格快照占满保存了过多高分辨率图片降低迭代保存频率或改输出 GIF这里我要专门提一句zip 包本身不是代码运行的一部分但它直接影响代码能不能顺利“复活”。我后来发新版本时学乖了每次打包前先写一个校验脚本生成所有文件的 SHA256 校验值发布后如果有人解压报错先让他核对校验值省去大量沟通成本。再来聊一个容易让人忽视的问题解压后的目录名不要带空格或中文冒号。某些科学计算库在解析路径时对特殊字符非常敏感你代码写得再好放在C:\桌面\最终版(2)\演化博弈代码2这种路径下也可能莫名报错。尽量把项目解压到一个不含中文和空格的纯英文路径下比如C:\work\evogame2能省很多事。7. 从 zip 到论文用这套代码能做什么如果你正在写课程论文或准备复现一篇顶刊的博弈论模拟结果这套 zip 里的代码应该能帮你快速起步。一个比较顺畅的流程是先在配置里定义好收益矩阵然后跑一遍复制者动态看相位图走势再用 Moran 过程验证有限群体下的随机效应最后把空间网格打开观察合作者在空间上的聚集模式。三种模型互相印证论文里可以写出很有说服力的结论段。我在实际使用中还发现把演化过程做成视频非常有利于答辩展示。代码里utils模块提供了保存逐帧图片的函数你跑完空间博弈后可以直接用 FFmpeg 把它们合成为视频ffmpeg -framerate 10 -i snapshot_%03d.png -pix_fmt yuv420p evolution.mp4这个命令生成的 mp4 可以在任何播放器上打开嵌进 PPT 也不会卡顿。视频里能直观看到合作者簇群的生长、分裂和扩散比单纯贴一张最终状态图直观太多。最后再多说一句关于 zip 本身的运维经验。经常折腾代码包的人建议准备一条常用命令在 Linux 或 macOS 下递归压缩目录时加上排除选项会更好zip -r 演化博弈代码2.zip 演化博弈代码2/ -x *__pycache__* -x *.git*这样可以避免把 Python 缓存和 git 历史文件卷进压缩包包体更小也不会给解压者带来冗余文件。另外压缩完成后用unzip -l看一眼文件列表确保关键脚本都在里面再发给别人。这个习惯救了我很多次至少避免过两三次把空目录或者缺少requirements.txt的半成品发出去的尴尬。本文还有配套的精品资源点击获取