diff src/matlab/chromagram_E.m @ 0:c52bc3e8d3ad tip

user: boblsturm branch 'default' added README.md added assets/.DS_Store added assets/playButton.jpg added assets/stopButton.png added assets/swapButton.jpg added data/.DS_Store added data/fiveoctaves.mp3 added data/glock2.wav added data/sinScale.mp3 added data/speech_female.mp3 added data/sweep.wav added nimfks.m.lnk added src/.DS_Store added src/matlab/.DS_Store added src/matlab/AnalysisCache.m added src/matlab/CSS.m added src/matlab/DataHash.m added src/matlab/ExistsInCache.m added src/matlab/KLDivCost.m added src/matlab/LoadFromCache.m added src/matlab/SA_B_NMF.m added src/matlab/SaveInCache.m added src/matlab/Sound.m added src/matlab/SynthesisCache.m added src/matlab/chromagram_E.m added src/matlab/chromagram_IF.m added src/matlab/chromagram_P.m added src/matlab/chromsynth.m added src/matlab/computeSTFTFeat.m added src/matlab/controller.m added src/matlab/decibelSliderReleaseCallback.m added src/matlab/drawClickCallBack.m added src/matlab/fft2chromamx.m added src/matlab/hz2octs.m added src/matlab/ifgram.m added src/matlab/ifptrack.m added src/matlab/istft.m added src/matlab/nimfks.fig added src/matlab/nimfks.m added src/matlab/nmfFn.m added src/matlab/nmf_beta.m added src/matlab/nmf_divergence.m added src/matlab/nmf_euclidean.m added src/matlab/prune_corpus.m added src/matlab/rot_kernel.m added src/matlab/templateAdditionResynth.m added src/matlab/templateDelCb.m added src/matlab/templateScrollCb.m
author boblsturm
date Sun, 18 Jun 2017 06:26:13 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/src/matlab/chromagram_E.m	Sun Jun 18 06:26:13 2017 -0400
@@ -0,0 +1,28 @@
+function Y = chromagram_E(d,sr,fftlen,nbin,f_ctr,f_sd)
+% Y = chromagram_E(d,sr,fftlen,nbin)
+%  Calculate a "chromagram" of the sound in d (at sampling rate sr)
+%  Use windows of fftlen points, hopped by ffthop points
+%  Divide the octave into nbin steps
+%  Weight with center frequency f_ctr (in Hz) and gaussian SD f_sd (in octaves)
+% 2006-09-26 dpwe@ee.columbia.edu
+
+if nargin < 3;   fftlen = 2048; end
+if nargin < 4;   nbin = 12; end
+if nargin < 5;   f_ctr = 1000; end
+if nargin < 6;   f_sd = 1; end
+
+fftwin = fftlen/2;
+ffthop = fftlen/4;  % always for this
+
+D = abs(specgram(d,fftlen,sr,fftwin,(fftwin-ffthop)));
+
+A0 = 27.5; % Hz
+A440 = 440; % Hz
+
+f_ctr_log = log(f_ctr/A0) / log(2);
+
+CM = fft2chromamx(fftlen, nbin, sr, A440, f_ctr_log, f_sd);
+% Chop extra dims
+CM = CM(:,1:(fftlen/2)+1);
+
+Y = CM*D;