wolffd@0: 
wolffd@0: 
wolffd@0: %seglengths
wolffd@0: %vols
wolffd@0: %periods
wolffd@0: %name
wolffd@0: 
wolffd@0: % 
wolffd@0: % datlen=size(data,1);
wolffd@0: % maxseglen=10;
wolffd@0: % clear segvol;
wolffd@0: % clear seglengths;
wolffd@0: % clear repetitions;
wolffd@0: % clear name;
wolffd@0: % 
wolffd@0: % filectr=0;
wolffd@0: % lineptr=15;
wolffd@0: % 
wolffd@0: % while lineptr <= size(textdata,1)
wolffd@0: %     if ~strcmp('',textdata{lineptr,1}) 
wolffd@0: %         path = textdata{lineptr,1};
wolffd@0: %     end
wolffd@0: %     if ~strcmp('',textdata{lineptr,2}) 
wolffd@0: %         filectr=filectr+1;
wolffd@0: %         name{filectr} = strcat(path,'\',textdata{lineptr,2});
wolffd@0: %         datapos=lineptr-14;
wolffd@0: %         
wolffd@0: %         segctr=0;
wolffd@0: %         while ~isnan(data(datapos,2))&& datapos < datlen
wolffd@0: %             segctr=segctr+1;
wolffd@0: %             seglengths(filectr,segctr)= data(datapos,2);
wolffd@0: %             if ~isnan(data(datapos,3))
wolffd@0: %                 repetitions(filectr,segctr)=data(datapos,3);
wolffd@0: %             end
wolffd@0: %             segvol(filectr,segctr)=data(datapos,4);
wolffd@0: %             datapos=datapos+1;
wolffd@0: %             lineptr = lineptr+1;
wolffd@0: %         end
wolffd@0: %     end
wolffd@0: %     lineptr = lineptr+1;
wolffd@0: % end
wolffd@0: 
wolffd@0: 
wolffd@0: compl_len=sum(seglengths,2);%ins secs
wolffd@0: compl_std=std(compl_len);
wolffd@0: compl_mean=mean(compl_len);
wolffd@0: 
wolffd@0: numsegs=seglengths>0;
wolffd@0: numsegs=sum(numsegs,2);
wolffd@0: meannumsegs=mean(numsegs)
wolffd@0: minnumsegs=min(numsegs)
wolffd@0: maxnumsegs=max(numsegs)
wolffd@0: 
wolffd@0: for i=1:maxnumsegs %segmentlänge, betrachtet werden jeweils alle segs 
wolffd@0:                     % mit ausreichender länge
wolffd@0:     valid_cols=find(numsegs>=i);
wolffd@0:     mean_seglen(i)=mean(seglengths(valid_cols,i));
wolffd@0:     min_seglen(i)=min(seglengths(valid_cols,i));
wolffd@0:     max_seglen(i)=max(seglengths(valid_cols,i));
wolffd@0:     std_seglen(i)=std(seglengths(valid_cols,i));%in s
wolffd@0: end
wolffd@0: 
wolffd@0: % valid_cols=find(numsegs==floor(meannumsegs));
wolffd@0: % for i=1:floor(mean(numsegs)) %segmentlänge, betrachtet werden jeweils alle segs 
wolffd@0: %                     % mit genau floor(mittlerer) länge
wolffd@0: %     mean_seglen2(i)=mean(seglengths(valid_cols,i));
wolffd@0: %     min_seglen2(i)=min(seglengths(valid_cols,i));
wolffd@0: %     max_seglen2(i)=max(seglengths(valid_cols,i));
wolffd@0: %     std_seglen2(i)=std(seglengths(valid_cols,i));%in s
wolffd@0: % end
wolffd@0: 
wolffd@0: %längen
wolffd@0: seglen_all=reshape(seglengths,115*6,1);
wolffd@0: goodlens=find(seglen_all>0);
wolffd@0: mean_seglen_all=mean(seglen_all(goodlens));
wolffd@0: std_seglen_all=std(seglen_all(goodlens));
wolffd@0: for i=1:maxnumsegs
wolffd@0:     segs(i)=sum(numsegs==i);
wolffd@0: end
wolffd@0: 
wolffd@0: period_segs=sum(repetitions>0,2);
wolffd@0: for i=1:maxnumsegs
wolffd@0:     havenumperiodsegs(i)=sum(period_segs==i);
wolffd@0: end
wolffd@0: segs_period=(repetitions>0);
wolffd@0: sum((period_segs==3)&(numsegs==4))
wolffd@0: sum((period_segs==2)&(numsegs==4))
wolffd@0: sum(((period_segs==4)&(numsegs==5)))
wolffd@0: sum(((period_segs==3)&(numsegs==5)))
wolffd@0: 
wolffd@0: better_half_ind=find((period_segs==3)&(numsegs==4));
wolffd@0: better_half_len=seglengths(better_half_ind,:);
wolffd@0: better_half_mean=mean(better_half_len,1);
wolffd@0: better_half_std=std(better_half_len);
wolffd@0: 
wolffd@0: %wiederholungen
wolffd@0: period_all=reshape(repetitions,115*6,1);
wolffd@0: goodreps=find(period_all>0);
wolffd@0: period_all=period_all(goodreps);
wolffd@0: peri_all_mean=mean(period_all);
wolffd@0: peri_all_std=std(period_all);
wolffd@0: 
wolffd@0: %frequenzen
wolffd@0: freqs = repetitions./(seglengths./1000);
wolffd@0: freqs_all=reshape(freqs,115*6,1);
wolffd@0: freqs_all=freqs_all(find(~isnan(freqs_all) & (freqs_all>0)));
wolffd@0: freqs_all_mean=mean(freqs_all);
wolffd@0: freqs_all_std=std(freqs_all);
wolffd@0: 
wolffd@0: 
wolffd@0: better_h_f=repetitions(better_half_ind,2:4)./(seglengths(better_half_ind,2:4)./1000);
wolffd@0: better_h_f_mean=mean(better_h_f);
wolffd@0: better_h_f_std=std(better_h_f);
wolffd@0: better_h_f_min=min(better_h_f);
wolffd@0: better_h_f_max=max(better_h_f);
wolffd@0: 
wolffd@0: better_half2_ind=find((period_segs==4)&(numsegs==5));
wolffd@0: 
wolffd@0: better_h2_f=repetitions(better_half2_ind,2:5)./(seglengths(better_half2_ind,2:5)./1000);
wolffd@0: better_h2_f_mean=mean(better_h2_f);
wolffd@0: better_h2_f_std=std(better_h2_f);
wolffd@0: better_h2_f_min=min(better_h2_f);
wolffd@0: better_h2_f_max=max(better_h2_f);
wolffd@0: 
wolffd@0: abweichung_f=better_h_f-repmat(better_h_f_mean,size(better_h_f,1),1);
wolffd@0: corrcoef(abweichung_f)
wolffd@0: abweichung2_f=better_h2_f-repmat(better_h2_f_mean,size(better_h2_f,1),1);
wolffd@0: corrcoef(abweichung2_f)