Мндя... да я не прав. Да и ответ у меня тогда не такой уж и простой будет... Правельно это сделать через расчет полюсов т.д. а там ведь функция бесселя у которой нули ваще какие-то не такие...
Вообщем есть такая идея: [b,a]=besself(5,0.2); % надо нормированную частоту, а то весь матлаб плющит randn('state',0); noise = randn(50000,1); x=filter(b,a,noise); h = lpc(x,15);
Вообщем применям линейное предсказание. Вообщем задачу можно решить быстро и с нужной точностью, но так сказать "в лоб". А и еще: h - это импульсная которую можно и в fdatool загнать
|