はてなブックマーク

この記事はPython Advent Calendarの17日目の記事です。去年の記事で書かなかったPythonを使った音響信号処理関連の内容について書きます。 目次 サンプリングレートを変換したい waveモジュールで24bitの音声ファイルを読みたい 群遅延をscipy.signal.group_delayで推定したい フレーム処理をnumpyのstrideトリックを使って行いたい(スペクトログラムを描画したい) TSP信号を生成したい 零点、極を可視化したい サンプリングレートを変換したい サンプリングレートを変換するためには、元のサンプリングレートと変換したいサンプリングレートが整数比であればアップ/ダウンサンプリングどちらか一方だけで済みますが、有理数比の場合はそれぞれのサンプリングレートの最小公倍数にアップサンプリングしてから、求めたいサンプリングレートにダウンサンプリングする

※2016/2/7 以下を追記 ・Dockerfileを画像を保存できるように修正 ・起動スクリプトを追記 ・Octave4.0でsubmit関数実行時にエラーが出ることに対する対策を追記 背景 CouseraのMachine Learningのコースの課題をこなすためにOctave開発環境が必要になるのですが、公式wikiによるとEl CapitanではOctaveはもはやvagrantを使って仮想マシンで動かすことが推奨されているようです。しかしOctaveを使うためだけに仮想マシンを一つ作るのも憚られるものがあるので、dockerで使えないかを試してみました。 以下は2016/2/1時点の手順です。 手順 以下のような手順で構築可能でした。4, 8あたりで引っかかりましたがなんとかこれで起動するようになりました。 Docker Toolboxをインストール docker-machin

Python その2 Advent Calendar 2015の13日目の記事です。 普段の仕事として主に音響信号処理のアルゴリズム開発やDSP実装などを行っているのですが、アルゴリズムを構築する際は最初にPythonを使ってアルゴリズムの検討を行い、その結果を踏まえてCやC++でPC上で動作するリアルタイムモデルを作成するという流れで行っています。音響信号処理に関するまとまった日本語記事は意外に存在しないため、ここではPythonで音響信号処理を行う方法およびこれまでに仕事等で使って来たノウハウなどを適当にまとめようと思います。 目次 Pythonで音響信号処理をするモチベーション オーディオファイルの読み書き リアルタイムにオーディオ処理を行いたい 周波数応答を表示したい デジタルフィルタを設計したい 零点、極と係数配列b, aを変換したい デジタルフィルタを時系列信号に適用したい 群