转载:《数学模型(第五版)》
数学建模:核军备竞赛
问题情境
在 20 世纪六七十年代的冷战时期,美苏两个核大国都声称为了保卫自己的安全, 而实行所谓核威慑战略,核军备竞赛不断升级.随着苏联的解体和冷战的结束,双方通过了一系列的核裁军协议.2010 年 4月8日,美国与俄罗斯领导人在捷克首都布拉格签署新的削减战略武器条约.根据这项新条约,美国和俄罗斯将在 7 年内将各自部署的战略核武器削减到不超过 1 550 枚,并把各自的战略核武器运载工具削减到不超过700 枚.
在什么情况下双方的核军备竞赛才不会无限扩张而存在暂时的平衡状态,处于这种平衡状态下双方拥有最少的核武器数量是多大,这个数量受哪些因素影响,当一方采取诸如加强防御、提高武器精度、发展多弹头导弹等措施时,平衡状态会发生什么变化?
模型假设
以双方的(战略)核导弹数量为对象,描述双方核军备的大小,假定双方采取如下同样的核威慑战略:
- 认为对方可能发起所谓第一次核打击,即倾其全部核导弹攻击己方的核导弹基地.
- 己方在经受第一次核打击后 ,应保存有足够的核导弹,给对方的工业、交通中心等目标以毁灭性的打击.
- 在任一方实施第一次核打击时,假定一枚核导弹只能攻击对方的一个核导弹基地,且摧毁这个基地的可能性是常数,它由一方的攻击精度和另一方的防御能力所决定.
图的模型
记 y = f ( x ) y=f(x) y=f(x)为甲方拥有 x x x枚核导弹时,乙方采取核威慑战略所需的最小核导弹数
x = g ( y ) x=g(y) x=g(y)为乙方拥有 y y y枚核导弹时,甲方采取核威慑战略所需的最小核导弹数
不妨让我们看看曲线
y
=
f
(
x
)
y=f(x)
y=f(x)应该具有什么性质。
当
x
=
0
x=0
x=0时
y
=
y
0
,
y
0
y=y_0,y_0
y=y0,y0是甲方在实施一次核打击后已经没有核导弹时,乙方为毁灭甲方的工业、交通中心等目标所需的核导弹数,以下简称乙方的威慑值;当
x
x
x增加时
y
y
y 应随之增加 ,并且由于甲方的一枚核导弹最多只能摧毁乙方的一个核导弹基地,所以
y
=
f
(
y= f(
y=f(
x
)
x)
x)不会超过直线
y = y 0 + x ( 1 ) y=y_{0}+x\ \ \ \ \ \ \ \ \ \ \ \ \ \ (1) y=y0+x (1)
这样,曲线
y
=
f
(
x
)
y=f(x)
y=f(x)应在图 1 所示的范围内,可以猜想它的大致形状如图 2.
曲线
x
=
g
(
y
)
x=g(y)
x=g(y)应有类似的性质
(
y
=
0
(y=0
(y=0 时
x
=
x
0
,
x
=
g
(
y
)
x=x_0,x=g(y)
x=x0,x=g(y)不超过直线
x
=
x
0
+
y
)
x=x_0+y)
x=x0+y),图 2中将两条曲线画在一起,可以知道它们会相交于一点,记交点为
P
(
x
m
,
y
m
)
P(x_m,y_m)
P(xm,ym),我们讨论
P
P
P 点的含义.
根据
y
=
f
(
x
)
y=f(x)
y=f(x)的定义,当
y
⩾
f
(
x
)
y\geqslant f(x)
y⩾f(x)时乙方是安全的( 在核威慑战略意义下),不妨称该区域为乙安全区,曲线
y
=
f
(
x
)
y=f(x)
y=f(x) 为(临界情况下的)乙安全线.类似地
,
x
≥
g
(
y
)
,x\geq g(y)
,x≥g(y)的区域为甲安全区
,
x
=
g
(
y
)
,x=g(y)
,x=g(y)为甲安全线.两个安全区的公共部分即为双方安全区,是核军备竞赛的稳定区域,而
P
P
P点的坐标
x
m
x_m
xm 和
y
m
y_m
ym 则为稳定状态下甲乙双方分别拥有的最小核导弹数,
P
P
P 点是平衡点.
平衡点怎样达到呢?不妨假定甲方最初只有
x
0
x_{_0}
x0枚导弹( 威慑值),乙方为了自己的安全至少要拥有
y
1
y_1
y1枚导弹,见图 2,而甲方为了安全需要将导弹数量增加到
x
1
x_{1}
x1,如此下去双方的导弹数量就会趋向
x
m
,
y
m
.
x_m,y_m.
xm,ym.
模型的精细化
为了研究
x
m
x_m
xm和
y
m
y_m
ym的大小与哪些因素有关,这些因素改变时平衡点如何变动,我们尝试寻求
y
=
f
(
x
)
y=f(x)
y=f(x)和
x
=
g
(
y
)
x=g(y)
x=g(y)的具体形式.
若
x
<
y
x<y
x<y,当甲方以全部
x
x
x枚核导弹攻击乙方的
y
y
y个核基地中的
x
x
x个时,记每个基地未被摧毁的概率为 s,以下简称乙方的残存率,则乙方(平均)有 sx个基地末被摧毁,且有
y
−
x
y-x
y−x个基地未被攻击,二者之和即为乙方经受第一次核打击后保存下来的核导弹数, 它应该就是图的模型中的威慑值
y
0
y_0
y0,即
y
0
=
s
x
+
y
−
x
y_0=sx+y-x
y0=sx+y−x,于是
y = y 0 + ( 1 − s ) x ( 2 ) y=y_0+\left(\begin{array}{c}1-s\\\end{array}\right)x\ \ \ \ \ \ \ \ \ (2) y=y0+(1−s)x (2)
由 0 < s < 1 0<s<1 0<s<1知直线(2)的斜率小于直线(1)的斜率,如图 3.
当 x = y x=y x=y时显然有 y 0 = s y y_0=sy y0=sy,所以有 y = y 0 s ( 3 ) y=\frac{y_0}{s}\ \ \ \ \ \ \ \ (3) y=sy0 (3)
若 y < x < 2 y y<x<2y y<x<2y,当甲方以全部 x x x枚核导弹攻击乙方的 y y y个核基地时 ,乙方的x-y 个将被攻击 2 次,其中 s 2 ( x − y s^2(x-y s2(x−y)个未被摧毁,且有 y − ( x − y ) = 2 y − x y-(x-y)=2y-x y−(x−y)=2y−x个被攻击 1次,其中 s ( 2 y − x ) s(2y-x) s(2y−x)个未被摧毁,二者之和即为图的模型中的 y 0 y_0 y0,即 y 0 = s 2 ( x − y ) + s ( 2 y − x ) y_0=s^2(x-y)+s(2y-x) y0=s2(x−y)+s(2y−x),于是
y = y 0 s ( 2 − s ) + 1 − s 2 − s x ( 4 ) y=\frac{y_0}{s\left(2-s\right)}+\frac{1-s}{2-s}x\ \ \ \ \ \ \ \ \ \ \ \ (4) y=s(2−s)y0+2−s1−sx (4)
直线(4)的斜率小于直线(2)的斜率,如图 3.
当
x
=
2
y
x=2y
x=2y时显然有
y
0
=
s
2
y
y_0=s^2y
y0=s2y,所以有
y = y 0 s 2 ( 5 ) y=\frac{y_0}{s^2}\ \ \ \ \ \ \ \ \ \ \ \ \ \ (5) y=s2y0 (5)
虽然上述过程可以继续下去,但是如果我们允许 x , y x,y x,y取连续值,考察 x = a y , a x=ay,a x=ay,a为大于零的任意实数 ,表示乙(临界)安全条件下甲乙双方导弹数量之比,那么由 x = y x=y x=y时的(3)式和 x = 2 y x=2y x=2y时的(5)式可以设想 y = f ( x ) y=f(x) y=f(x)的形式为
y = y 0 s a = y 0 s x / y , 0 < s < 1 ( 6 ) y=\frac{y_0}{s^a}=\frac{y_0}{s^{x/y}},0<s<1\ \ \ \ \ \ \ \ \ \ \ \ (6) y=say0=sx/yy0,0<s<1 (6)
它应该是图 3 中的光滑曲线,利用微积分的知识可以证明这是一条上凸的曲线.
x
=
g
(
y
)
x=g(y)
x=g(y)有类似的形式,曲线是向右凸的,当然,其中的 s 应为甲方的残存率.
由此可知,这样两条曲线
y
=
f
(
x
)
y=f(x)
y=f(x) 和
x
=
g
(
y
)
x=g(y)
x=g(y)必定相交,并且交点唯一.
进一步研究由(6)式表示的乙安全线 y = f ( x ) y=f(x) y=f(x)的性质:
- 若威慑值
y
0
y_0
y0 变大,则曲线整体上移 ,且变陡;若残存率 s 变大,则曲线变平.
甲安全线 x = g ( y ) x=g(y) x=g(y)有类似的性质.利用这些性质可以用上述模型解释核军备竞赛中平衡点 P ( x m , y m ) P(x_m,y_m) P(xm,ym)的变化
模型的解释
- 若甲方增加经费保护及疏散工业、交通中心等目标,则乙方的威慑值
y
0
y_{0}
y0将变大,而其他因素不变,那么乙安全线
y
=
f
(
x
)
y=f(x)
y=f(x)的上移会使平衡点变为
P
′
(
x
m
′
,
y
m
′
)
,
如图
P^\prime(x_m^{\prime},y_m^{\prime}),如图
P′(xm′,ym′),如图 4.显然有
x
m
′
>
x
m
,
y
m
′
>
y
m
x_m^{\prime}>x_m,y_m^{\prime}>y_m
xm′>xm,ym′>ym,说明虽然甲方的防御是被动的,但也会使双方的军备竞赛升级.
2.
若甲方将原来的固定核导弹基地改进为可移
动发射架,则乙安全线
y
=
f
(
x
)
不变(试说明其威慑值
y
0
、残存率
s
均不变),而甲方的残存率变大(威慑值
x
0
不变),于是甲安全线
x
=
g
(
y
)
向
y
轴靠近,平衡点
变为
P
′
(
x
m
′
,
y
m
′
)
,
如图
5.
显然有
x
m
′
<
x
m
,
y
m
′
<
y
m
,
说明
甲方的这种单独行为,会使双方的核导弹减少.
\begin{aligned}&2.\text{ 若甲方将原来的固定核导弹基地改进为可移}\\&\text{动发射架,则乙安全线}y=f(x)\text{不变(试说明其威慑值}\\&y_0\text{、残存率 }s\text{ 均不变),而甲方的残存率变大(威慑值}\\&x_0\text{ 不变),于是甲安全线 }x=g(y)\text{向 }y\text{ 轴靠近,平衡点}\\&\text{变为 }P^{\prime}(x_m^{\prime},y_m^{\prime}),\text{如图 }5.\text{显然有 }x_m^{\prime}<x_m,y_m^{\prime}<y_m,\text{说明}\\&\text{甲方的这种单独行为,会使双方的核导弹减少.}\end{aligned}
2. 若甲方将原来的固定核导弹基地改进为可移动发射架,则乙安全线y=f(x)不变(试说明其威慑值y0、残存率 s 均不变),而甲方的残存率变大(威慑值x0 不变),于是甲安全线 x=g(y)向 y 轴靠近,平衡点变为 P′(xm′,ym′),如图 5.显然有 xm′<xm,ym′<ym,说明甲方的这种单独行为,会使双方的核导弹减少.