Казват,че човек се учел докато е жив.Та сега ми се налага да уча Фурие
преобразувание защото ми се налага да навляза в дебрите на филтрацията,
по-точно да оценявам амплитуда на спектрални съставящи - главоболия и
недоспиване.Може ли някой от колегите в форума да обясни следните неща
защото,нито здрав сън,нито безсънни нощи нито Опенхайм,Шафер и Ненов не
помагат,нито даже бирата и водката
Ситуацията е следната:
1.Имаме задоволително крив сигнал представен със 256 отчета през равни
интервали от време с максимална честота в диапазона 30-70 кHz.
2.Отчетите са взети с честота 1 МHz,което според книгите е достатъчно.
Рзарядността на отчетите е 12 бита.
3.Подредили сме отчетите в масив x[n], n=0....255
4.Та значи на горната поредица ще и се прави Фурие анализ.Стана ясно
че се ползва алгоритъм за бързо преобразувание на Фурие.Като гледам
има 2 подхода - с децимация във времето и децимация в честотната област.
Разликата (според прочетеното) е тази,че в първия случай отчетите на
входния сигнал се разделята на четни-нечетни,подреждат се и се прави
така наречената пеперуда,като изходните отчети са подредени по номера,
докато при втория метод се ползват както са си подредени,но след FFT-то
трябва да се преподредят.По тази точка въпросите са:
-как се започва пеперудата - в смисъл за 256 точки може да се започне с
2 пеперуди по 128 точки,после 4 пеперуди по 64 точки,след това 8 по 32 и
така до края.Другия вариант - 2 128 точкови,които след това се делят на 64х4
точкови и така до края на преобразуванието.Кой от двата варианта е правилен
и как по дяволите се изчислява прословутия завъртащ множител в зависимост
от кой етап се намираме ?
- след това получаваме един масив от комплексни числа с размерност 256.
В него се стойностите Х(n).За за получим амплитудата на дадена спектрална
съставяща трябва да се работи с реални и имагинерни части(което не ми
харесва) а именно |X(n)|^2=X(n)re)^2+(X(n)im)^2.Това не може ли да се
избегне със използването на конволюция?А именно X(n).X(n)*
Ако номера с конволюцията е възможен не ми стана ясно какво се получава
ако преобразуванието е с децимация в времевата област?Ако е с децимация
в честотната област то тогава в масива на стойност n=0 трябва да е
амплитудата на основната честота,на n=1 - първи хармоник или бъркам?
Надявам се някой от колегите да отговори,за да мога да се измъкна от
кашата дето съм си я сътворил.Под ръка спектрален анализатор нямам,нито
осцилоскоп със възможност за Фурие анализ.Ако на някой също му е попадала
някоя хубава книга с картинки и блок схема на алгоритъм как се прави - да
сподели щото в книгите които имам нещо не мога да разгадая кое как е а
в нета сорсове те са предимно на Фортран.На С са малко и няма обяснения
кое какво е и не мога да се ориентирам.
).
