如何在Matlab中在正方形格子中数值模拟鲁棒拓扑边状态动画？

拓扑物理学通过承载无序免疫鲁棒边态传播，在光子学中实现了有趣的可能性。我想知道如何在Matlab中数值模拟方形晶格光子系统中的无序免疫边态传播，将边态传播显示为一个模拟电影。例如，参见本文。我的问题是，我们如何通过Matlab编码在电影模拟中数值模拟方形晶格中的散射自由边状态行为，如图2所示？ 提前谢谢