Dublin Core The Dublin Core metadata element set is common to all Omeka records, including items, files, and collections. For more information see, http://dublincore.org/documents/dces/. Title A name given to the resource note Description An account of the resource 日常のノート note 日常のノートの記録 note テキストの自由記述 <h1>ディジタル信号処理 - サンプリングと音の数学的表現</h1> <p><strong>概要</strong>: この資料では、ディジタル信号処理におけるサンプリングの基礎から、音の数学的表現までを体系的に学習できます。連続時間信号を離散時間信号に変換する過程と、音の要素（振幅・周波数・音色）の数学的記述について、具体例を用いて詳しく解説します。</p> <h2>目次</h2> <ul> <li><a href="#section1">1. サンプリングの基礎</a></li> <li><a href="#section2">2. サンプリングの例</a></li> <li><a href="#section3">3. 音の要素</a></li> <li><a href="#section4">4. 音色の数学的表現</a></li> <li><a href="#section5">5. まとめ</a></li> <li><a href="#section6">6. 参考文献・関連リンク</a></li> </ul> <h2>1. サンプリングの基礎</h2> <h3>1.1 基本概念</h3> <p>ディジタル信号処理では、<strong>連続時間信号</strong>を<strong>離散時間信号</strong>に変換する<strong>サンプリング処理</strong>が重要です。</p> <h4>連続時間信号</h4> <p>x_a(t) t: 時刻[sec]</p> <h4>サンプリング</h4> <p>t = nT (nは整数)</p> <h4>離散時間信号</h4> <p>x[n] = x_a(nT), n = -∞～＋∞</p> <p><strong>記号の説明:</strong></p> <ul> <li>n: 正規化された時刻のため無次元（1番目、2番目、3番目のような）</li> <li>T: サンプリング周期 [sec]</li> <li>F_s: サンプリング周波数 [Hz]</li> </ul> <h3>1.2 重要な関係式</h3> <p>$$T = \frac{1}{F_s}$$</p> <p>この関係式により、サンプリング周波数が決まれば、サンプリング周期も自動的に決まります。</p> <h2>2. サンプリングの例</h2> <p>離散時間信号と連続信号の関係を明確にするためには、<strong>標本化周期Tの情報</strong>が必要です。</p> <h3>2.1 例1: 1Hzの正弦波を4Hzでサンプリング</h3> <h4>連続時間信号</h4> <p>$$x_a(t) = \sin(2\pi \cdot 1 \cdot t)$$</p> <h4>サンプリング条件</h4> <ul> <li><strong>サンプリング周波数</strong>: F_s = 4 Hz</li> <li><strong>サンプリング周期</strong>: T = 1/F_s = 1/4 = 0.25 sec</li> </ul> <h4>離散時間信号</h4> <p>\begin{align} x[n] &= x_a(nT) \\ &= x_a(n \cdot 0.25) \\ &= \sin(2\pi \cdot 1 \cdot n \cdot 0.25) \\ &= \sin(0.5\pi n) \end{align}</p> <p><strong>ポイント</strong>: 4Hzでサンプリングするため、1秒間に4サンプルが得られます。</p> <h3>2.2 例2: ラ音（440Hz）のサンプリング</h3> <h4>連続時間信号</h4> <p>$$x_a(t) = \sin(2\pi \cdot 440 \cdot t)$$</p> <h4>サンプリング条件</h4> <ul> <li><strong>サンプリング周波数</strong>: F_s = 8 kHz = 8000 Hz</li> <li><strong>サンプリング周期</strong>: T = 1/8000 sec</li> </ul> <h4>離散時間信号</h4> <p>\begin{align} x[n] &= x_a(nT) \\ &= x_a(n/8000) \\ &= \sin(2\pi \cdot 440 \cdot n/8000) \\ &= \sin(0.11\pi n) \end{align}</p> <h2>3. 音の要素</h2> <p>音の要素についてまとめると、主に以下の3つがあります。</p> <h3>3.1 大きさ（振幅）</h3> <p>音の大きさは<strong>振幅</strong>によって決まります。</p> <p><strong>例:</strong></p> <ul> <li>x_a(t) = 0.05\sin(2\pi \cdot 440 \cdot t) （小さい音）</li> <li>x_a(t) = 3\sin(2\pi \cdot 440 \cdot t) （大きい音）</li> </ul> <p>振幅が大きいほど音も大きくなります。</p> <h3>3.2 高さ（周波数）</h3> <p>音の高さは<strong>周波数</strong>によって決まります。</p> <p><strong>例:</strong></p> <ul> <li>x_a(t) = \sin(2\pi \cdot 440 \cdot t) … ラの音（440Hz）</li> <li>x_a(t) = \sin(2\pi \cdot 880 \cdot t) … 1オクターブ高いラの音（880Hz）</li> </ul> <p>音楽では、1オクターブ高い音は元の周波数の2倍にするという慣例があります。</p> <h3>3.3 基本周波数（ピッチ）</h3> <h4>平均律について</h4> <ul> <li>音の高さ（ピッチ）は基本周波数で決まります</li> <li>平均律では、基準音（ラ、A4）は440Hzです</li> <li>1オクターブは12等分されており、n半音上の周波数fは次式で表されます：</li> </ul> <p>$$f = 440 \times 2^{n/12}$$</p> <p><strong>例:</strong></p> <ul> <li>n=0のとき: 440Hz</li> <li>n=12のとき: 880Hz（1オクターブ上）</li> </ul> <h2>4. 音色の数学的表現</h2> <h3>4.1 音色（ティンバー）の基礎</h3> <ul> <li>音色は<strong>波形の形状</strong>によって決まります</li> <li>これは数式で表すことができます</li> <li>一般的に、周期的な音の波形は<strong>フーリエ級数展開</strong>により、基本周波数とその整数倍の周波数（倍音）の和として表せます</li> </ul> <h3>4.2 音色の一般的な数式表現（フーリエ級数）</h3> <p>音色は基本周波数と倍音の組み合わせで決まります：</p> <p>$$x_a(t) = A_1\sin(2\pi ft) + A_2\sin(2\pi \cdot 2ft) + A_3\sin(2\pi \cdot 3ft) + \cdots$$</p> <h3>4.3 音色の例</h3> <h4>1. 純音（基本周波数のみ）</h4> <p>$$x_a(t) = \sin(2\pi \cdot 440 \cdot t)$$</p> <ul> <li>最もシンプルな音色</li> <li>倍音を含まない純粋な正弦波</li> </ul> <h4>2. 矩形波（奇数倍音のみ）</h4> <p>$$x_a(t) = \frac{4}{\pi}\left[\sin(2\pi \cdot 440 \cdot t) + \frac{1}{3}\sin(2\pi \cdot 1320 \cdot t) + \frac{1}{5}\sin(2\pi \cdot 2200 \cdot t) + \cdots\right]$$</p> <ul> <li>奇数倍音（3倍音、5倍音、7倍音...）のみを含む</li> <li>デジタル音源でよく使われる音色</li> </ul> <h4>3. 三角波（奇数倍音、振幅が1/n²で減衰）</h4> <p>$$x_a(t) = \frac{8}{\pi^2}\left[\sin(2\pi \cdot 440 \cdot t) - \frac{1}{9}\sin(2\pi \cdot 1320 \cdot t) + \frac{1}{25}\sin(2\pi \cdot 2200 \cdot t) - \cdots\right]$$</p> <ul> <li>矩形波より柔らかい音色</li> <li>高次倍音の影響が小さい</li> </ul> <h4>4. のこぎり波（全倍音を含む）</h4> <p>$$x_a(t) = \frac{2}{\pi}\left[\sin(2\pi \cdot 440 \cdot t) - \frac{1}{2}\sin(2\pi \cdot 880 \cdot t) + \frac{1}{3}\sin(2\pi \cdot 1320 \cdot t) - \cdots\right]$$</p> <ul> <li>全ての倍音を含む</li> <li>明るく豊かな音色</li> </ul> <h4>5. 楽器の音色例（バイオリン風）</h4> <p>$$x_a(t) = \sin(2\pi \cdot 440 \cdot t) + 0.5\sin(2\pi \cdot 880 \cdot t) + 0.3\sin(2\pi \cdot 1320 \cdot t) + 0.1\sin(2\pi \cdot 1760 \cdot t)$$</p> <ul> <li>基本周波数 + 2倍音 + 3倍音 + 4倍音</li> <li>各倍音の振幅比が音色を決定</li> </ul> <h3>4.4 音色の特徴</h3> <table border="1"> <tbody> <tr> <td><strong>基本周波数</strong></td> <td>音の高さを決定</td> </tr> <tr> <td><strong>倍音の種類</strong></td> <td>どの周波数成分を含むか</td> </tr> <tr> <td><strong>倍音の振幅比</strong></td> <td>各倍音の強さのバランス</td> </tr> <tr> <td><strong>包絡線</strong></td> <td>時間的な音量変化（ADSR）</td> </tr> </tbody> </table> <h2>5. まとめ</h2> <h3>5.1 重要なポイント</h3> <h4>サンプリング</h4> <ul> <li>連続時間信号 x_a(t) を離散時間信号 x[n] に変換</li> <li>サンプリング周期 T = 1/F_s が重要</li> <li>横軸が秒の場合、t = nT の関係</li> </ul> <h4>音の要素</h4> <ul> <li><strong>振幅</strong>: 音の大きさ</li> <li><strong>周波数</strong>: 音の高さ</li> <li><strong>音色</strong>: 波形の形状（倍音の組み合わせ）</li> </ul> <h4>音色の数学的表現</h4> <ul> <li>フーリエ級数展開による表現</li> <li>基本周波数と倍音の組み合わせ</li> <li>各倍音の振幅比が音色を決定</li> </ul> <h3>5.2 記号の統一</h3> <table border="1"> <tbody> <tr> <td>x_a(t)</td> <td>連続時間信号</td> <td>アナログ信号</td> </tr> <tr> <td>x[n]</td> <td>離散時間信号</td> <td>デジタル信号</td> </tr> <tr> <td>F_s</td> <td>サンプリング周波数</td> <td>サンプリング条件</td> </tr> <tr> <td>T</td> <td>サンプリング周期</td> <td>サンプリング条件</td> </tr> </tbody> </table> <h3>5.3 学習のポイント</h3> <ol> <li><strong>サンプリング定理</strong>の理解</li> <li><strong>フーリエ級数展開</strong>の応用</li> <li><strong>音の物理的性質</strong>と数学的表現の関係</li> <li><strong>実用的な計算例</strong>の習得</li> </ol> <h2>6. 参考文献・関連リンク</h2> <h3>参考文献</h3> <ul> <li>Oppenheim, A. V., & Schafer, R. W. (2010). <em>Discrete-Time Signal Processing</em> (3rd ed.). Prentice Hall.</li> <li>Proakis, J. G., & Manolakis, D. G. (2007). <em>Digital Signal Processing</em> (4th ed.). Prentice Hall.</li> <li>小野寺, 信. (2015). <em>ディジタル信号処理の基礎</em>. コロナ社.</li> <li>太田正哉. (2018). <em>例解ディジタル信号処理入門</em>. コロナ社.</li> </ul> <h3>関連リンク</h3> <ul> <li><a href="https://ja.wikipedia.org/wiki/%E3%83%95%E3%83%BC%E3%83%AA%E3%82%A8%E5%A4%89%E6%8F%9B" target="_blank" rel="noreferrer noopener">フーリエ変換の基礎 - Wikipedia</a></li> <li><a href="https://ja.wikipedia.org/wiki/%E3%82%B5%E3%83%B3%E3%83%97%E3%83%AA%E3%83%B3%E3%82%B0%E5%AE%9A%E7%90%86" target="_blank" rel="noreferrer noopener">サンプリング定理 - Wikipedia</a></li> <li><a href="https://www.asj.or.jp/" target="_blank" rel="noreferrer noopener">音響学の基礎 - 日本音響学会</a></li> </ul> <h3>実習用ツール</h3> <ul> <li><a href="https://www.mathworks.com/products/signal.html" target="_blank" rel="noreferrer noopener">MATLAB Signal Processing Toolbox</a></li> <li><a href="https://docs.scipy.org/doc/scipy/reference/signal.html" target="_blank" rel="noreferrer noopener">Python SciPy Signal Processing</a></li> <li><a href="https://octave.sourceforge.io/signal/" target="_blank" rel="noreferrer noopener">Octave Signal Processing</a></li> <li><a href="https://numpy.org/doc/stable/reference/routines.fft.html" target="_blank" rel="noreferrer noopener">Python NumPy FFT</a></li> </ul> Dublin Core The Dublin Core metadata element set is common to all Omeka records, including items, files, and collections. For more information see, http://dublincore.org/documents/dces/. Title A name given to the resource デジタル情報処理3 Creator An entity primarily responsible for making the resource unjuno Date Created Date of creation of the resource. 2025-10-06 Rights Information about rights held in and over the resource CC BY-NC 4.0