|
|
  |
fft в матлабе дает отрицательные значения для реальной части |
|
|
|
Dec 15 2011, 11:08
|
Частый гость
 
Группа: Участник
Сообщений: 136
Регистрация: 13-05-11
Пользователь №: 64 972

|
Цитата(petrov @ Dec 15 2011, 13:47)  Откуда следует что реальная часть должна быть положительной? Синус может иметь любую фазу, нарисуйте произвольный вектор на комплексной плоскости, он в том числе может иметь и отрицательную действительную часть. Вообще да, согласен. Только тогда остается главная проблема: почему на графике реальной части на месте частоты синуса кривая уходит то в минус, то в плюс, по сути получаются два пика вплотную друг к другу, отрицательный и положительный, один побольше, другой поменьше. У меня основная проблема эту хрень убрать В прикрепленном файле картинки для реальной и мнимой части, которые у меня получаются
Сообщение отредактировал RoadRunner - Dec 15 2011, 11:19
Эскизы прикрепленных изображений
|
|
|
|
|
Dec 15 2011, 11:18
|

Участник

Группа: Участник
Сообщений: 20
Регистрация: 28-02-06
Из: Минск
Пользователь №: 14 809

|
Цитата(RoadRunner @ Dec 15 2011, 14:08)  Вообще да, согласен. Только тогда остается главная проблема: почему на графике реальной части на месте частоты синуса кривая уходит то в минус, то в плюс, по сути получаются два пика вплотную друг к другу, отрицательный и положительный, один побольше, другой поменьше. У меня основная проблема эту хрень убрать Для того, чтобы не наблюдалось таких "соседних" пиков (большой мощности) необходимо, чтобы все частоты были кратны сетке FFT. Попробуйте подавать синус, чтобы поместилось 5, 10 или 20 полных периодов на длине FFT.
|
|
|
|
|
Dec 15 2011, 11:23
|
Частый гость
 
Группа: Участник
Сообщений: 136
Регистрация: 13-05-11
Пользователь №: 64 972

|
Цитата(Wasja @ Dec 15 2011, 14:18)  Для того, чтобы не наблюдалось таких "соседних" пиков (большой мощности) необходимо, чтобы все частоты были кратны сетке FFT. Попробуйте подавать синус, чтобы поместилось 5, 10 или 20 полных периодов на длине FFT. Этого и боялся..  Мне потом надо будет сигнал с АЦП туда загонять, неужто сидеть периоды выпиливать. Плохо пока представляю себе как это делать буду...
|
|
|
|
|
Dec 15 2011, 11:48
|

Участник

Группа: Участник
Сообщений: 20
Регистрация: 28-02-06
Из: Минск
Пользователь №: 14 809

|
Цитата(RoadRunner @ Dec 15 2011, 14:23)  Этого и боялся..  Мне потом надо будет сигнал с АЦП туда загонять, неужто сидеть периоды выпиливать. Плохо пока представляю себе как это делать буду... Зачем выпиливать (да и не получится это)? Все некратные частоты будут "размазыватся" между соседними спектральными линиями. А их обработка зависит от поставленной задачи...
|
|
|
|
|
Dec 15 2011, 12:56
|
Частый гость
 
Группа: Участник
Сообщений: 136
Регистрация: 13-05-11
Пользователь №: 64 972

|
Цитата(ivan219 @ Dec 15 2011, 15:12)  RoadRunner Могу ошибиться но по моему на картинке вы представили отдельно вещественную и мнимую часть спектра. Посчитайте модуль Sqrt(Sqr(Re) + Sqr(Im)) / N будет реальный спектр. угу, так и есть, вещественную я реальной называю от англ. "real" - буржуйский влияет  думается окно надобно накладывать, наверное уменьшит проблему Ханнинг помог)) В общем то и подбор целого числа периодов закономерно помог, спектр перестал сильно расплываться. Но когда окно наложишь, уже по барабану сколько периодов - и так хорошо получается. Большое спасибо за советы!
Сообщение отредактировал RoadRunner - Dec 15 2011, 12:51
Эскизы прикрепленных изображений
|
|
|
|
|
  |
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|