[x,Fs,bits] = wavread('');
x(6837 : 6840)=0;
%Для первого эллемента
d=length(x)
kk=24 %Количество интервалов
k=d/kk %Длина интервала
xx=x(1:d/kk);
shag=0.001; % шаг интегрирования (можно выбрать другой шаг)
x3=shag*cumtrapz(xx);
x2=shag*cumtrapz(x3);
x1=shag*cumtrapz(x2);
j=1;
k=k-1;
for i=1:1:k
    A(i,j)=1;
end;
j=2;
for i=1:1:k
    A(i,j)=x1(i);
end;
j=3;
for i=1:1:k
    A(i,j)=x1(i)^2;
end;
j=4;
for i=1:1:k
    A(i,j)=x1(i)*x2(i);
end;
j=5;
for i=1:1:k
    A(i,j)=x1(i)*x2(i)^2;
end;
j=6;
for i=1:1:k
    A(i,j)=x2(i)*x1(i)^2;
end;
j=7;
for i=1:1:k
    A(i,j)=x2(i)^2;
end;
j=8;
for i=1:1:k
    A(i,j)=x2(i);
end;
j=9;
for i=1:1:k
    A(i,j)=(x1(i)^2)*(x2(i)^2);
end;
cof=(A\xx((1:k))); %поиск коэффициентов

%запись полученных коэффициентов в массив C
for i=1:j
    C(i,1)=cof(i);
end

%Для второго эллемента
d=length(x) %Длина речевого сигнала
kk=24 %Количество интервалов
k=d/kk %Длина интервала
xx=x(d/kk:2*d/kk);
shag=0.001;
x3=shag*cumtrapz(xx);
x2=shag*cumtrapz(x3);
x1=shag*cumtrapz(x2);
j=1;
k=k-1;
for i=1:1:k
    A(i,j)=1;
end;
j=2;
for i=1:1:k
    A(i,j)=x1(i);
end;
j=3;
for i=1:1:k
    A(i,j)=x1(i)^2;
end;
j=4;
for i=1:1:k
    A(i,j)=x1(i)*x2(i);
end;
j=5;
for i=1:1:k
    A(i,j)=x1(i)*x2(i)^2;
end;
j=6;
for i=1:1:k
    A(i,j)=x2(i)*x1(i)^2;
end;
j=7;
for i=1:1:k
    A(i,j)=x2(i)^2;
end;
j=8;
for i=1:1:k
    A(i,j)=x2(i);
end;
j=9;
for i=1:1:k
    A(i,j)=(x1(i)^2)*(x2(i)^2);
end;

cof=(A\xx((1:k))); %поиск коэффициентов
%запись полученных коэффициентов в массив C
for i=1:j
    C(i,2)=cof(i);
end

%Для третьего эллемента
d=length(x)
kk=24
k=d/kk
xx=x(2*d/kk:3*d/kk);
shag=0.001;
x3=shag*cumtrapz(xx);
x2=shag*cumtrapz(x3);
x1=shag*cumtrapz(x2);
j=1;
k=k-1;
for i=1:1:k
    A(i,j)=1;
end;
j=2;
for i=1:1:k
    A(i,j)=x1(i);
end;
j=3;
for i=1:1:k
    A(i,j)=x1(i)^2;
end;
j=4;
for i=1:1:k
    A(i,j)=x1(i)*x2(i);
end;
j=5;
for i=1:1:k
    A(i,j)=x1(i)*x2(i)^2;
end;
j=6;
for i=1:1:k
    A(i,j)=x2(i)*x1(i)^2;
end;
j=7;
for i=1:1:k
    A(i,j)=x2(i)^2;
end;
j=8;
for i=1:1:k
    A(i,j)=x2(i);
end;
j=9;
for i=1:1:k
    A(i,j)=(x1(i)^2)*(x2(i)^2);
end;

cof=(A\xx((1:k))); %поиск коэффициентов

%запись полученных коэффициентов в массив C
for i=1:j
    C(i,3)=cof(i);
end

%Для последуюих элементов
for i = 4:kk
    d=length(x) %Длина речевого сигнала
    kk=24 %Количество интервалов
    k=d/kk %Длина интервала
    xx=x((i-1)*d/kk:i*d/kk);
    shag=0.001;
    x3=shag*cumtrapz(xx);
    x2=shag*cumtrapz(x3);
    x1=shag*cumtrapz(x2);
    j=1;
    k=k-1;
    for i=1:1:k
        A(i,j)=1;
    end;
    j=2;
    for i=1:1:k
        A(i,j)=x1(i);
    end;
    j=3;
    for i=1:1:k
        A(i,j)=x1(i)^2;
    end;
    j=4;
    for i=1:1:k
        A(i,j)=x1(i)*x2(i);
    end;
    j=5;
    for i=1:1:k
        A(i,j)=x1(i)*x2(i)^2;
    end;
    j=6;
    for i=1:1:k
        A(i,j)=x2(i)*x1(i)^2;
    end;
    j=7;
    for i=1:1:k
        A(i,j)=x2(i)^2;
    end;
    j=8;
    for i=1:1:k
        A(i,j)=x2(i);
    end;
    j=9;
    for i=1:1:k
        A(i,j)=(x1(i)^2)*(x2(i)^2);
    end;
    cof=(A\xx((1:k))); %поиск коэффициентов
    %запись полученных коэффициентов в массив C
    for i=1:j fff=4:kk
        C(i,fff)=cof(i);
    end
end