МНОГОМЕРНЫЕ СИГНАЛЫ И СИСТЕМЫ

Содержание: Введение. 18.1. Двумерные и многомерные сигналы. 18.2. Двумерные системы. Импульсный отклик систем. Двумерная свертка. 18.3. Частотные характеристики сигналов и систем. Частотный отклик системы. Импульсный отклик системы. Свойства двумерного преобразования Фурье. 18.4. Дискретизация двумерных сигналов. Прямоугольный растр дискретизации. Дискретные преобразования Фурье. Интерполяционный ряд восстановления двумерного сигнала. Произвольный растр дискретизации. Интерполяция дискретных сигналов. Прямоугольный и гексагональный растры дискретизации. 18.5. Частотный анализ многомерных сигналов. Периодические последовательности. Конечные последовательности. Многомерные последовательности. Литература.

18.1. Двумерные и многомерные сигналы [9].

Понятие многомерного сигнала. Многомерные сигналы представляют собой функции P независимых переменных при P>1. В общем случае, сигнал может быть непрерывным, дискретным или смешанным. Понятия непрерывности и дискретности аналогичны одномерным сигналам. Что касается смешанного сигнала, то это многомерный сигнал, который описывается функцией некоторого количества непрерывных и некоторого количества дискретных переменных. Пример смешанного двумерного сигнала: ансамбль непрерывных сигналов, изменяющихся во времени (t - вторая переменная), снимаемых с набора сейсмических приемников сейсмотрассы (номера датчиков - первая переменная).

Двумерный дискретный сигнал (цифровой массив) - это функция, определенная на совокупности пар числовых значений координат с определенным шагом дискретизации Dx и Dy. В общем случае, при различной физической размерности аргументов x и y, значения Dx и Dy не равны друг другу:

Элемент последовательности s_n,m представляет собой отсчет двумерной функции s в координатной точке (x=nDx,y=mDy), где значения x и y – независимые переменные (аргументы) функции. Для числовых массивов значения шага дискретизации по аргументам также могут приниматься равными 1 (независимо от размерности) и использоваться аргументация s(n,m) є s_n,m. Результаты геофизических съемок какого-либо одного геофизического параметра по поверхности земли относятся к двумерным функциям: дискретным - если это отсчеты в отдельных точках по определенной координатной сети (x,y), или смешанным - если это непрерывная регистрация данных по профилям (например - мощности экспозиционной дозы гамма излучения горных пород при аэросъемке). Но в настоящее время геофизические съемки относятся даже не к двумерным, а к многомерным функциям, так как регистрируется, как правило, сразу несколько физических параметров геологических сред. Так, например, при спектрометрической съемке естественной радиоактивности горных пород регистрируется содержание в горных породах урана, тория и калия, в гравиразведке - трехкоординатный вектор силы тяжести, и т.п. Если на какой-либо площади проведена съемка нескольких видов геофизики, то их результаты также могут рассматриваться в совокупности, как многомерная функция физических параметров данной геологической среды.

По определениям (18.1.1) двумерные функции и сигналы, равно как и многомерные, имеют бесконечную протяженность по координатам. На практике мы всегда имеем дело с конечными координатами наших данных. Учитывая это, будем считать, что значения наших сигналов за пределами определенных координат равны нулю.

Отметим некоторые двумерные последовательности (функции, сигналы), имеющие специальные названия.

Двумерный единичный импульс d(nDx,mDy) = d_n,m или единичный отсчет:

d_n,m = 1, при n = m = 0.

= 0, при n0, m0.

d_n,m = d_n d_m,

где d_n, d_m - одномерные единичные импульсы (импульсы Кронекера) по координатам n и m. Стилизованное графическое представление двумерного единичного импульса приведено на рис. 18.1.2.

Произвольное расположение двумерного единичного импульса по координатам n₁, m₁ соответственно записывается в виде: d((n-n₁)Dx,(m-m₁)Dy) = d_n-n1,m-m1. Попутно напомним, что математическая запись импульса Кронекера обозначает не единичный отсчет, а функцию, определяющую место положения единичного отсчета и нулевые значения по остальным координатам (аргументам).

Двумерный линейный импульс представляет собой последовательность единичных отсчетов по одной координате: s(n,m) = d(n) или s(n,m) = d(m). На рис. 18.1.3 приведены два двумерных линейных импульса, первый - по координате m = 0: s(n,m) = d(m), и второй импульс по координате n = 2: s(n,m) = d(n-2). Очевидно, что для P-мерных случаев точно таким же образом могут быть определены P-мерные единичные импульсы, P-мерные линейные импульсы, P-мерные площадные импульсы и т.д., хотя понятие импульса, заимствованное из теории одномерных сигналов, здесь несколько не к месту.

Двумерная единичная ступенька u(n,m), представленная на рис. 18.1.4, определяется выражением:

u(n,m) = 1, при n0 и m0,

u(n,m) = u(n) u(m),

Экспоненциальная последовательность: s(n,m) = aⁿb^m, -< n,m <, где а и b в общем случае комплексные числа. При а = exp(jw₁), b = exp(jw₂), |а|=1, |b|=1:

s(n,m) = exp(jnw₁+jmw₂) = cos(nw₁+mw₂)+jsin(nw₁+mw₂).

Разделимые последовательности. Разделимой называют последовательность, которую можно представить в виде произведения одномерных последовательностей. Так, для двумерной разделимой последовательности:

s(n,m) = s(n) s(m).

s(n,m) =s_{iЧ n}(n) s_{iЧ m}(m),

s(n,m) = s(n,i) d(m-i). (18.1.2)

Конечные последовательности. Важным классом сигналов являются последовательности конечной протяженности, для которых сигнал равен нулю вне определенной области, называемой опорной областью сигнала. На рис.18.1.5 условно представлена двумерная последовательность конечной протяженности, значения которой отличны от нулевых только внутри ограниченной прямоугольной области -3n2, -2m. Опорная область сигнала может быть произвольной формы и выходить за пределы сигнала, частично включая нулевые отсчеты. Отсчеты за пределами опорной области считаются равными нулю.

Периодические последовательности. Двумерные последовательности могут быть периодическими, регулярно повторяющимися в пространстве. Последовательность, удовлетворяющая условиям:

обладает периодичностью в двух направлениях, по n и по m. Значения М и N называют интервалами периодичности сигнала соответственно по координатам m и n (горизонтальными и вертикальными интервалами периодичности). Прямоугольная форма области периода (пример на рис.18.1.6) наиболее удобна при обработке данных, но не является единственно возможной.

Для двумерных последовательностей условия (18.1.3) могут рассматриваться как частный случай общих условий периодичности:

Системы осуществляют преобразование сигналов. Формализованная система - это оператор (операция) отображения входного сигнала на выходной: z(x,y) = Т[s(x,y)].

Используя базовые операции, любую двумерную последовательность можно разложить на сумму взвешенных двумерных единичных импульсов:

Обобщением скалярного умножения является пространственное маскирование:

Кроме линейных операций в системах используются также безынерционные нелинейные преобразования с независимым нелинейным воздействием на значения отсчетов входной последовательности. Пример операции - возведение в квадрат:

Линейные системы. Система считается линейной при выполнении двух условий:

1. Пропорциональное изменение входного сигнала вызывает пропорциональное изменение выходного сигнала.

2. Суммарный сигнал двух входных последовательностей дает суммарный сигнал двух соответствующих выходных последовательностей.

В выражении (18.2.1) значения s(i,j) можно рассматривать как скалярные множители для соответствующих единичных импульсов. Применяя оператор преобразования Т[.] к левой и правой части (18.2.1), получаем:

Инвариантность к сдвигу. Система инвариантна к сдвигу, если сдвиг входной последовательности приводит к такому же сдвигу выходной последовательности:

Т[s(n-N,m-M)] = z(n-N,m-M).

Линейность и инвариантность к сдвигу являются независимыми свойствами системы. Так, пространственное маскирование линейно, но не инвариантно к сдвигу, а безынерционные операторы нелинейны, но инвариантны к сдвигу.

В дальнейшем ограничимся рассмотрением только систем, широко распространенных при решении практических задач - линейных и инвариантных к сдвигу (ЛИС-системы).

Импульсный отклик на произвольно расположенный входной импульс, как следует из выражения (18.2.3), описывается выражением:

Для частного случая i = j = 0 имеем:

Используя принцип инвариантности к сдвигу, получим:

т.e. импульсный отклик на произвольно расположенный входной импульс равен сдвинутому импульсному отклику на входной импульс, расположенный в начале координат.

Двумерная дискретная свертка (18.2.5), является аналогом одномерной дискретной свертки. При замене переменных n-i = k, m-j = l, получим:

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

При обобщении этого выражения на многомерные системы, в векторной форме:

Разделимые системы. Если импульсный отклик системы может быть разделен:

z(n,m) = S_k h(k) S_l h(l) s(n-k,m-l), (18.2.7)

или: z(n,m) = S_k h(k) g(n-k,m), g(n-k,m) = S_l h(l) s(n,m-l).

Аналогичные разделимые системы могут существовать и в многомерном варианте.

Устойчивость системы. Интерес для практики представляют только устойчивые системы, обеспечивающие определенный конечный результат системной операции на конечные входные сигналы. Необходимым и достаточным условием устойчивости системы является абсолютная суммируемость ее импульсного отклика: S_k S_l |h(k,l)| < .

Специальные двумерные системы. На практике используются также системы с несколькими входами и/или выходами.

Допустим, система имеет i-входы и j-выходы, линейна и инвариантна к сдвигу по переменной t. Если на i-вход системы поступает одномерный единичный импульс d_i(t) при нулевых сигналах на остальных входах, то j-выходные сигналы будут импульсным откликом системы h_ij(t). При известном полном ансамбле значений h_ij для всех i-входов, для произвольной комбинации входных сигналов s_i(t) сигнал на j-выходе будет определяться выражением:

z_j(t) = S_i S_k h_ij(k) s_i(t-k). (18.2.8)

Частотный отклик системы. Допустим, что двумерная ЛИС-система имеет импульсный отклик h(kDx,lDy). Подадим на вход системы сигнал вида комплексной синусоиды:

s(n,m) = exp(jnDxЧ w_x+jmDyЧ w_y),

где w_x и w_y – значения частоты сигнала соответственно по координатам x и y. Принимая Dx = 1, Dy = 1 и выполняя двумерную свертку (18.2.5), получаем:

z(n,m) =h(k,l) exp[jЧ w_xЧ (n-k)+jЧ w_yЧ (m-l)] =

= exp(jnЧ w_x+jmЧ w_y) h(k,l) exp(-jkЧ w_x-jlЧ w_y) = H(w_x,w_y) exp(jnЧ w_x+jmЧ w_y).

Таким образом, выходной сигнал представляет собой комплексную синусоиду с теми же значениями частоты, что и у входного сигнала, с изменением амплитуды и фазы за счет комплексного множителя H(w_x,w_y), который носит название частотного отклика (частотной характеристики) системы. Для дискретных сигналов частотный отклик периодичен с периодом 2p по обеим частотным переменным:

H(w_x+2pk,w_y+2pl) = H(w_x,w_y).

Пример расчета частотного отклика системы.

H(w_x,w_y) =h(n,m)exp(-jnw_x-jmw_y) = 0.25+0.125[exp(-jw_x)+exp(jw_x)+exp(-w_y)+exp(jw_y)]+

+0.0625[exp(-jw_x-jw_y)+exp(-jw_x+jw_y)+exp(jw_x-jw_y)+exp(jw_x+jw_y)] = 0.25(1+cos w_x)(1+cos w_y).

Система является примером двумерного фильтра нижних частот. Частотный отклик системы на плоскости (w_x,w_y), приведенный на рис. 18.3.1, имеет осевую симметрию с коэффициентом передачи 1 в центре (w_x=0, w_y=0) со спадом до нуля при w_x=p и w_y=p.

h(k,l)= q(k)g(l) Ы Q(w_x)G(w_y)= H(w_x,w_y)

Q(w_x) = S_k q(k) exp(-jkЧ w_x). G(w_y) = S_l g(l) exp(-jlЧ w_y).

Импульсный отклик системы. Выражение (18.3.1), по существу, описывает разложение функции Н(w_x,w_y) в двумерный рад Фурье с коэффициентами разложения в виде отсчетов импульсного отклика h(k,l), т.е. прямое преобразование Фурье. Очевидно, что обратным преобразованием Фурье с интегрированием в пределах одного периода из частотного отклика H(w_x,w_y) можно получить импульсный отклик системы:

Пример расчета импульсного отклика фильтра.

Определить импульсный отклик идеального фильтра низких частот с прямоугольной частотной характеристикой вида: H(w_x,w_y) = 1 при |w_x|a<p, |w_y|b<p; H(w_x,w_y) = 0 в остальных случаях.

Импульсный отклик: h(k,l) = exp(jkw_x+jlw_y) dw_x dw_y.

Система разделима: h(k,l)= exp(jkw_x) dw_x exp(jlw_y) dw_y= Ч .

Определить импульсный отклик идеального кругового фильтра нижних частот:

H(w_x,w_y) = 1 при w_x²+w_y² <R²<p²; H(w_x,w_y) = 0 в остальных случаях.

Вычисления по круговой области целесообразно выполнять в полярных координатах: w =,

j = arctg(w_y/w_x), f = arctg(m/n), при этом выражение 18.3.2 перепишется в следующем виде:

h(n,m) = w exp[jw cos(f-j)] dj dw =

= w J_o(w) dw = J₁(R) /,

где J_o(…), J₁(…)- функции Бесселя 1-го рода 0-го и 1-го порядков соответственно.

S(w_x,w_y) = S_n S_m s(n,m) exp(-jnw_x-jmw_y). (18.3.3)

s(n,m) = S(w_x,w_y) exp(jnw_x+jmw_y) dw_xdw_y. (18.3.4)

z(n,m) = h(n,m) _** s(n,m) Ы H(w_x,w_y) S(w_x,w_y) = Z(w_x,w_y).

z(n,m) = c(n,m) s(n,m) Ы C(w_x,w_y) _** S(w_x,w_y) = Z(w_x,w_y).

аЧ s(n,m)+bЧ z(n,m) Ы aS(w_x,w_y)+bZ(w_x,w_y).

2) Пространственный сдвиг:

s(n-N,m-M) Ы S(w_x,w_y) exp(-jNw_x-jMw_y).

dS(w_x,w_y)/dw_x Ы -jn s(n,m),

dS(w_x,w_y)/dw_y Ы -jm s(n,m),

d²S(w_x,w_y)/(dw_x dw_y) Ы -nm s(n,m).

х*(n,m) Ы S*(-w_x,-w_y).

S(w_x,w_y) = S*(-w_x,-w_y).

Re [S(w_x,w_y)] = Re [S(-w_x,-w_y)].

Im [S(w_x,w_y)] = -Im [S(-w_x,-w_y)].

S_n S_m s(n,m) s*(n,m) = S(w_x,w_y) S*(w_x,w_y) dw_x dw_y.

S_n S_m |s(n,m)|² = |S(w_x,w_y)|² dw_x dw_y,

где левая часть уравнения представляет собой полную энергию дискретного сигнала s(n,m), a функция |S(w_n,w_m)|² - спектральную плотность энергии сигнала.

Прямоугольный растр дискретизации. Из способов обобщения одномерной периодической дискретизации на двумерный случай наиболее простым является периодическая дискретизация в прямоугольных координатах:

s(n,m) = s_a(nDx,mDy),

где Dx и Dy - горизонтальный и вертикальный интервалы дискретизации двумерного непрерывного сигнала s_a(x,y) с непрерывными координатами x и y. Ниже значения Dx и Dy, как и в одномерном случае, принимаются равными 1.

Интегральные преобразования Фурье аналоговых сигналов в непрерывной шкале частот W_x и W_y:

S_a(W_x,W_y) =s_a(x,y) exp(-jW_xx-jW_yy) dxdy. (18.4.1)

s_a(x,y) =S_a(W_x,W_y) exp(jW_xx+jW_yy₎ dW_xdW_y. (18.4.2)

S(k,l) =s(n,m) exp(-jn2pk/N-jm2pl/M), (18.4.3)

S(k,l) =exp(-jn2pk/N) s(n,m) exp(-jm2pl/M), (18.4.3')

s(n,m) =S(k,l) exp(-jn2pk/N-jm2pl/M). (18.4.4)

s(n,m) =exp(-jn2pk/N) S(k,l) exp(-jm2pl/M). (18.4.4')

Интерполяционный ряд восстановления двумерного сигнала. Если непрерывный сигнал s_a(x,y) является сигналом с ограниченным спектром, а периоды дискретизации выбраны достаточно малыми и спектры соседних периодов не перекрываются:

S_a(W_x,W_y) = 0 при |W_x|p/Dx, |W_y|p/Dx,

то, как и в одномерном случае, сигнал s_a(x,y) может быть восстановлен по дискретному сигналу с использованием двумерного аналога ряда Котельникова-Шеннона:

s_a(x,y) = S_n S_m s(n,m). (18.4.5)

Сигнал с неограниченным спектром также может быть дискретизирован, однако в этом случае имеет место наложение спектров в смежных периодах, при этом высокие частоты, большие частоты Найквиста, будут "маскироваться", как и в одномерном случае, под низкие частоты главного периода. Эффект "отражения" от границ периода дает еще более сложную картину вследствие интерференции частот, отраженных по разным координатам.

Произвольный растр дискретизации. Понятие прямоугольной дискретизации обобщается на произвольный растр дискретизации с линейно независимыми векторами v₁ = (v₁₁,v₂₁)^T и v₂ = (v₁₂,v₂₂)^T, где T - индекс транспонирования (рис. 18.4.1). Координаты двумерного периодического множества отсчетов на плоскости (x,y):

x = v₁₁n + v₁₂m,

y = v₂₁n + v₂₂m.

=

где = (x,y)^T, =(n,m)^T, =(v₁|v₂)- матрица дискретизации. Определитель матрицы не равен нулю, если вектора v₁ и v₂ линейно независимы. При дискретизации непрерывного сигнала s_a(x,y) матрицей формируется дискретный сигнал:

s() Ь s_a( ).

Двумерное интегральное преобразование Фурье непрерывного сигнала по непрерывному вектору = (W₁,W₂)^T:

S_a() =s_a() exp(-j^T) d, (18.4.6)

s_a() =S_a() exp(j^T) d, (18.4.7)

S() = S_n s() exp(-j^T), (18.4.8)

s() = S() exp(j^T) d. (18.4.9)

где: = (w_х,w_у)^T .

Выражение s() может быть получено дискретизацией выражения s_a() (18.4.7):

s() = s_a() =S_a() exp(j^T) d.

После подстановки в это выражение значения = ^T, получаем:

s() = S_a(/^T) exp(j^Т) d.

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

где - вектор целочисленных значений периодов дискретизированной функции по осям w_х и w_у. Сравнивая последнее выражение с выражением (18.4.9), получаем:

S() =S_a((-2p)/^T),

S(^T) = S_a(-), (18.4.11)

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

Как и в одномерном случае, интервалы дискретизации Dx и Dy определяют главный период двумерного спектра соответственно по осям w_x и w_y и частоты Найквиста: w_xN = p/Dx и w_yN = p/Dy. Спектр дискретного сигнала также является периодическим продолжением спектра аналогового сигнала. Для исключения искажений спектра (наложения спектров боковых периодов на главный период) предельные частоты сигнала должны быть меньше частот Найквиста. На рис. 18.4.2 приведен пример центральной части спектра дискретного сигнала при Dx=1 и Dy=1.

. (18.4.14)

Интерполяция дискретных сигналов. Для сигнала с ограниченным спектром изменением матрицы дискретизации можно подобрать матрицу периодичности таким образом, чтобы в правой части выражения (18.4.11) не было перекрытия спектров. Тогда для значений по точкам ^T области С главного периода спектра выражение (18.4.11) упрощается:

S(^T) = S_a() / |det |. (18.4.15)

S_a() = |det | S(^T) = |det | S(), О С. (18.4.16)

Из выражения (18.4.16) следует, что при корректной дискретизации непрерывной двумерной функции ее спектр с точностью до нормировочного множителя |det | может быть восстановлен по спектру дискретной функции. Соответственно, выполнив обратное преобразование Фурье левой и правой части равенства (18.4.16), получим уравнение восстановления непрерывной функции по ее дискретному варианту (многомерный аналог интерполяционного ряда Котельникова-Шеннона):

s_a() = s()exp(j^T (-)) d.

s_a() = s() f(-), (18.4.17)

где f(..) – интерполяционная функция:

f(-) = exp(j^T (-)) d. (18.4.18)

Все приведенные векторные уравнения могут быть обобщены на Р-мерные функции с заменой константы 4p² там, где она встречается, на (2p)^P.

Прямоугольный и гексагональный растры дискретизации. В принципе, сигнал с ограниченным спектром можно представить по различным растрам дискретизации. Выбор растра обычно производят из условия минимальной плотности отсчетов на плоскости, т.е. минимизацией величины |det|, при котором обеспечивается отсутствие наложений для частот анализируемых сигналов.

На практике для двумерных сигналов используют, как правило, только два варианта растров дискретизации - прямоугольный и гексагональный. Прямоугольному варианту соответствуют диагональные матрицы дискретизации и периодичности (18.4.13-14). Для гексагональной дискретизации, пример которой приведен на рис. 18.4.1, в частном случае при Dt = Dх каждый отсчет располагается на равном расстоянии от шести ближайших отсчетов, при этом матрицы дискретизации:

, .

Допустим, имеем сигнал с частотным спектром, ограниченным круговой областью частот w_r:

S_a(w_х,w_у) = 0 при w_х²+w_у² > w_r².

Круговая область частот вписывается без перекрытий в квадрат со стороной 2w_r или в шестиугольник со стороной 2w_r/. Матрицы дискретизации:

_пр =, det = p²/w_r²,

_гекс =, det = 2p²/(w_r²).

Поскольку плотность отсчетов пропорциональна 1/|det|, то отсюда следует, что для представления одного и того же сигнала гексагональный растр дискретизации требует на 13.4% меньше отсчетов по сравнению с прямоугольным. Эффективность "гексагональной" матрицы возрастает при увеличении размерности сигнала. Так, при 4-мерном сигнале для "гексагональной" матрицы требуется в 2 раза меньше отсчетов, чем для "прямоугольной".

18.5. Многомерный спектральный анализ [9].

Периодические последовательности. Двумерная последовательность s(n,m) прямоугольно периодична, если

s_п(n,m) = s(n,n+M) = s(n+N,m)

для всех (n,m) при целочисленных значениях N и M. Минимальные значения N и M, при которых выполняется данное равенство, называют горизонтальным и вертикальным периодами функции s_п(..), которыми ограничивается прямоугольная область R_N,M главного периода, содержащая NM независимых отсчетов: 0nN-1, 0mM-1.

Последовательность s_п(n,m) с периодами N и M можно представить в виде конечной суммы (ряда Фурье) комплексных синусоид с кратными частотами:

s_п(n,m) =(1/NM)S_п(k,l) exp(j2pnk/N+j2pml/M). (18.5.1)

S_п(k,l) =s_п(n,m) exp(-j2pnk/N-j2pml/M). (18.5.2)

Пример.

Разложить в ряд Фурье периодический сигнал: s_п(n,m) = d(n,m), 0n4, 0m2

(единичные импульсы с периодом: N = 5, M = 3).

S_п(k,l) =d(n,m) exp(-j2pnk/5-j2pml/3) = 1 для всех k и l, т.е. равномерная частотная характеристика в главном диапазоне. Соответственно сам сигнал может быть записан в виде двумерного ряда Фурье:

s_п(n,m) =(1/20)exp(j2pnk/5+j2pml/3).

Конечные последовательности. Если s(n,m) представляет собой последовательность конечной протяженности, имеющей опорную область R_N,M, то периодическую последовательность s_п(n,m) с главным периодом R_N,M можно сформировать периодическим продолжением s(n,m):

s_п(n,m) = s(n-aN,m-bM),

s(n, m) = s_п(n, m), (n, m) < R_N,M.

= 0, в остальных случаях.

Отсюда следует, что любой финитный сигнал может быть полностью определен своим периодическим продолжением и опорной областью.

Аналогично можно записать и для частотной области:

S_п(k, l) = S_a S_b S(k-aN, l-bM).

S(k, l) = S_п(k, l), 0kN, 0lM

= 0, в остальных случаях.

Отсюда значения s(n,m) и S(k,l) можно вычислить с использованием выражений (18.5.1-2) путем последовательности операций:

s(n,m) Ю s_п(n,m) Ы S_п(k,l) Ю S(k,l).

Практически это означает, что для получения ДПФ последовательности конечной протяженности достаточно из выражений (18.5.1-2) для рядов Фурье убрать знак периодичности, при этом следует помнить, что вычисление отсчетов s(..) вне опорной области приведет к вычислению значений не отсчетов s(..), а отсчетов s_п(..) периодического продолжения сигнала s(..).

Таким образом:

1. Дискретизация сигнала в пространственной области вызывает периодизацию частотного спектра сигнала.

2. Дискретизация частотного спектра сигнала вызывает периодизацию его пространственного представления.

3. Прямое и обратное ДПФ сигнала ограниченной протяженности автоматически означает периодизацию как его спектра, так и его пространственного представления.

4. Сигналы, ограниченные в пространстве, можно точно отобразить отсчетами их фурье-преобразования.

5. Частотное представление сигнала с ограниченным спектром обратным фурье-преобразованием может быть точно переведено в пространственную область.

6. Ограниченность как пространственного сигнала, так и его спектра является обязательным условием корректного ДПФ, т.к. в противном случае периодизация сигнала может привести к искажению его спектрального и пространственного представления.

Многомерные последовательности. Определение ДПФ для Р-мерной последовательности с опорной областью R_P = {: 0n_iN_i-1, i=1,2,3, ... ,P} производится введением диагональной матрицы значений N_i:

=,

при этом P-мерное ДПФ записывается в виде:

S() =s() exp(-j^T2p/). (18.5.3)

s() =S() exp(j^T2p/). (18.5.4)

Кратко рассмотрим особенности многомерных ДПФ (на примере двумерных последовательностей).

ДПФ суммы двух последовательностей с опорной областью на R_N,M равно сумме их ДПФ:

аs(n,m)+bz(n,m) Ы aS(k,l)+bZ(k,l),

но при этом все ДПФ должны быть одного размера и этот размер должен быть достаточным, чтобы включить всю опорную область суммарной последовательности аs(n,m)+bz(n,m). Практически это означает, что х(..) и z(..) должны иметь одну и ту же опорную область. Опорная область каждой последовательности при необходимости дополняется нулями.

Операция свертки двух функций в пространственной области отображается операцией умножения фурье-образов функций в частотной области, однако при этом линейная свертка полных пространственных сигналов при ее вычислении через ДПФ в силу периодического продолжения пространственных функций переходит в циклическую свертку (как и для одномерных сигналов). Результат свертки зависит от периодов N и М.

Допустим, что s(n,m) имеет опорную область R_P1,P2, a h(n,m) - R_Q1,Q2. Результат линейной свертки:

s(n,m) = S_k S_l h(k,l) s(n-k,m-l).

Опорная область последовательности s(n,m):

0nP1+Q1-1, 0mP2+Q2-1.

Следовательно, наложения периодов результата свертки не произойдет и циклическая свертка в главном частотном диапазоне будет равна линейной свертке при опорной области ДПФ:

NP1+Q1-1, MP2+Q2-1.

Пример.

Заданы последовательности (начало координат в нижнем левом углу): s(n,m) = , h(n,m) = .

Вычислить свертку z(n,m) = s(n,m) _** h(n,m).

ДПФ размера 2 х 2 для s(n,m). S(0,0) = s(0,0) + s(1,0) + s(0,1) + s(1,1) = 2+1+1+0 = 4

S(1,0) = s(0,0) - s(1,0) + s(0,1) - s(1,1) = 2 -1+1 -0 = 2

S(0,1) = s(0,0) + s(1,0) - s(0,1) - s(1,1) = 2+1 -1 -0 = 2

S(1,1) = s(0,0) - s(1,0) - s(0,1) + s(1,1) = 2 -1 -1+0 = 0

После аналогичного вычисления H(k,l) и перемножения S(k,l)= S(k,l) H(k,l):

S(k,l) = , H(k,l)= , Z(k,l)= .

После обратного ДПФ размера 2 х 2 получим результат циклической свертки: z(n,m) = .

Дополним опорные области s(.) и h(.) до размера 4 х 4 (для исключения искажения спектра, в принципе,

достаточен размер 3 х 3), и повторим вычисления:

s(n,m)= , h(n,m)= .

S(k,l)= , H(k,l)=.

S(k,l)= . s(n,m)= .

Сравнение данного результата с ДПФ размером 2 х 2 позволяет наглядно видеть эффект

цикличности свертки.

В настоящее время имеются разнообразные и весьма эффективные алгоритмы ДПФ. Для прямого вычисления P-мерного ДПФ требуется (N₁N_2...N_P)² операций умножения и сложения. Для многомерного ДПФ, как и для одномерного, существуют алгоритмы быстрых преобразований Фурье. Простейший из них в двумерном ДПФ - разбиение на строки и столбцы, который мы уже рассматривали. Аналогично, Р-мерное ДПФ может заменяться Р-операциями одномерных ДПФ, при этом общее количество операций умножения и сложения сокращается.

 

литература

9. Даджион Д., Мерсеро Р. Цифровая обработка многомерных сигналов. – М.: Мир, 1988. – 488 с.

Copyright ©2005 Davydov А.V.

к оглавлению

Знаете ли Вы, что такое мысленный эксперимент, gedanken experiment?
Это несуществующая практика, потусторонний опыт, воображение того, чего нет на самом деле. Мысленные эксперименты подобны снам наяву. Они рождают чудовищ. В отличие от физического эксперимента, который является опытной проверкой гипотез, "мысленный эксперимент" фокуснически подменяет экспериментальную проверку желаемыми, не проверенными на практике выводами, манипулируя логикообразными построениями, реально нарушающими саму логику путем использования недоказанных посылок в качестве доказанных, то есть путем подмены. Таким образом, основной задачей заявителей "мысленных экспериментов" является обман слушателя или читателя путем замены настоящего физического эксперимента его "куклой" - фиктивными рассуждениями под честное слово без самой физической проверки.
Заполнение физики воображаемыми, "мысленными экспериментами" привело к возникновению абсурдной сюрреалистической, спутанно-запутанной картины мира. Настоящий исследователь должен отличать такие "фантики" от настоящих ценностей.

Релятивисты и позитивисты утверждают, что "мысленный эксперимент" весьма полезный интрумент для проверки теорий (также возникающих в нашем уме) на непротиворечивость. В этом они обманывают людей, так как любая проверка может осуществляться только независимым от объекта проверки источником. Сам заявитель гипотезы не может быть проверкой своего же заявления, так как причина самого этого заявления есть отсутствие видимых для заявителя противоречий в заявлении.

Это мы видим на примере СТО и ОТО, превратившихся в своеобразный вид религии, управляющей наукой и общественным мнением. Никакое количество фактов, противоречащих им, не может преодолеть формулу Эйнштейна: "Если факт не соответствует теории - измените факт" (В другом варианте " - Факт не соответствует теории? - Тем хуже для факта").

Максимально, на что может претендовать "мысленный эксперимент" - это только на внутреннюю непротиворечивость гипотезы в рамках собственной, часто отнюдь не истинной логики заявителя. Соответсвие практике это не проверяет. Настоящая проверка может состояться только в действительном физическом эксперименте.

Эксперимент на то и эксперимент, что он есть не изощрение мысли, а проверка мысли. Непротиворечивая внутри себя мысль не может сама себя проверить. Это доказано Куртом Гёделем.

Понятие "мысленный эксперимент" придумано специально спекулянтами - релятивистами для шулерской подмены реальной проверки мысли на практике (эксперимента) своим "честным словом". Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМА Рыцари теории эфира

10.11.2021 - 12:37: ПЕРСОНАЛИИ - Personalias -> WHO IS WHO - КТО ЕСТЬ КТО - Карим_Хайдаров. 10.11.2021 - 12:36: СОВЕСТЬ - Conscience -> РАСЧЕЛОВЕЧИВАНИЕ ЧЕЛОВЕКА. КОМУ ЭТО НАДО? - Карим_Хайдаров. 10.11.2021 - 12:36: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от д.м.н. Александра Алексеевича Редько - Карим_Хайдаров. 10.11.2021 - 12:35: ЭКОЛОГИЯ - Ecology -> Биологическая безопасность населения - Карим_Хайдаров. 10.11.2021 - 12:34: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> Проблема государственного терроризма - Карим_Хайдаров. 10.11.2021 - 12:34: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> ПРАВОСУДИЯ.НЕТ - Карим_Хайдаров. 10.11.2021 - 12:34: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вадима Глогера, США - Карим_Хайдаров. 10.11.2021 - 09:18: НОВЫЕ ТЕХНОЛОГИИ - New Technologies -> Волновая генетика Петра Гаряева, 5G-контроль и управление - Карим_Хайдаров. 10.11.2021 - 09:18: ЭКОЛОГИЯ - Ecology -> ЭКОЛОГИЯ ДЛЯ ВСЕХ - Карим_Хайдаров. 10.11.2021 - 09:16: ЭКОЛОГИЯ - Ecology -> ПРОБЛЕМЫ МЕДИЦИНЫ - Карим_Хайдаров. 10.11.2021 - 09:15: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Екатерины Коваленко - Карим_Хайдаров. 10.11.2021 - 09:13: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вильгельма Варкентина - Карим_Хайдаров.