Одномерные символьные собственные функции лапласиана
Укажем одномерный лапласовский оператор.
In[1]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_81.png)
\[ScriptCapitalL] = -Laplacian[u[x], {x}];
Укажем однородные граничные условия Дирихле для собственных функций.
In[2]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_82.png)
\[ScriptCapitalB]1 = DirichletCondition[u[x] == 0, True];
Определим пять наименьших собственных значений и собственных функций.
In[3]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_83.png)
{vals, funs} =
DEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]1},
u[x], {x, 0, \[Pi]}, 5];
Проверим полученные собственные значения.
In[4]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_84.png)
vals
Out[4]=
![](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/O_37.png)
Проверим полученные собственные функции.
In[5]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_85.png)
funs
Out[5]=
![](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/O_38.png)
Визуализируем полученные собственные функции.
In[6]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_86.png)
Plot[Evaluate[funs + 2 Range[5]], {x, 0, \[Pi]}]
Out[6]=
![](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/O_39.png)
Укажем однородное граничное условие Неймана.
In[7]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_87.png)
\[ScriptCapitalB]2 = NeumannValue[0, True];
Найдем пять наименьших собственных значений и собственных функций.
In[8]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_88.png)
{vals, funs} =
DEigensystem[\[ScriptCapitalL] + \[ScriptCapitalB]2,
u[x], {x, 0, \[Pi]}, 5];
Проверим полученные собственные значения. Относительно условий Дирихле был добавлен нулевой режим.
In[9]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_89.png)
vals
Out[9]=
![](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/O_40.png)
Косинусы были заменены на синусы в полученных собственных функциях.
In[10]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_90.png)
funs
Out[10]=
![](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/O_41.png)
Визуализируем собственные функции.
In[11]:=
![Click for copyable input](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/In_91.png)
Plot[Evaluate[funs + 2 Range[5]], {x, 0, \[Pi]}]
Out[11]=
![](assets.ru/find-a-1d-laplacians-symbolic-eigenfunctions/O_42.png)