QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2389|回复: 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 − 链遗传。 根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因 型的分布。
    ! d) I6 M- L: e- N3 Y. g  i$ ^! |: \
    1 常染色体遗传模型9 n8 t7 W2 o- E4 W: x3 n+ x$ z0 [* P
    常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 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 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。& N# H( P4 b: M1 \! v
    ( v, D! `1 D. ~$ {7 V/ M7 @0 i3 C* R

    4 F& X7 E) w3 s$ R0 e$ X$ a6 ^/ q' _9 U# {2 @  R
    例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?- {- ]9 W/ C3 {  h
      i. ]- L* C' Z# h3 p* b
    (a)假设
    4 ]8 h0 Q$ ^5 L4 M5 |4 ]! P) e- y令n = 0,1,2,...。 / u& ]- l* i3 E/ Z. J1 R
    # d" U+ p, A: h3 V! R0 b5 H' t

    $ Q- C' F( B  n  Q( k3 u
    ; g: g! z. d+ d7 a(b)建模
    8 n/ O0 }# j. s; L6 E+ u3 L; l# k+ Z- c( g0 @

    1 q  N: m. }# s; R3 Y. t+ t
    8 G; y$ E, f, b8 M4 N/ x  l
    5 L+ K: i6 |2 G! o0 ?- L- a/ O  T
    ( l& ~7 T( S/ Z; M! I7 N编写如下 Matlab 程序:6 b4 W7 N5 @/ E4 M
    + W: v$ N9 s3 w( c
    syms n a0 b0 c0
    : s: B! i: f/ T/ vM=sym('[1,1/2,0;0,1/2,1;0,0,0]');
    1 t# U) K8 {2 P  J[p,lamda]=eig(M);3 [' Q" m5 C1 X" d$ z  O( y/ h
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];
    ' O. Y5 E( x; o3 ?! a, K, P0 t  mx=simple(x) & _/ J; E* i/ H$ ]0 {
    8 d; F% i9 A$ l
    & g. `9 Q! D; z) R
    2 }+ _* O4 B* W7 Z# [

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

    (c)模型的讨论

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

    0 z8 Q+ S/ c, [3 Z3 s  w
    4 F$ u1 N4 H7 q- S* g4 E: L3 |: ^: b
    ; I: D1 J% f) G2 J' `' a
    编写如下 Matlab 程序:7 }8 N, q$ A; o+ n+ q+ \

    ; `5 D) H2 Z6 U! Ssyms n a0 b0 c0
      C7 K- B8 }; YM=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');
    ' K! T0 E' Y% A0 H) N" V[p,lamda]=eig(M);
    $ g# i! O/ O7 K. a5 |9 X7 Bx=p*lamda.^n*p^(-1)*[a0;b0;c0];1 K" l/ a, Y% u: [& Z2 H9 g
    x=simple(x) 9 @7 T: s# ~( y  _) j! ?8 A
    ( o  C$ C, n  l
    1 m/ |; A* c, d3 i; I1 v: o$ P

    7 ?: c( ^& }2 A& G2 N  _2 l/ b* O( {, J4 C  B
    2 常染色体隐性病模型  J. Z$ q# n" [' q' u9 N+ u
    现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。0 ?0 Z3 @) e+ w1 y8 R! D6 X

    * w5 @( z2 _$ ?. o9 Y4 @(a)假设$ i! Z0 u1 B5 C& c3 D
    (i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。
    + r# y; V5 ?0 E6 c# P7 A. G7 o# g: |) E$ Y0 ^$ V

    + F- _" i$ R& ?- ?3 J& d+ M- e! |+ ~5 m; K& k% K/ A
    (b)建模
    5 s' r/ s8 o1 |/ k
    ! t$ O* y9 F  b. }5 P
    % v! ?& f; v5 A4 S4 j6 p+ a
    . ~8 n( c: O4 X/ z! c0 k6 Z) T: H2 q* O9 m6 n4 n- S& B
      Q7 P6 A$ R  D# m1 q; r& H
    (c)模型讨论7 X# Y5 R) `* |% q
    研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为. W1 m+ h2 J3 R) P( C6 G

    3 W, @1 d6 ^/ l3 ]2 Q$ I" V
    9 _5 p) A" M% y' Z5 D) R: F- I+ D! O4 D. q
    下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。+ b2 f$ f* S# J' |9 P2 H9 G5 S' w& M
    + ?! z0 a; K& w) W* B3 O1 {3 {& B% X
    3   X − 链遗传模型# h% q; u# ^9 C" V: |* l
    X − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。
    ! n9 O/ o" n0 w/ t$ F  A8 i1 f" w: b; o; S, e
    (a)假设
    . {. R# P4 E" e" e  r(i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。
    8 i& m+ R& a! T" o! {0 V8 W, Z$ }- A9 W0 b9 l
    (ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。* ^$ p4 n" ~+ w( t' v7 S
    ) e' m: X: r# a, K' W2 r
    5 ]8 A; e, B( h! V/ C3 c+ `

    2 G% d% J5 _0 }* S0 r# e) o0 p3 m5 J3 H3 Q! N1 i4 \

    ) L- |( f: Y$ D3 l( |* S' o! U" q/ o

    - l9 |5 z# ?& Y# X+ d/ ~8 s
    2 q9 ?& o* m5 T编写如下 Matlab 程序:
    " p5 A5 O, Q& Q, W) ?* V! M4 O; ~! T% _0 z
    syms n a0 b0 c0 d0 e0 f0
    * `( |0 ]+ i1 l  G/ p; qM=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;
    . u5 ?: u) D; }5 Z7 \  Y6 W 0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];
    2 |  n, _2 w- \4 XM=sym(M);
    , U" h, m3 Y: X& Z[p,lamda]=eig(M);
    - A  U6 m4 q9 n! R: X9 v0 K  Nx=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];
    $ J7 j. _: t: K; e) Q0 b( mx=simple(x) " K. R4 R+ Z, p+ K
    1 x7 k) ]0 \, t8 X6 E1 U- F
    3 e9 _# P+ r! ^" X. k- \) @: ~6 P
    由上述程序计算结果可以看出6 U1 J. ?0 e- A5 z5 r
    + l" ^! G6 ]+ ]: b5 I2 y
    & o4 Z6 ]! q6 g! r4 s) `- X8 h

    5 g7 Y5 [2 h& }% V- p0 @1 P$ T9 ]- V+ b0 N+ H. h0 ^/ @5 H8 G
    ————————————————( h% H& W8 v$ Y( F) p
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    , @: S. G( Z; [5 f; Z; z原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596
    9 f  a1 [2 y) s
    0 e- i/ N7 k( K$ O7 i* c! S3 {; j7 A* w1 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-7-28 21:55 , Processed in 0.340262 second(s), 51 queries .

    回顶部