4.3: Ecuaciones no homogéneas
- Page ID
- 118158
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Aquí consideramos el problema de valor inicial
\ begin {eqnarray}
\ label {inh1}
\ Box u&=&w (x, t)\\ mbox {on}\ x\ in\ mathbb {R} ^n,\ t\ in\ mathbb {R} ^1\
\ etiqueta {inh2}
u (x,0) &=&f (x)\\
\ etiqueta {xh3}
u_t (,0) &=&g (x),
\ end {eqnarray}
donde\(\Box u:=u_{tt}-c^2\triangle u\). Asumimos\(f\in C^3\),\(g\in C^2\) y\(w\in C^1\), que se dan.
Set\(u=u_1+u_2\), donde\(u_1\) es una solución de problema (\ ref {inh1}) - (\ ref {inh3}) con\(w:=0\) y\(u_2\) es la solución donde\(f=0\) y\(g=0\) en (\ ref {inh1}) - (\ ref {inh3}). Ya que tenemos soluciones explícitas\(u_1\) en los casos\(n=1\)\(n=3\),\(n=2\) y, queda por resolver
\ begin {eqnarray}
\ label {duhu2gl}
\ Box u&=&w (x, t)\\\ mbox {on}\ x\ in\ mathbb {R} ^n,\ t\ in\ mathbb {R} ^1\
\ etiqueta {duhu2in1}
u (x,0) &=&0\
\
uetiqueta {duhu2in2} _t (x,0) &=&0.
\ end {eqnarray}
El siguiente método se llama principio de Duhamel que puede considerarse como una generalización del método de variaciones de constantes en la teoría de ecuaciones diferenciales ordinarias.
Para resolver este problema, hacemos el ansatz
\ begin {ecuación}
\ label {duh1}
u (x, t) =\ int_0^t\ v (x, t, s)\ ds,
\ end {ecuación}
donde\(v\) es una función que satisface
\ begin {ecuación}
\ label {duh2}
\ Box v=0\\ mbox {para todos}\ s
\ end {ecuación}
y
\ begin {ecuación}
\ label {duh3}
v (x, s, s) =0.
\ end {ecuación}
De ansatz (\ ref {duh1}) y suposición (\ ref {duh3}) obtenemos
\ begin {eqnarray}
u_t&=&v (x, t, t) +\ int_0^t\ v_t (x, t, s)\ ds,\ nonumber\
\ label {duh4}
&=&\ int_0^t\ v_t (x, t, s).
\ end {eqnarray}
De ello se desprende\(u_t(x,0)=0\). La condición inicial\(u(x,t)=0\) se cumple por el ansatz (\ ref {duh1}). De (\ ref {duh4}) y ansatz (\ ref {duh1}) vemos que
\ begin {eqnarray*}
u_ {tt} &=&v_t (x, t, t) +\ int_0^t\ v_ {tt} (x, t, s)\ ds,\
\ triangle_x u&=&\ int_0^t\\ triangle_x v (x, t, s)\ ds.
\ end {eqnarray*}
Por lo tanto, dado que\(u\) es un ansatz para (\ ref {duhu2gl}) - (\ ref {duhu2in2}),
\ begin {eqnarray*}
u_ {tt} -c^2\ triangle_x u&=&v_t (x, t, t) +\ int_0^t (\ Box v) (x, t, s)\ ds\\
&=&w (x, t).
\ end {eqnarray*}
Así necesariamente\(v_t(x,t,t)=w(x,t)\), ver (\ ref {duh2}). Hemos visto que el ansatz proporciona una solución de (\ ref {duhu2gl}) - (\ ref {duhu2in2}) si por todos\(s\)
\ begin {ecuación}
\ etiqueta {duh5}
\ Caja v=0,\\ v (x, s, s) =0,\\ v_t (x, s, s) =w (x, s).
\ end {ecuación}
Dejar\(v^*(x,t,s)\) ser una solución de
\ begin {ecuación}
\ label {duh6}
\ Box v=0,\\ v (x,0, s) =0,\\ v_t (x,0, s) =w (x, s),
\ end {ecuación}
entonces
$$v (x, t, s) :=v^* (x, t-s, s)\]
es una solución de (\ ref {duh5}).
En el caso\(n=3\), donde\(v^*\) está dado por, ver Teorema 4.2,
$$v^* (x, t, s) =\ frac {1} {4\ pi c^2 t}\ int_ {\ parcial B_ {ct} (x)}\ w (\ xi, s)\ dS_\ xi.\]
Entonces
\ begin {eqnarray*}
v (x, t, s) &=&v^* (x, t-s, s)\\
&=&\ frac {1} {4\ pi c^2 (t-s)}\ int_ {\ parcial B_ {c (t-s)} (x)}\ w (\ xi, s)\ dS_\ xi.
\ end {eqnarray*}
de ansatz (\ ref {duh1}) sigue
\ begin {eqnarray*}
u (x, t) &=&\ int_0^t\ v (x, t, s)\ ds\\
&=&\ frac {1} {4\ pi c^2}\ int_0^t\\ int_ {\ parcial B_ {c (t-s)} (x)}\\ frac {w (\ xi, s)} t-s}\ dS_\ xi ds.
\ end {eqnarray*}
Cambio de variables por\(\tau=c(t-s)\) rendimientos
\ begin {eqnarray*}
u (x, t) &=&\ frac {1} {4\ pi c^2}\ int_0^ {ct}\\ int_ {\ int_ {\ parcial B_ {\ tau} (x)}\\ frac {w (\ xi, t-\ tau/c)} {\ tau}\ ds_\ xi d\ tau\\
&=& ac {1} {4\ pi c^2}\ int_ {B_ {ct} (x)}\\ frac {w (\ xi, t-r/c)} {r}\ d\ xi,
\ end {eqnarray*}
donde\(r=|x-\xi|\).
Fórmulas para los casos\(n=1\) y\(n=2\) seguir de fórmulas para la ecuación homogénea asociada con valores iniciales no homogéneos para estos casos.
Teorema 4.4. La solución de
$$\ Caja u=w (x, t),\\ u (x,0) =0,\\ u_t (x,0) =0, $$
donde\(w\in C^1\), viene dada por:
Caso\(n=3\):
$$u (x, t) =\ frac {1} {4\ pi c^2}\ int_ {B_ {ct} (x)}\\ frac {w (\ xi, t-r/c)} {r}\ d\ xi, $$
donde\(r=|x-\xi|\),\(x=(x_1,x_2,x_3)\),\(\xi=(\xi_1,\xi_2,\xi_3)\).
Caso\(n=2\):
$$u (x, t) =\ frac {1} {4\ pi c}\ int_0^t\\ izquierda (\ int_ {B_ {c (t-\ tau)} (x)}\\ frac {w (\ xi,\ tau)} {\ sqrt {c^2 (t-\ tau) ^2-r^2}}\ d\ xi\ derecha)\\ tau, $$
\(x=(x_1,x_2)\),\(\xi=(\xi_1,\xi_2)\).
Caso\(n=1\):
$$u (x, t) =\ frac {1} {2c}\ int_0^t\\ izquierda (\ int_ {x-c (t-\ tau)} ^ {x+c (t-\ tau)}\ w (\ xi,\ tau)\ d\ xi\ derecha)\ d\ tau. $$
OBSERVACIÓN. El integrando a la derecha en fórmula para\(n=3\) se llama potencial retardado. El integrando no se toma en\(t\), se toma en un momento anterior\(t-r/c\).