スキップしてメイン コンテンツに移動

Mathematica 9.0買ったドー

データ解析をやりたい。

 音声データをマルチマイクで加工して音源推定を行いたいのだ。いずれiPad用のアプリを開発するにせよまず、週末の限られた時間であれこれは試せないので効率的なプログラム言語が必要だ。

 パソコンがまだキラキラしている頃、大学出たての頃、MacでMathematicaが使えることを知ってどうしても欲しくて大枚をはたいてMathematica 2.2を買った。16万円はしたと思う。MacOS8.Xの頃ではなかったか。手続き型の関数を勉強したものだ。

 それから会社で使ってきたものの個人的には購入はできないできた。値段が高いから。MacもOSがOSXになり、MacOS9互換のRossetaもサポート停止しオールドマック(cube)をこれ専用にメインテナンスしてるのみ。

 一年に数回数式処理がしたくなるもののMaximaを使ったりしてみた。

 でも、やっぱり物足りない。効率的でない。

 20年ぶりにMathematicaを買った。HomeEdition for Mac 9.0だ。日本の代理店では7万円以上するが、Wolframのストアでは67000円だった。この差はなんだ？これを10年使い切れば元もとれよう。

 使ってみるとこれが凄い進化。基本の構造は変わってみたいだ。

 それでは早速使ってみよう。マニュアルを読めばすむことなのだが、備忘録として書き留めておく。というより思い出し。

 （１）音声データの取り込み
  • Directory[]　　　フォルダの確認
  • SetDirectory["~/Desktop"]　　作業ディレクトリをデスクトップに設定
  • data=Import["020A_120422_0812.WAV"]　音声データの取り込み。成功、簡単です。アウトプットに再生ボタンとかついているのでこれで聞こえた。なんとこの中にソナグラム（スペクトログラム）が既に表示されている。Mathematicaではスペクトログラム(Spectrogram)という。リストを与えれば表示するようだ。WAV形式のファイルを左右のチャンネルに分解しよう。
（２）音声データの理解
  • Import[]関数は強力だが自動でなんでもやりすぎる。ちょっと音声データをお勉強。まずフォーマットが"Sound"というのがある。
  • ここにサウンドの表現形式があった。
  • Short[InputForm[data]]で形式を確認できた。出力は、次の通りだった。Out[22]/Short=Sound[SampledSoundList[{{0., 0.,    0., << 1613276 >>, -1.19209*10^-7}, {0., << 1613279 >>}}, 96000]]
  • フムフム。
  • data[[1,1]]はSampledSoundListの形式で左右のチャンネルの音声データが入っている。
  • data[[1,2]]にサンプリングレートが入っている。
  • ListPlay[rightch, SampleRate -> rate]とすれば右の音声が再生できた。

（３）ソノグラム（スペクトログラム）の表示
  • Spectrogram[rightch]で表示できました。


  • Spectrogramの２つ目のオプションnは"長さn == 2^(Round[Log[2,(Sqrt[m])]] + 1)のパーティションを使う"とある。ここで、mはデータ長。計算すると2048だった。この数字を大きくしたり小さくするとフーリエ変換の変換粒度を変えられる。また、窓関数も変えられる。Optionで
    • SampleRate -> rateとすると横軸は時間縦軸は周波数で表示された。下の図は256ポイントでサンプリング。













  • 縦軸はダウンコンバートして調整すればいいのかな。これは専用ソフト(Amadeus Pro)にはまだかなわない感じだ。いやPlotRangeを使えばいいんだ。PlotRange ->{{0, 16}, {0, 12000}}とか。


いやー簡単です。ここまでで２時間くらい。まだまだ楽しめそう。

参考にしたページ

ラベル:

コメント

コメントを投稿

このブログの人気の投稿

瞬間日記のデータをday oneに移動させるのだ

瞬間日記をiPodを購入した２０１０年から使ってきたが、day oneが人気なので使ってみた。なかなかよかった。コマンドライン(CLI:Command Line Interface)から入力できるし、クラウドで同期できるのもいい。見た目もきれい。安定しているのもいい。 という訳で、瞬間日記からday oneにデータを引っ越しすることにした。 （このページで半角の ＆ｇｔ； と表示されたら半角の”＞”記号だと思ってください。HTMLの仕様のインプリバグでしょう。） 特定の日記アプリ間のデータの移動とは人生で一回しかないだろうから汎用性など考えずに間違えないように慎重にすることが大切だ。確かめて確かめて。このページはそんな備忘録である。また、瞬間日記のデータをday oneに移動させたという記事も見つけられなかったので誰かのお役に立てるかもしれない。でも自己責任でね。 瞬間日記は独自フォーマット(.dat3)や標準的な(.sqlite)、(.csv)でもMacにバックアップできる。使い勝手のいいcsvファイルは本文だけで写真情報が欠落している。瞬間日記のデータを引き上げるのは色々試したが、テキストばかりでなく写真も移動させたいの結局メールでバックアップすることにした。瞬間日記側で自分にひと月毎のバックアップをメールする。３年半分。ソチオリンピック開会式を見ながら作業でもれなく終了。 これは瞬間日記が３０枚しか１つのメールに写真を貼付できないので、ひと月分だけせっせとメールする。 届いたメールはOSX側のmail.appで開くが、~/Libraryに保存されているelmsでなく、単にテキストファイルでmail.app側の機能を使って出力することにする。 瞬間日記のデータを取り出し、day oneにエントリできるくらいの粒度のファイルにするまでの流れ： iPodの瞬間日記のひと月分に区切ったバックアップデータを自分のメールアドレスに送信する OSXのmail.appで受信したメール全部（件名"MomentDialy"で始まるメールを一つのメールボックスに束ねる。例えばMomentDiaryというメールボックスに全部入れる mail.app上のメールボックスMomentDiary内の全部のメールを選択する。つまり⌘+a。 ...
コメントを投稿
続きを読む

電解コンデンサの容量抜けを測ってみた

中古オーディオ屋でスーパーウーハーを買ってから7〜8年経つだろうか。 1989年に発売されたONKYOの SL-10 という機種で、現役で商品が店頭に並んでいた時にピアノ曲に深みが増すことにちょっとした感動を覚えて中古屋で見つけた時2万円で即購入したのだ。 スピーカのエッジのクズがエンクロージャのダクトから出てきて開けてみたら見事エッジ全体が下手っていた。これを奇麗にして、ウレタンエッジに交換することにした。写真1、と写真2がそのBefore/Afterだ。 写真1　SL-10のスピーカのエッジを交換してみた （エッジを取り除いた状態） 写真2　SL-10のエッジを交換したスピーカ また、SL-10が勝手に自己発振する様になっていたのでコンデンサの容量が抜けたと想像してこれらを交換することにした。これは電解コンデンサの容量が抜ける事でフィードバック回路の時定数が変わってしまいネガティブフィードバックがポジティブフィードバックに位相が回って知ったのでは予想した。 電解コンデンサを交換したSL-10のコンデンサ 交換した電解コンデンサ達（容量は抜けていなかった）  マルツ電波で発注しておいたオーディオ用のコンデンサを付け替えた。ただ、それだけでは詰まらないので交換したコンデンサを測ってみた。使った測定器は卓上テスタ（GBW 9000A）の付属機能。最大20μFまでしか測れなかったが、結論からすると交換した電解コンデンサの容量抜けは無かった。 測定例 定格（μF） 　　実測（μF） 10　　　　　　　10.19 10　　　　　　　10.00 4.7　　　　　　　4.78 3.3　　　　　　　3.27 1.0　　　　　　　1.01 0.47　　　　　　　0.46 というわけで直しても自己発振は直らないかも。明日組立ててみよう。
2 件のコメント
続きを読む

環境音があってこその野鳥の声、Woodieの野鳥録音の仕方

春雨にぬれそぼるの紅葉（越前市花筐公園にて） 街から自動車のエンジン音がなくなる日が来ることを考えていた。 野鳥の録音を趣味にしている人が増えてきているようだ。PCMレコーダが安価になり、また蒲谷鶴彦さんを初め松田道夫さん、上田秀雄さんなど先人の努力の賜物だと思っている。 私も１９７０年代の中学生時代のラジカセで生録（蝉時雨、虫の音とか）を初体験し、８０年代後半に始めたバードウオッチングで、鳥の声から鳥を識別できたらどんなに素敵だろうと識別勉強の目的で１９９０年代入ってMDにガンマイクを接続し録音を開始した。これは光ディスクに書き込む時のメカ音がかなり大きくてそれが録音されてしまい、それを避けるのに難儀した記憶がある。２０００年入ってsony初代のPCMソリッドレコーダPCM-D1を購入したりした。 鳥の声も野鳥の写真と同じで色々な楽しみ方があっていいと思う。まだ私の考えはまとまっていないが、好き嫌いの傾向はわかってきたので備忘録的に書き留めておく。 図鑑型：写真を多用したる野鳥図鑑では一種毎にその主を代表する個体の写真であって、識別ポイントがしっかり写っていることが重要であろう。そのためどのような環境下で撮影されたのかは二の次であり、複数の種が同じ写真に納まっていることはまずない。野鳥大鑑（蒲谷鶴雄）、野鳥の声（上田秀雄）は図鑑そのものだ。鳥の種類毎にその鳥の声しか入っていないようにしてある。川のせせらぎ、他の種の鳥の声など極力排除しようとしているし、そう音声の編集をしていると思う。つまり加工臭がかなりする。 ネイチャーフォト型：図鑑型に対して自然そのものの美と捉えたネイチャーフォトは鳥の生きている環境を重視しているようだ。写真から物語を紡ぎだせそうな秀作も多い。鳥の録音では日本野鳥紀行（蒲谷鶴彦）がそれだろうか。 最初は図鑑型の録り方を手本として、どうすれば一種類だけの鳥の声を録音できるだろうと試行錯誤してきた。その答えがパラボラ集音器やガンマイクだった。それでも、声ならばその狙った種以外にもどうしても他の種の声や、町の音、飛行機の音が入ってしまう。プロはすごいものだ。狙った鳥だけ浮かび上がらせる録音の技術も加工・編集の技術もさすがだと思う。 しかし今は、環境音がしっかり入っていてもそれはそれでいいと思える。電車や自動車の音が近...
1 件のコメント
続きを読む