   % debugging of vsp
imp=zeros(1,1024);
imp(1)=1;
for ii=1:72
ir(:,ii)=peakfilt(imp,ii*300+100,50,+20,44100);
end   
mn_plot(20*log10(abs(fft(ir))), [0 1200 -1 20]);
keyboard;
noi=rand(44100*5,1);
figure;
subplot(1,3,1);
specgram(conv(ir(:,4),noi));
subplot(1,3,2);
specgram(conv(ir(:,8),noi),1024,44100);
subplot(1,3,3);
specgram(conv(ir(:,12),noi),1024,44100);
keyboard;
out=vsp(noi,ir,ir,[0 45 90 45 0 -45 -90 -45 0]);
figure;
specgram(out(:,1));
