基于系统动力学的兰切斯特方程研究.docxVIP

基于系统动力学的兰切斯特方程研究 0 兰切斯特方程的应用 在第一次世界大战期间，英国的fw。为了使其对战争的描述更加准确, 人们对兰切斯特方程进行了种种改进和推广。但是, 以前的大多数文献在应用推广的兰切斯特方程时对参数进行了特殊设置, 又简化为最初的兰切斯特战争模型, 其结果没有达到推广的目的。本文先对兰切斯特方程进行了3种推广, 计算了所有参数对结果的影响;并将所得的几种推广结果相互结合, 使得模型对战争的描述更加全面。 1 基本作战模型 若用x,y代表作战双方, 假设x(t) 和y(t) 分别代表t时刻这两支部队的作战兵力, 兰切斯特方程为 {dx(t)/dt=-f(x,y)-αx(t)+u(t)dy(t)/dt=-g(x,y)-βy(t)+v(t)(1){dx(t)/dt=?f(x,y)?αx(t)+u(t)dy(t)/dt=?g(x,y)?βy(t)+v(t)(1) 式中:α0,β0表示非战斗减员系数;u(t) 和v(t) 分别表示作战双方的兵力补给。在不考虑自然损耗和兵力补给的条件下, 方程组 (1) 可以得到正规战、游击战和混合战3个不同的基本作战模型。 兰切斯特方程是以“接触性”的大规模作战样式为基础得来的, 上述的3个基本作战模型对战争的描述比较简单。为了更加准确地描述新的战争样式, 对该方程组进行一些改进和推广。 1.1 p2p、h1的比例 在考虑作战双方兵力补给和自然损耗时, 兰切斯特方程为 {dx(t)dt=-αx(t)+n1-E?G(t)y(t)1+Edy(t)dt=-βy(t)+n2-G(t)x(t)1+E(2)???dx(t)dt=?αx(t)+n1?E?G(t)y(t)1+Edy(t)dt=?βy(t)+n2?G(t)x(t)1+E(2) 式中:α和β为双方的非战斗减员系数;n1和n2分别为作战双方的兵力增援速率。方程最后一项表示的是一方对另一方随时间变化的杀伤效果, 其中G(t) 是描述战争激烈程度的参量,E=P2/P1为双方武器打击精度P2和P1的比值。这样, 方程组 (2) 的解可以写为 {x(t)=n2Λ2-n1βΛ1Λ2-αβ-De(μ-λ)t+μ+κΛ1Fe-(μ+λ)ty(t)=n1Λ1-n2αΛ1Λ2-αβ+De(μ-λ)t+μ+κΛ2Fe-(μ+λ)t(3)???x(t)=n2Λ2?n1βΛ1Λ2?αβ?De(μ?λ)t+μ+κΛ1Fe?(μ+λ)ty(t)=n1Λ1?n2αΛ1Λ2?αβ+De(μ?λ)t+μ+κΛ2Fe?(μ+λ)t(3) 式中: D=Λ1Λ22μ(μ+κ){(Ν2+n2α+μ+κΛ1Λ2-αβ)-μ+κΛ2(Ν1+n1β+μ-κΛ1Λ2-αβ)}F=Λ1Λ22μ(μ+κ){μ+κΛ1(Ν2+n2α-μ+κΛ1Λ2-αβ)+(Ν1+n1β-μ-κΛ1Λ2-αβ)}λ=(α+β)/2?κ=(β-α)/2?Λ1=G/(1+E)?Λ2=EG/(1+E)?μ=(κ2+Λ1Λ2)1/2。D=Λ1Λ22μ(μ+κ)??????????????????????????N2+n2α+μ+κΛ1Λ2?αβ??????μ+κΛ2?????N1+n1β+μ?κΛ1Λ2?αβ????????????????????F=Λ1Λ22μ(μ+κ)?????????????????????μ+κΛ1?????N2+n2α?μ+κΛ1Λ2?αβ?????+?????N1+n1β?μ?κΛ1Λ2?αβ??????????????????????????λ=(α+β)/2?κ=(β?α)/2?Λ1=G/(1+E)?Λ2=EG/(1+E)?μ=(κ2+Λ1Λ2)1/2。 方程 (2) 综合考虑了影响战争结局的几个主要因素, 并且涉及双方武器精度对作战效果的影响, 对作战过程的描述较为全面。根据实际的战斗情况, 实时调整方程里的时变参数, 可以得到较准确的战争结果预测。 1.2 明确趋势外力区中合成力量j的概率 在不考虑兵力补给、忽略非战斗减员的条件下, 描述多兵种混合作战的兰切斯特方程可写为 {dxi(t)dt=-n∑j=1αij?ijyj(t)?i=1???mdyj(t)dt=-m∑i=1βjiφjixi(t)?j=1???n(4) 式中:xi(t) 为x方第i类作战类型的兵力;yj(t) 为y方第j类作战类型的兵力;αij为y方第j类武器对x方第i类目标的损伤系数;βji为x方第i类武器对y方第j类目标的损伤系数;?ij为y方第j类武器向x方第i类目标射击的概率;φji为x方第i类武器向y方第j类目标射击的概率。 令兵力向量A= (x1,x2, …,xm,y1,y2, …,yn)T, 初始兵力向量A0= (x10,x20, …,xm0,y10,y20, …,yn0)T。分配系数矩阵定义为 Q