Home > melodyRecognition > pv2note.m

pv2note

PURPOSE ^

pv2note: Note segmentation on a give pitch vector

SYNOPSIS ^

function [note, pvCountInNote] = pv2note(pv, timeUnit, pitchTh, minNoteDuration, plotOpt)

DESCRIPTION ^

 pv2note: Note segmentation on a give pitch vector
    Usage: note = pv2note(pv, timeUnit, pitchTh, minNoteDuration, plotOpt)

    For example:
        pv=[49.21841 48.33041 48.33041 47.76274 48.04425 48.04425 47.76274 48.33041 48.04425 48.04425 48.04425 47.48574 47.48574 47.48574 47.48574 47.48574 47.76274 48.04425 48.04425 47.76274 47.76274 47.76274 47.21309 47.21309 47.21309 47.21309 47.21309 51.48682 51.48682 51.48682 51.48682 51.83658 51.83658 51.83658 51.83658 51.83658 54.92247 54.09792 54.09792 54.09792 54.09792 54.09792 54.09792 54.92247 54.92247 54.92247 54.92247 54.50529 54.50529 54.92247 54.50529 54.50529 54.09792 54.50529 55.34996 54.92247 54.92247 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 53.31086 53.31086 53.31086 53.31086 56.23796 55.78827 55.78827 56.23796 56.69965 56.69965 56.69965 57.17399 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 57.17399 57.17399 57.17399 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 59.76274 59.76274 59.76274 59.21309 59.21309 59.21309 58.68037 58.68037 58.68037 58.68037 54.09792 54.09792 54.09792 54.50529 54.50529 54.09792 54.09792 54.50529 54.92247 54.50529 54.50529 54.09792 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 50.47805 51.48682 51.83658 52.19354 52.55803 52.19354 52.55803 51.83658 51.83658 52.55803 52.93035 52.93035 52.93035 52.93035 52.93035 51.14399 51.14399 54.50529 53.31086 52.55803 52.19354 52.19354 52.19354 52.55803 52.93035 54.09792 54.50529 54.92247 55.78827 56.23796 56.23796 55.78827 55.34996 54.09792 54.09792 54.09792 51.48682 50.15444 50.15444 50.80782 51.14399 51.14399 51.14399 51.14399 52.19354 52.19354 51.83658 51.83658 51.83658 51.48682 51.48682 51.48682 51.83658 51.83658 51.48682 51.48682 51.48682 51.48682 51.48682 50.80782 50.80782 52.55803 51.48682 51.14399 50.80782 51.14399 51.48682 51.48682 51.48682 50.80782 50.80782 50.80782 48.91732 48.62138 48.33041 48.33041 48.33041 48.04425 48.91732 48.91732 48.91732 49.21841 49.21841 48.91732 48.62138 48.33041 48.33041 48.33041 49.83678 48.62138 48.62138 48.62138 48.62138 48.62138 48.91732 49.52484 49.52484 48.91732 48.62138 48.33041];
        timeUnit=256/8000;
        pitchTh=0.8;
        minNoteDuration=0.1;
        plotOpt=1;
        note=pv2note(pv, timeUnit, pitchTh, minNoteDuration, plotOpt);
        fprintf('Hit return to hear the pitch vector...'); pause; fprintf('\n');
        pvPlay(pv, timeUnit);
        fprintf('Hit return to hear the segmented note...'); pause; fprintf('\n');
        notePlay(note);

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SUBFUNCTIONS ^

SOURCE CODE ^

0001 function [note, pvCountInNote] = pv2note(pv, timeUnit, pitchTh, minNoteDuration, plotOpt)
0002 % pv2note: Note segmentation on a give pitch vector
0003 %    Usage: note = pv2note(pv, timeUnit, pitchTh, minNoteDuration, plotOpt)
0004 %
0005 %    For example:
0006 %        pv=[49.21841 48.33041 48.33041 47.76274 48.04425 48.04425 47.76274 48.33041 48.04425 48.04425 48.04425 47.48574 47.48574 47.48574 47.48574 47.48574 47.76274 48.04425 48.04425 47.76274 47.76274 47.76274 47.21309 47.21309 47.21309 47.21309 47.21309 51.48682 51.48682 51.48682 51.48682 51.83658 51.83658 51.83658 51.83658 51.83658 54.92247 54.09792 54.09792 54.09792 54.09792 54.09792 54.09792 54.92247 54.92247 54.92247 54.92247 54.50529 54.50529 54.92247 54.50529 54.50529 54.09792 54.50529 55.34996 54.92247 54.92247 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 53.31086 53.31086 53.31086 53.31086 56.23796 55.78827 55.78827 56.23796 56.69965 56.69965 56.69965 57.17399 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 57.17399 57.17399 57.17399 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 59.76274 59.76274 59.76274 59.21309 59.21309 59.21309 58.68037 58.68037 58.68037 58.68037 54.09792 54.09792 54.09792 54.50529 54.50529 54.09792 54.09792 54.50529 54.92247 54.50529 54.50529 54.09792 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 50.47805 51.48682 51.83658 52.19354 52.55803 52.19354 52.55803 51.83658 51.83658 52.55803 52.93035 52.93035 52.93035 52.93035 52.93035 51.14399 51.14399 54.50529 53.31086 52.55803 52.19354 52.19354 52.19354 52.55803 52.93035 54.09792 54.50529 54.92247 55.78827 56.23796 56.23796 55.78827 55.34996 54.09792 54.09792 54.09792 51.48682 50.15444 50.15444 50.80782 51.14399 51.14399 51.14399 51.14399 52.19354 52.19354 51.83658 51.83658 51.83658 51.48682 51.48682 51.48682 51.83658 51.83658 51.48682 51.48682 51.48682 51.48682 51.48682 50.80782 50.80782 52.55803 51.48682 51.14399 50.80782 51.14399 51.48682 51.48682 51.48682 50.80782 50.80782 50.80782 48.91732 48.62138 48.33041 48.33041 48.33041 48.04425 48.91732 48.91732 48.91732 49.21841 49.21841 48.91732 48.62138 48.33041 48.33041 48.33041 49.83678 48.62138 48.62138 48.62138 48.62138 48.62138 48.91732 49.52484 49.52484 48.91732 48.62138 48.33041];
0007 %        timeUnit=256/8000;
0008 %        pitchTh=0.8;
0009 %        minNoteDuration=0.1;
0010 %        plotOpt=1;
0011 %        note=pv2note(pv, timeUnit, pitchTh, minNoteDuration, plotOpt);
0012 %        fprintf('Hit return to hear the pitch vector...'); pause; fprintf('\n');
0013 %        pvPlay(pv, timeUnit);
0014 %        fprintf('Hit return to hear the segmented note...'); pause; fprintf('\n');
0015 %        notePlay(note);
0016 
0017 %    Roger Jang, 20060527, 20070531
0018 
0019 if nargin<1, selfdemo; return; end
0020 if nargin<2, timeUnit=256/8000; end
0021 if nargin<3, pitchTh=0.8; end
0022 if nargin<4, minNoteDuration=0.1; end
0023 if nargin<5, plotOpt=0; end
0024 
0025 note=[];
0026 pvCountInNote=[];
0027 pvInCurrentNote=pv(1);        % PV in the current note
0028 for i=2:length(pv)
0029     notePitch=median(pvInCurrentNote);
0030     noteDuration=length(pvInCurrentNote)*timeUnit;
0031     if abs(pv(i)-pv(i-1))>pitchTh & noteDuration>minNoteDuration        % If the pitch diff. is high and note duration is long...
0032         note=[note, notePitch, noteDuration];    % Create a new note
0033         pvCountInNote=[pvCountInNote, length(pvInCurrentNote)];        % Create a new note
0034         pvInCurrentNote=pv(i);
0035     else
0036         pvInCurrentNote=[pvInCurrentNote, pv(i)];
0037     end
0038 end
0039 % Process the last note
0040 note=[note, median(pvInCurrentNote), length(pvInCurrentNote)*timeUnit];        % Create the last note
0041 pvCountInNote=[pvCountInNote, length(pvInCurrentNote)];                % Create the last note
0042 
0043 if plotOpt
0044     plot((1:length(pv))*timeUnit, pv, '.-g');
0045     hold on; notePlot(note, 1); hold off
0046 end
0047 
0048 % ====== Self demo
0049 function selfdemo
0050 pv=[49.21841 48.33041 48.33041 47.76274 48.04425 48.04425 47.76274 48.33041 48.04425 48.04425 48.04425 47.48574 47.48574 47.48574 47.48574 47.48574 47.76274 48.04425 48.04425 47.76274 47.76274 47.76274 47.21309 47.21309 47.21309 47.21309 47.21309 51.48682 51.48682 51.48682 51.48682 51.83658 51.83658 51.83658 51.83658 51.83658 54.92247 54.09792 54.09792 54.09792 54.09792 54.09792 54.09792 54.92247 54.92247 54.92247 54.92247 54.50529 54.50529 54.92247 54.50529 54.50529 54.09792 54.50529 55.34996 54.92247 54.92247 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 54.50529 53.31086 53.31086 53.31086 53.31086 56.23796 55.78827 55.78827 56.23796 56.69965 56.69965 56.69965 57.17399 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 57.17399 57.17399 57.17399 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 56.69965 59.76274 59.76274 59.76274 59.21309 59.21309 59.21309 58.68037 58.68037 58.68037 58.68037 54.09792 54.09792 54.09792 54.50529 54.50529 54.09792 54.09792 54.50529 54.92247 54.50529 54.50529 54.09792 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 53.69992 50.47805 51.48682 51.83658 52.19354 52.55803 52.19354 52.55803 51.83658 51.83658 52.55803 52.93035 52.93035 52.93035 52.93035 52.93035 51.14399 51.14399 54.50529 53.31086 52.55803 52.19354 52.19354 52.19354 52.55803 52.93035 54.09792 54.50529 54.92247 55.78827 56.23796 56.23796 55.78827 55.34996 54.09792 54.09792 54.09792 51.48682 50.15444 50.15444 50.80782 51.14399 51.14399 51.14399 51.14399 52.19354 52.19354 51.83658 51.83658 51.83658 51.48682 51.48682 51.48682 51.83658 51.83658 51.48682 51.48682 51.48682 51.48682 51.48682 50.80782 50.80782 52.55803 51.48682 51.14399 50.80782 51.14399 51.48682 51.48682 51.48682 50.80782 50.80782 50.80782 48.91732 48.62138 48.33041 48.33041 48.33041 48.04425 48.91732 48.91732 48.91732 49.21841 49.21841 48.91732 48.62138 48.33041 48.33041 48.33041 49.83678 48.62138 48.62138 48.62138 48.62138 48.62138 48.91732 49.52484 49.52484 48.91732 48.62138 48.33041];
0051 timeUnit=256/8000;
0052 pitchTh=0.8;
0053 minNoteDuration=0.1;
0054 plotOpt=1;
0055 [note, pvCountInNote]=feval(mfilename, pv, timeUnit, pitchTh, minNoteDuration, plotOpt)
0056 fprintf('Hit return to hear the pitch vector...'); pause; fprintf('\n');
0057 pvPlay(pv, timeUnit);
0058 fprintf('Hit return to hear the segmented note...'); pause; fprintf('\n');
0059 notePlay(note, 1);

Generated on Tue 01-Jun-2010 09:49:37 by m2html © 2003