3
大学数学基礎解説
文献あり

どのような三角形の太鼓が一番良い音を鳴らすのか?

289
0

円盤、多角形、滑らかな境界を持つさまざまな形の太鼓を考えてみよう。
どのような形が最も良い音を響かせるだろうか?一般的に、「うなり」のある音は良い音とは感じにくい。太鼓と聞いて円盤を思い浮かべる人も多いかもしれないが、その円盤の対称性が「音のうなり」に関係している可能性がある。

まず、最も単純な図形のひとつである三角形について考えてみよう[図1]。
三角形の太鼓 三角形の太鼓

三角形の太鼓が最も良い音を鳴らすのは、どのような形の三角形のときだろうか?

本ノートでは、三角形の太鼓の中でも特に対称性の高い正三角形が、ある意味でうなりのない最も良い音を鳴らすことを数学的に証明することを試みる。

ラプラシアンのディリクレ固有値問題

太鼓の音は、ラプラシアンのディリクレ固有値問題によってモデル化される。倍音は各固有振動モードに対応しており、そのモードは膜に関連するラプラシアンのディリクレ固有値によって決まる。

ラプラシアンのディリクレ固有値問題は、有界領域 ΩR2 に対して次のように定式化される:
Δu=λuin Ω,u=0on Ω,
ここで u は振動モードに対応する固有関数であり、λ は対応する固有値である。2次元領域におけるラプラシアンΔは次式で定義される微分作用素である:
Δ:=2x2+2y2.
ラプラシアンの固有値は可算無限個存在し、非減少の実数列を形成する。第一固有値λ1は単純あり、λ1>0が成り立つ。すなわち固有値の列を次のように表すことができる:
0<λ1<λ2λ3.

うなりの周波数

物理的には、太鼓の倍音の周波数fnはディリクレ固有値λnに比例することが知られている。うなりの周波数は、2つの音の周波数差により決まる。具体的には、近い周波数f1f2 の音が同時に鳴ると、音波の干渉により音の振幅が周期的に変化する現象が発生する。これを「うなり」という。うなりの周波数fbeatは、次のように表される:
fbeat=|f1f2|.
f1,f2が近いほど、うなりはゆっくり発生し、周波数差が大きいと別々の音として認識される。

正三角形の太鼓

三角形の太鼓における第2、第3固有値λ2,λ3に注目しよう。正三角形の場合、
λ2=λ3=1129π2
が知られているMcCartin。このとき、うなりの周波数は
|λ3λ2|=0
であり、うなりは発生しない。したがってその意味で,正三角形の太鼓は「良い音」を鳴らすことが分かる。

わずかに歪んだ三角形の太鼓

では三角形をわずかに歪めるとどうなるか?実は歪んだ三角形においては|λ3λ2|>0となり、うなりが発生することが分かる。

実際に例えば、次の頂点を持つ三角形T(ε)を考える:
(0,0), (1,0), (1/2,3/2ε).
頂点(1/2,3/2ε)は、正三角形の頂点(1/2,3/2)からεだけずらされていることに注意しよう。ε=105のとき、Lehmann-Goerisch Methodを用いた固有値に対する精度保証付き数値計算で次の不等式が得られる:

122.8225λ2(T(ε))122.8231,  122.8232λ3(T(ε))122.8240.

このとき、
|λ3(T(ε))λ2(T(ε))|0.0001
であるため、うなりが発生することがわかる。

摂動を小さくするとどうなるか?

三角形Tの頂点の摂動εを限りなく0に近づけると、固有値λk (k=1,2,)は領域の摂動に対して連続に変化することから、ε0のときλ3,λ3は正三角形の重複固有値1129π2に収束する。したがってλ3λ20となり、うなりの周波数は0に収束する。

それでは任意の歪な三角形T(ε) (ε>0)に対して、うなりが存在することを証明できるだろうか?
具体的には次の予想を証明することはできるか?:

Conjecture 4.67 of "Shape Optimization and Spectral Theory" henrot

正三角形以外の三角形に対して、ディリクレ第2固有値λ2は単純である。

ここでλ2が単純であるとは、λ1<λ2<λ3が成り立つことを意味する。

小さな摂動に対する固有値の単純性の保証

Lehmann-Goerisch Methodを用いるような手法では、与えられたε>0についてλ2,λ3を分離できても、さらに小さいε>0に対して同様に分離できるかが常に問題となる。したがって、任意のε>0に対して固有値を分離することはできない。

プレプリントendoliuでは,「固有値の差分商」を精度保証付き数値計算で分離することによってλ2,λ3を分離した.例えば,具体的には以下の評価を得た:

(0,0),(1,0),(3/2,1/2ε)を頂点に持つ三角形をT(ε)とする.
任意のε(0,105]に対して,
λ2(T(ε))λ2(T(0))ε[89.5557,118.3308],
λ3(T(ε))λ3(T(0))ε[162.9557,196.4519].

正三角形T(0)に対してλ2(T(0))=λ3(T(0))であるから,この評価から任意のε(0,105]に対してλ2(T(ε))<λ3(T(ε))であることが分かる.

固有値の差分商公式

endoliuでは,重複固有値の方向微分公式に良く似た「固有値の差分商公式」を新しく導出し,固有値問題の精度保証付き計算の技術を用いて差分商公式を精度保証することにより,上記の評価を得た.

最後に差分商公式の主張について簡単に説明する.

  • パラメータ p=(x,y)(R2) に対して、頂点が (0,0),(1,0),(x,y) である三角形領域を Tp とする。固有値 λnp,,λNpp において重複している、つまり
    λnp==λNp(=:λ)
    であると仮定する。
    ここで、E を固有値 λ に対応する固有空間とする。

  • t>0,eR2 かつ e2=1 のとき、pt:=p+te(=:(x~,y~))p の摂動とする。記号を簡略化するために、Tt:=TptT0:=Tp、および λit:=λipt と表す。固有値 λnt,,λNt に対応する線形独立な固有関数をそれぞれ unt,,uNt とする。

  • また、St:R2R2を三角形領域 T0Tt に写す線形変換とする.

  • u~it:=uitSt(H01(T0))i=n,,N)として,E~t:=span{u~nt,,u~Nt}(H01(T0))
    と定める.

  • さらに、固有値の差分商を
    Dtλi:=(λitλ)/t(i=n,,N)
    と表し,次のように 2×2 行列 PtePe を定義する:Pte:=(St1StI)/t.

 

これらの設定のもと,以下の差分商公式が成り立つ:

固有値の差分商公式

dimE~t=dimE であると仮定する。E~t の基底を {ϕ~i}i=ni=NE の基底を {ϕi}i=ni=N とすると、次のように定義する:
Mt:=((Pteϕ~i,ϕj)T0),  Nt:=((ϕ~i,ϕj)T0)(i,j=n,,N).
このとき、差分商 Dtλi は次の行列固有値問題の (in+1) 番目の固有値となる:
Mtσ=μNtσ.
さらに、摂動された領域 Tt における固有関数 uit(i=n,,N) に対して、次が成り立つ:
u~it=uitSt=sniϕ~n++sNiϕ~N(E~t)(i=n,,N),
ここで、σi:=(sni,,sNi) は固有値Dtλiに対応する固有ベクトルとなる.

通常,近接固有値やその固有ベクトルの挙動を把握することは固有値問題の解析において非常に困難なことであるが,固有値の差分商に関する情報を用いることにより,近接固有値の挙動を追跡することができる.

参考文献

[1]
McCartin, Brian J., Eigenstructure of the equilateral triangle, Part I: The Dirichlet problem., Siam Review, 2003, 267-287
[2]
Henrot, Antoine., Shape optimization and spectral theory., De Gruyter Open, 2017
[3]
Ryoki Endo, Xuefeng Liu, Rigorous estimation for the difference quotients of multiple eigenvalues, arXiv:2305.14063, 2023
投稿日:20241027
更新日:20241027
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

「音の現象」を数学的に検証することに興味を持っています. 微分作用素の固有値問題に対する計算機援用証明について研究しています.

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中