instrucmentnote好乱啊，待整理

clear all;clc;close all;                                            
 load noteoctive
%=======================================================
strFileName={'violin'};%；flute；；piano1 ；trumpet piano2； trumpet0.5
 mkdir('parentFolder',strFileName{1});
 currPath = fileparts(mfilename('fullpath'));%[ST,I] = dbstack('-completenames')
 limitamp=0.5;         %violin0.5;flute0.5;piano0.6       piano1 0.4  piano2  1.5
 [s1,fs,nbits]=wavread(strFileName{1});               % s1=s1(440000:820000);
 %=======================================================
initialn=25;
findthreshold1;
finddiv_s;
savewavefiles;

B=[1,-0.9375];s=filter(B,1,s1);  %clear s1;
s=s/max(abs(s));      lenwinms=10;                      
lenwin=(lenwinms/1000)*fs;
noverlap=lenwin/3;nfft=lenwin; turesiglength=lenwin-noverlap;                    %
sigFramed = buffer(s,lenwin, noverlap,'nodelay'); %每一列代表一帧
%clear s;
%amp = sum(abs(sigFramed),1) ;%计算出信号的能量，此处用平均幅度代替

amp(find(amp<limitamp))=0;
amplogical=(amp~=0);
notei=amplogical;
difnotei=diff(notei);
Index1=find(difnotei==1);Index2=find(difnotei==-1);
lengthI=min(length(Index1),length(Index2));
Index1=Index1(1:lengthI);Index2=Index2(1:lengthI);
difff=Index2-Index1;
difIndex1=find(difff>100);%100difIndex2=find(difff<20);%win=hamming(lenwin); 
note_num=length(difIndex1);
%figure;subrow=ceil(note_num/5);
%mkdir('parentFolder1',strFileName{1});
for i=1:note_num
endpoint1(i)=turesiglength.*(Index1(difIndex1(i))-1)+1;
endpoint2(i)=turesiglength.*(Index2(difIndex1(i)));
% note{i}=s((turesiglength.*(Index1(difIndex1(i))-1)+1):turesiglength.*(Index2(difIndex1(i))));%%%%？？？？？？？？notei的值才是帧的下标？？？？？？？
note{i}=s(endpoint1(i):endpoint2(i));
% % tmp=[note{i}(1);note{i}(1:end-1)];%%zc
% % signs=(tmp.*note{i})<0;
% % signs=buffer(signs,lenwin, noverlap,'nodelay');
% % diffs=abs(tmp-note{i})>25*1e-4;%25*1e-4;3*1e-3; 求过零率时设定了一定的门限，此处为
% % diffs=buffer(diffs,lenwin, noverlap,'nodelay');
% % zcr=sum(signs.*diffs,1);III=find(zcr~=0,1,'first');JJJ=find(zcr~=0,1,'last');
% % onset1(i)=turesiglength.*(III-1)+1;onset2(i)=turesiglength.*(JJJ);
% % note{i}=note{i}(onset1(i):onset2(i));
% note{i}=note{i}((turesiglength.*(III-1)+1):turesiglength.*(JJJ));
% % % 
notecor=xcorr(note{i});
notecor=notecor(round(end/2):end);%notecor=notecor(length(notecor)/2:end);
d = diff(notecor);d1=[d(1);d(1:end-1)];
indmax = find(d1.*d<0 & d1>0);%%%局部极大值 
[kk,II]=max(notecor(indmax));
pitch(i)=fs./indmax(II);

dif_pitch=abs(pitch(i)-noteoctive.not_octij);
[notelabel_i(i),octivelabel_i(i)]=find(dif_pitch==min(min(dif_pitch)));
if i>1&(notelabel_i(i)==notelabel_i(i-1))&(octivelabel_i(i)==octivelabel_i(i-1))
    wavename{i}=[strFileName{1},'_',noteoctive.notelabel{notelabel_i(i)},noteoctive.octlabel{octivelabel_i(i)},'(2)'];
else
wavename{i}=[strFileName{1},'_',noteoctive.notelabel{notelabel_i(i)},noteoctive.octlabel{octivelabel_i(i)}];
end
% % N = sprintf('%6d', [strFileName{1},'_',noteoctive.notelabel{notelabel_i},noteoctive.octlabel{octivelabel_i}])
% % if i>1&(wavename{i}==wavename{i-1})
% %    wavename{i}=strcat(wavename{i},'(2)');
% % end
% if i<21
% subplot(4,5,i);plot(note{i});%%%画出分割结果
% % wavename=[strFileName{1},'_',noteoctive.notelabel{notelabel_i},noteoctive.octlabel{octivelabel_i}];
% title(wavename{i});
% else
%     if i==21
%         figure 
%     end
%     subplot(subrow-4,5,i-20);plot(note{i});%%%画出分割结果
% %wavename=[strFileName{1},'_',noteoctive.notelabel{notelabel_i},noteoctive.octlabel{octivelabel_i}];
% title(wavename{i});
% end
% wavwrite(note{i},fs,nbits,[currPath,'\parentFolder\',strFileName{1},'\',wavename{i},'.wav']); %
end

% figure;
% plot(s1);hold on;
% for i=1:note_num
% line([endpoint1(i) endpoint1(i)], [-1 1], 'Color', 'r');hold on;%endpoint1(i)
% line([endpoint2(i) endpoint2(i)], [-1 1], 'Color', 'b');hold on;%endpoint2(i)
% end
% title([strFileName{1},'result']);
% hold off;

%savewavefiles;
save ('noteesemble.mat','s1','fs','nbits','endpoint1',...
      'endpoint2','wavename','note_num','currPath','strFileName',...
      'lenwin','noverlap','turesiglength');%'onset1','onset2',
clear all;%clc;
load noteesemble.mat;
figure;subrow=ceil(note_num/5);
for i=1:note_num
note{i}=s1(endpoint1(i):endpoint2(i));
removenoise;
%note{i}=note{i}(onset1(i):onset2(i));
if i<21
subplot(4,5,i);plot(note{i});
% wavename=[strFileName{1},'_',noteoctive.notelabel{notelabel_i},noteoctive.octlabel{octivelabel_i}];
title(wavename{i});
else
if i==21
    figure 
end
subplot(subrow-4,5,i-20);plot(note{i});
%wavename=[strFileName{1},'_',noteoctive.notelabel{notelabel_i},noteoctive.octlabel{octivelabel_i}];
title(wavename{i});
end
wavwrite(note{i},fs,nbits,[currPath,'\parentFolder\',strFileName{1},'\',wavename{i},'.wav']); %
end
figure;
plot(s1);hold on;
for i=1:note_num
line([endpoint1(i) endpoint1(i)], [-1 1], 'Color', 'r');hold on;%endpoint1(i)
line([endpoint2(i) endpoint2(i)], [-1 1], 'Color', 'b');hold on;%endpoint2(i)
end
title([strFileName{1},' segment result']);
hold off;

可去掉removenoise，没用，

%note{i}\
absnote{i}=note{i}/max(abs(note{i}));                 %
sigFramednote = buffer(note{i},lenwin, noverlap,'nodelay'); 
amp = sum(abs(sigFramednote),1) ;
noteind=find(amp<(1/30*max(amp)));
difnoteind=diff(noteind);
maxdifnoteind=max(difnoteind);
if maxdifnoteind>100
onset=find(difnoteind==maxdifnoteind);
note{i}=note{i}((turesiglength.*(noteind(onset)-1)+1):turesiglength.*(noteind(onset+1)));
else
note{i}=note{i};
end