tomwalters@0
|
1 % Script to compare the AIM2006 output for four files (listed below)
|
tomwalters@0
|
2 % This scrips creates six figure windows, one for each stage of the AIM
|
tomwalters@0
|
3 % processing model and plots four panes in each window, one for each file
|
tomwalters@0
|
4
|
tomwalters@0
|
5 % (c) 2006-2008, University of Cambridge, Medical Research Council
|
tomwalters@0
|
6 % Tom Walters tcw24@cam.ac.uk
|
tomwalters@0
|
7 % http://www.pdn.cam.ac.uk/cnbh/aim2006
|
tomwalters@0
|
8 % $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
|
tomwalters@0
|
9 % $Revision: 585 $
|
tomwalters@0
|
10
|
tomwalters@0
|
11 filenames{1}='aa110p122s100t.wav';
|
tomwalters@0
|
12 filenames{2}='aa256p122s100t.wav';
|
tomwalters@0
|
13 filenames{3}='aa110p089s100t.wav';
|
tomwalters@0
|
14 filenames{4}='aa256p089s100t.wav';
|
tomwalters@0
|
15
|
tomwalters@0
|
16
|
tomwalters@0
|
17 analysis_start_time=0;
|
tomwalters@0
|
18 analysis_end_time=0.2;
|
tomwalters@0
|
19
|
tomwalters@0
|
20 display_start_time=0.13;
|
tomwalters@0
|
21 display_end_time=0.154;
|
tomwalters@0
|
22
|
tomwalters@0
|
23 bmm_scaling=0.25;
|
tomwalters@0
|
24 nap_scaling=1;
|
tomwalters@0
|
25 sai_scaling=0.25;
|
tomwalters@0
|
26 mellin_scaling=16;
|
tomwalters@0
|
27
|
tomwalters@0
|
28 ii=1;
|
tomwalters@0
|
29 %for ii=1:length(filenames)
|
tomwalters@0
|
30 name=filenames{ii};
|
tomwalters@0
|
31
|
tomwalters@0
|
32 % Find the pitch and scale of the input sound from the filename
|
tomwalters@0
|
33 inp=strfind(name,'p');
|
tomwalters@0
|
34 pitch=str2num(name(inp-3:inp-1));
|
tomwalters@0
|
35 ins=strfind(name,'s');
|
tomwalters@0
|
36 scale=str2num(name(ins-3:ins-1));
|
tomwalters@0
|
37
|
tomwalters@0
|
38 % Set the parameters for analysis of the sound
|
tomwalters@0
|
39 dcgc_options=setparams_dcgc_display(filenames{ii}, analysis_start_time, analysis_end_time);
|
tomwalters@0
|
40
|
tomwalters@0
|
41 %dcgc_output{ii}.filename=name;
|
tomwalters@0
|
42 %dcgc_output{ii}.scale=scale;
|
tomwalters@0
|
43 %dcgc_output{ii}.pitch=pitch;
|
tomwalters@0
|
44
|
tomwalters@0
|
45 % Analyse with aim2006
|
tomwalters@0
|
46 dcgc_output{ii}=aim(dcgc_options);
|
tomwalters@0
|
47
|
tomwalters@0
|
48 % dcgc_output{ii}.filename=name;
|
tomwalters@0
|
49 % dcgc_output{ii}.scale=scale;
|
tomwalters@0
|
50 % dcgc_output{ii}.pitch=pitch;
|
tomwalters@0
|
51
|
tomwalters@0
|
52 options.minimum_time=display_start_time;
|
tomwalters@0
|
53 options.maximum_time=display_end_time;
|
tomwalters@0
|
54
|
tomwalters@0
|
55 title_string=['Pitch:' num2str(pitch) ' Scale:' num2str(scale)];
|
tomwalters@0
|
56
|
tomwalters@0
|
57 figure(1)
|
tomwalters@0
|
58 set(gcf, 'Name', ['PCP: ' dcgc_output{ii}.info.current_pcp_module]);
|
tomwalters@0
|
59 %subplot(2,2,ii);
|
tomwalters@0
|
60 plot(dcgc_output{ii}.data.pcp, [display_start_time display_end_time]);
|
tomwalters@0
|
61
|
tomwalters@0
|
62 figure(2)
|
tomwalters@0
|
63 set(gcf, 'Name', ['BMM: ' dcgc_output{ii}.info.current_bmm_module]);
|
tomwalters@0
|
64 %subplot(2,2,ii);
|
tomwalters@0
|
65 plot_bmm(dcgc_output{ii}, options, bmm_scaling,title_string);
|
tomwalters@0
|
66
|
tomwalters@0
|
67 figure(3)
|
tomwalters@0
|
68 set(gcf, 'Name', ['NAP: ' dcgc_output{ii}.info.current_nap_module]);
|
tomwalters@0
|
69 %subplot(2,2,ii);
|
tomwalters@0
|
70 plot_nap(dcgc_output{ii},options, nap_scaling,title_string);
|
tomwalters@0
|
71
|
tomwalters@0
|
72 figure(4)
|
tomwalters@0
|
73 set(gcf, 'Name', ['SP: ' dcgc_output{ii}.info.current_strobes_module]);
|
tomwalters@0
|
74 %subplot(2,2,ii);
|
tomwalters@0
|
75 %plot(dcgc_output{ii}.data.nap);
|
tomwalters@0
|
76 %hold on;
|
tomwalters@0
|
77 plot_strobes(dcgc_output{ii}, options, nap_scaling,title_string);
|
tomwalters@0
|
78
|
tomwalters@0
|
79 max_frame=length(dcgc_output{ii}.data.sai);
|
tomwalters@0
|
80
|
tomwalters@0
|
81 figure(5)
|
tomwalters@0
|
82 set(gcf, 'Name', ['SAI: ' dcgc_output{ii}.info.current_sai_module]);
|
tomwalters@0
|
83 %subplot(2,2,ii);
|
tomwalters@0
|
84 plot_sai(dcgc_output{ii},max_frame, sai_scaling,title_string);
|
tomwalters@0
|
85
|
tomwalters@0
|
86 figure(6)
|
tomwalters@0
|
87 set(gcf, 'Name', ['Usermodule: ' dcgc_output{ii}.info.current_usermodule_module]);
|
tomwalters@0
|
88 %subplot(2,2,ii);
|
tomwalters@0
|
89 plot_mellin(dcgc_output{ii},max_frame, mellin_scaling, title_string);
|
tomwalters@0
|
90
|
tomwalters@0
|
91 %end |