求一维拉普拉斯算子的符号特征函数
设定一个一维拉普拉斯算子.
In[1]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_81.png)
\[ScriptCapitalL] = -Laplacian[u[x], {x}];
设定特征函数的齐次狄利克雷边界条件.
In[2]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_82.png)
\[ScriptCapitalB]1 = DirichletCondition[u[x] == 0, True];
求最小的五个特征值和特征函数.
In[3]:=
![Click for copyable input](assets.zh/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.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_84.png)
vals
Out[4]=
![](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/O_37.png)
查看特征函数.
In[5]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_85.png)
funs
Out[5]=
![](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/O_38.png)
可视化特征函数.
In[6]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_86.png)
Plot[Evaluate[funs + 2 Range[5]], {x, 0, \[Pi]}]
Out[6]=
![](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/O_39.png)
设定齐次诺伊曼边界条件.
In[7]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_87.png)
\[ScriptCapitalB]2 = NeumannValue[0, True];
求最小的五个特征值和特征函数.
In[8]:=
![Click for copyable input](assets.zh/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.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_89.png)
vals
Out[9]=
![](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/O_40.png)
特征函数中正弦取代了余弦.
In[10]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_90.png)
funs
Out[10]=
![](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/O_41.png)
可视化特征函数.
In[11]:=
![Click for copyable input](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/In_91.png)
Plot[Evaluate[funs + 2 Range[5]], {x, 0, \[Pi]}]
Out[11]=
![](assets.zh/find-a-1d-laplacians-symbolic-eigenfunctions/O_42.png)