% статистика по выбросам

tmp=load('train_spikes_2011.dat');
T11_1=datenum(tmp(:,1:6)); % момент переднего фронта
T11_2=datenum(tmp(:,1),tmp(:,2),tmp(:,3),tmp(:,7),tmp(:,8),tmp(:,9)); % момент заднего фронта
A11_1=tmp(:,10); A11_2=tmp(:,11); % передний и задний фронт

% средние для импульса
A11=mean([A11_1 A11_2]')';
T11=mean([T11_1 T11_2]')';
LT11=rem(T11+7/24,1)*24; % местное время в часах

tmp=load('train_spikes_2022.dat');
T22_1=datenum(tmp(:,1:6)); % момент переднего фронта
T22_2=datenum(tmp(:,1),tmp(:,2),tmp(:,3),tmp(:,7),tmp(:,8),tmp(:,9)); % момент заднего фронта
A22_1=tmp(:,10); A22_2=tmp(:,11); % передний и задний фронт

% средние для импульса
A22=mean([A22_1 A22_2]')';
T22=mean([T22_1 T22_2]')';
% оставляем 22 дня как для 2011
k=T22<datenum(2022,12,21);
A22=A22(k); T22=T22(k);
LT22=rem(T22+7/24,1)*24; % местное время в часах

figure(1)
VA=-10:0.5:10;
h11=histogram(A11,VA); hold on; h22=histogram(A22,VA); hold off
xlabel('Амплитуда выбросов, нТл')
ylabel('Количество выбросов')
legend('2011','2022')

figure(2)
subplot(211), plot(LT11,abs(A11),'o',LT22,abs(A22),'*','MarkerSize',3)
ylabel('Амплитуда выбросов, нТл'), legend('2011','2022')
set (gca,'Xlim',[0 24],'XTick',0:4:24)
VT=[0.5:0.5:23.5]';
N11=hist(LT11,VT)'; N22=hist(LT22,VT)';
subplot(212), bar(VT,[N11 N22])
set (gca,'Xlim',[0 24],'XTick',0:4:24)
ylabel('Количество выбросов'), legend('2011','2022')
xlabel('Местное время, час')

