clc; clear all; close all;
t=[0, 4, 8, 12, 16, 20, 24];   %ыяес глеяас поу цимомтаи летягсеис хеялойяасиас (ама 4 ыяес)
T=[5, 7, 10, 15, 13, 25, 5];   %есты ои летягсеис хеялойяасиым аисхгтгяа
fprintf('t: ');   fprintf('%4.1f ',t(1:end));
fprintf('\nT: '); fprintf('%4.1f ',T(1:end)); fprintf('\n');

N=length(T);                   %лецехос пимайа хеялойяасиым
%mean(T)
%std(T)

s=0;                           %лгдемислос ахяоислатос
for i=1:N                      %циа окес тис хеялойяасиес
    s=s+T(i);                  %ахяоисла циа еуяесг лесгс тилгс
end
meanT=s/N;                     %еуяесг лесгс тилгс хеялойяасиым

s=0;                           %лгдемислос ахяоислатос
for i=1:N                      %циа окес тис хеялойяасиес
    s=s+(T(i)-meanT)^2;        %ахяоисла циа еуяесг тупийгс апойкисгс хеялойяасиым
end
stdT=sqrt(s/(N-1));            %еуяесг тупийгс апойкисгс хеялойяасиым (Standard Deviation)

fprintf('\nлесг тилг: %4.2f\nтупийг апойкисг: %4.2f\n',meanT,stdT);

%амтийатастасг тилым стис хесеис тоу пимайа опоу упгявам айяаиес тилес

T1=T;                     %дглиоуяциа меоу пимайа поу ха апохгйеусеи то диояхылемо сгла (де хекоуле ма аккоиыхеи о аявийос пимайас T)
for i=1:N                 %циа окес тис тилес тоу пимайа
    if (abs(T(i))>3*stdT) %еуяесг айяаиас тилгс - ам лецакутеяг тоу 3пкасиоу тгс тупийгс апойкисгс
        if(i==1)
            T1(i)=T(i+1); %амтийатастасг ле 2О стоивеио тоу T
        elseif(i==N)
            T1(i)=T(i-1); %амтийатастасг ле пяотекеутаио стоивеио тоу T
        else              %ам хесг айяаиас тилгс диажояг тым айяым тилым тоу пимайа
                          %амтийатастасг тилгс опоу упгяве айяаиа тилг стом T
                          %ле лесг тилг пяогцоулемоу KAI еполемоу аяихлоу
            T1(i)=(T(i-1)+T(i+1))/2;
        end
    end
end

%ейтупысг 2 пимайым (аявийоу йаи диояхылемоу)
fprintf('\nаявийо сгла:     '); fprintf('%4.1f ',T(1:end));
fprintf('\nдиояхылемо сгла: '); fprintf('%4.1f ',T1(1:end));
fprintf('\n');

%сведиасг сглатым аявийоу йаи диояхылемоу
plot(t,T,'rs--','markersize',20);
hold on;
plot(t,T1,'b*--','linewidth',2);
grid on;
legend('аявийо сгла','диояхылемо сгла',2);
xlabel('ыяес');
ylabel('хеялойяасиа (╟C)');
title('ажаияесг айяаиым тилым сглатос');
