关于对《弹道导弹攻防对抗的建模与仿真》一书中部分章节模型的Python复现

Attack-Defense

Attack-Defense是对《弹道导弹攻防对抗的建模与仿真》（Attack-Defense Counterwork Modeling and Simulation of Ballistic Missile）一书中部分章节模型的Python复现。

本书作者：罗小明，出版社：国防工业出版社，文章内容属非涉密范畴，原书可在淘宝、京东、咸鱼等多家APP上购买查阅。

《弹道导弹攻防对抗的建模与仿真》主要是为弹道导弹攻防作战效能评估、体系效能评估、体系优化提供基础数学模型和计算方法，并给出了建模和仿真的基础知识和总体思路，力求反映或揭示卫星系统信息支援下弹道导弹攻防作战的一些客观规律，定量分析满足不同目标要求时所需要的各种弹道导弹数量，以及用于打击各种目标的弹道导弹的种类和数量。

若模型的建立存在问题或错误，请在Issues、 Blog(https://blog.csdn.net/weixin_40485904) 或WeChat(wwy18795980897) 上联系我，我将在第一时间进行整改维护，感谢！

1.战术弹道导弹武器系统作战效能评估的系统动力学模型

使用战术弹道导弹武器系统作战效能评估的系统动力学模型的Python复现.py脚本，该脚本是对书中第8章《战术弹道导弹武器系统作战效能评估的系统动力学模型》的复现。

系统的效能依赖于系统的结构，其效能的发挥以及增长又需要通过系统内、外部信息、物质流的交互来实现。因卫星信息支援下战术弹道导弹武器系统作战效能评估受到许多环节和因素的影响，且相互作用复杂，属于开放式反馈系统，要描述这样的系统需要多个状态变量，反映多个变量因果关系的模型必定是高阶次的非线性问题；多个决策者、多种功能又使得系统有多个反馈回路。这里主要结合作者观点，针对影响弹道导弹系统动力学模型建立DYNAMO方程式，根据用户输入，分析各类因素对系统动力学模型的影响情况，并绘制成图。

相关流程：

下图是把因果关系转换为DYNAMO语言的流程图，亦是对输入量和各因素之间相互关系所进行的量化。

此模型的难点在于DYNAMO方程式的建立，变量庞大，幂次较高，且相互间均有影响，如何妥善区分各变量之间的关系、变量与时间的关系、变量与常量之间的关系是解决模型的关键所在。其中，流位变量和流率变量计算难度大，书中方程式仅提供初始值和与其他变量间的关系，非计算机模型难以准确的计算出相应结果。此外，涉及比例关系的K系列值，文中均未提到具体范围，可能会对模型的结构造成影响。

仿真效果（K值不贴合模型，导致出现负数）：

2.有突防措施条件下多枚战略弹道导弹突防反导防御系统体系效能模型

使用有突防措施条件下多枚战略弹道导弹突防反导防御系统体系效能模型的Python复现.py脚本，该脚本是对书中第11章《有突防措施条件下多枚战略弹道导弹突防反导防御系统体系效能模型》的复现。

诱饵是一种战略弹道导弹最常使用的突防措施。这里主要结合作者观点，建立体系效能模型，定量研究多弹头和诱饵的影响下，多枚战略弹道导弹在反导防御系统采取不同的拦截方式和拦截策略情况下的突防效能，并从防御方的角度，建模分析了反导防御方要达到一定的防御效能所需的拦截弹数。在建模之前，给出以下假设条件：

(1)各枚战略弹道导弹弹头之间相互无影响。

(2)反导防御系统在某一段拦截时，若采取相同的拦截策略且在该段拦截弹未消耗完的条件下，采取相同突防措施的各枚战略弹道导弹弹头在该段的突防概率相同。

相关流程：

导弹进入助推段、上升段、中段和末段的模拟流程图如下图所示。

与第8章相比，公式幂次有所降低，但大量的公式、多变的因素、复杂的流程是脚本编写时最大的难点，无法使用函数、类等方法简化模型，且要考虑多种情况下子模型间的参数传递和数值干扰，确实是一个不小的挑战。此外，书中对中段路基拦截系统和中段海基拦截系统区分不大，为避免给使用者带来疑问或困扰，模型中提到的中段拦截系统均为路基。

仿真效果：

项目地址：

GitHub：https://github.com/wangwei39120157028/Attack-Defense

Gitee：https://gitee.com/wwy2018/Attack-Defense

欢迎start！！！