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

De MateWiki
Saltar a: navegación, buscar

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434
Línea 89: Línea 89:
 
N = [-T(2, :); T(1, :)];
 
N = [-T(2, :); T(1, :)];
 
% Representación de los campos tangente y normal
 
% Representación de los campos tangente y normal
plot(x, y, 'b-', 'LineWidth', 2); hold on;
+
hold on;
 +
plot(x, y, 'b-', 'LineWidth', 2);  
 
quiver(x, y, T(1, :), T(2, :), 0.5,'k');
 
quiver(x, y, T(1, :), T(2, :), 0.5,'k');
 
quiver(x, y, N(1, :), N(2, :), 0.5, 'r');
 
quiver(x, y, N(1, :), N(2, :), 0.5, 'r');
Línea 96: Línea 97:
 
xlabel('Eje x'); ylabel('Eje y');
 
xlabel('Eje x'); ylabel('Eje y');
 
grid on;
 
grid on;
 +
hold off
 
}}
 
}}
  

Revisión del 17:07 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