|
|
b =pi*alfa -произведение длительности на полосу частот. Обычно задают к-т alfa, начиная с 2. Он не обязательно целый.
tau - половина длины окна.
Учтите, t в формуле меняется от -tau до tau.
Есть приближения.
Приближения формулой:
W(n)=a0+a1*cos(fi(n))+a2*cos(2*fi(n))+a3*cos(3*fi(n)).
fi(n)=2pi*(n-(N-1)/2)/N.
a0=b0/c,a1=2*b1/c,a2=2*b2/c,a3=2*b3/c
c=b0+2*b1+2*b2+2*b3.
b[i]=[sh(pi*sqrt(alfa^2-i^2))]/sqrt(alfa^2-i^2)
alfa (см. начало) от 2 до 4. i - номер к-та b, целый и меньше равен alfa.
Все это описано в статье Харриса -Использование окон при гарм. анализе методом ДПФ.
Я ее (английский оригинал) несколько раз выкладывал на webfile.ru
А страницу или 2-3 назад на эту статью давал ссылку SM. Поищите.
Но прежде чем искать - акуратно сгенерируйте и подайте на вход Вашему 32-х р. Фурье 32-х разрядный синус частоты кратной 2pi/N.
Без окон. Все квантовательные и округлительные (если все правильно) должны быть (в разах) на уровне ~2^(-30).
|
|
E-mail: 
info@telesys.ru
