QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2375|回复: 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 − 链遗传。 根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因 型的分布。& u# Z" F' C; h  i* C& m' A' G5 L
    . r+ t7 Z+ Q& L( j! A
    1 常染色体遗传模型  ~7 r+ M+ G0 {. o; ?5 w
    常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 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 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。3 T$ x# t( A  X

    1 [1 a) U- t, c; g4 w& h
      k+ }. r' ]: R) x% C: l
    # N, D! `0 n6 W" _; s例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?  F6 n8 T' e( h+ E

      `2 @" k! e/ e  _" x(a)假设
    3 o- t$ p' D. X* s8 n$ `令n = 0,1,2,...。 0 Z; ~% ?) o0 y) B, z
    ( h- w$ d& A) q( v! }

    2 Q6 z! X5 a* K
    $ P6 A% ?/ E5 e/ F1 }7 q* X(b)建模
    * i6 t( G& h' E! S
    + {& A* B  V% q9 l% G& y  L1 O2 D- B9 O1 B
    " Z. p$ e# u. p. y: Y( `

    " ^9 |6 h* _7 T, n" R" o1 H/ X! ]0 V1 z
    编写如下 Matlab 程序:
    . v! e% A: y* ^# _
    ' u) v4 |; z( A+ B- jsyms n a0 b0 c0
    ) Y& }, A2 f- f2 V' \- xM=sym('[1,1/2,0;0,1/2,1;0,0,0]');
    3 K& o! T3 R  }: ?% z) C- z5 }[p,lamda]=eig(M);
      E6 J% g+ S$ ~  w$ cx=p*lamda.^n*p^(-1)*[a0;b0;c0];
    8 J' R8 n! B9 h; i: t5 Zx=simple(x) 8 q3 l+ R) \! b1 B% D
    # q6 ?! }3 T* n" d% s, T6 B

    9 U6 s! |7 v( u" R6 [' Y/ c  |. P2 l- c( J. M

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

    (c)模型的讨论

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

    ; C4 b, F; C, j! i2 A! n

    0 c0 b- T% r( [, t( ?" `2 Q7 s3 P% s; ]
    编写如下 Matlab 程序:- ?8 f3 n: r/ m8 R1 C8 G7 O
    0 Z7 x( j0 [' {+ q5 [% h/ f
    syms n a0 b0 c0
    , ], x; Y: L" ~# i4 fM=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');
    6 C7 |2 H* h4 s" V3 U- u& H[p,lamda]=eig(M);
    ; K) Y) }0 P, r) X) m- h, l: }x=p*lamda.^n*p^(-1)*[a0;b0;c0];
    9 Z( V# Q# u) J6 s+ _x=simple(x) ) k! ?. S0 x+ q6 S9 I: M
    : T3 i4 d, X, a; R# \
    + N! D8 h0 F4 {1 S
      H+ Y9 n5 Y+ |* B7 y8 d

    7 B2 Y$ Z* y  n! a2 常染色体隐性病模型
    ' b$ `* Q9 o& u- P1 }现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。
    . ~! _) o$ v, u% D3 q$ v& v1 X# s) [: I
    (a)假设- e6 L2 m* U; G% `
    (i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。
    ! u5 Q- b0 s9 H& L" x0 P1 p+ g9 w, ~8 b
    , h5 @8 y/ Z6 S5 E3 j2 W6 f

    ! S) E  ~9 F$ y0 ^0 }  m (b)建模5 P8 N$ X6 K1 i/ E* g1 v: ~

    ; C( @* I  H" @% D; e5 K3 n0 v; l
    ; \, R, n+ B9 i! i  h& a1 |* ?, m0 {+ P. ?# l

    7 R( S5 ^% c( C0 j* W3 @4 _& }* K! X- N: y4 ]( }6 @
    (c)模型讨论
    . l2 \) [) S$ U+ D! @研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
    ; c  z! b7 b/ `7 W* ~9 B1 n
    . S9 v; e; O* B/ ]: t( b4 c3 N' D3 s, [$ g" o

    & }8 ~% j2 q+ N下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。
    : K/ k* Z5 s# A, F6 D4 I
    % Q9 Q9 T+ ~' v# E0 Q3   X − 链遗传模型: P0 P  Z7 N$ E3 G/ t0 [
    X − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。) z/ [+ @, p+ ]0 I# y
    ! Y5 J8 Y, g8 ^- e& O. V
    (a)假设. Z) ^' x. ?  Y9 G( J, r+ m+ S
    (i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。# \- K. k8 D$ n( |& O0 Y# b: \) \3 T

    ( D: }# j7 e2 K* n$ k: W" L$ Q1 t8 c% T4 r(ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。6 G" O/ I1 }1 K; Y! m3 g

    3 y/ Y$ z1 J6 R. s3 ?# {/ Y% ^7 ^
    6 k$ Q# U9 L( t; Q# T( y. C9 ]7 d
    7 o& O- m8 ?$ V8 I! Z
    ( r# k! c& h2 y: w# P( b# K" o4 C) L
    ! O# r( D# ]3 j% S6 c# m5 p
    % {( J. E  c7 B, V$ `7 ^9 |: M& g- ]( N1 R8 y

    ! |5 O) ?4 J/ w; J3 }" ~编写如下 Matlab 程序:
    $ `/ F; e% R6 h" b1 b
    ; I6 j4 j# R9 Y9 X$ Xsyms n a0 b0 c0 d0 e0 f0$ b1 ]% r! V! W- t8 J1 Q1 S
    M=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;
    % t4 j7 R+ a1 | 0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];/ y& Y+ O: \) u3 j& z0 Z
    M=sym(M);! G1 @6 W- U7 t" V" J
    [p,lamda]=eig(M);
    4 ~2 |6 \; K8 q' E% [x=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];# Y' x1 z% y$ M( {  k6 N& n
    x=simple(x)
      I) \/ ^% \5 e4 [2 O1 F9 x. ?( |! ?% ?4 W5 W/ b/ M6 A
    * p1 @" G/ a: v1 t9 y; V- U* e8 j
    由上述程序计算结果可以看出9 Q0 G# m7 y+ L4 W5 y
    : O3 d! j4 d. A. E: F4 m

    ! x" ?' }0 i) J' i8 _' n' Y2 I! U3 J) K6 k8 K. r1 D
    6 Y% h+ I: C# w5 e) O8 w
    ————————————————9 _6 q- A* p, h  V8 ~# |% ^) o
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。3 S- M: X8 c; f
    原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596
    / z! H7 V  H0 s. u- m3 q1 }2 l8 g* `+ T9 T

    " E. i% Y& w1 e" x2 q+ h
    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-12 21:20 , Processed in 0.596552 second(s), 51 queries .

    回顶部