Diferencia entre revisiones de «Curvas de Bézier (Grupo 32)»

De MateWiki
Saltar a: navegación, buscar
(→‎Representación de la curvatura de la curva en función del parámetro t)
Línea 104: Línea 104:
 
[[Categoría:Teoría de Campos]]
 
[[Categoría:Teoría de Campos]]
 
[[Categoría:TC24/25]]
 
[[Categoría:TC24/25]]
 +
 +
En este apartado se adjunta una animación en Matlab que representa la dinámica de los vectores tangente, normal y de la circunferencia osculatriz mientras se recorre la curva de Bézier.
 +
Se observa que el vector tangente indica la dirección de la curva en cada punto; el vector normal, perpendicular al vector tangente apunta hacia el centro de la curva. Por otro lado, la circunferencia osculatriz es el círculo que se ajusta de la mejor manera posible a una curva en un punto específico, de manera que tiene la misma dirección (tangente) y curvatura que la curva en ese punto.
 +
La curvatura κ(t) de una curva parametrizada en el plano  x(t) e  y(t)  se calcula con la siguiente fórmula:
 +
k(t)=|x^' (t) y^'' (t)-y(t) x^'' (t)|/(〖〖(x〗^' (t)〗^2+〖y^' (t)〗^2 )^(3/2) )
 +
Donde:
 +
x(t)y y(t) son las componentes x e y  son las componentes de la curva en función del parámetro t.
 +
x'(t)y y'(t) son las primeras derivadas de las componentes x(t) e  y(t)
 +
x''(t)y y''(t) son las segundas derivadas de las componentes x(t) e y(t)
 +
La curvatura depende del cambio de la tangente a la curva en un punto, y el radio de curvatura. Se puede observar que los puntos con mayor radio en la circunferencia osculatriz, son los mismos que tienen menor curvatura. Por lo que se concluye que la relación entre la curvatura de una curva κ(t) y el radio de la circunferencia osculatriz R(t) es inversa.

Revisión del 18:58 5 dic 2024

Trabajo realizado por estudiantes
Título Derformación plana. Grupo 32
Asignatura Teoría de Campos
Curso 2024-25
Autores Nombres: Rocío Jamileth Ruiz Herrera, Mario Del Amo Domínguez, Diana Estefanía Sagal Tituaña, Jesús Gil Gutierrez y David Bretaña Blanco
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

Las cuervas de Bézier deben su nombre al ingeniero francés Pierre Bézier, quien las presentó en 1962 y posteriormente utilizó para el diseño de partes de automóviles para Renault. Actualmente, estas curvas son indispensables para la industria de la gráfica por computadora, el diseño industrial y la ingeniería, facilitando la creación de formas más fluidas y precisas.

Las curvas de Bézier de orden [math]n[/math] están definidas por los puntos de control [math]P_0,P_1,...,P_n[/math] y se pueden expresar mediante la siguiente fórmula:

[math] B(t)=\sum_{i=0}^n B_{i,n}(t) P_i\ [/math]

donde \(B_{i,n}(t)\) son los polinomios de Bernstein, dados por:

[math] B_{i,n}(t) = \binom{n}{i} t^i (1-t)^{n-i}\ [/math]

para \(t \in [0, 1]\), y donde \(\binom{n}{i}\) es el coeficiente binomial.

1 Representación de la curva de Bézier cúbica (n=3) junto con la curva poligonal que conecta los cuatro puntos coplanarios

Primero, seleccionamos cuatro puntos coplanarios cualesquiera. En este caso conseguimos los puntos coplanarios en el plano z=0 que posteriormente sustituiremos en la formula de Bézier:


[math] B(t) = (1-t)^3P_0 + 3t(1-t)^2P_1 + 3t^2(1-t)P_2 + t^3P_3, \; t \in [0,1] [/math]

Para mejor manejo de la fórmula la expresaremos por componentes en la base física cartesiana [x,y,z]:

[math] x(t) = (1 - t)^3 P_{0x} + 3(1 - t)^2 t P_{1x} + 3(1 - t) t^2 P_{2x} + t^3 P_{3x}[/math]

[math] y(t) = (1 - t)^3 P_{0y} + 3(1 - t)^2 t P_{1y} + 3(1 - t) t^2 P_{2y} + t^3 P_{3y}[/math]

[math] z(t) = (1 - t)^3 P_{0z} + 3(1 - t)^2 t P_{1z} + 3(1 - t) t^2 P_{2z} + t^3 P_{3z}.[/math]

Posteriormente elegimos los puntos coplanarios:
[math] P0=(0,0,0),P1=(1,2,0),P2=(4,3,0),P3=(3,0,0) [/math].

Tras sustituir los puntos en la fórmula obtenemos una gráfica en la que se puede observar la curva de Bézier de los cuatro puntos junto con la curva poligonal que une los puntos.

Curva de Bézier
%Puntos de control coplanarios
P0 = [0, 0,0];
P1 = [1, 2,0];
P2 = [4, 3,0];
P3 = [3, 0,0];
P = [P0; P1; P2; P3];
t = linspace(0, 1, 100);
% Polinomios de Bernstein
B0 = (1 - t).^3;
B1 = 3 * t .* (1 - t).^2;
B2 = 3 * t.^2 .* (1 - t);
B3 = t.^3;
% Coordenadas de la curva de Bézier
x = B0 * P0(1) + B1 * P1(1) + B2 * P2(1) + B3 * P3(1);
y = B0 * P0(2) + B1 * P1(2) + B2 * P2(2) + B3 * P3(2);
% Representación gráfica
plot(x, y, 'b', 'LineWidth', 2); hold on;
plot(P(:, 1), P(:, 2), 'ro--', 'LineWidth', 1.5); 
% Polígono de control
legend('Curva de Bézier', 'Curva poligonal');
title('Curva de Bézier Cúbica');
axis equal;
grid on;
xlabel('Eje x'); ylabel('Eje y');

2 Representación del campo tangente [math]T(t)[/math] y del campo normal [math]N(t)[/math] en varios puntos de la curva

A continuación calcularemos el campo tangente y y el vector normal de la curva de Bézier. Para obtener la tangente debemos derivar la función y dividirla por su módulo:

[math]\mathbf{T}(t)= \frac{\frac{dB(t)}{dt} }{\|\frac{dB(t)}{dt}\|}[/math]

Para mejor manejo lo podemos expresar en la base física cartesiana[x,y]

[math] dx(t)/dt = -3(1-t)^2 x_0 + 3(1-t)^2 x_1 - 6t(1-t) x_1 + 6t(1-t) x_2 - 3t^2 x_2 + 3t^2 x_3 [/math]
[math] dy(t)/dt = -3(1-t)^2 y_0 + 3(1-t)^2 y_1 - 6t(1-t) y_1 + 6t(1-t) y_2 - 3t^2 y_2 + 3t^2 y_3 [/math]

Como el vector normal es ortogonal al vector tangencial, el campo normal lo obtenemos rotando 90º el vector normal:

[math]\mathbf{N}(t) ={[ -T_y(t), T_x(t)]}[/math]

Siendo [math]T_x(t)[/math] y [math]T_y(t)[/math] las derivadas q componen el campo tangente en la base física cartesiana.

Campo tangente y Campo normal
%Puntos de control coplanarios
P0 = [0, 0,0];
P1 = [1, 2,0];
P2 = [4, 3,0];
P3 = [3, 0,0];
P = [P0; P1; P2; P3];
t = linspace(0, 1, 30);
% Polinomios de Bernstein
B0 = (1 - t).^3;
B1 = 3 * t .* (1 - t).^2;
B2 = 3 * t.^2 .* (1 - t);
B3 = t.^3;
% Coordenadas de la curva de Bézier
x = B0 * P0(1) + B1 * P1(1) + B2 * P2(1) + B3 * P3(1);
y = B0 * P0(2) + B1 * P1(2) + B2 * P2(2) + B3 * P3(2);
% Derivadas
dx = gradient(x, t);
dy = gradient(y, t);
% Campo tangente
T = [dx; dy] ./ sqrt(dx.^2 + dy.^2);
% Campo normal
N = [-T(2, :); T(1, :)];
% Representación de los campos tangente y normal
hold on;
plot(x, y, 'b-', 'LineWidth', 2); 
quiver(x, y, T(1, :), T(2, :), 0.5,'k');
quiver(x, y, N(1, :), N(2, :), 0.5, 'r');
legend('Curva de Bézier', 'Campo Tangente', 'Campo Normal');
title('Campo Tangente y Campo Normal');
xlabel('Eje x'); ylabel('Eje y');
grid on;
hold off


3 Representación de la curvatura de la curva en función del parámetro t

En este apartado se adjunta una animación en Matlab que representa la dinámica de los vectores tangente, normal y de la circunferencia osculatriz mientras se recorre la curva de Bézier. Se observa que el vector tangente indica la dirección de la curva en cada punto; el vector normal, perpendicular al vector tangente apunta hacia el centro de la curva. Por otro lado, la circunferencia osculatriz es el círculo que se ajusta de la mejor manera posible a una curva en un punto específico, de manera que tiene la misma dirección (tangente) y curvatura que la curva en ese punto. La curvatura κ(t) de una curva parametrizada en el plano x(t) e y(t) se calcula con la siguiente fórmula: 

k(t)=|x^' (t) y^ (t)-y(t) x^ (t)|/(〖〖(x〗^' (t)〗^2+〖y^' (t)〗^2 )^(3/2) )

Donde: x(t)y y(t) son las componentes x e y son las componentes de la curva en función del parámetro t. x'(t)y y'(t) son las primeras derivadas de las componentes x(t) e y(t) x(t)y y(t) son las segundas derivadas de las componentes x(t) e y(t) La curvatura depende del cambio de la tangente a la curva en un punto, y el radio de curvatura. Se puede observar que los puntos con mayor radio en la circunferencia osculatriz, son los mismos que tienen menor curvatura. Por lo que se concluye que la relación entre la curvatura de una curva κ(t) y el radio de la circunferencia osculatriz R(t) es inversa.

Obtenido de «https://mat.caminos.upm.es/w/index.php?title=Curvas_de_Bézier_(Grupo_32)&oldid=79152»