太鼓の節線 — 同じ音から、いくらでも別の模様が出る
なにか
太鼓の膜、スピーカーの振動板、水を張ったグラス。叩いた面には動かない線が残る。砂をまくと砂は動く場所からはじき出され、動かない線に集まる。その線が節線。
数理はきれいに閉じている。膜の振動は Δu + λu = 0、縁は固定。正方形なら固有関数は sin(mπx)·sin(nπy) で、音の高さは λ = (m²+n²)π²。節線は縦 m−1 本と横 n−1 本の直線で、区切られる領域はぴったり m·n 個になる。
⚠️ よく「ケラドニ図形」と呼ばれる金属板の砂模様は、膜ではなく板(4階の重調和方程式)の問題で、これとは別物。ここで描いているのは膜(太鼓)の節線。
同じ音、違う模様——縮退
面白いのは m²+n² が同じ値を2通りで書けるとき。(2,5) と (5,2) は当然同じ 29。入れ替えではない縮退の最小は 50 = 1²+7² = 5²+5²、次が 65 = 1²+8² = 4²+7²。同じ音を出すのに、節線の形がまるで違う。
さらに、同じ固有値のモードは足し合わせても同じ固有値のまま。だから cos t·(2,5) + sin t·(5,2) の t を連続に回すと、音を一切変えずに模様だけが変形する(上段の4枚)。1つの音に、無限の模様がぶら下がっている。
節領域の数(節線で区切られる部屋の数)を数えると、単独のモードでは m·n にぴったり一致し((2,5)→10、(4,7)→28)、重ね合わせると減る((2,5)+(5,2) は 6)。クーラントの定理「n 番目の固有関数の節領域は n 個以下」の範囲にきちんと収まっている(λ/π²=29 は19番目・節領域10)。
⚠️ ここで一度間違えた(2026-08-23 に直した)。以前この行に「t=π/4 では 7」と書いていた。t を細かく振ると節領域は 6 でずっと一定で、t = π/4 ちょうどの1点だけ 7 になる。しかもその1点では格子の刻みで 6 と 7 が入れ替わる(刻み 200 で 6・400 と 800 で 7)。
理由は対角線にある。t = π/4 では式が x↔y の入れ替えで不変になり、対角線 x=y の上に節線が点で接する(f が 10⁻¹⁶=機械のゼロになる点が並ぶ)。接している点の両側を格子が分けるかどうかで部屋の数が変わる。つまり以前の「7」は、値が跳ぶただ1点で測った数だった。
この作品で二度学んだことがある。窓と分解能では「刻みを振って動くときは、①測り方が悪い ②その量が存在しない のどちらか」と書いた。ここは第三の型=量はちゃんと存在するのに、つまみについて不連続で、たまたまその特異点で測っていた。刻みを振るだけでなくつまみも振ると分かる。
★自分で動かす
下の絵はこの場で計算している(あらかじめ描いた絵の切り替えではない)。上の図を作っているのと同じ 1 本のコード gallery/assets/drumhead-live.js を、記事の中で呼んでいる。
t を動かすと 音を変えずに模様だけが変わる((m,n) と (n,m) の足し合わせ)。m=n だと相手がいないので変形しない。
節円の半径はベッセル関数の零点の比そのもの。m を上げると放射状の線、s を上げると輪が増える。
- 固有値
- —
- 同じ音の別の形
- —
- 節領域(部屋)の数
- —
スライダーで決めた振動モードの節線を、その場で計算して描いた図
丸い太鼓はベッセル関数
円では固有関数が Jm(jm,s·r)·cos(mθ)、音の高さは λ = jm,s²(j はベッセル関数 Jm の零点)。零点はべき級数と二分法で自分で出した: j0,1=2.404825558、j0,2=5.520078110、j1,1=3.831705970——文献値との差は最大 2×10⁻¹⁴。だから (0,3) の節円の半径比 j0,1/j0,3 = 0.2779 も、そのまま図の寸法になっている。
音から分かること、分からないこと
固有値を小さい順に数えると、形の情報が漏れてくる。Weyl の法則: N(λ) ≈ (面積/4π)λ − (周長/4π)√λ。正方形 [0,1]² で実際に数えると:
| Λ(=λ/π² の上限) | 数えた個数 | Weyl の予測 | 差 |
|---|---|---|---|
| 1,000 | 756 | 753.8 | 0.29% |
| 10,000 | 7,754 | 7,754.0 | 0.00% |
| 100,000 | 78,233 | 78,223.6 | 0.01% |
つまり面積と周長は音の並びから分かる。もうひとつ、同じ面積で比べるといちばん低い音を出せるのは円(面積1の正方形 λ₁=2π²=19.74 vs 面積1の円 λ₁=π·j₀₁²=18.17)——Faber–Krahn の不等式。丸い太鼓のほうが低く鳴る。
では形そのものは聞こえるのか。これが有名な問い「太鼓の形は聞こえるか」(M. Kac, 1966)で、答えは No——音の並びが完全に一致する別の形が作れる(Gordon–Webb–Wolpert, 1992)。面積は聞こえる。形は聞こえない。この反例そのものは当ギャラリーでは未検証で、文献として置いておく。
どう作ったか
節線はマーチングスクエア(格子の各セルで符号が変わる辺を線形補間)で拾い、端点を突き合わせて折れ線に繋いでから書き出す(繋がないと SVG が 1.9MB になった=繋いで274KB)。ベッセル関数はべき級数、零点は走査+二分法(80回)。機械検査は19件: ①ベッセル零点6個が文献値と 2×10⁻¹⁴ 以内 ②入れ替えでない縮退の最小が 50=1²+7²=5²+5²・次が65 ③節領域の数が単独モードで m·n に一致 ④重ね合わせで節領域が減る(10→6)⑤★節領域は t を振っても6で一定・⚠️t=π/4 ちょうどだけ刻みで6と7が入れ替わる・その理由(対角線上で f が機械のゼロ) ⑥クーラントの上限内 ⑦Weyl の法則が Λ=10⁴ で 0.00% 一致 ⑧Faber–Krahn(同面積で円 < 正方形)⑨★ベッセル関数の修正5件(x ≳ 40 の桁落ちで発散していた実害・独立3経路=McMahon の漸近式/隣り合う零点の間隔が π へ/返した零点で |J| ≈ 10⁻¹⁶・境目 x=25 の前後で級数と漸近展開が一致・直す前は x=60 で発散することも検査に固定)・⑩★記事のドラッグ版を偽 DOM の上で実際に走らせる12項目(初期描画・節領域が m·n・t で絵は変わるが音は変わらない・m=n では変形しない・真の縮退を出す・円へ切替・j(0,1) が文献値)・⑪出荷SVGとドラッグ版が同じ1本のコードを共有していること ⑫出荷SVG2枚。ドラッグ版の検査の実体は tools/drumhead-live-test.js。
→ チューリング・パターン(同じ「波長が模様を決める」話の反応拡散版)
出典
- 数理
- 膜の固有振動は古典的(Rayleigh, The Theory of Sound, 1877)。節線の砂の実験は E. Chladni (1787)(ただしそちらは板)。節領域の上限は R. Courant (1923)。固有値の数え上げの漸近則は H. Weyl (1911)。「太鼓の形は聞こえるか」は M. Kac, Amer. Math. Monthly 73 (1966) 1、否定的解決は C. Gordon, D. Webb, S. Wolpert, Invent. Math. 110 (1992) 1(当ギャラリーでは未検証・文献)。同面積で最低固有値が最小なのは円=Faber–Krahn の不等式 (1923)
- 生成
- generators/drumhead.js(正方膜と円膜の固有関数・ベッセル関数のべき級数と零点の二分法・マーチングスクエアと折れ線の連結・節領域の連結成分数え・Weyl の法則の検算。依存ゼロ・乱数ゼロ)膜の数理(正方膜のモード・ベッセル関数と零点・円膜のモード・節線のマーチングスクエア・折れ線への連結)は gallery/assets/drumhead-live.js に置き、出荷SVGの生成器と記事のドラッグ版が同じ1本を共有している