nabbla (nabbla1) wrote,
nabbla
nabbla1

Category:

В поисках утраченной симметрии

Рассмотрим одномерный сигнал f(x) и его спектр F(w), а также двумерный сигнал f(x,y) и его пространственный спектр F(wx,wy), работаем в комплексных числах.

Как правило, симметрия в пространственной области будет так или иначе проявляться и в частотной:

Если f(x) действительный, то F(-w)=F*(w), т.е спектр в области отрицательных частот является избыточным, его можно восстановить из спектра в области положительных частот.

Если f(x,y) действительный, то F(-wx,-wy)=F*(wx,wy), т.е из 4 квадрантов независимыми будут только два.

Если f(x) обладает четной симметрией, f(-x)=f(x), то четной симметрией будет обладать и спектр.
Если f(x) ещё и действительный, то и спектр будет действительным.

Если f(x) обладает нечетной симметрией, f(-x)=-f(x), то нечетной симметрией будет обладать и спектр.
Если f(x) при этом является действительной функцией, то спектр будет чисто мнимым.

Если f(-x,y)=f(x,y), то F(-wx,wy)=F(wx,wy),
аналогично, если f(x,-y)=f(x,y), то F(wx,-wy)=F(wx,wy).
Наконец, если сигнал является действительным и обладает чётной симметрией по обеим осям, то спектр тоже будет являться действительным и чётным по обеим осям.



Эти (и похожие) свойства симметрии интуитивны и легко запоминаются. Действительно, если мы находим дифракцию на симметричном объекте, мы ожидаем той же симметрии и в картине дифракции, иначе получается абсурд - мы произведем поворот отверстия, переводящий его в себя, картина дифракции повернется и не совпадёт с исходной - противоречие.

Именно с этого и началась когда-то моя работа по расчету дифракции на сотовой бленде. Каково же было моё удивление, когда БПФ от идеально симметричной апертуры лишилось этой симметрии чуть менее чем полностью!

Возьмём предельно упрощенный случай, чтобы уловить суть проблемы, а именно: поле 4х4, посередине квадрат 2х2:

table1.png

Лучшего сложно пожелать:

  • Значения строго действительные

  • Чётная симметрия вокруг горизонтальной оси, сверху и снизу от неё равное количество отсчётов, по 8

  • Чётная симметрия вокруг вертикальной оси, слева и справа от неё равное количество отсчётов, по 8

  • Чётная центральная симметрия, что, в общем-то вытекает из предыдущих двух пунктов



Кажется, что и БПФ от этого изображения должен быть действительным и обладать четной симметрией по двум осям. Ага, щас:
table2.png
Симметрии нет НИ ОДНОЙ!


Но вспомним всё же, что начало отсчета у БПФ лежит в верхем левом углу, и переместим его поближе к центру, сделав циклический сдвиг на 2 (т.е на N/2) по обеим осям, того же самого можно добиться, если исходную "картинку" домножить на шахматную доску (-1)x+y:
table3.png

Что-то всё-таки проглядывает: если выкинуть нулевую строку и столбец и посмотреть только на область 3х3, то она хотя бы обладает свойством комплексной сопряженности, т.е F(2-wx,2-wy) = F*(2+wx,2+wy) - этим свойством должен обладать спектр любого действительного сигнала. Но даже и тут симметрия какая-то покоцанная: вместо поля 4х4 мы рассматриваем 3х3, сдвинутое вбок, и хорошо еще, остальные нули (хоть как-то проявилась симметрия исходного сигнала, а именно, нулевые высшие гармоники по обеим осям), но при произвольном действительном сигнале это вовсе не обязательно!


У меня есть объяснение, что же здесь случилось, в скором времени попробую "с божьей помощью" (а точнее, с помощью двух богов из разных пантеонов, это подсказка!) нарисовать красивую картинку, объясняющую, куда ж запрятались все симметрии и как их вернуть на место. А пока можете предложить свои объяснения, задача сформулировать это предельно "на пальцах" до сих пор остра, в неё упирается любая попытка обосновать преимущества уравновешенного троичного БПФ перед своим бинарным собратом.
Tags: математика, простые вещи, уравновешенное троичное БПФ
Subscribe

  • Так есть ли толк в ковариационной матрице?

    Задался этим вопросом применительно к своему прибору чуть более 2 недель назад. Рыл носом землю с попеременным успехом ( раз, два, три, четыре),…

  • Big Data, чтоб их ... (4)

    Наконец-то стряхнул пыль с компьютерной модели сближения, добавил в неё код, чтобы мы могли определить интересующие нас точки, и выписать…

  • Потёмкинская деревня - 2

    В ноябре 2020 года нужно было сделать скриншот несуществующей программы рабочего места под несуществующий прибор, чтобы добавить его в документацию.…

  • Post a new comment

    Error

    Anonymous comments are disabled in this journal

    default userpic

    Your IP address will be recorded 

  • 15 comments

  • Так есть ли толк в ковариационной матрице?

    Задался этим вопросом применительно к своему прибору чуть более 2 недель назад. Рыл носом землю с попеременным успехом ( раз, два, три, четыре),…

  • Big Data, чтоб их ... (4)

    Наконец-то стряхнул пыль с компьютерной модели сближения, добавил в неё код, чтобы мы могли определить интересующие нас точки, и выписать…

  • Потёмкинская деревня - 2

    В ноябре 2020 года нужно было сделать скриншот несуществующей программы рабочего места под несуществующий прибор, чтобы добавить его в документацию.…