Анализ спектров радиальной мощности

Question

Во-первых - мне нужно построить этот радиально усредненный спектр 1-D против волнового числа | k |. Но рано я вычислил | k | используя две матрицы, поэтому у меня есть матрица, представляющая | k |. Значит, мне тоже нужно измерить эту матрицу?

Второй - Если я удалю среднее значение из своих данных перед использованием функции FFT2, я ожидал увидеть только разницу в амплитуде значений, но на самом деле я отмечаю разницу в форме радиально усредненного спектра ,

Я пытаюсь использовать это для анализа данных магнитного поля. Спасибо за любую помощь !!

Answer 1

Объяснения второго вопроса следует из аддитивного свойства FFT:

FFT(c) = FFT(a+b) = FFT(a) + FFT(b) 

Если b этого смещения, а именно - среднее значения - что вы удалить из сигнала c, то спектр c-b будет равный вашему первоначальному спектру за вычетом спектра константы, b, но БПФ константы приводит к искренним колебаниям! Таким образом, удалив среднее значение сигнала, вы удаляете сбой sinc на нулевой частоте.