MATLAB Function Reference

curl

計算（compute）一個向量領域（vector field）的捲曲（curl）的速度（velocity）和角度（angular）的速度（velocity）。

Syntax

[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(...)

Description

[curlx,curly,curlz,cav] = curl(X,Y,Z,U,V,W) 計算三維向量領域U、V和W中和流向垂直（perpendicular）的捲曲角度（curl and angular）速度（單位是每一單位時間所有的弧度）。陣列X、Y和Z為U、V和W定義了座標（coordinates）而且這些陣列必須要是單調（monotonic）且三維格化的（3-D plaid）。（經由meshgrid所產生的就是。）

[curlx,curly,curlz,cav] = curl(U,V,W) 假設X、Y和Z是被下列的陳述式給決定：

[X Y Z] = meshgrid(1:n,1:m,1:p)

其中[m,n,p] = size(U).

[curlz,cav]= curl(X,Y,U,V) 計算二維（2-D）向量領域（vector field）U和V中與z垂直（perpendicular）的z部分（z-component）捲曲（curl）速度和角度（angular）速度（單位是每一時間單位所有的弧度。）。陣列X和Y定義在U和V的座標（coordinates），而且一定要是單調（monotonic）且是二維格化（2-D plaid）的（就好比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(...) 傳回捲曲度（curl）。

cav = curl(...) 只有傳回捲曲角度速度（curl angular velocity）。

Examples

這各範例使用上色的（colored）切片平面（slice planes）去展示（display）向量領域中在明確說明位置（location）的捲曲角度速度（curl angular velocity）。

load wind
cav = curl(x,y,z,u,v,w);
slice(x,y,z,cav,[90 134],[59],[0]); 
shading interp
daspect([1 1 1]); axis tight
colormap hot(16)
camlight

這個範例查看（views）在一個體積（volume）中的平面（planes）捲曲角度速度（curl angular velocity）同時繪出（plot）在同一個平面上的速度向量（velocity vector）（quiver）。

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

See Also

streamribbon, divergence

cumtrapz

customverctrl