我正在尝试绘制方向字段。我希望绘制一个归一化的矢量场。 %初始化矢量字段 F = @(X,Y)(2 X-9)(COS(Y))^ 2。。。; A = @(X,Y)的X-X + 1; B = @(X,Y)F(X,Y);
% plot the slope field
hold on
[x, y]=meshgrid(a:(b-a)/m:b , c:(d-c)/n:d);
quiver(x,y,A(x,y),B(x,y),1,'color','blue','linewidth',1.2)
[x,y] = ode45(f,[x0 b],y0);
plot(x,y,'linewidth',1,'color','black')
我搜索了论坛,他们建议使用以下代码
dy = f(x,y)
dx = ones(size(dy))
dyu = dy./sqrt(dx.^2+dy.^2)
dxu = dx./sqrt(dx.^2+dy.^2)
我已经尝试过但由于尺寸错误而无法正常工作。是否有人可以建议修改我的代码,以便我可以绘制标准化的矢量字段。
非常感谢提前!