Страница 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?
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 файла - вроде понятно, но как учесть переменный шаг? Спасибо большое за подсказку!
P.S. Хитрый пакет Origin как-то делает FFT 2-столбчатых данных, но как именно - не говорит. А я хочу прогу сделать, чтобы ASCII файлы по клику мышки открывались на просмотр, и тут же определялся период (если есть) понижения тока (дефектов).
Re: При снятии спектра с синусоиды - несущая частота не та
Добавлено: 02 фев 2022, 23:37
Andrey Petrzhik
Чтобы было нагляднее что же я хочу...
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 оператор даёт тот же результат).
Буду рад мудрым советам:)
Re: При снятии спектра с синусоиды - несущая частота не та
Добавлено: 03 фев 2022, 18:25
Andrey Petrzhik
Усп, прошу прощения... пики совпадают, фон только другой + шума больше. Ну и интенсивности пиков просели.
(как-то не аккуратно выделил видать пики в LV)
Осталось чуть-чуть, фильтр какой-нить запихнуть или усреднение?