我正在尝试在 MATLAB 中模拟不同的潜在涡流。我得到了流动函数和速度，它适用于一个单一的涡流。但我不知道如何改变漩涡在原点。有谁能够帮我？
[x,y] = meshgrid(-10:10);
Gamma = 5;%circulation
yc = 1; %aberration from the origin
xc = 1;
rc = sqrt((x-xc).^2+(y-yc).^2);
r = sqrt(x.^2 + y.^2);
u = - Gamma*(y-yc) / 2*pi*rc +2* Gamma* y / 2*pi*r -Gamma*(y-yc)/2*pi*rc;
v = Gamma * (x-xc) / 2*pi*rc-2* Gamma*x / 2*pi*r + Gamma*(x-xc)/2*pi*rc;
quiver(x,y,u,v);
我尝试将潜在涡流的势流函数转换为笛卡尔坐标，以模拟样方中的四个涡流：
[x,y] = meshgrid(-8:0.1:8);
Gamma = 5;
spacing = 0.1;
yc = 1;
xc = 1;
rc = sqrt((x+xc).^2+(y+yc).^2);
rc2 = sqrt((x+xc).^2+y.^2);
rc3 = sqrt(x.^2 +(y+yc).^2);
r = sqrt(x.^2 + y.^2);
u = Gamma/2*pi * (-(y+yc)./rc +(y+yc)./rc3 + y./rc2 - y./r);
v = Gamma/2*pi *(-(x+xc)./rc +(x+xc)./rc2 + x./rc3 - x./r);
quiver(x,y,u,v);
因此，我能够模拟不同的涡流，但现在涡流中间似乎有一个源。中间有不想要的来源的漩涡：
