Продолжение приложения б

begtime=datenum(2014,07,31,01,54,13);

ZTname='M=1.7,Az=237';

timeDo=50; timePosle=110;

z1=-125; z2=-90;

case '19'

filenameX1='2014212025635320_POLE__1_2.atr';

filenameY1='2014212025635320_POLE__1_3.atr';

filenameX2='2014212035635320_POLE__1_2.atr';

filenameY2='2014212035635320_POLE__1_3.atr';

filetime=datenum(2014,07,31,02,56,35.320);

begtime=datenum(2014,07,31,03,43,5);

ZTname='M=2.7,Az=217';

timeDo=50; timePosle=120;

z1=-125; z2=-90;

case '20'

filenameX1='2014213035956320_POLE__1_2.atr';

filenameY1='2014213035956320_POLE__1_3.atr';

filenameX2='2014213045956320_POLE__1_2.atr';

filenameY2='2014213045956320_POLE__1_3.atr';

filetime=datenum(2014,08,01,03,59,56.320);

begtime=datenum(2014,08,01,04,46,0);

ZTname='M=2.1,Az=262';

timeDo=50; timePosle=150;

z1=-125; z2=-75;

case '21'

filenameX1='2014213035956320_POLE__1_2.atr';

filenameY1='2014213035956320_POLE__1_3.atr';

filenameX2='2014213045956320_POLE__1_2.atr';

filenameY2='2014213045956320_POLE__1_3.atr';

Продолжение приложения б

filetime=datenum(2014,08,01,03,59,56.320);

begtime=datenum(2014,08,01,04,49,29);

ZTname='M=1.4,Az=266';

timeDo=30; timePosle=70;

z1=-125; z2=-90;

case '22'

filenameX1='2014213035956320_POLE__1_2.atr';

filenameY1='2014213035956320_POLE__1_3.atr';

filenameX2='2014213045956320_POLE__1_2.atr';

filenameY2='2014213045956320_POLE__1_3.atr';

filetime=datenum(2014,08,01,03,59,56.320);

begtime=datenum(2014,08,01,04,51,32);

ZTname='M=0.8,Az=262';

timeDo=30; timePosle=70;

z1=-125; z2=-100;

end

fmin=.3; dt=3;

delimeter=' ';

nStrok=9; L=86400; fs=200;

m=1.583e-6/2001.67*1000;

% load

data666 = importdata([p filenameX1], delimeter, nStrok);

dataX=data666.data;

data666 = importdata([p filenameX2], delimeter, nStrok);

dataX=[dataX; data666.data]*m;

data666 = importdata([p filenameY1], delimeter, nStrok);

dataY=data666.data;

data666 = importdata([p filenameY2], delimeter, nStrok);

dataY=[dataY; data666.data]*m;

Продолжение приложения б

cleardata666;

% filter

[b,a]=butter(3,fmin/fs*2,'high');

dataXf=filter(b,a,dataX-mean(dataX));

dataYf=filter(b,a,dataY-mean(dataY));

% cut

t1=round((begtime-timeDo/L-filetime)*L*fs);

t2=round((begtime+timePosle/L-filetime)*L*fs);

X=dataXf(t1:t2); Y=dataYf(t1:t2);

% spectrogram

sp01(X,fs,1000,900,0,z1,z2)

shadinginterp, ylim([1,30]);

print(gcf,'-dpng',[p2 nomer '_' ZTname '_X.png'],'-r300')

close(gcf);

sp01(Y,fs,1000,900,0,z1,z2)

shadinginterp, ylim([1,30]);

print(gcf,'-dpng',[p2 nomer '_' ZTname '_Y.png'],'-r300')

close(gcf);

function sp01(fname,fs,wl,overlap, meaning,zmin,zmax)

%load file

ifischar(fname)

delimeter=' ';

nStrok=9;

data666 = importdata(fname, delimeter, nStrok);

data1=data666.data;

clear data666;

else

data1=fname;

end

Продолжение приложения б

% computing

if meaning>1

data1=AntiTrendFast(data1,meaning);

end

[~,F,T,P]=spectrogram(data1,hann(wl),overlap,wl,fs);

F(1:3,:)=[];

P(1:3,:)=[];

%plot 1st

figure;

subplot('position',[0.04 0.75 0.94 0.22]);

set(gca,'fontSize',9)

plot((1:length(data1))./fs,data1); axis tight; grid on;

aaa=get(gca,'ylim');

aaa(1)=aaa(1)-0.02*(aaa(2)-aaa(1));

aaa(2)=aaa(2)+0.02*(aaa(2)-aaa(1));

Ylim(aaa);

%plot 2nd

subplot('position',[0.04 0.05 0.94 0.62]);

set(gca,'fontSize',9)

surf(T,F,10*log10(P),'edgecolor','none');

aa1=0; aa2=length(data1)/fs; aa3=F(1,1); aa22=F(size(F)); aa4=aa22(1);

axis ([aa1 aa2 aa3 aa4]);

set(gca,'yscale','log');

colorbar('east');

colormap(jet(4096));

set(gca,'clim',[zminzmax]);

holdon;