matlab,maple画常微分向量场.docVIP

作业 常微分方程向量场相关概念 1、常微分方程向量场定义 答： 设一阶微分方程 满足解的存在唯一性定理的条件。过中任一点 ，有且仅有一个解满足。称域为方程所定义的向量场。 常微分方程向量场性质 答： 性质1：解就是通过点的一条曲线（称为积分曲线），且就是该曲线上的点处的切线斜率，特别在切线斜率就是。 性质2：向量场可以用映射来表示，其中，为一个平面区域，表示由平面一点映射到一个二维向量。 性质3：若函数为上的连续函数，那么向量场也是连续的。 性质4：向量场的原函数不唯一，但是任意两个原函数之间只差一个常数。 3、利用向量场求常微分方程（组）近似解 答： 从几何上看，方程的一个解就是位于它所确定的向量场中的一条曲线，该曲线所经过的每一点都与向量场在这一点的方向相切。形象的说，解就是始终沿着向量场中的方向行进的曲线，因此，求方程满足初始值的解，就是求通过点的这样的一条曲线。 4、利用向量场研究常微分方程定性理论 答： 向量场对于求解微分方程的近似解和研究微分方程的几何性质极为重要，因为，可根据向量场的走向来近似求积分曲线，同时也可根据向量场本身的性质来研究解的性质。 李雅普诺夫稳定性 考虑方程和，现用向量场判断方程的李雅普诺夫稳定性。分别绘制其向量场如下图： 从该图中的向量场方向可以看出，其所有解都渐进稳定于直线。 从该图中的向量场方向可以看出，当初值条件满足时，其解都渐进稳定于某竖直的直线；当初值条件满足时，其解都渐进稳定于直线。 奇点 考虑方程和，现用向量场判断方程的奇点。分别绘制其向量场如下图： 从该图中的向量场方向可以明显观察出，点是这两个常微分方程的零点。 极限环 考虑方程 现用向量场判断方程的极限环。绘制其向量场如下图： 从该图中的向量场方向可以看出，在圆内，有一部分初值点使得其解都渐进稳定于圆，而另一部分趋近于点；在圆外，有一部分初值点使得其解都渐进稳定于圆，而另一部分趋近于无穷远处。 二、基于绘制常微分方程向量场 1、基于绘制常微分方程向量场的思路。 解：向量场计算与与绘制步骤如下： 给定常微分方程~以及区域； 给定x轴分段数（如30），为保证所绘制的向量场中相邻的向量上下左右的距离相等，令； 令x、y的取值分别为； 求由微分方程右边~求在每个点处的斜率k； 由斜率k求对应的方向向量，其中，； 用数据和调用函数，设置颜色等参数绘制向量场。 常微分方程边值问题的解曲线，绘制在上述向量场中的步骤如下： 给定初值点，用‘o’绘制初值点； 调用函数数值求解该边值问题； 调用函数绘制解曲线。 利用绘制常微分方程向量场的程序设计。 解： function [x,y,vx,vy]=vectFieldData(fxy,bx,by,snx) %向量场数据生成通用函数 if nargin4, snx = 30; else, snx = snx*30; end %x轴分段数 sny = floor(snx*(by(2)-by(1)) /(bx(2)-bx(1))); %y轴分段数 hx = (bx(2)-bx(1)) /snx; %x步长 hy = (by(2)-by(1)) /sny; %y步长 x = bx(1) :hx :bx(2); y = by(1) :hy :by(2); [x,y] = meshgrid(x,y); %形成二维矩形数据 k = fxy(x,y); %求斜率 vx = 1./(1+k.^2); %向量横坐标cos a=g(tan a) vy = vx.*k; %向量纵坐标sin a=h(cos a) end 例如方程，其向量场用蓝色表示。当给定初值条件和时的解曲线用红色表示。 绘制程序如下。 %绘制向量场 clc,clear,clf; fxy=@(x,y) y.*(.5-x)+1; %定义函数 f(x,y) bx=[0 6]; %x区间 by=[0 3]; %y区间 [x,y,vx,vy]=vectFieldData(fxy,bx,by); %生成向量场数据 quiver(x,y,vx,vy,color,b) %绘制向量场 hold on %数值求解过定点的常微分方程 y0=[0.2 %给定初值条件 1.8]; plot(zeros(length(y0),1),y0,o) %绘制初值点 [sx,sy] = ode45(fxy, bx, y0); %数值求解常微分方程边值 plot(sx,sy,r,lineWidth,2); %绘制解曲线clc,clear,clf; title(dy/dx= - 4y+2y^2); fxy=@(x,y) -4.*y+2.*y.*y; %定义函数 f(x,y) bx=[0 3]; %x区间 by=[0 3.2]; %y区间 y0=[2+