- •Сигналы и линейные системы
- •Тема 13. Многомерные сигналы и системы.
- •Содержание
- •Введение.
- •13.1. Двумерные и многомерные сигналы [9].
- •13.2. Двумерные системы.
- •13.3. Частотные характеристики сигналов и систем.
- •13.4. Дискретизация двумерных сигналов [9].
- •13.5. Многомерный спектральный анализ [9].
- •Литература
13.1. Двумерные и многомерные сигналы [9].
Рис.
13.1.1.
В общем случае, двумерный непрерывный сигнал представляет собой функцию, значения которой зависят от двух независимых переменных (аргументов, координат):
s(x,y) = sin(x2+y2), -< x,y <(13.1.1)
График функции (в пределах одного периода) приведен на рис. 13.1.1.
Двумерный дискретный сигнал (цифровой массив) - это функция, определенная на совокупности пар числовых значений координат с определенным шагом дискретизации x и y. В общем случае, при различной физической размерности аргументов x и y, значения x и y не равны друг другу:
sn,m = s(nx,my), -< n,m <. (13.1.1')
Элемент последовательности sn,m представляет собой отсчет двумерной функции s в координатной точке (x=nx,y=my), где значения x и y – независимые переменные (аргументы) функции. Для числовых массивов значения шага дискретизации по аргументам также могут приниматься равными 1 (независимо от размерности) и использоваться аргументация s(n,m) sn,m. Результаты геофизических съемок какого-либо одного геофизического параметра по поверхности земли относятся к двумерным функциям: дискретным - если это отсчеты в отдельных точках по определенной координатной сети (x,y), или смешанным - если это непрерывная регистрация данных по профилям (например - мощности экспозиционной дозы гамма излучения горных пород при аэросъемке). Но в настоящее время геофизические съемки относятся даже не к двумерным, а к многомерным функциям, так как регистрируется, как правило, сразу несколько физических параметров геологических сред. Так, например, при спектрометрической съемке естественной радиоактивности горных пород регистрируется содержание в горных породах урана, тория и калия, в гравиразведке - трехкоординатный вектор силы тяжести, и т.п. Если на какой-либо площади проведена съемка нескольких видов геофизики, то их результаты также могут рассматриваться в совокупности, как многомерная функция физических параметров данной геологической среды.
По определениям (13.1.1) двумерные функции и сигналы, равно как и многомерные, имеют бесконечную протяженность по координатам. На практике мы всегда имеем дело с конечными координатами наших данных. Учитывая это, будем считать, что значения наших сигналов за пределами определенных координат равны нулю.
Отметим некоторые двумерные последовательности (функции, сигналы), имеющие специальные названия.
Рис.
13.1.2.
n,m = 1, при n = m = 0.
= 0, при n0, m0.
n,m = n m,
где n, m - одномерные единичные импульсы (импульсы Кронекера) по координатам n и m. Стилизованное графическое представление двумерного единичного импульса приведено на рис. 13.1.2.
Произвольное расположение двумерного единичного импульса по координатам n1, m1 соответственно записывается в виде: ((n-n1)x,(m-m1)y) = n-n1,m-m1. Попутно напомним, что математическая запись импульса Кронекера обозначает не единичный отсчет, а функцию, определяющую место положения единичного отсчета и нулевые значения по остальным координатам (аргументам).
Рис.
13.1.3.
На рис. 13.1.3 приведены два двумерных линейных импульса, первый - по координате m = 0: s(n,m) = (m), и второй импульс по координате n = 2: s(n,m) = (n-2).
Очевидно, что для P-мерных случаев точно таким же образом могут быть определены P-мерные единичные импульсы, P-мерные линейные импульсы, P-мерные площадные импульсы и т.д., хотя понятие импульса, заимствованное из теории одномерных сигналов, здесь несколько не к месту.
Рис.
13.1.4.
u(n,m) = 1, при n0 и m0,
= 0, в остальных случаях.
u(n,m) = u(n) u(m),
где u(j) представляют собой единичные ступеньки соответственно по координатам n и m: u(j)=1 при j0, u(j)=0 при j<0. Двумерная единичная ступенька отлична от нуля в одном квадранте (n,m)- плоскости.
Экспоненциальная последовательность: s(n,m) = anbm, -< n,m <, где а и b в общем случае комплексные числа. При а = exp(j1), b = exp(j2), |а|=1, |b|=1:
s(n,m) = exp(jn1+jm2) = cos(n1+m2)+jsin(n1+m2).
Экспоненциальные последовательности, как и в одномерном случае, являются собственными функциями двумерных линейных систем, инвариантных к сдвигу.
Разделимые последовательности. Разделимой называют последовательность, которую можно представить в виде произведения одномерных последовательностей. Так, для двумерной разделимой последовательности:
s(n,m) = s(n) s(m).
Разделение возможно для немногих практических сигналов. Однако любое двумерное множество с конечным числом ненулевых отсчетов разлагается на конечную сумму разделимых последовательностей:
s(n,m) =sin(n) sim(m),
где N- число ненулевых строк или столбцов массива. В крайнем случае, для этого достаточно выразить s(n,m) в виде суммы отдельных строк:
s(n,m) = s(n,i) (m-i). (13.1.2)
Рис.
13.1.5.
Рис.
13.1.6.
s(n,m+M) = s(n,m),
s(n+N,m) = s(n,m), (13.1.3)
обладает периодичностью в двух направлениях, по n и по m. Значения М и N называют интервалами периодичности сигнала соответственно по координатам m и n (горизонтальными и вертикальными интервалами периодичности). Прямоугольная форма области периода (пример на рис.13.1.6) наиболее удобна при обработке данных, но не является единственно возможной.
Для двумерных последовательностей условия (13.1.3) могут рассматриваться как частный случай общих условий периодичности:
s(n+N1, m+M1) = s(n,m), (13.1.4)
s(n+N2, m+M2) = s(n,m),
D = N1M2 - N2M1 0.
Рис.
13.1.7.
Понятие периодичности можно обобщить на многомерные сигналы. P-мерный сигнал s() будет представлять собой P-мерную периодическую последовательность, если существует P линейно независимых P-мерных целочисленных N-векторов периодичности, с которыми выполняется условие:
s() = s(+), i = 1,2,3, ... ,P.
Столбцы векторов Ni образуют матрицу периодичности N размером P х P. Векторы периодичности матрицы линейно независимы при наличии у матрицы ненулевого определителя. Абсолютное значение определителя равно числу отсчетов в периоде. Последовательность s() прямоугольно периодична для случаев диагональной матрицы N. Если функция s() периодична с матрицей периодичности N, то для любого целочисленного вектора Р имеет место s(+) = s(), и матрица PN также будет матрицей периодичности для s(). Отсюда следует, что любая многомерная периодическая последовательность имеет не единственную матрицу периодичности.