Diferencia entre revisiones de «Ecuación de ondas»

Revisión del 17:44 23 may 2024

1 Introducción

2 Modelización de los desplazamientos transversales

Para modelar el comportamiento de los desplazamientos transversales de la cuerda vibrante, utilizamos la ecuación de ondas en una dimensión. Dado que la cuerda está fija en los extremos y tiene una densidad [math] d [/math] y tensión constante [math] \tau_0 [/math] de manera que la velocidad de propagación es [math] c = \tau_0/d = 1 [/math], la ecuación de ondas se simplifica.

La ecuación de ondas en una dimensión para los desplazamientos transversales [math] u(x,t) [/math] de la cuerda es:

Dado que [math] c = 1 [/math], la ecuación se reduce a:

Esta ecuación diferencial parcial (EDP) debe acompañarse de las condiciones de contorno y las condiciones iniciales para estar completamente especificada.

2.1 Condiciones de contorno

Dado que la cuerda está fija en los extremos, las condiciones de contorno son:

2.2 Condiciones iniciales

Las condiciones iniciales especifican la posición inicial de la cuerda [math] u_0(x) [/math] y su velocidad inicial o impulso [math] u_1(x) [/math]:

2.3 Sistema de EDP

Juntando la ecuación de ondas con las condiciones de contorno e iniciales, el sistema completo que modela el comportamiento de los desplazamientos transversales de la cuerda es:

Este sistema describe completamente la evolución temporal de los desplazamientos transversales de una cuerda vibrante con los extremos fijos, dada su posición e impulso iniciales.

2.4 Resolución del sistema por separación de variables en términos de los coeficientes de Fourier de los datos iniciales

Primero, planteamos la solución de la ecuación de onda mediante separación de variables. Asumimos que la solución \( u(x,t) \) se puede escribir como un producto de dos funciones independientes, una que depende solo de \( x \) y otra que depende solo de \( t \):

Sustituimos esta forma en la ecuación de onda:

Esto se simplifica a:

Dividimos ambos lados por [math]X(x) T(t) [/math]:

Dado que el lado izquierdo depende solo de \( t \) y el lado derecho solo de \( x \), ambos deben ser iguales a una constante que llamaremos \(-\lambda\). Esto nos lleva a dos ecuaciones ordinarias:

Para satisfacer las condiciones de contorno \( u(0, t) = 0 \) y \( u(1, t) = 0 \), tenemos:

La ecuación para \( X(x) \) tiene la solución general:

Aplicando las condiciones de contorno, obtenemos que \( B = 0 \) y que \( \sqrt{\lambda} \) debe ser un múltiplo de \( \pi \):

Por lo tanto, \( \lambda = (n \pi)^2 \) y las funciones propias son:

La ecuación para \( T(t) \) se convierte en:

Con la solución general:

Combinando \( X(x) \) y \( T(t) \), la solución general para \( u(x,t) \) es:

Donde \( A_n \) y \( B_n \) son coeficientes a determinar a partir de las condiciones iniciales.

1. 1. 1. Determinación de los coeficientes \( A_n \) y \( B_n \)

Las condiciones iniciales son:

Para \( t = 0 \):

Usamos la ortogonalidad de las funciones seno para determinar \( A_n \):

Para la derivada temporal:

Por lo tanto:

Finalmente, la solución de la ecuación de onda en términos de los coeficientes de Fourier de los datos iniciales es:

Esta expresión proporciona la solución completa de la ecuación de onda en función de las condiciones iniciales \( u_0(x) \) y \( u_1(x) \).

2.4.1 Ejemplo de solución periódica en tiempo

En esta parte del documento para comprender lo anteriormente explicado de manera óptima se va a proceder a analizar una serie de ejemplo. El primero de ellos se trata de la representación gráfica de una solución periódica en tiempo. Esta viene dada suponiendo que los datos iniciales son [math]u_0(x)=e^{-100(x-\frac{1}{2})^2}[/math] y [math]u_1(x)=0[/math]. Es por ello que la representación de la solución en el intervalo t ∈ [0, 2] es la siguiente:

Para la obtención de esta representación se ha creado el siguiente código en Matlab:

% 
clear
close all
clc

% Definimos nuestros datos iniciales:
u_0=@(x) (exp(-100.*(x-1/2).^2));
u_1=@(x) (0);

% Definimos nuestro intervalo de tiempo t:
t=0:0.001:2;

% Definimos nuestro intervalo x:
x=0:0.001:1;

% Definimos el número de términos de la serie a dibujar:
n=50;

% Calculamos los términos de fourier asociados a u_0:
ak_0=zeros(n,1);
bk_0=zeros(n,1);
ck_0=zeros(n,1);

f_sen_0=@(x,k) u_0(x).*sin(k*pi*x);
f_cos_0=@(x,k) u_0(x).*cos(k*pi*x);
f_medio_0=@(x,k) u_0(x).*1/2;

for k=1:n
    % Vector con los valores de f_sin:
    y_sen_0=f_sen_0(x,k);
            
    % Vector con los valores de f_cos:
    y_cos_0=f_cos_0(x,k);
            
    % Vector con los valores de f_sin:
    y_medio_0=f_medio_0(x);

    ak_0(k)=trapz(x,y_sen_0);
    bk_0(k)=trapz(x,y_cos_0);
    ck_0(k)=trapz(x,y_medio_0);
end

% Calculamos los términos de fourier asociados a u_1:
ak_1=zeros(n,1);
bk_1=zeros(n,1);
ck_1=zeros(n,1);

f_sen_1=@(x,k) u_1(x).*sin(k*pi*x);
f_cos_1=@(x,k) u_1(x).*cos(k*pi*x);
f_medio_1=@(x,k) u_1(x).*1/2;

for k=1:n
    % Vector con los valores de f_sin:
    y_sen_1=f_sen_1(x,k);
            
    % Vector con los valores de f_cos:
    y_cos_1=f_cos_1(x,k);
            
    % Vector con los valores de f_sin:
    y_medio_1=f_medio_1(x);

    ak_1(k)=trapz(x,y_sen_1);
    bk_1(k)=trapz(x,y_cos_1);
    %ck_1(k)=trapz(x,y_medio_1);
end

% Sustituimos en la solución calculada por separación de variables:
u=@(xx,tt) 0;
for i=1:n
    u=@(xx,tt) u(xx,tt) + (ak_0(i).*cos(i*pi*tt)+bk_1(i).*sin(i*pi*tt)).*sin(i*pi*xx);
end

% Graficamos la solución:
[X,T]=meshgrid(x,t);
surf(X,T,u(X,T),'EdgeColor','flat')
title('Solución ecuación de ondas')
xlabel('x')
ylabel('t')

Como se observa en la imagen anterior la función muestra un comportamiento que se repite regularmente con el tiempo. A este tipo de funciones se les denomina periódicas en tiempo. Esto se puede observar de manera óptima para su comprensión en la siguiente gráfica, en la cual se observa el comportamiento explicado a anteriormente.

centro

2.4.2 Ejemplo de la onda viajera

En este caso se va