UR1=4;  nR1=-1;
UR2=1;  nR2=-1;
for i=1:n
    if nR1<0 
        if result(i)>UR1  nR1=i; end
    else
        if nR2<0 && result(i)<UR2; nR2=i; end
    end
end
plot(t,oscil,t,result, t(nR1), oscil(nR1),'bo',t(nR2), oscil(nR2),'ro');

clear all;clc;close all;
%t=0:0.01:5*pi;
%oscil=10*sin(5*t);
load('oscil.mat', 'oscil');
n=length(oscil);
t=1:n;
for ii=1:230
    oscil(ii)=oscil(ii)*0.3;
end
for ii=n-130:n
    oscil(ii)=oscil(ii)*0.05;
end
window = 125;
if(mod(window,2)==0)
window=window+1;
end
hw=(window-1)/2; %размах окна влево и вправо от текущей позиции
 
n=length(oscil);
result=zeros(n,1);
result(1)=oscil(1); %первый элемент берем из исходного массива oscil как есть
 
for i=2:n %организовываем цикл по числу элементов
init_sum = 0;
if(i<=hw) %если индекс меньше половины окна, мы находимся в начале массива, 
%нужно брать окно меньшего размера
k1=1; %в качестве начала окна берем первый элемент
k2=2*i-1; %конец окна 
z=k2; %текущий размер окна
elseif (i+hw>n) %если индекс+половина окна больше n - мы приближаемся к концу массива и размер окна
    %также нужно уменьшать
    k1=i-n+i; %начало окна 
    k2=n; %конец окна - последний элемент массива
    z=k2-k1; %размер окна
else %если первые два условия не выполняются, мы в середине массива
    k1=i-hw;
    k2=i+hw;
    z=window;
end
 
for j=k1:k2 %организуем цикл от начала до конца окна
init_sum=init_sum+abs(oscil(j)); %складываем все элементы
 
end
result(i)=init_sum/(z); %и делим на текущий размер окна
end
UR1=2.5;  nR1=-1;
UR2=0.04;  nR2=-1;
UR3=0.04;  nR3=-1;
for i=1:n
    if nR1<0 
        if result(i)>UR1  nR1=i; end
    else
        if nR2<0 && result(i)<UR2  nR2=i; end
    end
end
for i=n-1:-1:1
    if nR3<0
        XX=result(i);
        if result(i)>UR3  nR3=i; end;
    end
end
%plot(result);
plot(t,oscil,t,result, t(nR1), oscil(nR1),'bo',t(nR2), oscil(nR2),'ro',t(nR3), oscil(nR3),'ro');