QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2376|回复: 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 − 链遗传。 根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因 型的分布。
    % ]& m: \- U' j, Q+ C" ~
    - E. G/ n; ^! ?( R+ o& {! N/ X2 P1 常染色体遗传模型2 H$ F, U& r; g
    常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 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 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。) T* T# C0 g% @; @# W' I( H

    ' ?$ D( M0 X( w* S( ~& V& j/ a# e
    $ u" `7 M  O" t, {  a( ~8 P5 J" H/ g2 {! k& }$ a
    例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?
    ; r- e/ e# r/ R6 b1 {& H& r* c6 s: u
    # \1 a6 F0 r1 ~; j3 _+ r(a)假设0 H+ y$ j: g' M7 E* S$ T
    令n = 0,1,2,...。
    + O# _5 d7 w, A  F# H3 U" r7 e; I( S; |3 w. F/ z4 P

    0 h; N# ?4 V& f5 L1 p2 V0 x7 V5 H0 g# p3 D! i5 P" M
    (b)建模5 D$ A% D9 q* O  F5 {/ P" `  N0 ]3 ~
    & r) B) T9 Q+ k+ g9 L  Y. E

    , m0 x4 }, {+ o. X8 A; K: q' |, J: I
    9 Z0 L! F3 ~; N+ A
    9 }' s, {( M" W2 B9 f& i- a  _9 b: b4 N! n1 A; }, R6 ~0 P' l' p3 Y& l
    编写如下 Matlab 程序:
    9 h; {% j+ s# m0 i
    5 W% n' Q1 Q' z. R- fsyms n a0 b0 c08 J! [$ a8 f: E4 `0 ^) ^( G2 ]
    M=sym('[1,1/2,0;0,1/2,1;0,0,0]');$ q9 a+ x0 J! B
    [p,lamda]=eig(M);3 g8 G% a! w. _/ ~
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];
    3 X& i5 f  I8 `# A" M! ux=simple(x) . g/ ^& s, S' R0 {8 Y; R

    7 O9 ?- ?, }  s) j$ G; e3 s7 }0 \0 `3 N# H2 T2 H

    ( b  H. C+ `) v" n( S/ O2 j5 X' @

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

    (c)模型的讨论

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


    9 t/ m9 m9 X  t) p2 ~" d) i- j
    0 U! D7 t% R+ I/ @  V6 W
    9 `2 Q% ?6 E/ T( s编写如下 Matlab 程序:
    0 ?* z1 w$ g0 ~$ ^, T  J# S, C# I3 z
    + ~/ N6 y& U# ~2 N" g2 c7 e3 [+ J: ?syms n a0 b0 c04 o3 U& M; L$ i1 G; D4 n
    M=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');" W% R+ o! D$ c- Q0 @
    [p,lamda]=eig(M);& @* s4 X- k3 u: o
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];3 g7 @4 L( M6 w
    x=simple(x) % E6 f) c& W  }$ A7 r4 o

    : C, w/ [. R9 h, }- s$ i2 @3 x4 \& ^) W

    ' N" M. b$ w: Z: ]) ~& z9 u& p. z
    2 常染色体隐性病模型
    6 I# T8 F) b, y: M; f; ~现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。
    7 V' k* ^6 U6 [0 d+ r" P  `/ z8 f3 g& H/ N, H6 X% y+ _
    (a)假设
    # _7 g" N7 d1 |(i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。' a' c: x& Z/ @8 @' A& b

    : D* \. |9 A8 p% Z: ?: Y- e
    9 S$ }$ q: X+ K7 p, c" a
    6 O+ A3 D. m9 U+ [0 y. P6 h9 k# ` (b)建模- u/ j8 Z+ V. g/ ?3 L/ q
    8 @, W9 g  P7 B6 Y# f; f/ `
    / B. j8 p3 L4 ~! V4 \* X, N

    3 Q. d# M2 @; S% r# ?2 P
    " N3 w# B. L! t! |8 u& D
    & b* o& l. _: f! ~' f(c)模型讨论
    1 |5 h3 u0 a1 E7 s) w研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
    ! j  p) F2 z' [1 T7 t1 G2 a" X+ R& u# B% f; I* z! a0 B
    ; d" K3 W: ?* S' R, t% a" E( y

    / G) |! b/ l/ L- n; j下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。# h/ C; o/ {0 n% w
    " H; E3 O  F: s, k" c
    3   X − 链遗传模型
    4 u4 E, s" a" f% t0 \0 H) GX − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。
    9 @) Y. J% \! V
    4 v5 Y% B) E& X(a)假设5 K" @! A0 R2 A( \! a, |4 S* |
    (i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。
    , e9 P: V2 ?0 b/ x. k( s, f2 s
    - g) o: M$ B) d8 B, B3 q) W, ^. M& [/ r(ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。
    ) w5 w& n$ i1 V# @: ?; J5 P! a( C& u" z

    ' P8 |" ~# C& a  v2 x
    1 b0 u! i$ k* T9 L, ^9 [8 g" b# y, F4 Q0 [* b7 c

    9 S; Q( e; p3 C" q
    . Q3 O) [8 \( E/ }' x( u; N2 t- B. c: e
    5 Q4 k$ d: @* l4 `
    编写如下 Matlab 程序:
    * a+ [* w# a5 t  W* ]1 l/ a* c- p
    7 x1 _0 W! h, ?. O8 Y- r$ T7 A/ l$ s" Qsyms n a0 b0 c0 d0 e0 f0
    ) W+ }* z; e; B; L; z8 g% P# BM=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;/ ]8 e& ^$ E) s- Y- R* E
    0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];7 U/ x4 d' b8 f" u4 i& P5 Y. B  E9 C
    M=sym(M);8 h4 U' s* f1 q+ \/ \: N
    [p,lamda]=eig(M);8 p2 e7 {0 s) ~, V
    x=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];
    ) K# {) H9 C, F$ y- \7 b; Sx=simple(x) - a: ?- n1 q+ ^2 }2 U4 e" |2 ?

    ! @3 S* l% u) y1 `: N" l6 @. f% m+ z3 [: {; P
    由上述程序计算结果可以看出
    8 t, L# `% k3 @
    7 \" U( L  d+ a. ~( m+ Q
    2 x0 [( r4 E0 s$ K; u6 l
    4 b/ E' X1 \) q9 A- t4 S
    1 V: c$ q$ K! a3 ~# U————————————————, D  u6 e7 t% P" s3 h! o2 p$ U
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。, V1 F, Z$ H0 C6 R2 G
    原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596
    # ^# K. ^5 S; L' c, V8 Q& {4 }
    $ S8 S# k8 v' k- P+ P; a7 ?6 {, ?2 W; r, E( ^8 X# B9 m$ ^6 u
    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 11:45 , Processed in 0.422152 second(s), 51 queries .

    回顶部