计算一个 L-形区域的特征函数

设定一个 L-形区域. .

In[1]:=

L = Polygon[{{1, 0}, {2, 0}, {2, 2}, {0, 2}, {0, 1}, {1, 1}}];

设定拉普拉斯算子.

In[2]:=

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

设定狄利克雷边界条件.

In[3]:=

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

计算 L-形区域上的特征函数.

In[4]:=

{vals, funs} = NDEigensystem[{\[ScriptCapitalL], \[ScriptCapitalB]}, u[x, y], {x, y} \[Element] L, 6];

查看特征值.

In[5]:=

vals

Out[5]=

可视化特征函数.

In[6]:=

Plot3D[#, {x, y} \[Element] L, PlotPoints -> 75, Mesh -> None, PlotStyle -> Directive[Orange, Specularity[White, 30]], BoxRatios -> {1, 1, 0.8}] & /@ funs

Out[6]=