デジタル情報処理3
note Item Type Metadata
note
ディジタル信号処理 - サンプリングと音の数学的表現
概要: この資料では、ディジタル信号処理におけるサンプリングの基礎から、音の数学的表現までを体系的に学習できます。連続時間信号を離散時間信号に変換する過程と、音の要素（振幅・周波数・音色）の数学的記述について、具体例を用いて詳しく解説します。
目次
1. サンプリングの基礎
1.1 基本概念
ディジタル信号処理では、連続時間信号を離散時間信号に変換するサンプリング処理が重要です。
連続時間信号
x_a(t) t: 時刻[sec]
サンプリング
t = nT (nは整数)
離散時間信号
x[n] = x_a(nT), n = -∞～＋∞
記号の説明:
- n: 正規化された時刻のため無次元（1番目、2番目、3番目のような）
- T: サンプリング周期 [sec]
- F_s: サンプリング周波数 [Hz]
1.2 重要な関係式
$$T = \frac{1}{F_s}$$
この関係式により、サンプリング周波数が決まれば、サンプリング周期も自動的に決まります。
2. サンプリングの例
離散時間信号と連続信号の関係を明確にするためには、標本化周期Tの情報が必要です。
2.1 例1: 1Hzの正弦波を4Hzでサンプリング
連続時間信号
$$x_a(t) = \sin(2\pi \cdot 1 \cdot t)$$
サンプリング条件
- サンプリング周波数: F_s = 4 Hz
- サンプリング周期: T = 1/F_s = 1/4 = 0.25 sec
離散時間信号
\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}
ポイント: 4Hzでサンプリングするため、1秒間に4サンプルが得られます。
2.2 例2: ラ音（440Hz）のサンプリング
連続時間信号
$$x_a(t) = \sin(2\pi \cdot 440 \cdot t)$$
サンプリング条件
- サンプリング周波数: F_s = 8 kHz = 8000 Hz
- サンプリング周期: T = 1/8000 sec
離散時間信号
\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}
3. 音の要素
音の要素についてまとめると、主に以下の3つがあります。
3.1 大きさ（振幅）
音の大きさは振幅によって決まります。
例:
- x_a(t) = 0.05\sin(2\pi \cdot 440 \cdot t) （小さい音）
- x_a(t) = 3\sin(2\pi \cdot 440 \cdot t) （大きい音）
振幅が大きいほど音も大きくなります。
3.2 高さ（周波数）
音の高さは周波数によって決まります。
例:
- x_a(t) = \sin(2\pi \cdot 440 \cdot t) … ラの音（440Hz）
- x_a(t) = \sin(2\pi \cdot 880 \cdot t) … 1オクターブ高いラの音（880Hz）
音楽では、1オクターブ高い音は元の周波数の2倍にするという慣例があります。
3.3 基本周波数（ピッチ）
平均律について
- 音の高さ（ピッチ）は基本周波数で決まります
- 平均律では、基準音（ラ、A4）は440Hzです
- 1オクターブは12等分されており、n半音上の周波数fは次式で表されます：
$$f = 440 \times 2^{n/12}$$
例:
- n=0のとき: 440Hz
- n=12のとき: 880Hz（1オクターブ上）
4. 音色の数学的表現
4.1 音色（ティンバー）の基礎
- 音色は波形の形状によって決まります
- これは数式で表すことができます
- 一般的に、周期的な音の波形はフーリエ級数展開により、基本周波数とその整数倍の周波数（倍音）の和として表せます
4.2 音色の一般的な数式表現（フーリエ級数）
音色は基本周波数と倍音の組み合わせで決まります：
$$x_a(t) = A_1\sin(2\pi ft) + A_2\sin(2\pi \cdot 2ft) + A_3\sin(2\pi \cdot 3ft) + \cdots$$
4.3 音色の例
1. 純音（基本周波数のみ）
$$x_a(t) = \sin(2\pi \cdot 440 \cdot t)$$
- 最もシンプルな音色
- 倍音を含まない純粋な正弦波
2. 矩形波（奇数倍音のみ）
$$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]$$
- 奇数倍音（3倍音、5倍音、7倍音...）のみを含む
- デジタル音源でよく使われる音色
3. 三角波（奇数倍音、振幅が1/n²で減衰）
$$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]$$
- 矩形波より柔らかい音色
- 高次倍音の影響が小さい
4. のこぎり波（全倍音を含む）
$$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]$$
- 全ての倍音を含む
- 明るく豊かな音色
5. 楽器の音色例（バイオリン風）
$$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)$$
- 基本周波数 + 2倍音 + 3倍音 + 4倍音
- 各倍音の振幅比が音色を決定
4.4 音色の特徴
| 基本周波数 | 音の高さを決定 |
| 倍音の種類 | どの周波数成分を含むか |
| 倍音の振幅比 | 各倍音の強さのバランス |
| 包絡線 | 時間的な音量変化（ADSR） |
5. まとめ
5.1 重要なポイント
サンプリング
- 連続時間信号 x_a(t) を離散時間信号 x[n] に変換
- サンプリング周期 T = 1/F_s が重要
- 横軸が秒の場合、t = nT の関係
音の要素
- 振幅: 音の大きさ
- 周波数: 音の高さ
- 音色: 波形の形状（倍音の組み合わせ）
音色の数学的表現
- フーリエ級数展開による表現
- 基本周波数と倍音の組み合わせ
- 各倍音の振幅比が音色を決定
5.2 記号の統一
| x_a(t) | 連続時間信号 | アナログ信号 |
| x[n] | 離散時間信号 | デジタル信号 |
| F_s | サンプリング周波数 | サンプリング条件 |
| T | サンプリング周期 | サンプリング条件 |
5.3 学習のポイント
- サンプリング定理の理解
- フーリエ級数展開の応用
- 音の物理的性質と数学的表現の関係
- 実用的な計算例の習得
6. 参考文献・関連リンク
参考文献
- Oppenheim, A. V., & Schafer, R. W. (2010). Discrete-Time Signal Processing (3rd ed.). Prentice Hall.
- Proakis, J. G., & Manolakis, D. G. (2007). Digital Signal Processing (4th ed.). Prentice Hall.
- 小野寺, 信. (2015). ディジタル信号処理の基礎. コロナ社.
- 太田正哉. (2018). 例解 ディジタル信号処理入門. コロナ社.
コメント