1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
| figure ;
[...]
Esub = 0.2; % subplot thick
Ysub = 0.1 ;
for ii=1:length(filecodecanal)
% ouverture fichiers et données
dataName{ii}=[datadir fileroot filecodecanal{ii} '.' solNumb fileappend]
[DAT,I]=rdmseed2(dataName{ii});
dataRaw = cat(1,DAT(:).d);
timeRaw = cat(1,DAT(:).t);
data{ii,:}=dataRaw(:);
Head{ii}=I;
time{ii}=timeRaw;
% calcul spectre des données
[f{ii},s{ii}] = spectre(data{ii,:},Fs(ii));
sub1=subplot('Position',[.05 Ysub .45 Esub]) ;
set(sub1,'YTick',[],'XTick',[]);
plot(time{ii},data{ii,:},'k')
incSub=incSub+2;
hold on
sub2=subplot('Position',[.55 Ysub .45 Esub]) ;
set(sub2,'YTick',[],'XTick',[]);
semilogy(f{ii},s{ii})
Ysub = Ysub + (0.6/length(filecodecanal))
end |
Partager