数学实验全部答案..docVIP

问题一:鹿群增长模型 问题二:养老保险问题 问题三: 金融公司的支付基金流动 问题四:保险金问题 摘要: 本篇实验报告主要是针对实验十:简单的鹿群增长问题而建立的模型。并且将此模型的求解方法,运用到其他的类似的模型当中。对该模型的求解,运用差分方程组和线性代数的有关知识,通过用matlab编程,实现对矩阵的特征值和特征向量的自动求解。以及将已知矩阵进行对角化。并且用该模型的建模思想和求解方法,对课后的四个实验任务,分别进行了模型的建立和求解。具体的四个实验任务如下: 鹿群增长模型的建立,算法编程以及程序的可行性验证; 养老保险问题模型的建立与求解; 金融公司支付基金的流动模型的建立与求解; 人寿保险计划模型的建立与求解; 针对这几个实验任务, 我分别建立了不同的数学模型,运用Matlab编程进行求解。通过书上给出的实际数据进行了算法的可行性检验,并且通过实际数据给出了该模型的优略性评价。 问题一:鹿群增长模型 问题重述: 假设在一个自然生态地区生长着一群鹿,在一段时间内鹿群的增长受资源制约的因素较小。这里所说的资源包括:有限的食物﹑空间﹑水等。试建立一个简单的鹿群增长模型,并以适当的数据给出结果。 给出数据一:x0=0.8 ,y0=1 ,a1=0.3 ,a2=1.5 ,b1=0.62 ,b2=0.75 ,s=0.8; 数据二: x0=2.8 ,y0=3.4 ,a1=0.4 ,a2=1.8 ,b1=0.61 ,b2=0.72 ,s=0.7; 情况下的结果. 模型假设: (1)只考虑母鹿,并将其分为两组,一岁以下为幼鹿组,其余的为成年组; (2)不考虑饱和状态,即在所考虑的时间段内,种群的增长基本上是不受自然资源的制约; (3)鹿的生育数与鹿的总数成正比。 符号说明: : 第年幼鹿的数量; : 第年成年鹿的数量; : 幼鹿的生育率; : 成年鹿的的生育率; : 幼鹿的存活率; :成年鹿的存活率; : 系数矩阵; :矩阵A的特征值; :矩阵A的特征值; : 开始时幼鹿的数量; :开始时成年鹿的数量; : 刚出生的幼鹿在哺乳期的存活率; 模型的建立: 问题分析: 根据鹿群数量增长的关系模型,建立幼鹿和成年鹿的数量关系式(观测时间取为一年),建立如下的线性差分方程组: (1 ) 令 , 代入方程(1)中,可以得到: (2) 问题转化为对(2)进行求解。 模型的求解原理: 利用线性代数的知识,我们可以将矩阵A进行对角化,即找到矩阵P, 使得为对角矩阵,并且对角线上的元素为A的特征值,即有: 从而有: ; 令 可得: 即可得到: 通过矩阵的方法可以得到第一年到第n年的幼鹿和成年鹿的数量分别为: 根据该公式可以通过编程直接给出第一年到第n年的幼鹿和成年鹿的数量。下面将给出编程的程序。 实验程序: function [x,y]=nianshu(n) % 根据已知数据x0,y0,a1,a2,b1,b2,s进行求解 % 其中x0是初始的幼鹿的数目,y0是初始的成年鹿的数目; % a1是幼鹿的生育率,a2是成年鹿的生育率 % b1是幼鹿的存活率,b2是成年鹿的存活率 % s是刚出生的幼鹿在哺乳期的存活率 % n是年数; % 数据一:x0=0.8 ,y0=1 ,a1=0.3 ,a2=1.5 ,b1=0.62 ,b2=0.75 ,s=0.8 % 数据二:x0=2.8 ,y0=3.4 ,a1=0.4 ,a2=1.8 ,b1=0.61 ,b2=0.72 ,s=0.7 x0=0.8 ;y0=1 ;a1=0.3 ;a2=1.5; b1=0.62 ;b2=0.75 ; s=0.8 % x0=2.8 ;y0=3.4;a1=0.4 ;a2=1.8 ;b1=0.61;b2=0.72;s=0.7 c=-(s*a1+b2); d=s*(b2*a1-a2*b1); m1=[-c+sqrt(c^2-4*d)]/2 m2=[-c-sqrt(c^2-4*d)]/2 a=[(m1-b2) (m2-b2);b1 b1]; cc=(inv(a)*[x0;y0]);4 c1=cc(1); c2=cc(2); for i=1:n x(i)=c1*(m1-b2)*(m1^i)+c2*(m2-b2)*(m2^i); y(i)=c1*b1*(m1^i)+c2*b1*(m2^i); end x=x y=y 运行结果: 对数据一: x0=0.8 ,y0=1 ,a1=0.3 ,a2=1.5 ,b1=0.62 ,b2=0.75 ,s=0.8 ,n=6进行处理后得到的结果为: 年份 幼鹿的数量