自然进化与气候影响模型.docVIP

PAGE 1 承 诺 书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。 我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 西安电子科技大学 参赛队员 (打印并签名) ：1. 李栋 2. 王晓辉 3. 叶鑫林 指导教师或指导教师组负责人 (打印并签名)： 冯海林 周杰 日期： 2010年 月 日 赛区评阅编号（由赛区组委会评阅前进行编号）： 2009高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）： 目录 TOC \o 1-3 \h \z \u 自然进化与气候影响模型 2 一 问题重述 3 二 问题分析 4 三 模型假设与符号说明 5 3.1模型假设 5 四 模型建立 6 4.1问题一模型 6 4.2问题二模型 6 五 模型求解 7 5.1一个初始值情况 7 5.2不同初始值情况 8 5.3 模型解释 11 六 模型误差和模型改进 12 七 参考文献 13 自然进化与气候影响模型 摘要 模型采用矩阵计算和孟德尔遗传定律结合的方法，通过计算机迭代，计算出基因频率的稳定比例以及基因型和表现型的比例。分析温度升高0.5℃的问题二可知，其实质和问题一一样，只需计算过程中改变存活率这个参数。 针对转移概率矩阵的元素可变、状态转移存在非线性函数关系这一情况，模型中放弃马尔科夫链这一经典模型。而是采用了迭代的计算机数值求解。迭代的过程是这样的： step1:由第k代的基因频率U(k)算出下一代初始基因型结构P0(k+1)； step2:将P0(k+1)和存活率L相乘，得出第k+1代的基因型结构P(k+1)； step3:通过矩阵乘法，算出第k+1代的基因频率U(k+1)； step4:重复step1-step3,直到第300次迭代。 具体的，step1中的算法用到了孟德尔遗传定律，可以通过多项式乘法简单算出。其它两步的算法中，各个需要的矩阵可以根据线性关系求出。之所以要300次循环，是考虑到稳定需要的时间和计算机运算时间两方面。实际计算中，效果很好。得到稳定的基因型结构后，按照表现型将基因型分类，求出表现型的比例。另外，模型中设立了多组初始值。设立了均匀变化的10组，通过多次试验，验证稳定性和稳定点的唯一性。 问题二中，只是改变了step2中的存活率L，算法和问题一一样。计算出两问都有唯一的稳定点，但是两问中的稳定点数值上不相等。说明环境的改变影响了这个昆虫种群的进化。 关键词 遗传定律 矩阵运算 迭代  一 问题重述 昆虫学家在对某种昆虫的研究发现，昆虫身体颜色特征由常染色体基因决定；决定该颜色的基因有三种：记为A，B和C，A和B均为显性基因，C为隐性基因。具有AA和AC基因的昆虫表现为白色，BB和BC基因的表现为灰色，CC表现为透明，而AB表现为黑白相间。 在正常的年份时（年平均气温正常），具有不同基因的子代个体成活率有一些差异，比如AA基因的子代个体有63%的个体能成年并繁殖后代， 对于BB基因这一数据为65.3%，CC为68.6%，AC为73.5%，BC为71.6%，AB为74.2%。 （1）建立适当