Определение собственных функций для L-образной области
Укажем область с L-образной формой.
In[1]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/In_38.png)
L = Polygon[{{1, 0}, {2, 0}, {2, 2}, {0, 2}, {0, 1}, {1, 1}}];
Укажем оператор Лапласа.
In[2]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/In_39.png)
\[ScriptCapitalL] = Laplacian[u[x, y], {x, y}];
Укажем граничное условие Дирихле.
In[3]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/In_40.png)
\[ScriptCapitalB] = DirichletCondition[u[x, y] == 0., True];
Определим собственные функций для L-образной области.
In[4]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/In_41.png)
{vals, funs} =
NDEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]},
u[x, y], {x, y} \[Element] L, 6];
Рассмотрим полученные собственные значения.
In[5]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/In_42.png)
vals
Out[5]=
![](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/O_18.png)
Визуализируем собственные функции.
In[6]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/In_43.png)
Plot3D[#, {x, y} \[Element] L, PlotPoints -> 75, Mesh -> None,
PlotStyle -> Directive[Orange, Specularity[White, 30]],
BoxRatios -> {1, 1, 0.8}] & /@ funs
Out[6]=
![](assets.ru/compute-eigenfunctions-in-an-l-shaped-region/O_19.png)