function [mjd_mean,y_mean,y_std,y_n]=mean_of_day(mjd,y);
%
% Вычисление среднего за сутки значениЯ компоненты полЯ
%        [mjd_mean,y_mean,y_std,y_n]=mean_of_day(mjd,y);
%
% mjd - модифицированнаЯ юлианскаЯ дата
%   y - измеренные значениЯ

kd=find(diff(fix(mjd))>0);    % ищем индексы смены суток
kd=kd(:); kd=[kd;length(mjd)];% добавлЯем последние сутки
n=length(kd);            % кол-во суток с измерениЯми

mjd_mean=zeros(n,1); y_mean=zeros(n,1); y_std=zeros(n,1); y_n=zeros(n,1);
for k=1:n,
    kk=fix(mjd)==fix(mjd(kd(k)));
    mjd_mean(k)=mean(mjd(kk));        % mjd длЯ вычислениЯ статистики
    y_mean(k)=mean(y(kk));
    y_std(k)=std(y(kk));
    y_n(k)=sum((kk));
end
