北师MATLAB课件实验课五.pptVIP

数学模型实验课（五） 矩阵特征值，特征向量 常微分方程求解 问题1： 兔子的繁殖III 如果一对兔子每一个月可以生一对兔子，并且兔子在出生二个月以后就具有繁殖后代的能力，三个月后就离开群体。由一对兔子开始，一年可以繁殖成多少对兔子？求这个种群的稳定分布： 1 运用数值模拟方法。 2 运用矩阵特征值和特征向量方法。 问题2 用人口增长的Logistic模型 dy/dt=r(1-y/K)y 拟合美国人口数据： 1790 1800 1810 1820 1830 1840 1850 3.9 5.3 7.2 9.6 12.9 17.1 23.2 1860 1870 1880 1890 1900 1910 1920 31.4 38.6 50.2 62.9 76.0 92.0 106.5 1930 1940 1950 1960 1970 1980 1990 2000 123 132 151 179 204 227 251 281 注：模型的解 y=k/[1+(k/y(0)-1)*exp(-r*t)]是非线性函数 先拟合线性模型 un=r-m*yn 其中 un= (yn+1-yn)/yn得到r和k=r/m的近似值， 以此近似值为参数的初值拟合非线性函数 y=k/[1+(k/y(0)-1)*exp(-r*t)] 问题3 利用Matlab求解常微分方程： dp/dt=(3-6*p(t))/ (2000+2*t), p(0)=0.001. 解析解、数值解。 一. 矩阵的特征值和特征向量 [V,D]=eig（A） A:n×n二维数组，待求特征值和特征向量的矩阵。 V:n×n二维数组，列向量为特征向量 D:对角矩阵，矩阵A的特征值 例1 a=[1,-3,2;4,4,-1;6,3,5];[v,d]=eig(a) 例2 b=[0,1,1;1,0,0;0,1,0];[m,n]=eig(b) 差分方程数值模拟 模型 II a0(n) = a1(n-1) a1(n) = a0(n-1) + a1(n-1) a0(1)=1 a1(1) = 0 模拟幼兔成兔数量和比例随时间的变化 Matlab程序 a=[0 1;1 1]; x=[1 0]; for k=2:12 y=a*x(:,k-1); x=[x y]; end zz=repmat(sum(x),[2 1]); z=x./zz; t=1:12; plot(t,x(1,:),‘r*’,t,x(2,:),‘b*’),grid %幼兔成兔数量随时间变化 plot(t,z(1,:),‘r*’,t,z(2,:),‘b*’),grid %幼兔成兔比例随时间变化 二. 常微分方程求解 1.先求解表达式，再求数值解 Syms x r = dsolve(eq1,eq2,..., ’cond1,cond2,..., v) eq：微分方程 Cond：条件 v：独立变量 2. 二、三阶龙格库塔法求常微分方程的数值解 函数语句：y关于t的二、三阶微分方程 [tout,yout]=ode23(ypfun,[to,tfinal],yo): %调用ode23.m文件解微分方程 参数说明： ypfun:输入参数，被积分的函数名。 to：输入参数，t的初始值 tfinal：输入参数，t的终值。 yo：输入参数，y的初始值。 3. 其他的求解语句 ode45 ode113 ode15s ode 23s ode23t ode23tb 其他的输出语句 odeplot odeprint odephas2 odephas3 例1. 求解 x’=rx(1-x/k)， 0 ≤ t ≤ 20, x(0)=1 建立M文件 logic. M function y=logic(t,x);r=0.3;k=20; y=r.*x.*(1-x./k); 操作 x0=1; [t,x]=ode23(‘logic’,[0,20],x0); plot(t,x) [t;x] syms x k r S=dsolve(‘Dx=r*x*(1-x/k)’, ‘x(0)=0.01’); r =0.3; k =8; s=subs(S);%将参数带入函数表达式 t=0:40;ss=subs(s,’t’,t);％将自变量代入函数表达式 plot(t,ss),grid ss 美国人口数据 年份 1790 1800