北京工业大学数学建模选修第二次作业.docxVIP

数学建模第二次作业 微分方程稳定性分析 dx 由上式得x=0y=0的实根x0=0,yo=0为微分方程的平衡点,记作P0=x0,yo.易得 dx 由上式得-x=02y=0的实根x0=0,yo=0为微分方程的平衡点,记P0=x0,yo.易得 dx 由上式得y=0-2x=0的实根x0=0,yo=0为微分方程的平衡点 x y lim lim 所以P0点不稳定 dx 由上式得-x=0-2y=0的实根x0=0,yo=0为微分方程的平衡点,记作P0=x0,yo.易得 营养平衡问题 1、模型建立 营养变化速度的数学描述为: dN 易得N=RK+Ce-Kt.平衡点N0=R 种群增长问题 1、模型建立 病菌的数目为N,单位成员的增长率为r1,处于周界表面的那些病菌由于寒冷而受到损伤,它们死亡的数量与N12成比例,其比例系数为r dN 解得N 由上式得d 解得N 画出Nt图像,即微分方程的解族 单种群开发模型 1、模型建立 考虑单种群开发方程: dx 令fx=dxdt=r1-xNx-Ex=0,得到两个平衡点x0=N1-Er,x1=0.又知fx0=E-r0,f Compertz模型 1、模型建立 记时刻t渔场中鱼量为x（t），关于x（t）的自然增长和人工捕捞作如下假设： 在无捕捞条件下x（t）服从compertz模型即dx 同时h（x）=Ex，记F（x）=f（x）-h（x），即 dx 这里不需要解方程已得到x（t）的动态变化过程，只希望知道渔场的稳定量和保持稳定的条件，即时间t足够长以后渔场鱼量的趋向，并由此确定最大持续产量，令F（x）=0得到平衡点X0=Ne-Er，由F（x）的一阶导函数可知，F’（ X0） =EX0 作图：y=f（x）与y=Ex交点p的横坐标就是稳定点X0。则X0*= 综上所述，将强度控制在E*，或者说渔场鱼量保持在最大鱼量的1/e 种族相互依存的模型 1、模型建立 微分方程 d 令 f 解得平衡点p0(0,0),p1( ? ? A= P=-tr(A) q= det(A) 对于p0, 因为P0，不稳定 对于p1 A= P=-3r1(- 对于p2 A= P=-r2(r 对于p3 A= P=-{(r1 q={( 蝴蝶效应与混沌解 首先建立描述微分方程组的外部函数(函数名：lorenz.m)，内容如下： function xdot = lorenz(t,x) xdot = [-8/3*x(1)+x(2)*x(3); -10*x(2)+10*x(3); -x(1)*x(2)+28*x(2)-x(3)]; 再调用ode45()求解，输入： [t,x]=ode45(@lorenz,[0 100],[0 0 1e-10]); plot3(,x(:,2),x(:,1)); plot3(,x(:,2),x(:,3)); plot3(,x(:,3),x(:,1)); 得到的图形依次为：x2~x1， （2） σ σ σ 药物动力学模型 1)模型建立 d 2)模型求解 利用matlab求解方程，程序如下： syms o1 o2 o3; [y1 y2]=dsolve(Dy1=-y1*o1,Dy2=y1*o1-y2*o2,y1(0)=o3,y2(0)=0,t) 得 y1 =o3*exp(-o1*t) y2 =-(-o1+o2)/(o1-o2)^2*o1*o3*exp(-o2*t)-1/(o1-o2)*o1*exp(-o1*t)*o3 创建数据集： xdata=[1;2;3;4;6;8;10;12;16]; ydata=[0.7;1.2;1.4;1.4;1.1;0.8;0.6;0.5;0.4]; 调用cftool工具箱，创建自定义拟合函数：y2 =-(-o1+o2)/(o1-o2)^2*o1*o3*exp(-o2*x)-1/(o1-o2)*o1*exp(-o1*x)*o3 得出拟合结果： General model: f(x) = -(-o1+o2)/(o1-o2)^2*o1*o3*exp(-o2*x)-1/(o1-o2)*o1*exp(-o1*x)*o3 Coefficients (with 95% confidence bounds): o1 = 0.497 (0.1562, 0.8378) o2 = 0.1571 (0.06348, 0.2507) o3 = 2.273 (1.273, 3.274) Goodness of fit: SSE: 0.05594 R-square: 0.9526 Adjusted R-square: 0.9368 RMSE: 0.09