第五次实验.pdfVIP

清华大学《数学实验》实验报告  实验四 常微分方程数值解 姓名：芦琛璘 班级：化33 学号：2013011934 实验目的： 1．学习线性代数方程组的两类数值解法：直接法和迭代法；  2．通过数据拟合讨论超定线性方程组的最小二乘解；  3 ．学习线性方程组数值解法的 Matlab  实现，并勇气解决实际问题。  实验内容： 【问题1】 种群的数量因繁殖而增加，因自然死亡而减少，对于人工饲养的 种群 （比如家畜）而言，为了保证稳定的收获，各个年龄的种群数量应维持不变。 种群因雌性个体的繁殖而改变，为方便起见以下种群数量均指其中的雌性。种群 年龄记做k=1,2,…,n 当年年龄k 的种群数量记做 ，繁殖率记做 ，（每个 雌性个体一年繁殖的数量），自然存活率记做 ，( 1 ， ，为一年的死 亡率)，收获量记做 ，，则来年年龄k 的种群数量 应为要求各个年龄的种群数 量每年维持不变就是要使 ~ , k1,2 ,3… （1） 若已知，给定收获量建立求各年龄的稳定种群数量的模型（用矩阵、 向量表示）。 （2） 设n=5, b1=b2=b5=0, b3=5, b4=3, s1=s4=0.4, s2=s3=0.6，如要求 h1～h5 为500，400，200，100，100，求x1～x5. （3） 要使h1～h5 均为500，如何达到？ 【模型建立】 （1）记年龄为k 的种群数量为 ，欲使各年龄的种群数量稳定，则必须使 ,，由已知给出的模型得： ~ 1 k1,2,3…n 其中， 为种群繁殖率, 为自然存活率, 为收获量,所以方程组可以写为 Axb 其中： 1 0 b11 2 3 ⋯ 2 1 2 1 1 0 0 0 3 A = 0 2 1 0 0 , x . , b . ⋮ 0 3 ⋱ ⋮ . . 0 0 0 ⋯ 1 . .