MATLAB dibuja un diagrama de fase / campo de dirección / campo vectorial de ecuación diferencial

Para ecuaciones diferenciales

$$y^{\prime }=1-y^2$$

Dibuja su diagrama de fase, donde $X\in [-3,3]$ ,$y\in [-.\; 3..3]$ , tamaño de paso$h=0.2$ , flecha longitud$C=0.01$

Entonces el código MATLAB es

clc,clear,close all
c=0.01;
x_0=-3:0.2:3;
y_0=-3:0.2:3;
[x,y]=meshgrid(x_0,y_0);
d=sqrt(1+(1-y.^2).^2);
u=c./d;
v=c*(1-y.^2)./d;
quiver(x,y,u,v);
xlim([-3,3])
ylim([-3,3])
%		下面几行代码用来绘制初值为 y(0)=1.6 的数值解图像
%		[X,Y]=ode45(@(x,y) 1-y.^2,[-3,3],[0,1.6]);
%		hold on
%		plot(X,Y(:,2),'-r',"LineWidth",2)
%		hold off
%%%%%%%%%% 以下代码保存成 ode1.m 文件 %%%%%%%%%%
%		function dy= ode1(~,y)
%    		dy = 1:1-y(2).^2;
%    	end

---

9 de febrero de 2021 20:35:54