QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2371|回复: 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 − 链遗传。 根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因 型的分布。
    3 @% U: w) Z$ ^8 N8 A4 ?- Z1 m8 y: X# @/ d" h1 ]$ ?* |  q" O
    1 常染色体遗传模型
    6 k/ k/ T0 {$ b  X1 ~常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 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 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。. s! }: q6 T6 }, q! C1 W

    2 K6 ?7 L5 y0 r! V2 P+ J" {- J
    6 J, S9 s  W8 t, B5 |- d' I2 _" J! e8 d0 E4 m
    例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?
    ' {. F4 U7 i' c. N) o5 }
    % V4 m5 z7 X9 W(a)假设$ F4 b+ F' L* T1 U0 s: U
    令n = 0,1,2,...。
    ( F' ?$ v3 k5 p
    1 x5 P( t1 l9 O: S7 E# F: L; k/ m4 V, D3 ~! w  {
    * r! a( N9 }) B, I/ S6 B$ z
    (b)建模
    8 o4 }* J) f1 z1 W/ s( z
    9 k* l; Z  v/ B7 B+ C! T' z* X$ ]4 N8 t0 Q
    / M0 T7 r1 d; V* O
    , T' ^2 D* G* t5 |9 ~: j

    2 x3 h0 B; \0 h% }! W3 {% S编写如下 Matlab 程序:5 q( @( F$ l* I. u4 A0 \
    6 W: K' i: S; m/ k- J  V. i
    syms n a0 b0 c0
    . v" K7 S! T' eM=sym('[1,1/2,0;0,1/2,1;0,0,0]');3 `' G, Y; S/ m$ I* t
    [p,lamda]=eig(M);/ r# O! f$ y3 b) V/ ]+ X4 Z8 Q
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];# |/ X* p0 W' ^, p2 Q
    x=simple(x) & f2 i% h% ]( z( `) `$ N2 J
    ' @3 \* g  L! Q' T6 m

    ! N1 {/ @' X/ {/ R( P) ]% a2 z& D  q. S! Y- ^& g

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

    (c)模型的讨论

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


    ( D1 y5 S) z1 q3 I" Q. e) F1 e$ ~1 d! K% Q( u

    8 D6 Y, ]0 L& ]4 N5 t, Z1 _5 ]/ O& c编写如下 Matlab 程序:. O. O" S3 e+ O+ Q
    9 m  T3 ?  ^% p8 W7 x2 r
    syms n a0 b0 c0
    , L4 u4 Q# }& k) d# y7 }4 PM=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');
    2 ~/ c! c* |$ X* s[p,lamda]=eig(M);. O1 _5 A& l6 m
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];
    ( F, L2 @$ Z3 c* Z1 ~x=simple(x)
    - c; x" e+ O4 h( M2 v" g7 [, n( \1 o# O7 V, s6 M, z
    0 d1 r8 O+ _+ n; I, B, Y3 n
    5 _% ~1 [9 t+ w' T  j1 ^6 |

    ; s* C, c7 f/ @0 e: T! c0 c2 常染色体隐性病模型. \: ]* ~5 T/ b: q2 T5 z
    现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。( O$ Y, _3 ?/ ~1 J2 z6 ^
    ' p# c$ @/ k  I2 ^1 G
    (a)假设
    ' c; M- ^  B  s" e8 V; ~* j/ E(i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。
    2 [& A! G$ c+ |0 X2 p' a' o3 u5 f1 b& g, U' d
    % }% f$ p( ], H3 F: {
    : s1 p( _2 q+ p6 w+ ?) T, A9 D0 w
    (b)建模
    1 e3 S+ D# C& E5 E% ~6 @  v
    8 G5 W# ?7 V$ p) d% n9 @0 y5 U2 ]6 ]4 W/ p  ~0 u

    7 d9 G7 @9 G9 @/ k
    * i6 A: p; Q* L. m! R; }0 a( f- J1 k) ^  G' k2 P& F) X
    (c)模型讨论
    ( \; Z# `; }1 L/ H2 L& }4 i研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
    ) T8 r: f2 {' E( R+ v
    9 c. r# B+ O0 b- o5 R+ m; F
    3 o* N$ \6 F0 p, @# G3 I. n- Z2 m" U0 q, L7 D
    下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。
    4 d6 M% ^0 b2 ]2 ^1 E. ^, O" U: t8 ?; b& {2 d+ I7 s
    3   X − 链遗传模型5 u# i+ b  F( H" L6 Q5 c
    X − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。
    1 D, `5 T# W( `) F
    ! i7 s$ T) k# l/ V4 b5 ]) ](a)假设) i4 _  _5 ]4 w; |& h) b3 P
    (i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。3 W5 g" O$ N! n3 o1 c
    % Y  S* g3 k4 p- X4 A- s0 ]
    (ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。
    : H, t7 ]8 ~# i  i$ Q5 b& ^4 o* n( l  {
    % |, G. j( p  N( B

    / Y. t7 V1 Z1 z6 s# Y; h- M% G0 Z! X% B: U* }8 P4 x

    , ~6 x  a5 ]5 r" S& x
    5 V! Q* w: J" Q+ q" {% F
      A1 i6 w: }1 W2 K# g% ~" k+ F% j, H6 ]) p; a/ }' `* ?9 ^: r
    编写如下 Matlab 程序:
    3 D0 r( i$ g& c: t' ~+ ]4 t7 P: U3 l$ V7 B3 V' V' l: h
    syms n a0 b0 c0 d0 e0 f0
      B7 ~' J. W% C4 q$ FM=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;
    ! c; F9 o7 s# H) p* _% w& _, a, ]/ _ 0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];& {  p- s0 x) K; C* `* o. \
    M=sym(M);& [7 y) |$ N" N# n3 r7 ~
    [p,lamda]=eig(M);
    : G( B9 p- p" ?) a# r' Hx=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];
    ) `( d  I9 q5 e' ^6 A3 }x=simple(x)
    , R6 [* h1 `, n9 n( o2 V# O, K. k7 U/ h" c" F& [

    3 r  u! ]5 {! N+ [9 \由上述程序计算结果可以看出
    4 a. F  _+ S' D" D  k0 p$ h6 e) q  z, r7 r. [; e; y: _

      A, }; C6 f2 }; `: T
    5 {3 ~* @+ n4 L# n. }. \7 L# s+ Y  @* A/ R: c1 k
    ————————————————
    9 X+ Y: |, c7 C2 {9 b5 D& t3 ~版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。) a4 p3 g5 ~9 `3 o4 D) e, \# r
    原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596* A: {- _( \7 |9 t
    / i* V! w: _) p! z% r5 B0 q: Q: K. R

    " P: R: U9 z3 Y0 K4 q: Q" _% e/ ]7 A
    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-11 05:14 , Processed in 0.399104 second(s), 50 queries .

    回顶部