常微分方程模型及分析——Lanchester战斗模型

赭红色的锆石

已于 2024-07-23 20:10:39 修改

阅读量436

点赞数 6

文章标签：算法数学建模

于 2024-07-23 20:07:26 首次发布

本文链接：https://blog.csdn.net/Zircon_Xray/article/details/140644685

版权

来源是西电2024数学建模双创课中《常微分方程模型及分析》的作业。

题目：在第一次世界大战中，甲、乙两方战役的胜负，往往由双方兵力多少及战斗力强弱两个因素决定，其中，兵力因战斗及非战斗减员而减少，因增援而增加;战斗力与射击次数及命中率有关。请按照正规战争和游击战争两类情况应用常微分方程进行建模，分析并预测战役结局。

分析：根据Lanchester战斗模型，记甲、乙两方兵力分别为关于时间 $t$ 的函数 $x\left( t \right)$ 和 $y\left( t \right)$ ，并且连续可导。显然有“战斗力的变化率＝后勤补给率－自然损失率（非战斗减员）－对方的杀伤率（战斗减员）”。

常规战

常规战模型及其解

讨论常规战下的模型

$\left\{ \begin{gathered} x'(t) = P(t) - ax - by \hfill \\ y'(t) = Q(t) - cx - ey \hfill \\ \end{gathered} \right.$

其中a,e分别为甲和乙的自然损失率，其中b,c分别为乙和甲的杀伤率，P(t),Q(t)分别为甲和乙的后勤补给率。为简化，不妨假设P(t)=p,Q(t)=q分别为一常数。

下面简要描述一下求解以下常微分非自治线性方程组的过程：

$\left\{ {\begin{array}{*{20}{l}} {x'(t) = p - ax - by} \\ {y'(t) = q - cx - ey} \end{array}} \right.$

先求解对应的齐次方程组

$\left\{ {\begin{array}{*{20}{l}} {x'(t) = - ax - by} \\ {y'(t) = - cx - ey} \end{array}} \right.$

计算特征方程：

$\left| {\begin{array}{*{20}{c}} {a + \lambda }&b \\ c&{e + \lambda } \end{array}} \right| = 0$

展开行列式，得到特征方程：

${\lambda ^2} + (a + e)\lambda + (ae - bc) = 0$

使用求根公式求解二次方程：

${\lambda _{1,2}} = \frac{{ - (a + e) \pm \sqrt {{{(a + e)}^2} - 4(ae - bc)} }}{2}$

根据特征根，齐次方程组的通解为：

${X_h}(t) = {C_1}{e^{{\lambda _1}t}}\left[ {\begin{array}{*{20}{c}} 1 \\ {\frac{c}{{{\lambda _1} + a}}} \end{array}} \right] + {C_2}{e^{{\lambda _2}t}}\left[ {\begin{array}{*{20}{c}} 1 \\ {\frac{c}{{{\lambda _2} + a}}} \end{array}} \right]$

其中 $C_1$ 和 $C_2$ 是任意常数。

接着寻找非齐次方程组的特解，假设特解的形式为常数向量：

${X_p}(t) = \left[ {\begin{array}{*{20}{c}} A \\ B \end{array}} \right]$

代入非齐次方程组，得到：

$\left[ {\begin{array}{*{20}{c}} p \\ q \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {p - aA - bB} \\ {q - cA - eB} \end{array}} \right]$

解这个线性方程组，得到特解：

$\begin{gathered} A = \frac{{pb - qe}}{{ab + be + ce}} \hfill \\ B = \frac{{pc - qa}}{{ac + bc + ce}} \hfill \\ \end{gathered}$

构造非齐次方程组的通解，非齐次方程组的通解是齐次解与特解之和：

$X(t) = {X_h}(t) + {X_p}(t)$

因此，原方程组的通解为：

$X(t) = {C_1}{e^{{\lambda _1}t}}\left[ {\begin{array}{*{20}{c}} 1 \\ {\frac{c}{{{\lambda _1} + a}}} \end{array}} \right] + {C_2}{e^{{\lambda _2}t}}\left[ {\begin{array}{*{20}{c}} 1 \\ {\frac{c}{{{\lambda _2} + a}}} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {\frac{{pb - qe}}{{ab + be + ce}}} \\ {\frac{{pc - qa}}{{ac + bc + ce}}} \end{array}} \right]$

常规战的平方率

为分析射击次数及命中率对战局的影响，采取简化模型

$\left\{ \begin{gathered} x'(t) = - by \hfill \\ y'(t) = - cx \hfill \\ \end{gathered} \right.$

两式相除，得

$\frac{{dy}}{{dx}} = \frac{{cx}}{{by}}$

得

$b{y^2} - c{x^2} = by_0^2 - cx_0^2 = K$

有

${\left( {\frac{{{y_0}}}{{{x_0}}}} \right)^2}\begin{array}{*{20}{c}} {} \\ {} \end{array}\begin{array}{*{20}{c}} > \\ < \end{array}\begin{array}{*{20}{c}} {} \\ {} \end{array}\frac{c}{b}$

记射击率为 $r$ ，命中率为 $p$ ，则有

$\frac{c}{b} = \frac{{{r_x}{p_x}}}{{{r_y}{p_y}}}$

综合所有因素考虑，可以得到结论：常规战胜负取决于开战前力量(人数)对比，且此比值平方放大，具体如图：

游击战

游击战模型

$\left\{ \begin{gathered} x'(t) = P(t) - ax - gxy \hfill \\ y'(t) = Q(t) - ey - hxy \hfill \\ \end{gathered} \right.$

其中a,e分别为甲和乙的自然损失率，其中g,h分别为乙和甲的杀伤率，P(t),Q(t)分别为甲和乙的后勤补给率。

与常规战模型比较，最显著的特征是游击战的战斗伤亡率和自身部队战斗力正相关。

游击战的线性率

同样，只考虑战斗减员，简化模型为

$\left\{ \begin{gathered} x'(t) = - gxy \hfill \\ y'(t) = - hxy \hfill \\ \end{gathered} \right.$

相除，交叉相乘得

$g\left( {y - y_0^{}} \right) = h\left( {x - x_0^{}} \right)$

得

$gy - hx = gy_0^{} - hx_0^{} = L$

进一步有

$\frac{{{y_0}}}{{{x_0}}}\begin{array}{*{20}{c}} {} \\ {} \end{array}\begin{array}{*{20}{c}} > \\ < \end{array}\begin{array}{*{20}{c}} {} \\ {} \end{array}\frac{h}{g}$

其中 $h = {r_x}\frac{{{S_{rx}}}}{{{S_y}}},g = {r_y}\frac{{{S_{ry}}}}{{{S_x}}}$ ，里面 $r$ 为射击率、 ${S_r}$ 一次射击的有效面积、 $S$ 为游击队员的活动面积，即：

$\frac{{{S_y}{y_0}}}{{{S_x}{x_0}}}\begin{array}{*{20}{c}} {} \\ {} \end{array}\begin{array}{*{20}{c}} > \\ < \end{array}\begin{array}{*{20}{c}} {} \\ {} \end{array}\frac{{{r_x}{S_{rx}}}}{{{r_y}{S_{ry}}}}$

结论：战前力量对比与队员活动面积对比同样重要，如图：

其他：游击队VS常规部队

游击队VS常规部队模型

$\left\{ \begin{gathered} x'(t) = P(t) - ax - gxy \hfill \\ y'(t) = Q(t) - cx - ey \hfill \\ \end{gathered} \right.$

其中a,e分别为甲和乙的自然损失率，其中g,c分别为乙和甲的杀伤率，P(t),Q(t)分别为甲和乙的后勤补给率。

游击队VS常规部队的抛物律

简化模型为

$\left\{ \begin{gathered} x'(t) = - gxy \hfill \\ y'(t) = - bx \hfill \\ \end{gathered} \right.$

同理推得

其中 $b = {r_x}{p_x},g = {r_y}\frac{{{S_{ry}}}}{{{S_x}}}$ ，里面 $r$ 为射击率、 ${S_r}$ 一次射击的有效面积、 $S$ 为游击队员的活动面积，即：

${\left( {\frac{{{y_0}}}{{{x_0}}}} \right)^2}\begin{array}{*{20}{c}} {} \\ {} \end{array}\begin{array}{*{20}{c}} > \\ < \end{array}\begin{array}{*{20}{c}} {} \\ {} \end{array}2\frac{{{r_x}}}{{{r_y}}}\frac{{{S_x}{p_x}}}{{{S_{ry}}}}\frac{1}{{{x_0}}}$

结论：该模型适合以弱胜强。如图所示：

总结

常规战胜负取决于开战前力量(人数)对比，且此比值平方放大，具有集中优势兵力(如三大战役)；游击战战前力量对比与队员活动面积对比同样重要；当游击队与常规部队交战时，更有利于以弱胜强。

赭红色的锆石

关注

6
点赞
踩
12

收藏

觉得还不错? 一键收藏
打赏
0
评论
常微分方程模型及分析——Lanchester战斗模型

在第一次世界大战中，甲、乙两方战役的胜负，往往由双方兵力多少及战斗力强弱两个因素决定，其中，兵力因战斗及非战斗减员而减少，因增援而增加;战斗力与射击次数及命中率有关。请按照正规战争和游击战争两类情况应用常微分方程进行建模，分析并预测战役结局。
复制链接

扫一扫