%%%% example: load Virgo_SeismicData file and make one sample graph (figure %%%% 3 of VIR-0390A-15) %%%% load('VIR-0390A-15.mat') figure(3) n1=1; % VSR4 data, EW %n1=4; % 1 yr data, EW loglog(VirgoSeismicPL(n1).f(:,7),VirgoSeismicPL(n1).spectra(:,1),'k') hold on loglog(VirgoSeismicPL(n1).f(:,6),VirgoSeismicPL(n1).spectra(:,3),'b') loglog(VirgoSeismicPL(n1).f(:,5),VirgoSeismicPL(n1).spectra(:,4),'c') loglog(VirgoSeismicPL(n1).f(:,4),VirgoSeismicPL(n1).spectra(:,5),'g') loglog(VirgoSeismicPL(n1).f(:,3),VirgoSeismicPL(n1).spectra(:,6),'y') loglog(VirgoSeismicPL(n1).f(:,3),VirgoSeismicPL(n1).spectra(:,7),'m') loglog(VirgoSeismicPL(n1).f(:,3),VirgoSeismicPL(n1).spectra(:,9),'r') grid on xlabel('frequency [Hz]') ylabel(VirgoSeismicPL(n1).units) hold off legend('1-st percentile','10-th percentile','25-th percentile','50-th percentile','75-th percentile','90-th percentile','99-th percentile') title(['V1:',VirgoSeismicPL(1).sensor.description]) hold off axis([0.05 150 1e-12 5e-5]) print( gcf, '-dpng', '-r95','Virgo_CEB_percentiles.png');