-
Notifications
You must be signed in to change notification settings - Fork 1
/
acontrast_demo.m
49 lines (37 loc) · 1.35 KB
/
acontrast_demo.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
% visualization parameters
font_size=12; % self-explanatory
disp_band=[0 12e3]; % spectrogram band to display
clipping=0; % clipping for spectrogram
imscale=15; % image scaling parameter
regression_timescale=.0075; % sets the smoothness of the derivative calculation (in s)
pow_scale='log'; % use log power weighting (try 'lin' for linear)
pow_weight=1; % enable power weighting
% Load the sample sound
[signal, fs] = wavread('finchdoublet.wav');
signal_t=[1:length(signal)]/fs;
tic
[feature_matrix,raw_features,feature_names,consensus,f,t]=acontrast_contour(signal,fs,...
'regression_timescale',regression_timescale,'pow_weight',pow_weight,'pow_scale',pow_scale);
[envelope_score,envelope,sonogram,sono_f,sono_t]=acontrast_envelope(signal,fs,'regression_timescale',regression_timescale);
contour_score=sum(abs(feature_matrix));
toc
%% Visualization
figure
h(1) = subplot(4,1,1:2);
imagesc(sono_t,sono_f,max(log(abs(sonogram)+imscale),clipping));
ylim([disp_band])
ylabel('Frequency [Hz]','FontSize', font_size);
set(gca,'XTick',[]);
axis xy;
title('Sonogram','FontSize', font_size);
colormap(hot)
h(2) = subplot(4,1,3);
plot(signal_t,contour_score);
ylabel('Cont.','FontSize', font_size);
set(gca,'XTick',[]);
axis xy;
h(3)=subplot(4,1,4);
plot(signal_t,abs(envelope_score));
ylabel('Env.','FontSize',font_size);
xlabel('Time (in s)','FontSize',font_size);
linkaxes(h,'x')