Страница 1 из 1

При снятии спектра с синусоиды - несущая частота не та

Добавлено: 05 дек 2019, 21:59
dakishi
Добрый вечер!
Можете подсказать в чём здесь дело ? Я создаю сигнал (синусоиду с частотой 0.01 Hz - несущей частотой, например - см. рисунок), разлагаю в Фурье спектр при помощи FFT Power Spectrum and PSD и основная мощность в спектре, в итоге, приходится на 10 Hz (max index(ex) на рисунке) - как это возможно ? С использованием инструмента из той же палитры, выделяющего мнимую и реальную часть (второй рисунок) получается и того хуже - шум.

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 06 дек 2019, 01:08
Kosist
Функция возвращает кластер из 3-х составляющих: начальную частоту f0, шаг частоты df и мощность. Поэтому найдя индекс максимальной мощности N, соответсвующую частоту Вы найдете как f = f0 + N*df. Вы так определили, что это 10 Hz?
Frequency.png

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 10 дек 2019, 20:27
dakishi
Kosist писал(а):Функция возвращает кластер из 3-х составляющих: начальную частоту f0, шаг частоты df и мощность. Поэтому найдя индекс максимальной мощности N, соответсвующую частоту Вы найдете как f = f0 + N*df. Вы так определили, что это 10 Hz?
Работает, но когда частота выше порога, идет неправильный расчёт. В данном случае реальная частота 3, а по FFT 4, причина, очевидно, в шуме на FFT - можете обьяснить, почему это происходит ?

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 10 дек 2019, 21:08
zxc_pavel
Он вообще отрицательный у Вас, это не спектр, а действ. часть комплексного представления.
Поставьте или магнитуду или корень из суммы квадратов re и im, как на картинке

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 02 фев 2022, 23:11
Andrey Petrzhik
Добрый день.

А подскажите пожалуйста что делать если данные для FFT - это не сигнал, т.е. условный синус с шумами, а набор точек X-Y двумя столбиками, где X - это координата, а Y - измеряемая величина? Причём шаг по X немного плавает.

Т.е. как загнать для FFT функцию с неким постоянным шагом dX из 1 колонки ASCII файла - вроде понятно, но как учесть переменный шаг? Спасибо большое за подсказку!
fft-magnitude.jpg
P.S. Хитрый пакет Origin как-то делает FFT 2-столбчатых данных, но как именно - не говорит. А я хочу прогу сделать, чтобы ASCII файлы по клику мышки открывались на просмотр, и тут же определялся период (если есть) понижения тока (дефектов).

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 02 фев 2022, 23:37
Andrey Petrzhik
Чтобы было нагляднее что же я хочу...
исходник.jpg
fft-magnitude_origin.jpg
финал.jpg
ASCII файлик тоже приложил (в rar, а то великоват), вдруг у кого будет время поиграться...

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 03 фев 2022, 10:03
zxc_pavel
Можно сплайном делать выравнивание по шагу

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 03 фев 2022, 11:10
Andrey Petrzhik
О! Спасибо большое! А можно для полного счастья в Labview 2012 сохранить?
Да, и ещё небольшой шаг... а как мне 1/f построить по оси x?
Если бы строился XY graph, я просто массив X-ов беру и 1/x делаю, а тут как?
К сожалению пока не вкурил чё там в кластере, может всё и очевидно конечно...

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 03 фев 2022, 11:28
zxc_pavel
Там в кластере, который на waveform graph идет 3 элемента: x0, dx, массив игреков. Конечно, можно сделать 1/dx.
Сохранил в 12

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 03 фев 2022, 13:29
Andrey Petrzhik
Что-то не придумал как покрасивее сделать 1/f из f... но вроде правильно должно быть.
Картинка уже близка к нужной, но периоды немного "уплыли", см. например пик 77.30, выделенный на финальном графике, да и видны все пики менее явно чем в Origin.
Где-то слыхал, что магнитуда - это не амплитуда, но в labview не пойму что лучше выбрать (Real-Im оператор даёт тот же результат).
Буду рад мудрым советам:)
1f.jpg
spline_res.jpg

Re: При снятии спектра с синусоиды - несущая частота не та

Добавлено: 03 фев 2022, 18:25
Andrey Petrzhik
Усп, прошу прощения... пики совпадают, фон только другой + шума больше. Ну и интенсивности пиков просели.
(как-то не аккуратно выделил видать пики в LV)
LVvsOrigin.jpg
Осталось чуть-чуть, фильтр какой-нить запихнуть или усреднение?