1Dラプラス演算子の記号的固有関数を求める
1Dのラプラス演算子を指定する.
In[1]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_81.png)
\[ScriptCapitalL] = -Laplacian[u[x], {x}];
固有関数に対する同次ディリクレ境界条件を指定する.
In[2]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_82.png)
\[ScriptCapitalB]1 = DirichletCondition[u[x] == 0, True];
固有値と固有関数を小さいものから5個求める.
In[3]:=
![Click for copyable input](assets.ja/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.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_84.png)
vals
Out[4]=
![](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/O_37.png)
固有関数を調べる.
In[5]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_85.png)
funs
Out[5]=
![](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/O_38.png)
固有関数を可視化する.
In[6]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_86.png)
Plot[Evaluate[funs + 2 Range[5]], {x, 0, \[Pi]}]
Out[6]=
![](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/O_39.png)
同次ノイマン(Neumann)境界条件を指定する.
In[7]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_87.png)
\[ScriptCapitalB]2 = NeumannValue[0, True];
固有値と固有関数を小さいものから5個求める.
In[8]:=
![Click for copyable input](assets.ja/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.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_89.png)
vals
Out[9]=
![](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/O_40.png)
また,固有関数で正弦の代りに余弦が使われている.
In[10]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_90.png)
funs
Out[10]=
![](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/O_41.png)
固有関数を可視化する.
In[11]:=
![Click for copyable input](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/In_91.png)
Plot[Evaluate[funs + 2 Range[5]], {x, 0, \[Pi]}]
Out[11]=
![](assets.ja/find-a-1d-laplacians-symbolic-eigenfunctions/O_42.png)