Wolfram言語

微分固有系

車の音響固有モードを解析する

Miniの切断面の近似に対する音響の固有値と固有関数を計算する.

切断面の画像をインポートする.

In[1]:=
Click for copyable input
img = Import[ "http://upload.wikimedia.org/wikipedia/commons/d/d3/Mini_cross_\ section.jpg"]
Out[1]=

マスクツールを使って境界のグラフィックスを生成する.

In[2]:=
Click for copyable input
boundary = \!\(\* GraphicsBox[ {GrayLevel[1.], {AbsoluteThickness[1], LineBox[CompressedData[" 1:eJxdkmtI0wEUxZdKWCaazxQlHyQLRckXQmQH+2sWCs5HivlhSYaKjxbF8rU0 0U2n8ndzsyznllvbSkrRAidGPjPMkmI0bRk+FlqZm8+cJRXkhD4cLj8OF+65 HM/0/PgMMwqF4vtHf+eGPMHP+YQAI1FXm8O8ROippM+RZbdQwGJaJa+IkFuR PX+OK4WZnvh+1L0FfBubH8nxcnQN2TblTonwOmtLGapWwImhMKzSRZj+1TH+ sl0JSWm43WJhM2SF7OFM3UP4ExzZYIgQbg0B28v2SkjfqvoNVkLwGVGqwHQF Prkk0iw3BDh2kSkefyMH60HstdHEBkTqJlvv2MuR8ozRulRCmtisjJFoZVWH HOPq2GagHNTysPMSXg1SDZ19A49bccPrcg7dkUSed5BPcVErKL1m40lfuVCU HNakbksR10d+LDNyMOujp8ZsyiCO/bly6Ew1NlLyay2GZIArm+/DrsSySxRl ZFQEbd7S8RYLNiJmInhKYzP0jpsVRmYV0kati4JKxejKoYbsK+Eg/5ECE+/E mFx/RewnuDi9NaRJGpQgO6No8ayBh6DCdHV16T1YsuJdC+oEJs7iegZ/8BDC o/GFjc5vl3f2/WfiZheOCNCr0TYI7CS43+ZQfEHLQ/DTlTg1XWzKR6RmVIX3 3MZMwMnKfjmJ2Rhhvpa6yzv+8HuxuattLdzWLz0f824y3eumUhPlp+6a/qFq 2+vr2XEH1K3U+NAEDtbqMzfSappwfSB5zvkmG9bmBxdotHroadV7IlPY//KT cGdqHCZ4JJTd7QmyKR5Yn6f0jiE86ER2HtPd9WBFO3XqCT6i5688qfUjYama FQzk8VEecICQOpNIyllb/vKtEf/38zd0gE9/ "]]}}, Background->GrayLevel[0.], ContentSelectable->True, ImagePadding->{{0., 0.}, {0., 0.}}, ImageSize->{273.3999999999956, 170.4}, PlotRange->{{0, 3910.}, {0, 2535.}}, PlotRangePadding->Automatic]\);

グラフィックスを離散化する.

In[3]:=
Click for copyable input
bdr = BoundaryDiscretizeGraphics[boundary]
Out[3]=

切断面の固有値と固有関数を6個計算する.

In[4]:=
Click for copyable input
{vals, funs} = NDEigensystem[{-Laplacian[u[x, y], {x, y}]}, u[x, y], {x, y} \[Element] bdr, 6];

固有値を調べる.

In[5]:=
Click for copyable input
vals
Out[5]=

車の切断面の2つ目の固有関数を可視化する.

完全なWolfram言語入力を表示する
In[6]:=
Click for copyable input
Show[img, ContourPlot[funs[[2]], {x, y} \[Element] bdr, Axes -> None, Frame -> None, AspectRatio -> Automatic, ColorFunction -> Function[f, {Opacity[0.75], ColorData["TemperatureMap"][f]}]], ImageSize -> Automatic]
Out[6]=

関連する例

de en es fr ko pt-br ru zh