Помогите мне понять функцию FFT (Matlab)
Я полагаю, что вы можете просто установить переменную окружения так, чтобы она смотрела по нужному пути. Поэтому, если у вас есть сценарий оболочки, который устанавливает некоторые переменные, скажем, setEnv.sh, который вы можете добавить:
export CORE_PEER_FILESYSTEMPATH=<Your ledger's path>
Затем вы можете запустить этот файл и создать его исходную точку, например: . ./setEnv.sh [ 114]
Если вы не хотите использовать сценарии оболочки, вы можете просто сделать это в терминале:
source CORE_PEER_FILESYSTEMPATH=<Your ledger's path>
5 ответов
1) Помимо отрицательных частот, который является минимальной частотой, обеспеченной функцией FFT? Это - нуль?
fft(y) возвращает вектор с 0-th к (N-1)-th образцы DFT y, где y (t) должен думаться, как определено на 0... N-1 (следовательно, 'периодическое повторение' y (t) может считаться периодическим сигналом, определенным по Z).
Первый образец fft(y) соответствует частоте 0. Преобразование Фурье реального, дискретного времени, периодические сигналы имеют также дискретный домен, и это является периодическим и Эрмитовым (см. ниже). Следовательно, преобразование для отрицательных частот является сопряженными из соответствующих образцов для положительных частот.
Например, если Вы интерпретируете (периодическое повторение) y как периодический вещественный сигнал, определенный по Z (выбирающий период == 1), затем домен fft(y) должен быть интерпретирован как N равномерно расположенные точки 0, 2π/N... 2π (N-1)/N. Образцы преобразования при отрицательных частотах-π...-π/N являются спряганием образцов на частотах π... π/N, и равны образцам на частотах π... 2π (N-1)/N.
2) Если это - нуль, как мы выводим нуль на печать на логарифмическом масштабе?
Если Вы хотите потянуть своего рода График Боде, можно вывести преобразование на печать только для положительных частот, игнорируя образцы, соответствующие самым низким частотам (в особенности 0).
3) Результат всегда симметричен? Или это просто, кажется, симметрично?
Это имеет Эрмитовую симметрию если y реально: Ее действительная часть симметрична, ее мнимая часть антисимметрична. Иными словами, его амплитуда симметрична и его антисимметричная фаза.
4) Если я использую брюшной пресс (fft (y)) для сравнения 2 сигналов, я могу потерять некоторую точность?
Если Вы имеете в виду abs(fft(x - y)), это в порядке, и можно использовать его для понимания частотное распределение различия (или ошибка, если x является оценкой y). Если Вы имеете в виду abs(fft(x)) - abs(fft(y)) (???) Вы теряете, по крайней мере, информацию о фазе.
Ну, если Вы хотите понять Быстрое преобразование Фурье, Вы хотите вернуться к основам и понять сам DFT. Но, это не то, что Вы спросили, таким образом, я просто предложу, чтобы Вы сделали это в свободное время :)
Но, в ответе на Ваши вопросы:
- Да, (за исключением отрицательных сторон, поскольку Вы сказали), это - нуль. Диапазон 0 к (N-1) для входа N-точки.
- В MATLAB? Я не уверен, что понимаю, что Ваш вопрос - выводит нулевые значения на печать, поскольку Вы были бы любое другое значение... Хотя, как справедливо указано duffymo, нет никакого естественного журнала нуля.
- Это чрезвычайно подобно sinc (кардинал синуса) функция. Это не обязательно будет симметрично, все же.
- Вы не потеряете точности, у Вас просто будет ответ величины (но я предполагаю, что Вы уже знали это).
Консалтинг "С числовыми Рецептами в C", говорится в Главе 12 по "Быстрому преобразованию Фурье":
Частотные диапазоны от отрицательного ФК до положительного ФК, где ФК является критической частотой Nyquist, которая равна 1 / (2*delta), где дельта является интервалом дискретизации. Таким образом, частоты могут, конечно, быть отрицательными.
Вы не можете вывести на печать что-то, что не существует. Нет никакого естественного журнала нуля. Вы или выведете частоту на печать как ось X или выберете диапазон, который не включает нуль для Вашей оси полужурнала.
Присутствие или отсутствие симметрии в частотном диапазоне зависят от природы функции во временном интервале. У Вас может быть график в частотной области, которая не симметрична об оси y.
Я не думаю, что, принимая абсолютное значение как этот является хорошей идеей. Вы захотите читать гораздо больше о свертке, исправлении и обработке сигналов для сравнения двух сигналов.
- результат fft может быть 0. уже отвеченный другими людьми.
- для графического изображения 0 частот прием должен установить его на очень маленькое положительное число (я использую exp (-15) с этой целью).
- уже отвеченный другими людьми.
- если Вы только интересуетесь величиной, да можно сделать это. это применимо, скажем, во многих проблемах обработки изображений.
Половина Вашего вопроса:
3) Результаты операции FFT зависят от природы сигнала; следовательно нет ничего требующего, что это быть симметричным, хотя, если это - Вы, может получить еще некоторую информацию о свойствах сигнала
4) Это сравнит величины пары сигналов, но те, которые равным, не делают никакой гарантии, что FFTs идентичны (не забывайте о фазе). Это может, однако, быть достаточно в Ваших целях, но необходимо быть уверены в этом.