首先，注意矩阵乘积*和元素乘.*之间的区别。如果两个操作数是标量（单个值），它们是相同的，但是对于矩阵输入，它们做的事情非常不同。因为你在矩阵t，x和y上操作，你应该更新你的函数：

f = @(t,x,y) (0.02*sin(2*pi*t).*(x.^2 - 2*x) + x);  
dfdx = @(t,x,y) (0.02*sin(2*pi*t).*(2*x - 2) + 1);
dxdt = @(t,x,y) -pi*sin(pi*f(t,x,y)).*cos(pi*y);
dydt = @(t,x,y) pi*cos(pi*f(t,x,y)).*sin(pi*y).*dfdx(t,x,y);

字符串
每次*或^或/的两个操作数都是向量或矩阵时，您都希望使用逐元素运算符（.*，.^，./）。
我还分离出了dxdt和dydt，因为表示整个x和y值数组的向量输出很复杂。
接下来，在调用d时，必须传入计算的t、x和y。此外，由于我们不想绘制3D矢量图，因此我们将选择单个t，而不是它们的数组（如果您更喜欢绘制3D图，则应在meshgrid调用中包含t）：

t = 0;
[x,y] = meshgrid(0:.1:2,0:0.1:1);
u = dxdt(t,x,y);
v = dydt(t,x,y);

型
现在我们应该能够绘制：

quiver(x,y,u,v);

型
对于多个t值，您可以绘制图形并设置动画：

for t = linspace(0,4,100)
   u = dxdt(t,x,y);
   v = dydt(t,x,y);
   quiver(x,y,u,v);
   pause(0.5)
end

型
（虽然不幸的是，我没有看到t的任何变化？）

有两件事你应该知道：
1.请使用元素操作，例如.*和.^，而不是矩阵或向量操作，例如*和^;
1.应将d解叠加为两个单独的偏导数，即dxdt和dydt

代码

f = @(t,x,y) (0.02.*sin(2.*pi.*t).*(x.^2 - 2.*x) + x);  
dfdx = @(t,x,y) (0.02.*sin(2.*pi.*t).*(2.*x - 2) + 1);
dxdt = @(t,x,y) -pi.*sin(pi.*f(t,x,y)).*cos(pi.*y); 
dydt = @(t,x,y) pi.*cos(pi.*f(t,x,y)).*sin(pi.*y).*dfdx(t,x,y);

[x,y] = meshgrid(0:.1:2,0:0.1:1);
t = linspace(0, 4, 100);
arrayfun(@(p) {figure;quiver(x,y,dxdt(p,x,y),dydt(p,x,y))}, t, UniformOutput=false);

字符串