Sunday, November 19, 2006

可能现在blond是比以前少了是因为以前可能由于交通和观念的原因

blond只和blond结婚(或近似这样),现在结婚虽然不是random,但比以前random很多,我做了一点back of the envelope calculation,在情况下,blond的数量是大于0的。

suppose 我们一开始有两个人口,纯非blond的是p_1,纯blond的是p_2,没有杂非blond。之前他们之间的结婚是一定和自己一样的人结婚,但之后是random的。

suppose在任何时候t有三类人:纯blond的比例是A_t,杂非blond的比例是B_t,纯非blond的比例是C_t,我们知道A_t和C_t的最大限制是p_2/(p_1+p_2)<1,p_1/(p_1+p_2)<1。

那么他们evlove的规则是如下:

A_{t+1}=(A_t*B_t)/2+(B_t^2)/4+A_t^2 %a blond with a hybrid non-blond has 1/2 chance to get a blond child, 2 hybrid nonblonds have 1/4 chance for a blond child, 2 blonds definitely have a blond child

B_{t+1}=(B_t^2)/2+(A_t*B_t)/2+A_t*C_t+(B_t*C_t)/2

C_{t+1}=(C_t^2)+(B_t*C_t)/2+(B_t^2)/4

可能有一点不完整(check了一次了,对A_{t+1}check了几次,这个最重要的式子估计无错吧,其他两个大家将就一点吧,看在我又粗又短的手指来打了这么多奇形怪状符号份上)。

我暂时无时间研究有无及如何得到equilibrium(因为这是back-of-envelope),但估计可以做simulation看看分布如何。无论如何,我们知道A_{t+1}不会为0,因为C_t<1,所以B_t和A_t是不会同时为0的。因此A_{t+1}是不会为0的。

大家放心吧,羊牛可能会变成珍稀动物,但不会绝种的。。。

所以马拉松还得继续跑。。。

0 Comments:

Post a Comment

<< Home