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

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

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

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


person Wanderson Ferreira    schedule 04.08.2013    source источник
comment
Возможно, вы захотите опубликовать это на сайте обмена стеками по физике (или математике?). Я не совсем уверен в первом пункте, но может помочь немного больше информации, например |k| не только функция радиуса?   -  person Buck Thorn    schedule 05.08.2013
comment
Более общий вопрос: какую симметрию вы ожидаете (или наблюдаете) в своем изображении? Является ли операция, которая дает вам |k| создание матрицы с такой же симметрией?   -  person Buck Thorn    schedule 05.08.2013


Ответы (1)


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

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

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

person Buck Thorn    schedule 04.08.2013