How do you make a Moving Quiver Plot ?

3 次查看（过去 30 天）

James McGinley 2019-4-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/458411-how-do-you-make-a-moving-quiver-plot

编辑： James McGinley 2019-4-25

I am trying to use a quiver plot to plot the velocity field. The exact solution should be able to move to a certain location, however when I attempt to create a code to show this, my graph just changes colors. Any advice on how to fix this?

clear
clc
close all
%Start with all given variables
%Number of grid points
N = 31;
%This is in both x and y direction
Nx = N;
Ny = N;
%Viscousity is 1e -1
mu = 1e-1;
% The Domain is 2pi
L=2*pi;
% The domain is the same in both x and y
Lx = L;
Ly = L;
dx = Lx/(Nx-1);
dy = Ly/(Ny-1);
NXX = 33; NYY = 33;
XLin = linspace(0,Lx,NXX);
YLin = linspace(0,Ly,NYY);
% Create N equally spaced values between 0 and L in both directions
xlin = linspace(0,Lx,Nx);
ylin = linspace(0,Ly,Ny);
[x,y] = meshgrid(xlin,ylin);
x=x';
y=y';
u_exact = zeros(31,31);
v_exact = zeros(31,31);
T = 0;
for n = 1:1000
    
    u_exact = cos(x)*sin(y)*exp(-2*mu*T);
    v_exact = -sin(x)*cos(y)*exp(-2*mu*T);
    figure(3)
    hold on
    quiver(x,y,u_exact,v_exact)
    pause(0.01)
    title('The Exact Value of Taylor-Green vortex')
    T = T+ 0.01;
    
end