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 が動かなくなったおそらく，10.9.2 にアップデートした時点で動かなくなったのだと思うついでに，MacPortsのselfupdateも通らない有様 ... いろいろなブログをはじめとした情報を参考に試したが，結局のところ，Octave…

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

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

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

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

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

タイプグリッドありグリッドなし 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…