QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2344|回复: 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 − 链遗传。 根据亲体基因遗传给后代的方式,建立模型,利用这些模型可以逐代研究一个总体基因 型的分布。
    . W4 z7 O  ~4 m8 V
    " D, j# R; A% j7 y1 常染色体遗传模型
    8 n3 _/ M, {" ?常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 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 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。9 R$ @* h# T; `) ~0 C* u5 t3 Y1 U6 `
    ! p3 \7 m8 C$ A0 v9 a: [- h" F7 P
    2 L0 A, O+ K- E6 L3 u% j7 T
    4 l" M/ ^( D! g: e9 Z; E
    例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?
    ! C% J. b" w7 T
    # i5 ^$ M& j( h( h& [(a)假设
    / s# J; B1 q3 b/ I! y令n = 0,1,2,...。 * n. G7 _6 d/ w/ p4 C

    , g- G/ K; [7 \1 s. V9 k; h  `! d+ B0 a
    : \- j2 R( M% `, Z; U" O6 X4 c2 e+ U7 m
    (b)建模
    # J6 ~2 u; h2 B$ b0 \2 [: a
    6 ~2 C4 Z  @+ w  V5 B; D, v! }. f6 V0 E1 Z) B; u+ x

    2 y$ E3 Y* S' [3 Y. l9 u5 \) F2 O, O
    2 G! y1 z# `( ^2 O
    编写如下 Matlab 程序:
    ) ~; a/ z# D6 c$ {/ ^( R- a% Y2 l6 v, O' H
    syms n a0 b0 c0. l6 k" Z" U7 E  q. ?5 g* P
    M=sym('[1,1/2,0;0,1/2,1;0,0,0]');* R- Y; A! n2 x% J* E1 w6 `/ o! q% F
    [p,lamda]=eig(M);
    , W! n7 d8 t- rx=p*lamda.^n*p^(-1)*[a0;b0;c0];
    8 N4 D+ X3 z+ d) _) A& ~x=simple(x) " w) s2 N: \2 k5 v
    & ~1 \) ^+ E+ c) ], Q
    & c! [# s/ W: g# v0 k2 a0 K

    / u+ U; B9 y' B7 g: _+ {1 `' Y

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

    (c)模型的讨论

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


    1 t3 z6 r4 X! `( a& K2 ^: C" P0 k8 o6 p) Q5 W
    & K4 ]% {6 `4 S" S- D) h; C- r
    编写如下 Matlab 程序:
    1 M! d$ }0 F8 R& R* W' R  q: D3 m! D$ y* ^' Q) _5 E- ^8 ]
    syms n a0 b0 c0& G# k  }3 _5 j5 Q  q" q. e8 k8 g* J* Z
    M=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');1 \4 s$ v$ m  ]3 W5 {5 N9 l  I
    [p,lamda]=eig(M);* Y: p% J+ d9 g9 v) g
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];, p: O5 v# A7 R0 w/ Z8 x/ j! C- q
    x=simple(x) ! }; e4 _" \9 H* W' E! [

    6 P' g! K0 T1 q8 ^. W" _* d
    - {4 J* G& Z( u+ A$ {+ @* _( ]/ U: e/ a* m$ C7 ]6 z

    ) u/ F' W( O/ g+ e: o  E  e2 常染色体隐性病模型1 y, W$ q! T* c
    现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。
    6 N+ Z+ K; N6 \4 v, F1 ~
    ) B' ~+ I. x7 [! {% \(a)假设1 N& P; r3 R5 @8 ^
    (i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。# n- m! a2 b. _9 b; G

    6 B; i$ \( ^" W/ A, X
    8 x$ N, P1 ~+ J5 [; a! o+ m+ C3 I6 w3 Q
    (b)建模- H+ c! ?$ P1 E0 \# A7 n' U# ?
    , {' {4 @2 A" M" |& F

    " i/ c( G* T/ C  Q3 s
    0 E5 C' `8 [; k# p9 U7 a/ y% i( L5 H! D# f
    , G5 x% [9 G1 S0 x( Q' w  |" q! }
    (c)模型讨论
    & ^/ h, M# v- q# c: l, S( [研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
    % ^( \  Q" z: u% V1 Q: X' g
    ; G* {1 \- |9 s( L3 `/ K; `
    , ?- r1 w* a8 [) h3 T/ z: S! f3 |/ ]8 e
    下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。
    9 N; f0 e2 d; v- B6 ~
    1 z% H. x, r# m8 U% a7 n3   X − 链遗传模型
    3 F; v9 }. g' p+ \8 \X − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。; a# i4 o4 n( D
    0 S/ `  G8 y- Y" Z& R+ B( c/ o
    (a)假设( A7 g, `! ^* ]7 [
    (i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。% i" B) k# C$ O, {/ R

    ) Z- |, l$ Z3 F) V5 P(ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。
    & U% i. \  s- c  n+ B$ y" ~
    # I: R" A0 l' f  B& k4 D* H0 @; T. j7 w4 z# _- {) L. @6 \) e
    9 u8 |7 x0 e- a1 @" t; J" Q$ n& w
    - }/ [, X( W& A! P1 Z
    5 v- Z6 d& c* z; F& q! l
    ; S8 G" ~1 m( i- Q: G7 F) \
    " }# \* ~8 v. E- }

    + c  h6 l/ I8 E& g. s编写如下 Matlab 程序:
    6 K8 Q/ |! p+ n5 p# a" ]
    0 R% W/ Y/ N9 Zsyms n a0 b0 c0 d0 e0 f0% p" o7 a6 L2 {/ V3 k3 [' D, f
    M=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;4 D, p/ [& L- L2 g$ r
    0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];8 k' T9 a9 z6 x5 h
    M=sym(M);: e' `! K4 K$ [% q  v$ w- H
    [p,lamda]=eig(M);
    8 \) t& w; K( g1 n$ tx=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];
    3 D$ k% ^9 X( Yx=simple(x) # ~% Z) x/ F5 ^. V6 x7 h- w( F6 Y1 D( }+ x
    & T7 s  ^1 s7 F& a8 z( \; F: u

    1 M# f& G' S: t+ a4 ^% y由上述程序计算结果可以看出- v* g  e0 u" y& _- U. ~

    5 @& K+ l6 W7 m, l; z+ }9 e! q  u! I9 r! k+ k
    , h+ R9 r* h, D1 X0 S  w; M6 j
    " W! [( s# n# O1 n  |4 b7 M
    ————————————————
    2 H6 p) X* I' o+ g& t- _版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。; B3 O3 l$ Z, Z
    原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596
    + J# w$ S* K- ~8 W4 w' I3 r; q1 A$ x( A

    ( q* g$ K5 c9 x; o: L& V4 D
    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-4-12 14:06 , Processed in 0.396477 second(s), 51 queries .

    回顶部