Chris@3: 
Chris@3: load MIDI_FB_ellip_pitch_60_96_22050_Q25.mat;
Chris@3: 
Chris@3: dirac = zeros(50000, 1);
Chris@3: dirac(1) = 1.0;
Chris@3: 
Chris@3: delays = zeros(120, 1);
Chris@3: 
Chris@3: for n = 21:120
Chris@3:     f = filter(h(n).b, h(n).a, dirac);
Chris@14:     [~,pos] = max(f(2:10000));
Chris@14:     [~,neg] = max(-f(2:10000));
Chris@14:     delays(n) = 2 + (pos + neg) / 2;
Chris@3: end
Chris@14: 
Chris@14: csvwrite('delays.csv', round(delays));