作业帮 > 综合 > 作业

matlab画平面向量场

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/04/29 05:41:33
matlab画平面向量场
微分方程y'=y(1-y):
y(0)=0.2
给定区域为0
x1=linspace(0,6,50);
y1=linspace(0,2,30);
[xx,yy]=meshgrid(x1,y1);
dy=yy.*(1-yy);
dx=ones(size(xx));
dx=dx./(sqrt(dx.^2+dy.^2)+eps);
dy=dy./(sqrt(dx.^2+dy.^2)+eps);
quiver(xx,yy,dx,dy)
axis([0,6,0,2])
[x1,yy1]=ode23('fequ',x1,0.2);
[x1,yy2]=ode23('fequ',x1,1.8);
hold on
plot(x1,yy1,x1,yy2);
streamline(xx,yy,dx,dy,0,0.5);
streamline(xx,yy,dx,dy,0,1.5);
streamline(xx,yy,dx,dy,1,0.5);
streamline(xx,yy,dx,dy,1,1.5);
function s=fequ(tt,xx)
s=xx.*(1-xx);