matlab矢量模拟,MATLAB画状态方程的向量图

在matlab中画向量图用到的函数是：quiver

函数 通过help 命令可以得到

quiver函数的常用方法，在本文中，只介绍关于状态方程的向量图的用法：

(1)如给定的形式如：x'=1-x^2-y^2;y'=2xy;

matlab 程序如下：

%对微分方程组：x'=1-x.^2-y.^2,y'=2xy;绘制向量场

x0=-3:.5:3;%选定范围

y0=-3:.5:3;

[x y]=meshgrid(x0,y0);%设定绘制点

d=sqrt((1-x.^2-y.^2).^2+(2*x.*y).^2);

u=(1-x.^2-y.^2)./d;

v=(2*x.*y)./d;%u v是相应点(x,y)的指示方向

hold on

quiver(x,y,u,v,0.1,'b');%绘画向量场

(2)如给定形式是：x'=1-x^2;

这时只有一个变量x,我们就添加t作为第二个变量，得到的方程：t'=1;x'=1-x^2;

matlab程序如下：

clear

c=0.3;

x0=-5:.5:5;

y0=-5:.5:5;

[x,y]=meshgrid(x0,y0);

d=sqrt(1^2+(1-y.^2).^2);

u=c./d;

v=c*(1-y.^2)./d;

hold on

quiver(x,y,u,v);