QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2347|回复: 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 z5 k  k# y' n! @7 [7 b0 [  m6 b

    0 u1 Z% P( |; P5 _* N1 常染色体遗传模型
    9 i$ d" C1 V1 h- i. [常染色体遗传中,后代从每个亲体的基因对中各继承一个基因,形成自己的基因对, 基因对也称为基因型。如果我们所考虑的遗传特征是由两个基因 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 的可能性相等。下面给出双亲 体基因型的所有可能的结合,以及其后代形成每种基因型的概率,如下表所示。, E8 q& V; O! m0 B- v2 {; S+ j
    : ^# z( o  ?( v
    * I) _4 c1 Z2 V* |* K" Z
    4 s- B" [  b  O. {/ a. {" f
    例 5 农场的植物园中某种植物的基因型为 AA, Aa 和 aa 。农场计划采用 AA型的 植物与每种基因型植物相结合的方案培育植物后代。那么经过若干年后,这种植物的任 一代的三种基因型分布如何?
    1 d' {3 O9 r2 l6 V4 C( o. a8 g& f1 H4 q- H5 R
    (a)假设
    : x2 G$ I+ |# p' e4 P. _+ b令n = 0,1,2,...。 3 M  i9 k5 Q/ \3 |4 j

    + Z( t4 \: U3 T3 u0 i9 b/ A
    - c# G) O8 [& p0 M" p/ U
    5 h: l( [" }% r% f- w4 V2 t(b)建模
    " O) d% M( o# M
    9 r3 O1 f: l; D# t( [8 l' F! h  l: S. H) X

    & C$ O- F7 @( `) w3 W8 w0 w2 R9 k) l  Q  e; ^$ U5 z' ]2 l& J( l

    9 }. y' h7 d' X/ ~编写如下 Matlab 程序:' d7 h7 [& T$ B% q
    ! P( _6 B" p2 l
    syms n a0 b0 c0
    9 e6 L% k1 ^4 ]& d4 `& }4 ?M=sym('[1,1/2,0;0,1/2,1;0,0,0]');
    ' \; O# I1 M4 `[p,lamda]=eig(M);
    6 A* }! n  L9 v5 ?* D9 P4 x0 L9 ax=p*lamda.^n*p^(-1)*[a0;b0;c0];* H3 b( [" L$ K  V5 J, V
    x=simple(x)
    + v2 m* S/ D; B1 L0 T# n% _
    / y& q& F; k2 y) g+ s
    1 \1 C. y+ T, `. A4 d8 t1 O8 Z' Z
    : W- ~. E2 P  a: s

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

    (c)模型的讨论

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


    " S1 ?0 b. G4 d; p) k5 w6 Y5 F0 L. V2 ~! G" }# y: U( C" Q$ D

    * Y5 M: j0 N1 [' Y( U6 T: u9 r编写如下 Matlab 程序:
    1 ^& |" F8 x0 m2 e, x8 ]& B. m" \0 C4 Z! O( {, F: B7 t8 V
    syms n a0 b0 c0
    - z6 T6 P4 D$ V' r! R$ w3 |* KM=sym('[1,1/4,0;0,1/2,0;0,1/4,1]');
    0 k4 e& S1 _/ X; Z; V3 w  x  t( S4 n[p,lamda]=eig(M);+ U9 m( f3 c/ ]/ W% `
    x=p*lamda.^n*p^(-1)*[a0;b0;c0];
    9 W+ @1 d( B% G& O% Xx=simple(x)
    ; x3 `, b' B$ t
    - s& z* {& Z6 u: U- m0 ?5 ]% e- R( ~; k) N# q% ]8 U( v

    $ _& S2 ^3 n! J6 Z, M7 E+ N5 E) W8 q( X; b3 ^" @: z/ ]! v( v
    2 常染色体隐性病模型" C# g3 d; h% a5 l8 a. T
    现在世界上已经发现的遗传病有将近 4000 种。在一般情况下,遗传病与特殊的种 族、部落及群体有关。例如,遗传病库利氏贫血症的患者以居住在地中海沿岸为多,镰 状网性贫血症一般流行在黑人中,家族黑蒙性白痴症则流行在东欧犹太人中间。患者经 常未到成年就痛苦地死去,而他们的父母则是疾病的病源。假若我们能识别这些疾病的 隐性患者,并且规定两个隐性患者不能结合(因为两个隐性患者结合,他们的后代就可 能成为显性患者),那么未来的儿童,虽然有可能是隐性患者,但决不会出现显性特征, 不会受到疾病的折磨。现在,我们考虑在控制结合的情况下,如何确定后代中隐性患者 的概率。, d" ~  [% I- y- ]' f7 M
    ( F$ N' {9 u! _
    (a)假设; B/ W6 G' X/ ]: A& C
    (i)常染色体遗传的正常基因记为 A ,不正常基因记为a ,并以 AA, Aa,aa 分别 表示正常人,隐性患者,显性患者的基因型。
    - n+ J/ ~4 ?, `# N4 a5 l" Y9 A
    + B; {/ x) h: J& x7 v
    2 }0 l2 o  M0 m6 q$ {/ _
    (b)建模
    . p) n7 E  u, t9 A* q' C  N
    4 D: ?! V* Q% m& T% e# @
    ! p! F" q) x- h+ w
    & b$ V. j8 B9 G; K" f- }
    7 g0 P1 }( |1 M& P% r. e2 ]4 O0 k- ]0 n- ~7 \7 B  T4 x4 ~5 r% F
    (c)模型讨论
    ' `6 W# \0 m7 T" C4 z% p研究在随机结合的情况下,隐性患者的变化是很有意思的,但随机结合导致了非线 性化问题,超出了本章范围,然而用其它技巧,在随机结合的情况下可以把(24)式改写为
      o. B9 A5 c4 }1 v7 I- T, `% J2 c5 h7 S0 n% @

      w8 e0 {. r, S7 a" ]' h3 @+ ?) B' o8 K; s5 x5 u6 V# ]. ?) {
    下面给出数值的例子: 某地区有 10%的黑人是镰状网性贫血症隐性患者,如果控制结合,根据(24)式 可知下一代(大约 27 年)的隐性患者将减少到 5%;如果随机结合,根据(25)式, 可以预言下一代人中有 9.5%是隐性患者,并且可计算出大约每出生 400 个黑人孩子, 其中有一个是显性患者。
    , y5 W9 g7 h0 R' m) N( c
    5 @* ?- L6 P& x1 ?7 t' V7 a" _3   X − 链遗传模型# |. d2 N+ D) h, D
    X − 链遗传是指雄性具有一个基因 A 或a ,雌性具有两个基因 AA,或 Aa ,或 aa 。 其遗传规律是雄性后代以相等概率得到母体两个基因中的一个,雌性后代从父体中得到 一个基因,并从母体的两个基因中等可能地得到一个。下面,研究与 X − 链遗传有关 的近亲繁殖过程。
    * d( G) [2 S7 h7 G1 `6 V. t, x  a' q6 Q8 o& S) r; U2 s. G: P) \- M5 Z% l
    (a)假设( D; @% y4 v6 m6 d1 m: F# W
    (i)从一对雌雄结合开始,在它们的后代中,任选雌雄各一个成配偶,然后在它 们产生的后代中任选两个结成配偶。如此继续下去。7 ~/ E+ H6 E) Z6 {
    2 I, _9 L2 L& E$ ]- }
    (ii)父体与母体的基因型组成同胞对,同胞对的形式有 (A, AA) , (A, Aa) , (A,aa) , (a, AA) ,(a, Aa) ,(a,aa) 六种。初始一对雌雄的同胞对,是这六种类型 中的任一种,其后代的基因型如下表所示。4 `4 W+ j( m- i' D6 A4 K

    + n; e" W4 r3 a0 ^- q( }& L! c2 w! d1 M6 B+ P0 f0 ?
    4 Y+ {$ @7 W% f: p

    1 I1 u/ c, N2 o& S! L  i; ]. C2 l# u1 k
    $ ~6 l4 V( w' t8 E9 S  {# ?
    # e+ k9 W' h' u) L: L

    - R% X0 }5 E1 W. R0 b编写如下 Matlab 程序:- c" l# l) `1 v- R0 p! `
      ~/ b. l# x! \# R
    syms n a0 b0 c0 d0 e0 f0) `# L2 q. U' H1 c$ w! ?7 M: R
    M=[1 1/4 0 0 0 0;0 1/4 0 1 1/4 0;0 0 0 0 1/4 0;8 L! i4 n* Q. c3 l4 R
    0 1/4 0 0 0 0;0 1/4 1 0 1/4 0;0 0 0 0 1/4 1];
    4 Q0 j. V  a# g% K  S/ C# CM=sym(M);
    * M. h4 t% p! D0 l) O/ K( D, R& N[p,lamda]=eig(M);
    9 x- m% @/ Q- e4 D; k. ]% G/ G+ ~& Ix=p*lamda.^n*p^(-1)*[a0;b0;c0;d0;e0;f0];) \. j# P) c( V! g3 N1 Y9 r5 k
    x=simple(x)
    9 |" ~; z" V# b3 V  V# _
    , o1 ~- [3 x7 d3 e  k6 ]+ N- I: L( M& a$ i$ ^, X9 D- M
    由上述程序计算结果可以看出
    * X: q( M$ a) \3 y, m+ X% V& G" J/ A2 i6 ~1 m
    # B6 q- `" }) e$ y) g' ^/ g/ N% g
    , ], `, `$ Q4 v' [$ @

    3 o3 ]) k; g: ]2 I+ O————————————————  j( w# U. J3 b: |- K  W
    版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
    * T8 K1 P- [; {. F原文链接:https://blog.csdn.net/qq_29831163/article/details/89646596
    9 R# t0 D& d  w3 a  E9 D* a7 T; B0 b' t8 P& Y4 l  t; z, i
    , n; T, W; {4 G/ L
    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-14 16:42 , Processed in 0.402773 second(s), 51 queries .

    回顶部