建模培训微分方程部分.pptxVIP

微分方程篇;结束;;;;;;;;;;;;;前面的讨论引导我们考虑两条连续可微的曲线，有如图9-6所画的那样，因为曲线y=f(x)表示国家y需要的最少导弹数，那么曲线上方的区域代表国家y满意的导弹水平，与此类似，曲线x=g(y)右方的区域代表国家x满意的导弹水平。因此图9-6中最黑的区域表示对两个国家来说都满意的导弹水平。曲线y=f(x)和x=g(y)的交点表示双方都满意的最低水平。为了弄明白的确如此，假定国家y拥有y0枚导弹而且探察到国家x有x0枚导弹国家为了适应它的目标将必须增添充分多的导弹以达到图9-6中得坐标1.;反过来，国家x将必须添充分多的导弹以达到坐标2. 这个过程将会继续直到双方同时感到满足，请留意，最黑区域中任何点都会使两个国家满意，而且，其中存在的许多点是非常可能在实际中出现的。图9-6中得交点表示两个国家为适应其目标所需要的最低力量水平;交点的惟一性;;;;;;;;练习.管理不可再生资源：能源危机;用微分方程建模;;;;;;;;;;;;;;;;;; 数值近似方法;欧拉法：用一段段的切线来逼近解曲线;微分方程组建模;一阶线性方程;恰当方程;伯努利方程;常系数齐次线性方程;;;一阶线性微分方程组;;微分方程稳定性理论简介;;;;;;;;;;a;;注意;若初始点存在一个邻域，局部稳定；任意初始点， 全局稳定。对线性方程局部和全局是等价的。; 微分方程数值解;; 2．梯形法;; 龙格－库塔法;微分方程的解析解; 结 果：u = tg(t-c);例 2求微分方程 的通解，并验证。;例 3求微分方程 在初值条件 下的特解，并画出解函数的图形。;例4 求微分方程组 在初值条件 下的特解，并画出解函数的图形。;微分方程的数值解;;Matlab提供的ODE求解器; 1、在解n个未知函数的方程组时，x0和x均为n维向量，m-文件中的待解方程组应以x的分量形式写成.;[x,y]=ode23(fun,[0,0.5],1);fun=inline(-2*y+2*x^2+2*x,x,y); [x,y]=ode23(fun,[0,0.5],1);图中，y1的图形为实线，y2的图形为“*”线，y3的图形为“+”线.;解 令 y1=x，y2=y1’;2、取t0=0，tf=3000，输入命令： [T,Y]=ode15s(‘cxd,[0 3000],[2 0]); plot(T,Y(:,1),-);例1 马尔萨斯模型,以1994 年我国人口为12亿为初值，求解常微分方程 ;常微分方程组初值问题;捕食者与被捕食者问题 ;创建MATLAB的函数文件;------兔子数量; ------狐狸数量;抛射曲线实验，假设阻力与速度成正比。在微分方程中增加阻力项 ;2008电影《集结号》展现出视听震撼的战争场面，92式山炮，炮弹初速: 198米/秒,最大射程:2788米 利用实验程序确定阻力系数 k;实验数据:;嫦娥一号轨道数据实验;假设五个轨道上最大速度从10.3(公里/秒)逐步增加到10.9(公里/秒) 10.3，10.45，10.6，10.75，10.9;转换为一阶微分方程组;右端函数的函数文件 function z=orbit(t,y) GM=3.986005e05; z(1,:)=y(2); z(2,:)=-GM*y(1)./((y(1).^2+y(3).^2).^(3/2)); z(3,:)=y(4); z(4,:)=-GM*y(3)./((y(1).^2+y(3).^2).^(3/2));;function [Vmax,H]=orbitlab(v,h,T) T0=T*60*60; Y0=[-(6378+h),v*cos(-pi/2),0,v*sin(-pi/2)]; [T,Y]=ode23(orbit,[0,T0],Y0); x=Y(:,1);y=Y(:,3); vx=Y(:,2);vy=Y(:,4); V=sqrt(vx.^2+vy.^2); Vmax=max(V); H=max(x); plot(x,y,[0,-(6378+h)],[0,0],ro);实验数据 (h=200km)