Определение собственных функций закрепленной мембраны
Вычислим первые шесть собственных функций для круглой мембраны с зажатыми краями.
Укажем оператор Лапласа.
In[1]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/In_33.png)
\[ScriptCapitalL] = -Laplacian[u[x, y], {x, y}];
Укажем граничное условие Дирихле.
In[2]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/In_34.png)
\[ScriptCapitalB] = DirichletCondition[u[x, y] == 0, True];
Найдем шесть наименьших собственных значений и собственных функций.
In[3]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/In_35.png)
{vals, funs} =
NDEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]},
u[x, y], {x, y} \[Element] Disk[], 6];
Рассмотрим полученные собственные значения.
In[4]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/In_36.png)
vals
Out[4]=
![](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/O_16.png)
Визуализируем собственные функции.
In[5]:=
![Click for copyable input](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/In_37.png)
Table[Plot3D[funs[[i]], {x, y} \[Element] Disk[], PlotRange -> All,
PlotLabel -> vals[[i]], PlotTheme -> "Minimal"], {i, Length[vals]}]
Out[5]=
![](assets.ru/compute-eigenfunctions-for-a-clamped-membrane/O_17.png)