annotate demo/demo_aim2006_dcgc_one.m @ 4:537f939baef0 tip

various bug fixes and changed copyright message
author Stefan Bleeck <bleeck@gmail.com>
date Tue, 16 Aug 2011 14:37:17 +0100
parents 74dedb26614d
children
rev   line source
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