matlabmaple画常微分向量场

1、作业1、 常微分方程向量场相关概念1、常微分方程向量场定义答： 设一阶微分方程 满足解的存在唯一性定理的条件。过中任一点 ，有且仅有一个解满足。称域为方程所定义的向量场。2、 常微分方程向量场性质答：性质1：解就是通过点的一条曲线（称为积分曲线），且就是该曲线上的点处的切线斜率，特别在切线斜率就是。性质2：向量场可以用映射来表示，其中，为一个平面区域，表示由平面一点映射到一个二维向量。性质3：若函数为上的连续函数，那么向量场也是连续的。性质4：向量场的原函数不唯一，但是任意两个原函数之间只差一个常数。3、利用向量场求常微分方程（组）近似解答： 从几何上看，方程的一个解就是位于它所确定的向量场中

2、的一条曲线，该曲线所经过的每一点都与向量场在这一点的方向相切。形象的说，解就是始终沿着向量场中的方向行进的曲线，因此，求方程满足初始值的解，就是求通过点的这样的一条曲线。4、利用向量场研究常微分方程定性理论答： 向量场对于求解微分方程的近似解和研究微分方程的几何性质极为重要，因为，可根据向量场的走向来近似求积分曲线，同时也可根据向量场本身的性质来研究解的性质。(1). 李雅普诺夫稳定性考虑方程和，现用向量场判断方程的李雅普诺夫稳定性。分别绘制其向量场如下图：从该图中的向量场方向可以看出，其所有解都渐进稳定于直线。从该图中的向量场方向可以看出，当初值条件满足时，其解都渐进稳定于某竖直的直线；当初

3、值条件满足时，其解都渐进稳定于直线。(2). 奇点考虑方程和，现用向量场判断方程的奇点。分别绘制其向量场如下图：从该图中的向量场方向可以明显观察出，点是这两个常微分方程的零点。(3). 极限环考虑方程现用向量场判断方程的极限环。绘制其向量场如下图：从该图中的向量场方向可以看出，在圆内，有一部分初值点使得其解都渐进稳定于圆，而另一部分趋近于点；在圆外，有一部分初值点使得其解都渐进稳定于圆，而另一部分趋近于无穷远处。二、基于绘制常微分方程向量场1、基于绘制常微分方程向量场的思路。解：向量场计算与与绘制步骤如下：(1). 给定常微分方程以及区域；(2). 给定x轴分段数（如30），为保证所绘制的向量

4、场中相邻的向量上下左右的距离相等，令；(3). 令x、y的取值分别为；(4). 求由微分方程右边求在每个点处的斜率k；(5). 由斜率k求对应的方向向量，其中，；(6). 用数据和调用函数，设置颜色等参数绘制向量场。常微分方程边值问题的解曲线，绘制在上述向量场中的步骤如下：(7). 给定初值点，用o绘制初值点；(8). 调用函数数值求解该边值问题；(9). 调用函数绘制解曲线。2、 利用绘制常微分方程向量场的程序设计。解：function x,y,vx,vy=vectFieldData(fxy,bx,by,snx)%向量场数据生成通用函数if nargin4, snx = 30; else,

5、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)end3、 利用绘制常

6、微分方程向量场举例。简介：MATLAB是matrix和laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。和Mathematica、Maple并称为三大数学软件。它在数学类科技应

7、用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数

8、库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。发展历程：20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。1984年由Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。例如方程，其向量场用蓝色表示。当给定初值条件和时的解曲线用红色表示。绘制程序如下。%绘制向量场clc,clear,clf;fxy=(x,y) y.*(.5-x

9、)+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); %绘制解曲线4.（2）的程序如下：clc,clear,clf;ti

10、tle(dy/dx= - 4y+2y2);fxy=(x,y) -4.*y+2.*y.*y; %定义函数 f(x,y)bx=0 3; %x区间by=0 3.2; %y区间y0=2+.00025 %给定初值条件 2-.000005 2-.004 1;x,y,vx,vy=vectFieldData(fxy,bx,by); %生成向量场数据quiver(x,y,vx,vy,color,b) %绘制向量场%绘制解曲线hold on%数值求解过定点的常微分方程plot(zeros(length(y0),1),y0,o) %绘制初值点sx,sy = ode45(fxy, bx(1),2, y0(1); %数

11、值求解常微分方程边值plot(sx,sy,r,lineWidth,2); %绘制解曲线sx,sy = ode45(fxy, bx, y0(2:end); %数值求解常微分方程边值plot(sx,sy,r,lineWidth,2); %绘制解曲线line(2.03 2.03,2.3 3.2,LineWidth,2,Color,0 0 0);line(0.5 3.1,-0.03 -0.03,LineWidth,2,Color,0 0 0);(1). 例如方程，其向量场用蓝色表示。当给定初值条件和时的解曲线用红色表示。绘制程序如下。%绘制向量场clc,clear,clf;hold on title(

12、dy/dx= y / (x+sqrt(x2+y2);fxy=(x,y) y ./ (x+sqrt(x.*x+y.*y); %定义函数 f(x,y)b=1bx=-.5*b b; %x区间by= 0 b; %y区间y0=.3 %给定初值条件 .0000001;x,y,vx,vy=vectFieldData(fxy,bx,by,1); %生成向量场数据quiver(x,y,vx,vy,color,b) %绘制向量场%绘制解曲线hold on%数值求解过定点的常微分方程bx0=-0.1 1;plot(bx0(1)*ones(length(y0),1),y0,o) %绘制初值点sx,sy = ode45

13、(fxy, bx0, y0); %数值求解常微分方程边值plot(sx,sy,r,lineWidth,2); %绘制解曲线(2). 例如方程，其向量场用蓝色表示。当给定初值条件时的解曲线用红色表示。绘制程序如下。%绘制向量场clc,clear,clf;hold on title(dy/dx=cos(y)+2*x);fxy=(x,y) (cos(y)+2*x); %定义函数 f(x,y)b=2bx=-b b; %x区间by=-.2*b b; %y区间y0=.3 ; %给定初值条件x,y,vx,vy=vectFieldData(fxy,bx,by,1); %生成向量场数据quiver(x,y,vx

14、,vy,color,b) %绘制向量场%绘制解曲线hold on%数值求解过定点的常微分方程bx0=-1.2 1;plot(bx0(1)*ones(length(y0),1),y0,o) %绘制初值点sx,sy = ode45(fxy, bx0, y0); %数值求解常微分方程边值plot(sx,sy,r,lineWidth,2); %绘制解曲线(3). 例如方程，其向量场用蓝色表示。当给定初值条件时的解曲线用红色表示。绘制程序如下。%绘制向量场clc,clear,clf;hold on title(dy/dx=cos(y)*exp(x);fxy=(x,y) (cos(y).*exp(x);

15、%定义函数 f(x,y)b=2.7bx=-b b; %x区间by=-.2*b b; %y区间y0=-.45 ; %给定初值条件x,y,vx,vy=vectFieldData(fxy,bx,by,1); %生成向量场数据quiver(x,y,vx,vy,color,b) %绘制向量场%绘制解曲线hold on%数值求解过定点的常微分方程bx0=-1.2 2;plot(bx0(1)*ones(length(y0),1),y0,o) %绘制初值点sx,sy = ode45(fxy, bx0, y0); %数值求解常微分方程边值plot(sx,sy,r,lineWidth,2); %绘制解曲线三、基于

16、绘制常微分方程向量场1、基于绘制常微分方程向量场的思路。解：用计算各点的斜率的方法在网格点上画出向量场的方向可以得到向量场。2、 利用绘制常微分方程向量场的程序设计。解：利用函数或者函数可以绘制向量场。3、 利用绘制常微分方程向量场举例。Maple简介：Maple是目前世界上最为通用的数学和工程计算软件之一，在数学和科学领域享有盛誉，有“数学家的软件”之称。Maple起源于1988年，由加拿大安大略滑铁卢的一家公司，Waterloo Maple Inc.（亦称Maplesoft枫软）进行开发和商业销售。最新版是Maple 18。它的主要竞争者是Mathematica。现在共有Personal（

17、个人版）、Professional（专业版）、Academic（学术版）、Government（政府版）和Student（学生版）五个版本。Maple系统内置高级技术解决建模和仿真中的数学问题，包括世界上最强大的符号计算、无限精度数值计算、创新的互联网连接、强大的4GL语言等，内置超过5000个计算命令，数学和分析功能覆盖几乎所有的数学分支，如微积分、微分方程、特殊函数、线性代数、图像声音处理、统计、动力系统、等。Maple不仅仅提供编程工具，更重要的是提供数学知识。Maple是教授、研究员、科学家、工程师、学生们必备的科学计算工具，从简单的数字计算到高度复杂的非线性问题，Maple都可以帮助

18、您快速、高效地解决问题。用户通过Maple产品可以在单一的环境中完成多领域物理系统建模和仿真、符号计算、数值计算、程序设计、技术文件、报告演示、算法开发、外部程序连接等功能，满足各个层次用户的需要，从高中学生到高级研究人员。解：在区域内画出方程的向量场。程序：方法一：with(plots);de := diff(y(x), x) = 2*x*y(x);dfieldplot(de, y(x), x = -2 . 2, y = -2 . 2, color = red, thickness = 2);方法二：with(DEtools);de := diff(y(x), x) = 2*x*y(x);D

19、Eplot(de, y(x), x = -2 . 2, y = -2 . 2, color = red, thickness = 2);得到的效果是一样的。例子2.：在区域内画出方程的向量场。with(plots);de := diff(y(x), x) = -8 x;dfieldplot(de, y(x), x = -2 . 2, y = -2 . 2, color = green, thickness = 2);例子3：画向量场V=y,cos x,z的图形V:=vector(y,cos(x),z);Fieldplot3d(V,x=-2.2,y=-2.2,z=0.4)Matlab画图和maple画图的区别Mathematica、matlab和maple都是三款主流、功能十分强大的数学软件。他们各有千秋，在画图方面，mathmatical效果要比另外两款更好一点，maple画图功能以及符号计算功能要比matlab强大，它集成了世界上最新的算法，所以不管是运行时间，还是效果都要比matlab好。而matlab数值计算、仿真模拟，工程计算要比maple强大，总而言之，只能说各有千秋，综合着使用最好，要用哪块，哪个软件这方面功能更强大就用哪个。总结：世界上最流行的三大数学软件有matlab，mathametica,maple.MATLAB是一款容易使用、可以由多种操作系统