wave_form.nb

音色の周波数解析
Mathematica による Wave File の解析

2003/2/26 平野拓一（東京工業大学）
2003/2/26 Takuichi Hirano (Tokyo Institute of Technology)

Wave ファイルの読み込み

In[1]:=

Out[1]=

$c:\\hira2\\public_html\\hobby\\edu\\sonic_wave\\sound_analysis\\sound_instruments$

In[2]:=

In[3]:=

音を出して確認

In[4]:=

$ListPlay[WaveFormPiano, SampleRate -> 22050, PlayRange -> {-2^16, 2^16}]$

[Graphics:HTMLFiles/index_13.gif]

Out[4]=

時間波形を描く

In[5]:=

$ListPlot[WaveFormPiano, PlotRange -> {-2^16, 2^16}, PlotJoined -> True, PlotStyle -&g ... gt; None, FrameLabel -> {"Time", "Amplitude", "", ""}]$

[Graphics:HTMLFiles/index_16.gif]

Out[5]=

１周期を切り取る

[入力パラメータ]
smprate: wave ファイルのサンプリングレート
freq: 解析する音の周波数（うまく１周期を切り出すために指定）
n1: 切り出しを開始するリスト要素の開始番号

In[6]:=

$smprate = 22050 ; freq = (27.5 * 2^3) * (2^1/12)^3 ; t = 1/freq ; n1 = 20000 ; n = Floor[(t/1) ... t; None, FrameLabel -> {"Time", "Amplitude", "", ""}]$

[Graphics:HTMLFiles/index_19.gif]

Out[13]=

離散フーリエ変換してスペクトルを描く

In[14]:=

$General :: spell1 : スペル間違いの可能性があります．新規シンボル\" FWaveFormPiano2 \"はすでにあるシンボル\" WaveFormPiano2 \"に似ています．$

[Graphics:HTMLFiles/index_23.gif]

[Graphics:HTMLFiles/index_24.gif]

切り出した１周期をいくつも並べて周期波形を作り、音を聴いてみる

In[21]:=

$freq = 440 * (2^1/12)^3 ; t = 1/freq * 1 ; dur = 1.0 ; n = Ceiling[dur/t] ; WaveFormPiano3 = F ... , {i, 1, n}]] ; ListPlay[WaveFormPiano3, SampleRate -> 22050, PlayRange -> {-2^16, 2^16}] ;$

[Graphics:HTMLFiles/index_26.gif]

Converted by Mathematica (February 28, 2003)

音色の周波数解析Mathematica による Wave File の解析

Wave ファイルの読み込み

時間波形を描く

１周期を切り取る

離散フーリエ変換してスペクトルを描く

切り出した１周期をいくつも並べて周期波形を作り、音を聴いてみる

音色の周波数解析
Mathematica による Wave File の解析