matlab求差分方子程模型.pptVIP

一阶线性常系数差分方程 Florida沙丘鹤属于濒危物种，它在较好自然环境下，年均增长率仅为1.94%，而在中等和较差环境下年均增长率分别为 -3.24% 和 -3.82%，如果在某自然保护区内开始有100只鹤，建立描述其数量变化规律的模，并作数值计算。 Matlab实现 Xn ＝ (1+r)^n X0 1. 在指令窗口输入： X0=100; r=0.0194; n=20; Xn=(1+r)^n*X0 执行后得到： Xn =146.8563 2. 建立M脚本文件 在M文件编辑器中输入： X0=100; r=0.0194; n=20; Xn=(1+r)^n*X0 单击工具栏中的保存按钮，文件名sqh1.m保存以后选择Debug:run菜单，则会在指令窗口输出： Xn =146.8563 3.建立M函数文件 格式： function 输出变量＝函数名称（输入变量） 例如：输入 function Xn=sqh(n,r) X0=100; Xn=(1+r)^n*X0; end 单击工具栏中的保存按钮，文件名sqh2.m保存以后在指令窗口调用sqh2函数,输入： Xn＝sqh2(20,0.0194) 输出同上。 function Xn=sqh(n,r,X0) Xn=(1+r)^n*X0; end 单击工具栏中的保存按钮，文件名sqh3.m 保存以后在指令窗口调用sqh3函数,输入： Xn＝sqh3(20,0.0194,100) 输出同上。 Matlab实现 Xk+1 ＝(1+r)Xk function x=sqh(n,r) x(1)=100; for k=1:n x(k+1)=(1+r)*x(k); end 单击工具栏中的保存按钮，文件名sqh4.m保存以后在指令窗口调用sqh4函数,输入： xn＝sqh4(20,0.0194) 输出同上。 xn = Columns 1 through 6 100.0000 101.9400 103.9176 105.9336 107.9888 110.0837 Columns 7 through 12 112.2194 114.3964 116.6157 118.8780 121.1843 123.5353 Columns 13 through 18 125.9318 128.3749 130.8654 133.4042 135.9922 138.6305 Columns 19 through 21 141.3199 144.0615 146.8563 利用plot 绘图观察数量变化趋势 k=0:20； y1= sqh4(20,0.0194); plot(k,y1) 得到Figure1 在Edit菜单下Copy Figure，可以将图粘贴到Word文档。 人工孵化是挽救濒危物种的措施之一，如果每年孵化a只鹤放入保护区，观察在自然条件下沙丘鹤的数量如何变化 Xk+1=（1+r）Xk +a 设在中等自然条件下每年孵化5只鹤放入保护区. k=0:20; %一个行向量 y=rgsqh(20,-0.0324,5); plot(k,y) 每年放养多少只会保证物种数量的稳定？ 令Xk=Xk+1 得a=－rX0 当X0＝100，r=－0.0324时， a=3.24 当a=3， a=4时，如图 高阶线性常系数差分方程——一年生植物的繁殖 一年生植物春季发芽，夏天开花，秋季产种，没有腐烂，风干，被人为掠取的那些种子可以活过冬天，其中一部分能在第2年春季发芽，然后开花，产种，其中的另一部分虽未能发芽，但如又能活过一个冬天，则其中一部分可在第三年春季发芽，然后开花，产种，如此继续，一年生植物只能活1年，而近似的认为，种子最多可以活过两个冬天，试建立数学模型研究这种植物数量变化的规律，及它能一直繁殖下去的条件。 记一