読者です 読者をやめる 読者になる 読者になる

shingoushori's dialy

音信号処理を専ら研究していた元博士後期課程の学生によるメモ

整数の逆数、その2進数表現

1から16までの整数の逆数を、2進数で表現してみた. <ソースコード> digits=16; max=2^(digits-1); for hoge=1:16 disp([num2str(hoge,'%02d') ' ' dec2bin(hoge,digits) ' ' dec2bin(round(max/hoge),digits)]); end <算出結果> 01 0000000000000001 100000…

Mac OS X 10.9.2 Mavericks で Octave が動かなくなった

Mac OS X 10.9.2 Mavericks で Octave が動かなくなったおそらく,10.9.2 にアップデートした時点で動かなくなったのだと思うついでに,MacPortsのselfupdateも通らない有様 ... いろいろなブログをはじめとした情報を参考に試したが,結局のところ,Octave…

MIDIファイルを使った自動作曲の試み,いきなりつまずく

Mathematicaで,MIDIファイルを使った自動作曲を試みた. しかし,Mathematica内でMIDIをwavとして出力する,効率のよい方法がなかった. とりあえず,正弦波で合成するものを書いたが,処理が悲劇的に遅い. たった1秒,4つの正弦波を合算するだけのものを…

MATLAB / Octaveで,IIRフィルタの極・零点の位置を算出 -> z平面上にプロット

z平面上にプロットするのであれば, zplane関数で簡単にできてしまいます.零点-極プロット - MATLAB zplane - MathWorks 日本zplane関数から,極・零点の位置を抜き出す方法がわからなかった.そこで,算出しプロットしようと思い立った次第です. 今回は,…

MATLAB / Octaveで,スペクトログラムやスペクトル平均を算出 Ver.2

Ver.1 から以下の点について改良した. <改良点> データの先頭と終端まで,フレームの重みを均一化. スペクトルとスペクトル平均のプロットに,縦横軸の表記を追加. MATLABのエディタで開いたとき,階層化されるようコメントを整理. <ソースコード> 諸元 …

MATLAB / Octaveで,epsファイルをカラーで書き出す

画像ファイルの書き出しには,print関数かsaveas関数を用いることができる.両関数とも,MATLABとOctaveの両方で用いれる.ところが,print関数でepsファイルをカラーで書き出すとき,それもMATLABとOctaveを行き来するようなときは,以下の順で引数を入力す…

MATLAB / Octaveで,スペクトログラムやスペクトル平均を算出 Ver.1

(2013.08.20 追記) 改良版 Ver.2をつくりました -> Ver.2 spectrogram関数を用いてしまえば,一発でプロットできてしまいます. 短時間フーリエ変換を使用したスペクトログラム - MATLAB spectrogram しかし,窓関数やゼロ詰めといった諸元を調整するには,…

MATLAB / Octaveで,二次元ヒストグラムを作図するとき,グリッド線を消す

タイプグリッドありグリッドなし surf pcolor <ソースコード> % set sample table table=((1:16)-8)'*((1:16)-8); % surf version figure;surf(table); saveas(gcf,'surf_edge_on.jpg'); figure;surf(table,'EdgeColor','none'); saveas(gcf,'surf_edge_off.…

テスト信号の準備

<概要> 動作チェックに欲しいwavファイルを,MATLAB/Octaveで制作.今回は,エネルギーが一定のものを用意.dB値の基準を,複数のプログラム間で揃える必要があるため. <制作するデータ> constant minimum bit : 全サンプルが,正の最小値 constant full bi…