战争模型

战战战争争争模模模型型型Abstract:介绍战争中的不同的战争模式，并对战争的胜负进行预测Keywords:游击战、阵地战及混合模型1Introduction在第一次世界大战中，Lanchester就提出了几个预测战争结局的数学模型，游击战、正规战争（阵地战）或者混合战争。后来人们在此基础上对这些模型做了改进和进一步的解释，用来分析历史上的一些著名的战争，得到较好的模拟效果。Lanchester所提出的模型是非常简单的，只考虑双方兵力的多少和战斗力的强弱。兵力的减员包括战斗减员和非战斗减员，又由后备力量的增援而增加。战斗力即杀伤对方的能力，与射击率、射击命中率及战争的类型有关系。当然这些模型没有考虑战争中的士气、政治、经济等相关的因素。但是通过这一类型的模型的分析，可以帮助我们有效得理解毛泽东所提出的游击战争在中国开展的必要性，如何用空间换取时间，最后并取得战争的胜利。一般的战争模型：用x(t);y(t)表示甲乙双方在t时刻所对应的人数。每一方的战斗减员率都取决于双方的兵力和战斗力，分别用f(x;y);g(x;y)表示，而每一方的非战斗减员率(疾病、逃兵等相关因素)与本方的兵力成正比。每一方的增援是固定的函数u(t);v(t)x′=f(x;y)x+u(t);0y′=g(x;y)y+v(t);0(1.1)针对于不同的战争类型，f;g有不同的表达形式。根据不同的形式，分析对战争结局的不同影响。正规战争模型：这一表达式中，我们只需分析甲方的战斗减员率f(x;y).对于甲方的士兵而言，处于乙方的每个士兵的监视和杀伤范围之内，一旦甲方的某个士兵被杀伤，乙方的火力就会立刻集中在其余士兵身上，也就是说甲方的战斗减员与乙方的兵力有关。可以假设f=ay,a表示乙方平均每个士兵对甲方士兵的杀伤率，称为乙方的战斗乙方系数。这里a可以有进一步的分解a=rypy,其中ry是射击率（单位时间的射击次数），py是每次射击的命中率。1这样，就可以根据以上假设给出战斗模型的表达式：x′=rypyyx+u(t);0y′=rxpxxy+v(t);0(1.2)做如下假设，忽略战斗中的非战斗减员，则上述模型可以写成如下形式：x′=rypyy=ay;0y′=rxpxx=bx;0x(0)=x0;y(0)=y0(1.3)对于该方程，因为这是一个初始值问题，可以应用Matlab中的ode45来进行求解。但是在写该程序的时候，必须对该方程进行预处理：把上式写成如下形式：x′y′=0rypyrxpx0xy(1.4)对应的Matlab代码：innity=12;options=odeset('RelTol',1e-4,'AbsTol',[1e-41e-41e-5]);[T,Y]=ode45(@rigid,[012],[x0y0],options);functiondy=rigid(t,y)dy=zeros(2,1);dy(1)=-rypyy(1);dy(2)=-rxpxy(2);这是一个初始值问题，这样可以计算出方程的解。而同样对这样的方程不直接进行求解，可以直接判断双方的胜负：dydx=bxay=rxpxxrypyy(1.5)很容易计算出上面方程的解，ay2bx2=k(1.6)而根据方程的初始条件，有ay20bx20=k(1.7)根据以上表达式可以知道，如果k0,这说明会在某一点，当x的值为零时，y的值大于零。这也说明当一方的兵力为零的时候，另外一方还有兵力存在，这说明有兵力存在的一方2胜利。同理可知，当k0时，另外一方获胜。而k=0时，双方战平。根据以上分析可以知道，根据初始条件所确定的k值，可以确定双方的胜负。y20x20ba=rxpxrypy(1.8)这说明双方的兵力之比y0x0以平方的关系影响着战争的结局。很容易可以看出，如果初始兵力之比增加为原来的两倍，则影响战争结局的能力将增加到原来的4倍。所以由于这个原因正规战争模型也称为平方律模型。游游游击击击战战战争争争模模模型型型双方都用游击部队作战。这种模型和前面的模型是有所不同的，甲方士兵在乙方士兵看不到的某个区域sx内活动。乙方士兵不是向甲方士兵开火，而是向这个隐蔽区域射击，而且不知道杀伤情况。这时甲方士兵的伤亡率不仅和乙方的兵力有关，而且随着甲方兵力的增加而增加。因为在一个有限的区域内，己方士兵越多，被杀伤的就越多。这样可以假设f=cxy。而原来的正面战争模型中甲方的伤亡率决定于乙方的人数。这里的战斗系数可以表示为c=rypy=rysrysx,这里ry仍然为射击率，而命中率等于乙方一次射击的有效面积sry成正比，与甲方的活动面积sx成反比。类似地有g=dxy;d=rxpx=rxsrxsy这样上述方程就可以写成如下形式：x′=rysrysxxy=cxy;0y′=rxsrxsyxy=dxy;0x(0)=x0;y(0)=y0(1.9)类似于上面的战争模型，同样我们可以写出下列的相平面上的解：cydx=m;cy0dx0=m(1.10)与前一模型的不同之处是该模型的相轨线是直线。这时战争的胜负仍然可以根据m值的正负来决定。进一步可以知道y0x0dc=rxsrxsxrysrysy(1.11)初始的兵力之比以线性关系影响战争的结局。进一步可以得到混混混合合合战战战争争争模模模型型型：3x′=rysrysxxy=cxy;0y′=rxpxx=bx;0x(0)=x0;y(0)=y0(1.12)可以得到对应的相轨线：cy22bx=n;cy202bx0=n;(1.13)这时相轨线是抛物线。并且乙方取得胜利的条件是：y20x202bcx0=2rxpxsxrysryx0(1.14)通过这个表达式，很容易可以得到在战争中如何处于弱势的游击一方以空间换取对对方的战斗力的消弱的。可以作如下假设：设甲方初始兵力为x0=100,命中率为px=0:1,火力rx是乙方火力ry的一半，活动区域面积sx=0:1平方公里，而乙方的射击的有效面积是sry=1平方米，那么乙方胜利的条件是：y20x202bcx0=20:10:110621100=100(1.15)根据上式可以知道，两者之比必须大于10，也就是说初始乙方必须10倍于甲方的兵力。硫硫硫磺磺磺岛岛岛战战战役役役：二次世界大战中，经典的战役有硫磺岛战役，根据美军的战地记录，发现模型结果与实际情况吻合很好。硫磺岛战役时日军重要的空军基地。美军在1945年2月19日开始进攻，战斗持续一个多月，日军21500人全部被俘或者阵亡。在战斗中，美军投入兵力73000人，伤亡20265人。战争进行到28天时美军宣布占领该岛，实际战斗到36天才停止。用A(t);J(t)分别表示美军和日军的人数dAdt=aJ(t)+u(t);a0dJdt=bA(t);b0A(0)=0;J(0)=21500(1.16)4而美军的战地记录给出的增援为u(t)=54000;0t16000;2t313000;5t60;others(1.17)对以上方程进行积分，则可得：A(t)=A(0)at∑=1J()+t∑=1u()J(t)=J(0)bt∑=0A():(1.18)这里时间t的取值最大为36，而且J(36)=0,并且根据∑36=1A()=2037000,可以计算出b的值为0.0106.把这个系数带入进(1.18)可以依次算出J(t);t=1;2;


;36.然后跟据(1.18)估计a的值。令t=36，得a=∑36=1u()A(36)∑36=1J()(1.19)分子上是美军的总的伤亡人数，为20265人，分母可由已经得到的J(t)得到，为372500人，这样可以计算出系数a=0:0544,把这个结果带入到（1.18）可得：A(t)=0:0544t∑=1J()+t∑=1u()(1.20)根据这一表达式就可以计算出美军每天的理论值。5