流函数 \mathrm{\varphi（x,y,t）=1-tanh[(y-B(t)cos(k(x-ct)))/((1+k^2 (B(t))^2 sin^2 (k(x-ct))))^(1/2) ] \mathrm{B(t)=}\mathrm{B}_\mathrm{0}\mathrm{+\varepsilon cos(wt+\theta)} B_0=1.2、 c=0.12、 k=0.84 、 w=0.4、 \varepsilon=0.3 、 \theta=\pi/2 \mathrm{U(x,y,t)=-}\frac{\mathrm{\partial\varphi}}{\mathrm{\partialy}} V(x,y,t)=\frac{\partial\varphi}{\partialx} （x, y）所在的位置生成的海流流场matlab代码

以下是生成该海流流场的 MATLAB 代码：

% 海流流场参数
B0 = 1.2;
epsilon = 0.3;
w = 0.4;
theta = pi/2;
c = 0.12;
k = 0.84;

% 海流流场计算函数
phi = @(x,y,t) 1-tanh((y-B(t)*cos(k*(x-c*t)))./sqrt(1+k^2*(B(t)*sin(k*(x-c*t))).^2));
u = @(x,y,t) epsilon*k*B(t)*sin(k*(x-c*t)).*sech((y-B(t)*cos(k*(x-c*t)))./sqrt(1+k^2*(B(t)*sin(k*(x-c*t))).^2)).^2;
v = @(x,y,t) -epsilon*k*B(t)*sin(k*(x-c*t)).*tanh((y-B(t)*cos(k*(x-c*t)))./sqrt(1+k^2*(B(t)*sin(k*(x-c*t))).^2)).*sqrt(1+k^2*(B(t)*sin(k*(x-c*t))).^2);

% 海流流场绘制
x = linspace(0,10,100);
y = linspace(-2,2,100);
[X,Y] = meshgrid(x,y);
t = 0;
Z = phi(X,Y,t);
quiver(X,Y,u(X,Y,t),v(X,Y,t));

其中，`phi`、`u` 和 `v` 分别表示流函数、x 方向流速和 y 方向流速的计算函数，`quiver` 函数用于绘制流场。在代码中，我们选择了 x 轴范围为 0 到 10，y 轴范围为 -2 到 2，并在时间 t=0 时绘制了该流场。你可以根据需要进行修改。