解说数学建模ppt.pptVIP

模型及其求解 设种群按年龄等间隔的分成n个年龄组，记i=1,2,···，n,时段记作k=0,1,2···,且年龄组区间与时段长度相等(若5岁为一个年龄组，则5年为一个时段)。以雌性个体为研究对象 记在时段k第i年龄组的数量为xi(k);第i年龄组的繁殖率为bi，表示每个个体在一个时段内繁殖的数量；第i年龄组死亡率为di，表示一个时段内死亡数与总数的比，si=1-di是存活率。 * 注意：第k时段的第i年龄组活过来的，是第k+1时段的第i+1年龄组 Xi+1(k+1)=sixi(k) i=1，2,···,n-1, k=0,1,···· 各年龄组在第k时段繁殖的数量和是第k+1时段的第1年龄组 X1(k+1)= k=0,1,···· 记在时段k种群各年龄组的数量为 X(k)=[x1(k),x2(k),····,xn(k)]’ * 这样，有x(k+1)=Lx(k),k=0,1,···· 给定在0时段，各年龄组的初始数量x(0) 就可以预测任意时段k,各年龄组的数量 设一种群分成5个年龄组， 繁殖率b1=0,b2=0.2,,b3=1.8,b4=0.8,b5=0.2 存活率s1=0.5,s2=0.8,s3=0.8,s4=0.1 各年龄组现有数量都是100只， 用matlab计算x(k) * b=[0,0.2,1.8,0.8,0.2]; s=diag([0.5,0.8,0.8,0.1]); L=[b;s,zeros(4,1)]; x(:,1)=100*ones(5,1); n=30; for k=1:n x(:,k+1)=L*x(:,k); end round(x) k=0:30; subplot(1,2,1),plot(k,x),grid * 将x(k)归一化后的向量记做x’(k),称为种群按年龄组的分布向量，即各年龄组在k时段在数量上占总数的百分比。 y=diag(1./sum(x)) ;% sum(x)对列求和 Z=x*y Subplot(1,2,2),plot(k,z),grid 结果分析：时间充分长以后，种群按年龄组的分布x’(k)趋向稳定。 * * * * * * * * * * * * * 用Matlab求解差分方程问题 一阶线性常系数差分方程 高阶线性常系数差分方程 线性常系数差分方程组 * 差分方程是在离散时段上描述现实世界中变化过程的数学模型 例1、 某种货币1年期存款的年利率是r ，现存入M元，问年后的本金与利息之和是多少？ Xk+1=(1+r)xk , k = 0 , 1 , 2 ····· 以k=0时x0=M代入，递推n次可得n年后本息为 * 污水处理厂每天可将处理池的污水浓度降低一个固定比例q，问多长时间才能将污水浓度降低一半？ 记第k天的污水浓度为ck,则第k+1天的污水浓度为 ck+1=(1-q)ck,k=0,1,2,···· 从k=0开始递推n次得 以cn=c0/2代入即求解。 * 一阶线性常系数差分方程 濒危物种的自然演变和人工孵化 问题 Florida沙丘鹤属于濒危物种，它在较好自然环境下，年均增长率仅为1.94%，而在中等和较差环境下年均增长率分别为 -3.24% 和 -3.82%，如果在某自然保护区内开始有100只鹤，建立描述其数量变化规律的模型，并作 数值计算。 * 模型建立 记第k年沙丘鹤的数量为xk,年均增长率为r，则第k+1年鹤的数量为 xk+1=(1+r)xk k=0,1,2······ 已知x0=100, 在较好，中等和较差的自然环境下 r=0.0194, -0.0324,和-0.0382 我们利用Matlab编程，递推20年后观察沙丘鹤的数量变化情况 * Matlab实现 首先建立一个关于变量n ，r的函数 function x=sqh(n,r) a=1+r; x=100; for k=1:n x(k+1)=a*x(k); end * 在command窗口里调用sqh函数 k=(0:20); y1=sqh(20,0.0194); y2=sqh(20,-0.0324); y3=sqh(20,-0.0382); round([k,y1,y2,y3]) * 利用plot 绘图观察数量变化趋势 可以用不同线型和颜色绘图 r g b c m y k w 分别表示 红绿兰兰绿洋红黄黑白色 ： + o * .