MATLAB 色盲具有遗传性,由两种基因c和s的遗传规律决定,男性只有一个基因c或s;女性有两个基因cc,cs或ss,当某人具有基因c或cc时则呈色盲表征。基因遗传关系是:男孩等概率地继承母亲两个基因中的一个;女孩继承父亲的那个基因,并等概率的继承母亲的一个基因。由此可看出,当母亲是色盲时男孩一定色盲,女孩却不一定。用Markov链模型研究非常极端的近亲结婚情况下的色盲遗传,即同一对父母的后代婚配。父母的基因组合共有6种类型,形成Markov链模型的6种状态,问哪些是吸收状态?如果父亲非色盲而母亲为色盲,问平均经过多少代其后代就会变成全为色盲或者全不为色盲的状态,变成这两种状态的概率各为多大?

时间: 2024-03-03 09:51:01 浏览: 26
可以将该问题建模为马尔可夫链,假设状态S为当前家族中出现色盲的人数,S = {0, 1, 2},则该马尔可夫链的状态转移矩阵为: ``` P = [1/2, 1/2, 0; 1/4, 1/2, 1/4; 0, 1/2, 1/2]; ``` 其中,P(i,j)表示从状态i转移到状态j的概率。 根据定义,吸收状态是指无法转移出去的状态。在本问题中,当S为0或2时,都无法转移到其他状态,因此0和2是吸收状态。 接下来,可以使用马尔可夫链的理论来回答问题。对于问题1,吸收状态即为0和2,因此这两个状态是吸收状态。 对于问题2,可以使用瞬时状态分布来计算经过多少代后家族中的人全部为色盲或全部不为色盲的状态。假设当前状态为S0,则经过n步后的状态分布为: ``` π(n) = π(0) * P^n ``` 其中,π(0)为初始状态分布,即初始时只有父母两人,父亲非色盲,母亲为色盲,因此: ``` π(0) = [0, 1, 0] ``` 计算得到转移矩阵的特征值和特征向量: ``` [V, D] = eig(P); ``` 其中D为特征值矩阵,V为特征向量矩阵。由于P是不可约矩阵,因此特征值都为实数。将特征向量按照特征值从大到小排列,则有: ``` λ1 = 1, λ2 = 1/4, λ3 = 0 v1 = [1; 1; 1], v2 = [-1; 0; 1], v3 = [1; -2; 1] ``` 由于P是不可约的,因此λ1 = 1是P的最大特征值,对应的特征向量v1是正常化后的平稳分布。因此,当n趋向于无穷大时,π(n)趋向于平稳分布。计算得到平稳分布为: ``` π∞ = [4/13, 5/13, 4/13] ``` 因此,当n趋向于无穷大时,家族中全部为色盲或全部不为色盲的概率都趋向于4/13,即平稳分布中对应的概率。由于是近亲结婚情况,这个结果并不奇怪。 至于平均经过多少代后家族中的人全部为色盲或全部不为色盲,可以计算状态0和状态2的期望吸收时间: ``` T0 = (V(1,:)*inv(V)*[1; 0; 0])/D(1,1) T2 = (V(3,:)*inv(V)*[0; 0; 1])/D(3,3) ``` 其中,V(1,:)和V(3,:)分别是特征向量矩阵中对应特征值为1和0的行向量,[1; 0; 0]和[0; 0; 1]分别是状态0和状态2的指示向量。计算得到T0 ≈ 1.8462,T2 ≈ 2.3077。因此,平均经过约2代后,家族中的人全部为色盲或全部不为色盲。 需要注意的是,这个结果只是一个平均值,具体情况还会受到随机性的影响。

相关推荐

最新推荐

recommend-type

基于遗传算法的MATLAB16阵元天线的优化.doc

利用Matlab编制一个遗传算法或粒子群算法程序,并实现对间距为半波长均匀直线阵综合,指标如下: 阵元数:16元 副瓣电平: 增益:>11dB 要求撰写设计报告,内容包括:所采用的算法基本原理,目标函数的设计,各个...
recommend-type

双层规划模型的遗传算法求解的Matlab源码-双层规划模型的遗传算法求解的Matlab源码.doc

双层规划模型的遗传算法求解的Matlab源码-双层规划模型的遗传算法求解的Matlab源码.doc 非常实用,值得一看
recommend-type

遗传算法及其Matlab程序

***~~~***~~~***~~~遗传算法的matlab程序,侵权删除。
recommend-type

装箱问题遗传算法MATLAB实现.doc

装箱问题遗传算法MATLAB实现.doc,这份文档介绍了装箱问题遗传算法MATLAB实现,装箱问题遗传算法MATLAB实现.doc
recommend-type

三个遗传算法matlab程序实例

遗传算法(Genetic Algorithm)是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型,是一种通过模拟自然进化过程搜索最优解的方法,它最初由美国Michigan大学J.Holland教授于1975年首先提出来的...
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

MATLAB遗传算法自动优化指南:解放算法调优,提升效率

![MATLAB遗传算法自动优化指南:解放算法调优,提升效率](https://help-static-aliyun-doc.aliyuncs.com/assets/img/zh-CN/8487939061/p208348.png) # 1. MATLAB遗传算法概述** 遗传算法是一种受生物进化启发的优化算法,它模拟了自然选择和遗传的过程。在MATLAB中,遗传算法工具箱提供了丰富的函数和类,用于创建和运行遗传算法。 **1.1 遗传算法的基本原理** 遗传算法的工作原理如下: - **初始化:**创建由随机个体组成的初始种群。 - **评估:**根据目标函数计算每个个体的适应度。 -
recommend-type

failed to register layer: ApplyLayer exit status 1 stdout: stderr: archive/tar: invalid tar header

很抱歉听到这个问题。这个错误通常是由于Docker镜像损坏或Docker存储库中缺少所需的图层引起的。您可以尝试以下几种方法来解决这个问题: 1. 清理Docker镜像和容器:使用命令“docker system prune”清理不再使用的Docker对象。 2. 检查Docker存储库:确保您正在使用可靠的Docker存储库,并检查是否存在所有必需的图层。 3. 重新下载Docker镜像:如果您确定问题是由于损坏的Docker镜像引起的,则可以尝试重新下载Docker镜像。 4. 更新Docker版本:如果您使用的是旧版Docker,则可能会出现此问题。尝试更新到最新版本的Docke
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。