curl
向量场的旋度和角速度
语法
说明
示例
输入参数
输出参数
详细信息
算法
curl
使用有限差分计算其定义中的偏导数。对于内部数据点,使用中心差分计算偏导数。对于沿边缘的数据点,使用单侧(正向)差分计算偏导数。
例如,假设有一个二维向量场 F,它由位于 X
和 Y
位置的矩阵 Fx
和 Fy
表示,大小为 m
×n
。位置是由 [X,Y] = meshgrid(x,y)
创建的二维网格,其中 x
是长度为 n
的向量,y
是长度为 m
的向量。然后 curl
计算偏导数 ∂Fy / ∂x 和 ∂Fx / ∂y,如下所示
dFy_dx(:,i) = (Fy(:,i+1) - Fy(:,i-1))/(x(i+1) - x(i-1))
和dFx_dy(j,:) = (Fx(j+1,:) - Fx(j-1,:))/(y(j+1) - y(j-1))
(对于内部数据点。)
dFy_dx(:,1) = (Fy(:,2) - Fy(:,1))/(x(2) - x(1))
和dFy_dx(:,n) = (Fy(:,n) - Fy(:,n-1))/(x(n) - x(n-1))
(对于左边缘和右边缘的数据点。)
dFx_dy(1,:) = (Fx(2,:) - Fx(1,:))/(y(2) - y(1))
和dFx_dy(m,:) = (Fx(m,:) - Fx(m-1,:))/(y(m) - y(m-1))
(对于上边缘和下边缘的数据点。)
向量场的数值旋度等于 curlz = dFy_dx - dFx_dy
,角速度为 cav = 0.5*curlz
。
扩展功能
版本历史记录
在 R2006a 之前推出