Genere una expansión de función propia
Calcule la expansión de función propia de la función con respecto a la base proporcionada por el operador laplaciano con condiciones de límite de Dirichlet en el intervalo
.
In[1]:=
![Click for copyable input](assets.es/generate-an-eigenfunction-expansion/In_109.png)
basis = DEigensystem[{-Laplacian[u[x], {x}],
DirichletCondition[u[x] == 0, True]}, u[x], {x, 0, \[Pi]}, 6,
Method -> "Normalize"][[2]]
Out[1]=
![](assets.es/generate-an-eigenfunction-expansion/O_54.png)
Calcule los coeficientes de Fourier de la función .
In[2]:=
![Click for copyable input](assets.es/generate-an-eigenfunction-expansion/In_110.png)
f[x_] := E^(-x) x^2 (\[Pi] - x) Sin[4 x]
In[3]:=
![Click for copyable input](assets.es/generate-an-eigenfunction-expansion/In_111.png)
coeffs = (Table[Integrate[f[x] basis[[i]], {x, 0, Pi}], {i, 6}] //
FullSimplify);
Defina como la
suma parcial de la función.
In[4]:=
![Click for copyable input](assets.es/generate-an-eigenfunction-expansion/In_112.png)
eigexp[x_, n_] := Sum[coeffs[[i]] basis[[i]], {i, n}]
In[5]:=
![Click for copyable input](assets.es/generate-an-eigenfunction-expansion/In_113.png)
eigexp[x, 3] // N
Out[5]=
![](assets.es/generate-an-eigenfunction-expansion/O_55.png)
Compare la función con su expansión de función propia de distintos valores de .
In[6]:=
![Click for copyable input](assets.es/generate-an-eigenfunction-expansion/In_114.png)
Table[Plot[{f[x], eigexp[x, i]} // Evaluate, {x, 0, Pi}], {i, 3, 6}]
Out[6]=
![](assets.es/generate-an-eigenfunction-expansion/O_56.png)