一个简单的恋爱动力学模型

非线性动力学 数学模型 恋爱


Part I.

有一个简单的模型能够描述恋爱中的小情侣之间的拉扯游戏,并且能够很好地帮我们复习线非性动力学中的主要知识,我觉得很有趣,下面是具体介绍:

有一对男女,我们就叫他们静希草十郎苍崎青子好了。

文件:静希草十郎 苍崎青子

不妨设:

C(t)tQ(t)t

值越大表示对对方的态度越亲密,正值表示亲密,负值表示疏远。

为了简单起见,二人的恋爱模型可以表示为

(1)C˙=aQQ˙=bC

常数a、b均为正值。

用形象的理解就是,在由C、Q张成的平面空间中的每一个点(c,q)上都分配了一个向量:

(aQ,bC)

我们在4个象限中分别画一下向量的朝向,不难看出,表示二人亲密状态的曲线轨迹应该是一个圆环,也就是说:二人的亲密关系始终在

(+,+)(+,)(,)(,+)(+,+)...

四种状态中循环,如下图中所示:

image-20241003174633944

可以用通解法进一步求解,(1)式改写为:

(2)ddt(CQ)=(0ab0)(CQ)

求解参数矩阵的特征方程:

(3)(0ab0)ν=λν

(4)det(λabλ)=0λ1=iab,λ2=iab

将特征值带回(3)式,特征向量为

(5)ν1=(iab1),ν2=(iab1)

通解为:

(6)(CQ)=c1eiabt(iab1)+c2eiabt(iab1)

用欧拉公式进一步化简,得

(7){C=bω(c1+c2)sinωtQ=(c1+c2)cosωtω=ab,c1c2

这是一个椭圆的参数方程,与之前的分析一致。

因此,两人的关系会一直处于一个循环的状态,不会有任何实质性进展,但好消息是两人会在1/4的时间内同时持亲密态度,并且两人至少不会彼此莫不关心:也就是到达系统稳定点(0,0)。

综上,笨蛋草十郎永远不会得逞,谁让青子是蘑菇我老婆呢?

Part II.

在大多数恋爱关系中,我们采取的策略显然会更加复杂,于是我们可以建立一个稍微更为普适的模型:

(8)ddt(xy)=(abcd)(xy)

特征方程为

(9)det(aλbcdλ)=0

特征值和特征向量的求解略显复杂,但我们可以可以通过刻画平面中的特征值和特征方向来观察系统的演化,展开行列式有:

(10)λ2τλ+Δ=0τ=tr(A)=a+dΔ=det(A)=adbc

特征值与矩阵的迹和行列式值的关系为

(11)λ=τ±τ24Δ2,τ=λ1+λ2,Δ=λ1λ2

通解为:

(12)x(t)=c1eλ1tx1+c2eλ2tx2

可以分情况讨论

Δ<0

这种情况下特征方向为两条直线,且一条直线上的解指数增长,另一条直线上的解指数减小,这种情况,平面中心的不动点称为鞍点(saddle point),系统的相图可以刻画出两人恋爱状态的演化轨迹:

image-20241003181236454

这种情况下,两人恋爱的最终结果可能是走向第一象限的热恋状态,也有可能进入第三象限的对抗状态,具体哪一种则取决于两人的初始感觉。

综上所述,我们可以用一张图来表示所有不动点的情况,我们可以借由不动点的类型轻易地画出系统的大致相图,并判断系统的最终状态:

image-20241003185603548