第5讲+差分方程模型.pptVIP

模型及其求解 记一棵植物春季产种的平均数为c,种子能活过一个冬天的(1岁种子)比例为b,活过一个冬天没有发芽又活过一个冬天的（2岁种子）比例仍为b,1岁种子发芽率a1，2岁种子发芽率a2。 设c,a1,a2固定，b是变量，考察能一直繁殖的条件 记第k年植物数量为Xk，显然Xk与Xk-1,Xk-2有关，由 Xk-1决定的部分是 a1bcXk-1,由Xk-2决定的部分是 a2b(1-a1)bcXk-2 Xk= a1bcXk-1 + a2b(1-a1)bcXk-2 Xk= a1bcXk-1 + a2b(1-a1)bcXk-2 实际上，就是Xk= pXk-1 + qXk-2 我们需要知道x0,a1,a2,c, 考察b不同时，种子繁殖的情况。在这里假设 X0=100,a1=0.5,a2=0.25,c=10,b=0.18~0.20 这样可以用matlab计算了 Xk= a1bcXk-1 + a2b(1-a1)bcXk-2 function x=zwfz(x0,n,b) c=10;a1=0.5;a2=0.25; p=a1*b*c;q=a2*b*(1-a1)*b*c; x1=x0; x2=p*x1; for k=3:n x(k)=p*x(k-1)+q*x(k-2); end K=(3:2