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

Простейшие вопросы в области инженерной разработки
Ответить
dakishi
assistant
assistant
Сообщения: 120
Зарегистрирован: 05 сен 2019, 21:01
Версия LabVIEW: 2019
Контактная информация:

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

Сообщение dakishi »

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

Activity Gold
expert
expert
Сообщения: 1236
Зарегистрирован: 21 фев 2011, 23:44
Награды: 2
Версия LabVIEW: 2013-2020
Благодарил (а): 23 раза
Поблагодарили: 30 раз
Контактная информация:

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

Сообщение Kosist »

Функция возвращает кластер из 3-х составляющих: начальную частоту f0, шаг частоты df и мощность. Поэтому найдя индекс максимальной мощности N, соответсвующую частоту Вы найдете как f = f0 + N*df. Вы так определили, что это 10 Hz?
Frequency.png
Мы делили апельсин - много наших полегло...
dakishi
assistant
assistant
Сообщения: 120
Зарегистрирован: 05 сен 2019, 21:01
Версия LabVIEW: 2019
Контактная информация:

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

Сообщение dakishi »

Kosist писал(а):Функция возвращает кластер из 3-х составляющих: начальную частоту f0, шаг частоты df и мощность. Поэтому найдя индекс максимальной мощности N, соответсвующую частоту Вы найдете как f = f0 + N*df. Вы так определили, что это 10 Hz?
Работает, но когда частота выше порога, идет неправильный расчёт. В данном случае реальная частота 3, а по FFT 4, причина, очевидно, в шуме на FFT - можете обьяснить, почему это происходит ?
Вложения
Problem 10.1.JPG
Аватара пользователя
zxc_pavel
adviser
adviser
Сообщения: 242
Зарегистрирован: 02 июн 2016, 12:07
Версия LabVIEW: 15
Благодарил (а): 9 раз
Поблагодарили: 23 раза
Контактная информация:

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

Сообщение zxc_pavel »

Он вообще отрицательный у Вас, это не спектр, а действ. часть комплексного представления.
Поставьте или магнитуду или корень из суммы квадратов re и im, как на картинке
Вложения
1012.png
Andrey Petrzhik
user
user
Сообщения: 81
Зарегистрирован: 17 ноя 2015, 11:57
Версия LabVIEW: 12
Благодарил (а): 4 раза
Контактная информация:

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

Сообщение Andrey Petrzhik »

Добрый день.

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

Т.е. как загнать для FFT функцию с неким постоянным шагом dX из 1 колонки ASCII файла - вроде понятно, но как учесть переменный шаг? Спасибо большое за подсказку!
fft-magnitude.jpg
P.S. Хитрый пакет Origin как-то делает FFT 2-столбчатых данных, но как именно - не говорит. А я хочу прогу сделать, чтобы ASCII файлы по клику мышки открывались на просмотр, и тут же определялся период (если есть) понижения тока (дефектов).
Andrey Petrzhik
user
user
Сообщения: 81
Зарегистрирован: 17 ноя 2015, 11:57
Версия LabVIEW: 12
Благодарил (а): 4 раза
Контактная информация:

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

Сообщение Andrey Petrzhik »

Чтобы было нагляднее что же я хочу...
исходник.jpg
fft-magnitude_origin.jpg
финал.jpg
ASCII файлик тоже приложил (в rar, а то великоват), вдруг у кого будет время поиграться...
Вложения
#2577_13_322_Ag05__NC.rar
(297.15 КБ) 34 скачивания
Аватара пользователя
zxc_pavel
adviser
adviser
Сообщения: 242
Зарегистрирован: 02 июн 2016, 12:07
Версия LabVIEW: 15
Благодарил (а): 9 раз
Поблагодарили: 23 раза
Контактная информация:

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

Сообщение zxc_pavel »

Можно сплайном делать выравнивание по шагу
Вложения
0302.png
0302-2.png
0302.vi
(30.88 КБ) 37 скачиваний
Andrey Petrzhik
user
user
Сообщения: 81
Зарегистрирован: 17 ноя 2015, 11:57
Версия LabVIEW: 12
Благодарил (а): 4 раза
Контактная информация:

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

Сообщение Andrey Petrzhik »

О! Спасибо большое! А можно для полного счастья в Labview 2012 сохранить?
Да, и ещё небольшой шаг... а как мне 1/f построить по оси x?
Если бы строился XY graph, я просто массив X-ов беру и 1/x делаю, а тут как?
К сожалению пока не вкурил чё там в кластере, может всё и очевидно конечно...
Аватара пользователя
zxc_pavel
adviser
adviser
Сообщения: 242
Зарегистрирован: 02 июн 2016, 12:07
Версия LabVIEW: 15
Благодарил (а): 9 раз
Поблагодарили: 23 раза
Контактная информация:

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

Сообщение zxc_pavel »

Там в кластере, который на waveform graph идет 3 элемента: x0, dx, массив игреков. Конечно, можно сделать 1/dx.
Сохранил в 12
Вложения
0302-v12.vi
(22.74 КБ) 35 скачиваний
Andrey Petrzhik
user
user
Сообщения: 81
Зарегистрирован: 17 ноя 2015, 11:57
Версия LabVIEW: 12
Благодарил (а): 4 раза
Контактная информация:

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

Сообщение Andrey Petrzhik »

Что-то не придумал как покрасивее сделать 1/f из f... но вроде правильно должно быть.
Картинка уже близка к нужной, но периоды немного "уплыли", см. например пик 77.30, выделенный на финальном графике, да и видны все пики менее явно чем в Origin.
Где-то слыхал, что магнитуда - это не амплитуда, но в labview не пойму что лучше выбрать (Real-Im оператор даёт тот же результат).
Буду рад мудрым советам:)
1f.jpg
spline_res.jpg
Вложения
0302-v12AMP.vi
(34.25 КБ) 38 скачиваний
Andrey Petrzhik
user
user
Сообщения: 81
Зарегистрирован: 17 ноя 2015, 11:57
Версия LabVIEW: 12
Благодарил (а): 4 раза
Контактная информация:

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

Сообщение Andrey Petrzhik »

Усп, прошу прощения... пики совпадают, фон только другой + шума больше. Ну и интенсивности пиков просели.
(как-то не аккуратно выделил видать пики в LV)
LVvsOrigin.jpg
Осталось чуть-чуть, фильтр какой-нить запихнуть или усреднение?
Ответить
  • Похожие темы
    Ответы
    Просмотры
    Последнее сообщение

Вернуться в «Для чайников»