function mas_c=LZInWindow(y,window,st,b);
% оконная процедура вычисления меры сложности Леппля-Зива
% Параметры: 
%y - входной ряд
%window- длина окна
%st- количество состояний (система счисления при кодировании)
%    st mожет быть 2,3,или 4
%b- значение нормального уровня прибыльности (для классификации состояний больше 2)
% Алгоритм работы следующий: 
%Для каждого окна:
%1) вычисляются прибыльности
%2) ряд прибыльностей преобразовуются в последовательность битов (в случае st=2) 
%   или большего кол-ва кодов {0,1,2} в случае st=3 и  {0,1,2,3} в случае st=4)
%3) Вычисляется алгоритмическая сложность с помощью процедуры Лемпля-Зива
% Классификация (пункт 2) осуществляется кодированием значения
% волатильности соответствующим кодом. В случае 2-х битов (st=2) 0 -
% отрицательная прибыльность, 1- положительная. В случае 3-х символов,
% уровень b задает порог нормальной прибыльности. Код 0 присваивается, если 
%прибыльность меньше -b (резко отрицательна), код 1 - если прибыльность 
%в диапазоне -b...b, и 2 в случае положительной прибыльности ret>b. 
% В случае 4-х состояний существует 3 границы: {-b,0,b} и соответственно 4
% состояния.
%
% Алгоритмическая сложность вычисляется следующим образом: Алгоритму доступны 
% только% 2 операции: 1) добавления (вставка) нового бита в существующую
% последовательность и 2) копирование уже сформированной
% последовательности. Алгоритмическая сложность - это количество операций
% 1) и 2), необходимых для того, чтобы сформировать заданную
% последовательность. Для случайной последовательности длиной н алгоритмическая
% сложность может быть вычислена за формулой c_rand=n/log{2}(n).
% Относительная алгоритмическая сложность данной последовательности
% вычисляется как отношение полученой сложности к случайной:
% С=с/с_rand, 
% где с - алгоритмическая сложность (кол-во действий 1) и 2)), 
% с_rand - алгоритмическая сложность случайной последовательности. 
n=length(y);
y=(y(2:n)-y(1:n-1))./y(2:n);
y=y/std(y);
n=length(y);
mas_c=[];

for i=1:n-window
    ywind=y(i:i+window);
    c=LZMeasureComplexity(ywind,st,b);
    mas_c=[mas_c,c];
end
figure;
plot(mas_c);