Определение символических собственных значений
Укажем одномерный лапласовский оператор.
In[1]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_92.png)
\[ScriptCapitalL] = -Laplacian[u[x], {x}];
Зададим однородное граничное условие Дирихле.
In[2]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_93.png)
\[ScriptCapitalB] = DirichletCondition[u[x] == 0, True];
Найдем выражения для 5 самых маленьких собственных значений на отрезке .
In[3]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_94.png)
DEigenvalues[{\[ScriptCapitalL], \[ScriptCapitalB]}, u[x], {x, a, b},
5]
Out[3]=
![](assets.ru/compute-symbolic-eigenvalues/O_43.png)
Укажем оператор Эйри.
In[4]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_95.png)
\[ScriptCapitalL] = -Laplacian[u[x], {x}] + x u[x];
Найдем 5 наименьших собственных значений и соответствующих собственных функций.
In[5]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_96.png)
{vals, funs} =
DEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]},
u[x], {x, 0, 1}, 5];
Собственные значения представляют собой корни трансцендентного уравнения.
In[6]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_97.png)
vals[[1]] // TraditionalForm
Out[6]//TraditionalForm=
![](assets.ru/compute-symbolic-eigenvalues/O_44.png)
С точностью определим трансцендентальное собственное значение.
In[7]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_98.png)
N[vals[[1]], 500] // TraditionalForm
Out[7]//TraditionalForm=
![](assets.ru/compute-symbolic-eigenvalues/O_45.png)
Визуализируем собственные функции.
In[8]:=
![Click for copyable input](assets.ru/compute-symbolic-eigenvalues/In_99.png)
Plot[Evaluate[funs + Range[5]], {x, 0, 1}, ImageSize -> Medium,
PlotTheme -> {"Business", "Bare"}, AspectRatio -> 1]
Out[8]=
![](assets.ru/compute-symbolic-eigenvalues/O_46.png)