QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2377|回复: 0
打印 上一主题 下一主题

[建模教程] 差分方程模型(四):遗传模型

[复制链接]
字体大小: 正常 放大
浅夏110 实名认证       

542

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2020-11-14 17:15
  • 签到天数: 74 天

    [LV.6]常住居民II

    邮箱绑定达人

    群组2019美赛冲刺课程

    群组站长地区赛培训

    群组2019考研数学 桃子老师

    群组2018教师培训(呼伦贝

    群组2019考研数学 站长系列

    跳转到指定楼层
    1#
    发表于 2020-6-9 11:11 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    随着人类的进化,人们为了揭示生命的奥妙,越来越重视遗传学的研究,特别是遗 传特征的逐代传播,引起人们更多的注意。无论是人,还是动植物都会将本身的特征遗 传给下一代,这主要是因为后代继承了双亲的基因,形成自己的基因对,基因对将确定 后代所表现的特征。下面,我们来研究两种类型的遗传:常染色体遗传和 x − 链遗传。 根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因 型的分布。
    6 ~2 c; W6 b- l# y8 I1 @* m5 D, [- M
    1 常染色体遗传模型
    3 p' @0 c4 o# G4 T. y8 g! F! u8 x: y常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 A 和 a 控制的,那么 就有三种基因对,记为 AA, Aa,aa 。例如,金鱼草由两个遗传基因决定花的颜色,基 因型是 AA的金鱼草开红花, Aa 型的开粉红色花,而 aa 型的开白花。又如人类眼睛 的颜色也是通过常染色体遗传控制的。基因型是 AA或 Aa 的人,眼睛为棕色,基因型 是aa 的人,眼睛为蓝色。这里因为 AA和 Aa 都表示了同一外部特征,我们认为基因 A 支配基因a ,也可以认为基因a 对于 A 来说是隐性的。当一个亲体的基因型为 Aa ,而 另一个亲体的基因型是aa 时,那么后代可以从aa 型中得到基因a ,从 Aa 型中或得到 基因 A ,或得到基因a 。这样,后代基因型为 Aa 或 aa 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。; l# W  \; T0 A6 f: [

    - z6 J/ f" O9 y$ h( O4 b6 U5 L, B% b( P. I2 y2 Z

    # |: K9 W$ J3 W0 `% K* P9 d例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?+ v: p6 |) b6 U2 r: `' o& g
    - Q% i+ l, e- B# B+ v! L* A2 [
    (a)假设) U: Z& v6 _% F" b: v
    令n = 0,1,2,...。
    6 P7 P, ]( e: w' I, b( _3 L1 S/ D8 r) `6 [3 i/ e
    # |. D5 w3 s: l7 P- Y0 C2 ^
    ( h1 g9 G# q6 P7 k9 L2 I" @- U" i
    (b)建模
      W+ q# D5 C6 g+ O4 y. |6 F% k! z( ]8 U4 }) Z

    % j. p. K+ j" L0 y. V
    0 `+ P/ V8 f, x3 B. V' t4 G! b% E( ]+ v# R

    . b3 F3 f  M0 }' a2 }# g编写如下 Matlab 程序:6 c* o- m0 W3 {4 u

    ) N% e* c9 D7 A( F% `5 b, dsyms n a0 b0 c0
    5 L8 b4 [( O- E; ]% S& M/ VM=sym('[1,1/2,0;0,1/2,1;0,0,0]');5 N* F6 t! b( d$ `/ Y+ u* W, Y% Z
    [p,lamda]=eig(M);9 ^5 ]) C* }: H% j5 j" ?% }# T$ w
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];2 A" @  F7 a% p4 B4 Z
    x=simple(x)
    2 i& w5 n& q+ c/ F- M
    ! E- z+ B: x3 a% g; i1 o, N; g. U2 b" q
    % r. c0 Q4 O5 r3 b+ q* {

    即在极限的情况下,培育的植物都是 AA型。

    (c)模型的讨论

    若在上述问题中,不选用基因 AA型的植物与每一植物结合,而是将具有相同基因 型植物相结合,那么后代具有三种基因型的概率如下表所示。

    . z6 q( D3 S: R0 `- c

    ) B5 p. C5 {/ `9 |. I0 F2 d9 J4 `
    编写如下 Matlab 程序:, S% X& J' M4 k6 _& [0 Z# w
    + V- i/ {2 z( X3 ^7 v( U1 e
    syms n a0 b0 c0
    / `5 B9 [) z3 V. Y9 g: @3 U. rM=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');6 o4 \3 _/ O+ T4 |1 }' N
    [p,lamda]=eig(M);
    0 ^, o& O& `, z2 tx=p*lamda.^n*p^(-1)*[a0;b0;c0];  i2 G: K' r4 U) _. o# K
    x=simple(x)
    $ G# G& m8 z7 ?! ?4 @( n% o, z) t# k3 q0 A# _% s

    3 F/ d9 ^% y' h& d9 C# ~
    ; I2 p" A5 ]- p" `# g, p; A$ r) I" W% Q7 h9 u# @
    2 常染色体隐性病模型
    ( g* \* \& Q/ C现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。# r1 D" H9 e9 r# V- C; {

    * _; {3 t/ y7 {% V(a)假设. w7 `3 @0 f4 Y. ?( `5 Y
    (i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。& P1 C. O$ x$ ]% u0 \: l, p) q

    8 |5 {& Y1 E' v! B+ m2 O
    / n) d( F$ q# f" q5 N  r
    7 u! z, K9 s) G! q (b)建模2 h) k1 S! ~3 W. l# j- e
    - B0 {2 P: C# J# B, d

    6 ?# V5 s: Z+ t6 A4 v2 [$ b" |
    # }9 e! L* a" I4 Y: o
    / S1 V! N' [- |( Z* [1 q/ T
    ( N( W2 c- u8 I(c)模型讨论) F% j1 t% y0 A4 S
    研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
    1 o. G+ K) N4 l* K9 j9 ?( Q
    4 F) m( y$ M* G$ C, E* \- m0 ?0 x4 u& M7 d+ S

    5 e0 S4 Y2 i; b: o5 Q7 Z下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。
    ' o7 c; F+ g, ^" o
    ' J+ w9 X0 A4 ^3   X − 链遗传模型
    . k, L3 [7 Z+ t) ]X − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。
    5 L9 h8 c8 b. z0 n; w4 F5 W; Z1 }
    $ [: M( K7 w5 O! R(a)假设
    # u. N; j5 g2 s& A(i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。
    6 ~/ Y: {9 t6 e7 \( W. c
    : G5 ~5 m9 e0 g! O. u, y6 _(ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。
    ; v% {1 s( l9 N5 ~/ \# O
    ! }9 e  z, [6 {1 _0 v0 _  l# p3 s: E# {3 f, C4 x7 Q, m

    7 a; Z) {8 {5 Z0 d/ [+ D" p+ x7 v" h; g

    . @# h  m) b5 X- \" x' J
    ) R9 Q4 z5 s2 r, U2 H, a9 `6 @' x- q- T' }
    : ?0 S  L% b3 I
    编写如下 Matlab 程序:* m  A% l/ d& ]( p" ?0 {

    3 C) [) N: t( Y# |syms n a0 b0 c0 d0 e0 f02 p- b6 M. H* d+ b4 d& z
    M=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;; Q" t& \/ J) S/ ~& M; v) [
    0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];
    : l/ V" R7 a+ aM=sym(M);
      n% _/ s8 P" y7 T[p,lamda]=eig(M);
    / w. J$ T2 c+ P+ j' D! }8 [; bx=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];
    / @* {$ m, A4 \  D, T0 o( [8 E: ]9 O6 px=simple(x)
      r1 r7 O, k* H# m& N0 f) B* M3 S7 ~5 U6 W; b, _5 x
    5 h! \1 ~7 p$ S6 q
    由上述程序计算结果可以看出
    5 v8 A- q9 w. }; K
    9 X4 J1 n. U2 w+ b2 ~0 L7 _) b) g. X2 ^
    ) K' h8 ^0 u1 V) l7 G2 [# R

    + z! j; d% f7 X$ L7 o+ H————————————————
    / B! D1 l$ c  t# g1 M版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    + Z! \: |, \* c: I. t原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596
    # C) `4 _  Q  |2 l$ C; |* t8 F3 B) m* O5 M5 i7 }+ N0 D

    # R! k; A  M0 C. w) g& v! c# \
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-14 14:04 , Processed in 0.356818 second(s), 51 queries .

    回顶部