[an error occurred while processing this directive]
|
такая реализация Дурбина (матлаб-код)?
R(1..) - АКФ, R(1) - энергия.
a(1..) - выходные коэффициенты фильтра.
E = R(1);
for i = 0:NDurbin-1acc=R(i+2);
for j = 0:i-1
acc = acc - a(j+1)*R(i-j-1+2);
endpk=-acc/E;
a(i+1)=-pk;
E = acc*pk+Err;
tmp=a;
for j = 0:i-1
a(j+1) = a(j+1)+pk*tmp(i-j-1+1);
end
end;
E-mail: info@telesys.ru