Wolfram 语言™

创建球内的拉普拉斯算子特征函数图集

定义一个三维拉普拉斯算子.

In[1]:=

\[ScriptCapitalL] = -Laplacian[u[x, y, z], {x, y, z}];

设定齐次狄利克雷边界条件.

In[2]:=

\[ScriptCapitalB] = DirichletCondition[u[x, y, z] == 0, True];

求球内最小的十六个特征值和特征函数.

In[3]:=

\[CapitalOmega] = Ball[{0, 0, 0}, 2];

In[4]:=

{vals, funs} = DEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]}, u[x, y, z], {x, y, z} \[Element] \[CapitalOmega], 16];

特征值以 BesselJZero 的形式给出.

In[5]:=

vals[[1]] // TraditionalForm

Out[5]//TraditionalForm=

生成特征函数图集.

In[6]:=

Grid[Partition[ ParallelTable[ DensityPlot3D[ funs[[i]] // N // Evaluate, {x, y, z} \[Element] \[CapitalOmega], Boxed -> False, Axes -> False, ColorFunction -> Hue, Method -> {"ShrinkWrap" -> True}, ImageSize -> 125], {i, 16}], 4]]

Out[6]=