Campus Livre
curl
Rotacional e velocidade angular de um campo vetorial.
Categoria: Graphics > 2-D and 3-D Plots
Sintaxe
cz = curl(X, Y, U, V)Componente z do rotacional de um campo 2-D.
[cz, cav] = curl(X, Y, U, V)Também retorna a velocidade angular.
[curlx, curly, curlz, cav] = curl(X, Y, Z, U, V, W)Componentes do rotacional e velocidade angular de um campo 3-D.
Descrição
O rotacional mede quanto um campo gira em torno de cada ponto. É estimado por diferenças finitas na grade. Para a rotação de corpo rígido (U, V) = (-y, x) o rotacional em z vale 2 em todo ponto, e para um campo sem rotação, como (x, y), vale 0.
A velocidade angular é metade do rotacional.
Argumentos de entrada
X, Y, Z- Coordenadas da grade, como de
meshgrid. U, V, W- Componentes do campo vetorial, com o tamanho da grade.
Argumentos de saída
cz- Componente z do rotacional,
dV/dx - dU/dy. cav- Velocidade angular: metade do rotacional.
Exemplo
Rotação contra campo radial
[X, Y] = meshgrid(1:4);
cz = curl(X, Y, -Y, X)
cz0 = curl(X, Y, X, Y)Saída
cz = 4×4
2 2 2 2
2 2 2 2
2 2 2 2
2 2 2 2
cz0 = 4×4
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0