Вычислите завихрение и скорость вращения векторного поля
[curlx,curly,curlz,cav] = curl(X,Y,Z,U,V,W)
[curlx,curly,curlz,cav] =
curl(U,V,W)
[curlz,cav]= curl(X,Y,U,V)
[curlz,cav]= curl(U,V)
[curlx,curly,curlz] = curl(...)
[curlx,curly] = curl(...)
cav = curl(...)
[curlx,curly,curlz,cav] = curl(X,Y,Z,U,V,W) вычисляет завихрение (curlx, curly, curlz) и скорость вращения (cav) перпендикуляр к потоку (в радианах на единицу измерения времени) 3-D векторного поля UVW.
Массивы XY, и Z, которые задают координаты для UV, и W, должно быть монотонным, но не должны быть расположенными равными интервалами. XY, и Z должен иметь то же число элементов, как будто произведенный meshgrid.
[curlx,curly,curlz,cav] =
curl(U,V,W) принимает XY, и Z определяются выражением
[X Y Z] = meshgrid(1:n,1:m,1:p)
где [m,n,p] = size(U).
[curlz,cav]= curl(X,Y,U,V) вычисляет вихревой z- компонент и перпендикуляр скорости вращения к z (в радианах на единицу измерения времени) 2D векторного поля U, и V.
Массивы X и Y, которые задают координаты для U и V, должно быть монотонным, но не должны быть расположенными равными интервалами. X и Y должен иметь то же число элементов, как будто произведенный meshgrid.
[curlz,cav]= curl(U,V) принимает X и Y определяются выражением
[X Y] = meshgrid(1:n,1:m)
где [m,n] = size(U).
[curlx,curly,curlz] = curl(...), [curlx,curly] = curl(...) возвращает только завихрение.
cav = curl(...) возвращает только вихревую скорость вращения.
Используйте окрашенные плоскости разбиения, чтобы отобразить вихревую скорость вращения в заданных местоположениях в векторном поле.
load wind cav = curl(x,y,z,u,v,w); h = slice(x,y,z,cav,[90 134],59,0); shading interp daspect([1 1 1]); axis tight colormap(hot(16)); camlight set([h(1),h(2)],'ambientstrength',.6);
Отобразите вихревую скорость вращения в одной плоскости объема. Затем постройте векторы скорости в той же плоскости.
load wind k = 4; x = x(:,:,k); y = y(:,:,k); u = u(:,:,k); v = v(:,:,k); cav = curl(x,y,u,v); pcolor(x,y,cav); shading interp hold on quiver(x,y,u,v,'y'); hold off colormap('copper');
