curl
벡터장의 회전(Curl)과 각속도(Angular Velocity)
구문
설명
예제
입력 인수
출력 인수
세부 정보
알고리즘
curl
는 유한 차분을 사용하여 정의에서 편도함수를 계산합니다. 내부 데이터 점의 경우 편도함수는 중심 차분을 사용하여 계산됩니다. 가장자리에 있는 데이터 점의 경우 편도함수는 단방향(전향) 차분을 사용하여 계산됩니다.
예를 들어, 위치 X
와 Y
에서 행렬 Fx
와 Fy
로 표현되는 크기 m
×n
의 2차원 벡터장 F를 가정하겠습니다. 위치는 [X,Y] = meshgrid(x,y)
로 만들어진 2차원 그리드입니다. 여기서 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 이전에 개발됨