1.求微分方程的解析解,并画出图形,.docxVIP

1．求微分方程的解析解 , 并画出图形， y’= y + 2 x, y(0) = 1,0 x1; 把所有项移到方程左边并除以 x ，得 y - y - 2 = 0 x x dsolve(Dy*(1/x)-y*(1/x)-2=0,y(0)=1,x) ans = 3*exp(x) - 2*x - 2 作图： fplot(3*exp(x)-2*x-2,[0,1,1,4.5]) 4.5 4 3.5 3 2.5 2 1.5 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 图 1 微分方程 y’= y + 2 x 的图形解 2．求微分方程的数值解 , 并画出图形， y y cos x 0, y(0) 1, y (0) 0 先把原方程化为一阶方程组。 令 y1 = y ， y2 = y ，则有 ì = y2 ? y1 ? = - y1 cosx ?íy2 ? (0) 1, y2 (0) = 0 ?? y1 首先建立 M-文件函数 : function f=jie(x,y) f=[y(2);-y(1)*cos （x ）]; 计算 : [x,y]=ode23(jie,[0,2*pi],[1,0]) x = 0 0.0001 0.0005 0.0025 0.0125 0.0625 0.1541 0.2788 0.4345 0.6266 0.8806 1.1500 1.4125 1.6535 1.8848 2.1178 2.3597 2.6056 2.8572 3.1416 y = 1.0000 0 1.0000 -0.0001 1.0000 -0.0005 1.0000 -0.0025 0.9999 -0.0125 0.9980 -0.0624 0.9882 -0.1529 0.9616 -0.2717 0.9085 -0.4083 0.8158 -0.5514 0.6569 -0.6884 0.4592 -0.7686 0.2529 -0.7963 0.0606 -0.7987 -0.1242 -0.8009 -0.3129 -0.8228 -0.5185 -0.8855 -0.7499 -1.0086 -1.0276 -1.2127 -1.4177 -1.5535 作图： y1=y(:,1); y2=y(:,2); plot(x,y1,x,y2,r),gtext(y1),gtext(y2) 5 0 -5 -10 0 1 2 3 4 5 6 7 图 2 微分方程 y y cos x 0, y(0) 1, y (0) 0 的图形解 3．两种相似的群体之间为了争夺有限的同一种食物来源和生活空间而进行生存竞争时，往往是竞争力较弱的种群灭亡， 而竞争力较强的种群达到环境容许的最大数量。 假设有甲、乙两个生物种群，当它们各自生存于一个自然环境中，均服从 Logistic 规律。 （1） x1 (t), x2 (t) 是两个种群的数量； （2） r1 , r2 是它们的固有增长率； （3） n1, n2 是它们的最大容量； （4） m2 (m1 ) 为种群乙 ( 甲 ) 占据甲 ( 乙) 的位置的数量，并且 m2 x2 ; m1 x1. 计算 x1 (t ), x2 (t ) , 画出图形及相轨迹图。解释其解变化过程。 根据给定的具体数据，得到以下方程组 ì x1 + 2 ? dx1 ? = x1(1- ) ? x1 (0) = 10 ? 100 ídt ? x2 + 0.5 x2 (0) = 10 ?dx = x2 (1- ? ) ? dt 100 ? 编制程序如下： 1、 建立 m文件 jingzheng.m 如下： function dx=jingzheng(t,x) dx=zeros(2,1); dx(1)=x(1)*(1-(x(1)+2)/100); dx(2)=x(2)*(1-(x(2)+0.5)/100); 2、 在命令窗口执行如下程序： [t,x]=ode45(jingzheng,0:0.1:15,[10,10]) plot(t,x(:,1),-,t,x(:,2),*),grid 100 90 80 70 60 50 40 30 20 10 0 5 10 15 图 3 种群数量与时间的关系图（ 1） 3、 画出相轨迹图： plot(x(:,1),x(:,2)) 100 90 80 70 60 50 40 30 20 10 20 30 40 50 60 70 80 90 100 10 图 4 相轨迹图 从图 3 可以看出，在开始阶段，两个种群的数量都随着时间增长；随着时间的推移，大约经过 6 个单位时间，种群数量达到饱和；当时间足够长时，两者都将趋向于环境的最大容量。同时也看出，两