clear all;
EnCoder;%Задание изначальной матрицы
%Работа декодера сводится к следующему. Изначально задается матрица,
%в которой выбирается начальная ветвь, ошибка которой =0, ошибки остальных
%представляются равными настолько большому числу, что должны отпастьв
%процессе отбора. Далее работа разбивается на бабочки, разбивка происходит
%в ручную, как и задание параметров для каждой бабочки. Параметрами
%являются следующие X1,X2 - входные состояния, Х3-Х4 - выходные состояния
%бабочки. A- переход в бабочке между Х1-Х3 И Х2=Х4, B- переход между двумя
%другими. Далее считается суммарная метрика и выбирается минимальная,
%приходящая в узлы Х3 и Х4.
%Если приходящая в узел Х3 метрика из Х1 меньше
%чем из Х2, то берем путь (1-19), по которому до этого пришли в Х1 и
%присваиваем строке (2-20), означающей путь в Х3,а первый элемент
%собственно само значение этого состояния ( задается для всех).
%Использование временных матриц обусловленно тем, что нельзя изменять
%метрику пути, приходящую в Х3 сразу, так как в следующей бабочке эта
%информация будет использоваться для других вычислений

AA_CB=[0,0];  %Вектора переходов
BD_DC=[0,1];
BC_DD=[1,0];
AB_CA=[1,1];

%Задание исходной матрицы путей
    Path_Matrix(1:4,1:20)=0;
    Path_Matrix(1:4)=1:4;
    T_Path(1:4,1:20)=0; %временная матрица
%Задание матрицы ошибок
    Err_Matrix(1:4)=(100);
    Err_Matrix(1)=(0); %задание начальной точки (1)
    T_err(1:4)=0; %временный вектор ошибок

%----Функция выполнения операции с одной бабочкой

function Buttefly(X1,X2,X3,X4,A,B,Y)
    X1_X3= sum(abs(A-Y)); %вычисление ошибок, приходяших в точки
    X1_X4= sum(abs(B-Y));
    X2_X4=X1_X3;
    X2_X3=X1_X4;

    T_Path(1:4,1)=1:4; %В каждый узел приходит один ветвь, значит
            %окончание каждому окончанию пути соответствует одна ветвь
    if (X1_X3<X2_X3) %выбор пути для Х3
        ERR_T(X3)=X1_X3; %Вес пути присваиваем новый
        T_Path(X3,2:20)=Path_Matrix(X1,1:19);
        % элементы со второго получают значение пути, пришедшего в Х1
    else
        ERR_T(X3)=X2_X3;
        T_Path(X3,2:20)=Path_Matrix(X2,1:19);
    end

    if (X1_X4<X2_X4)
        ERR_T(X4)=X1_X4; %Вес пути присваиваем новый
        T_Path(X4,2:20)=Path_Matrix(X1,1:19);
     else
        ERR_T(X4)=X2_X4;
        T_Path(X4,2:20)=Path_Matrix(X2,1:19);
    end
end

function main()
L=15;
D=5;
for i=1:5
y=MagicPeopleVodooPeople(i,L,D);
end

function y=MagicPeopleVodooPeople(i,L,D)
y=i*(L+D)
end

function main()
L=15;
D=5;
for i=1:5
y=MagicPeopleVodooPeople(i);
end

function y=MagicPeopleVodooPeople(i)
y=i*(L+D)
end

function main()
global L D
L=15;
D=5;
for i=1:5
y=MagicPeopleVodooPeople(i);
end

function y=MagicPeopleVodooPeople(i)
global L D
y=i*(L+D)
end

function main()
A=1:3:60;
B=sin(A);
for i=1:2:20
    T(i)=sin(A(i)^2);
   %Еще какой то код....
    M(i)=B(i)^2+A(i);
end
disp(T);
disp(M);

function main()
A=1:3:60;
B=sin(A);
T=sin(A.^2);
   %Еще какой то код....
M=B.^2+A;

disp(T);
disp(M);