Masassiah Blog

読書,商品紹介,大相撲,DIY など

MATLAB

文字抽出

MATLABを使って文字抽出を行ってみました。文字抽出というより、文字とその他の部分で画像を2値化するプログラムです。 以前のボリュームデータの可視化で使ったisosurfaceを使って実現しました。isosurfaceは3次元でしか使えない(?)ので画像データ(もちろん…

音響FDTD2次元

以前に書いた音響FDTD1次元バージョンを拡張して、同じくMATLABで音響FDTD(Finite Difference Time Domain)2次元バージョンのプログラムを作ってみた。 音圧と粒子速度に関する運動方程式と連続の式を、空間、時間領域で差分方程式に展開して、逐次計算して…

CIP法

はじめに CIP法のプログラムをMATLABで作ってみました.CIP法は高精度数値解法で,移流方程式を解くのに優れていると言われております.とりあえず1次元のプログラムです.方向分離解法を使えば,1次元のCIP法だけで多次元化することもできます(M型CIP法やC…

音響FDTD1次元バージョン

MATLABで音響FDTDプログラムを作成しました.とりあえず1次元バージョンです. *自由音場における音圧p,粒子速度uに関する偏微分方程式を差分化して,逐一音圧と粒子速度を交互に計算していきます. *密度と体積弾性率は室温20度の値を用いています. *Mur…

HRIRをwavファイルに畳み込み

HRIR(Head-Related Impulse Response,頭部インパルス応答)をwavファイル(1チャンネルのドライソース信号)に畳み込むMATLABプログラム. load LHRIR load RHRIR [wn,fs,bits]=wavread('hoge1.wav'); voiceLR(:,1)=conv(LHRIR,wn); voiceLR(:,2)=conv(RHRIR,wn…

ボリュームデータの可視化

MATLABを使ってボリュームデータを可視化する方法を示します. %===parameter===% th=0.5; SIGMA=10; x=0:30;y=0:30;z=0:30; [X,Y,Z]=meshgrid(x,y,z); Vol=exp(-((X-15).^2+(Y-15).^2+(Z-15).^2)/2/SIGMA/SIGMA); Vol=Vol+0.25*rand(size(Vol)); %===3D ima…