|
|
Про 256 точек. Вам что, обязательно нужно чтобы в спектре было именно 256 точек? А в фактически интересующей, при действительном сигнале, первой половине спектра 128 точек? Никак нельзя иметь 512 во всем и 256 в первой половине? Или 1024 и 512? С окнами, конечно.
Понятно, что при Fdiskr=2*60*100=12 КГц, сигнале 30 Гц, 256 точечной выборке и БПФ спектр будет хреновенький, т.е. основная гармоника размажется. И Вы набираете 400 точечную выборку и, видимо заранее зная вх. частоту, передискретизируете сигнал, предварительно его проинтерполировав. Приводите к 256 точечной выборке, но на другой частоте дискретизации. Или я что-то не понял?
Но при 512, тем более 1024 точечной выборке на той же Fdiskr и окнах будет существенно лучше и без передискретизации. Не похоже, по крайней мере из Вашего объяснения, чтобы нужен был скользящий спектр. Скачущий возможно.
1024 т. Фурье на 218x считается порядка миллисекунды, в зависимости от частоты процессора. Времени, при 12 кгц дискретизации, вагон.
Вполне можно обеспечить скачки в ~12 отсчетов при 1024 т. БПФ.
Или примерно так: при высоких частотах сигнала берете 256 т. выборку,
при низких выборку большего размера, а в спектре, лучше предварительно умноженном на ФНЧ, выбрасываете середину, т.е. область на границе полож. и отрицательных частот, так, чтобы суммарное число оставшихся частот в положит. и отрицательной областях была равна 256.
Ну а если все описанное не подходит, то, в дополнение к сделанному Вами, несколько советов.
1. Посмотрите книжку Марпла, название вроде "Цифровой спектральный анализ". Она где-то в Итернете выложена. Ссылки периодически появляются в этой и микроконтроллерной конф. Я их не знаю, к сожалению. Там разные методы анализа рассмотрены, одно из их назначений - спектр. анализ коротких (относительно) временных рядов, Ваш случай с 256 т. выборкой, 12 КГц дискретизацией и 30 Гц сигналом.
Эти методы не панацея, да и не быстрые, надо сказать. Требуют для хорошей работы максимум информации о сигнале. Но может Вам подойдет что либо. Если найдете ссылку, то там же и другая литература по ЦОС.
В частности Гольденберг,Матюшкин,Поляк. В числе прочего и передискретизация рассматривается.
2. Еще с аппроксимацией (интерполяцией) и передискретизацией. Посмотрите cic фильтры. Начальную информацию можно найти в описание микросхем цифровых приемников, например ad6634,ad6620. Ну и у других фирм, конечно. Плюс поиск в интернете. Но у них есть одно неприятное звено, интегратор, реализуемый как БИХ фильтр. Для Вашей задачи, по крайней мере как я ее понял, он может и не подойти.
E-mail: 
info@telesys.ru
