СТАЦИОНАРНЫЕ ЛИНЕЙНЫЕ СИСТЕМЫ

Входной сигнал системы может представлять собой m - мерный вектор (m входных сигналов), а выходной сигнал n - мерный вектор, при этом система будет иметь m входов и n выходов. Пример такой системы в геофизике: трехканальный гамма-спектрометр, на три входа решающего блока которого поступают сигналы от калиевого, радиевого и ториевого каналов спектрометра, а на три выхода выводятся сигналы содержаний калия, урана и тория, при этом системный оператор реализует алгоритм решения системы трех линейных уравнений с тремя неизвестными.

Для определения системы необходимо задать характер, тип и области допустимых величин входных и выходных сигналов. Как правило, системы выполняются на сигналы одного типа по входу/выходу и подразделяются на системы непрерывного времени (аналоговые или дискретные сигналы на входе и выходе) и цифровые системы. Совокупность системного оператора Т и пространства сигналов образует математическую модель системы.

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

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

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

Графическое отображение операций (цифровая форма) приведено на рис. 11.1.1.

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

Для систем, с размерностью 2 и более существует также еще одна базовая операция, которая называется операцией пространственного маскирования, которая может рассматриваться как обобщение скалярного умножения. Так, для двумерных систем:

где u(x,y) – двумерный входной сигнал, c(x,y) – пространственная маска постоянных (весовых) коэффициентов. Пространственное маскирование представляет собой поэлементное произведение значений сигнала с коэффициентами маски.

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

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

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

Другой важной особенностью ЛИС - систем является то, что любые их комбинации также являются ЛИС - системами, а любую сложную ЛИС - систему можно разложить на комбинации простых систем. Так, например, при последовательном (каскадном) соединении систем, когда выходной сигнал одной системы служит входным сигналом для второй и т.д., образуемая система в целом также является ЛИС - системой, если линейны и инвариантны к сдвигу все системы, в нее входящие, при этом по отношению к общей системной операции преобразования порядок соединения входящих в нее систем значения не имеет.

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

При обработке данных на ЭВМ ограничение по каузальности системного оператора снимается. В программном распоряжении системы могут находиться как "прошлые", так и "будущие" значения входных отсчетов, при этом уравнение (11.1.3) будет иметь вид:

Представим, что на одной полоске бумаги выписаны по порядку сверху вниз значения данных x(k). На второй полоске бумаги находятся записанные в обратном порядке значения коэффициентов системы b_n. Для вычисления y(k) располагаем вторую полоску против первой таким образом, чтобы значение b₀ совпало со значением x(k), перемножаем все значения b_n с расположенными против них значениями x(k-n) и суммируем результаты перемножения. Результат суммирования является выходным значением сигнала y(k). Сдвигаем окно системы - полоску коэффициентов b_k, на один отсчет последовательности x(k) вниз (по порядку возрастания номеров k) или массив x(k) сдвигаем на отсчет вверх и вычисляем аналогично следующее значение, и т.д.

Описанный процесс свертки в вещественной области массива данных x(k) с нерекурсивным оператором системы b_n (массивом весовых коэффициентов системы) обычно называют нерекурсивной цифровой фильтрацией данных, а саму систему, если она выполняет только данную операцию, нерекурсивным цифровым фильтром (НЦФ).

Уравнение НЦФ: y_k =b_n x_k-n, b_n = 0,2. Начальные условия - нулевые.

Входной сигнал: x_k = {0,0,0,0,10,0,10,0,0,0,0}. Выходной сигнал: y_k = {0,0,2,2,4, 4, 4,2,2,0,0}.

Заметим: сумма коэффициентов b_n сглаживающего НЦФ равна 1, при этом сумма значений выходного сигнала равна сумме значений входного сигнала. Сигнал "размазался" по аргументу.

Рекурсивные цифровые системы. Системы, которые описываются полным разностным уравнением (11.1.2), принято называть рекурсивными цифровыми системами (РЦС) или рекурсивными цифровыми фильтрами (РЦФ), так как в вычислении текущих значений выходного сигнала участвует не только входной сигнал, но и значения выходного сигнала, вычисленные в предшествующих циклах расчетов. С учетом последнего фактора рекурсивные системы называют системами с обратной связью. Пример рекурсивной системы приведен на рис. 11.1.3.

Полное окно рекурсивной системы состоит из двух составляющих: нерекурсивной части b_n,, аналогичной окну нерекурсивной системы и ограниченной в работе текущими и "прошлыми" значениями входного сигнала (при реализации на ЭВМ возможно использование и “будущих” отсчетов сигнала), и рекурсивной части a_m, которая работает только с "прошлыми", ранее вычисленными значениями выходного сигнала. Техника вычислений для РЦС приведена на рис. 11.1.4.

Пример. Уравнение РЦС: y_k= b_ox_k+a₁y_k-1, при b_o=a₁=0.5, y_-1=0. Входной сигнал: x_k = {0, 10, 0, 0, 0,....}.

у_o = 0,5x_o + 0,5y_-1 = 0 y₁ = 0,5x₁ + 0,5y_o = 5 y₂ = 0,5x₂ + 0,5y₁ = 2.5 y₃ = 0,5x₃ + 0,5y₂ = 1.25

y₄ = 0,5x₄ + 0,5y₃ = 0.625 y₅ = 0,5x₅ + 0,5y₄ = 0.3125 и т.д.

Выходной сигнал: y_k = {0, 5, 2.5, 1.25, 0.625, 0.3125, 0.15625,...}

Из примера можно видеть, что реакция РЦС на конечный входной сигнал, в принципе, может иметь бесконечную длительность, в отличие от реакции НЦС, которая всегда ограничена количеством членов b_k (окном системы).

Пример. Уравнение РЦС: y_k = b_ox_k+a₁y_k-1+a₂y_k-2, при b_o=a₁=a₂=1. x_k = {0, 1, 0, 0, 0,....}, y_-1 = y_-2 = 0.

y₄ = x₄+y₃+y₂ = 3. y₅ = x₅+y₄+y₃ = 5. y₆ = x₆+y₅+y₄ = 8. y₇ = x₇+y₆+y₅ = 13. и т.д.

Заметим: сумма коэффициентов обратной связи a_i больше 1 и выходной сигнал идет "в разнос".

Уравнение РЦФ: y_k=(x_k+x_k-1)/2+y_k-1, начальные условия - нулевые.

Входной сигнал: x_k={0,0,2,2,8,0,4,4,0,..}.

Выполните фильтрацию. Контроль: y_k= {0,0,1,3,8,12,14,18,20,20,...}.

Стационарные и нестационарные системы. Система считается стационарной и имеет постоянные параметры, если ее свойства (математический алгоритм оператора преобразования) в пределах заданной точности не зависят от входного и выходного сигналов и не изменяются ни во времени, ни от каких-либо других внешних факторов. Математически это означает задание системы уравнениями типа (11.1.1-2) с постоянными значениями коэффициентов a_j и b_i и реакция системы на какое-либо воздействие не зависит от времени (координат) его приложения. В противном случае система является нестационарной или параметрической (системой с переменными параметрами). Среди последних большое значение имеют так называемые адаптивные системы обработки данных. В этих системах производится, например, оценивание определенных параметров входных и выходных сигналов, по результатам сравнения которых осуществляется подстройка параметров преобразования (переходной характеристики системы) таким образом, чтобы обеспечить оптимальные по производительности условия обработки сигналов или минимизировать погрешность обработки.

Импульсный отклик системы. По определению, импульсными характеристиками систем (второй широко используемый термин - импульсный отклик систем) называются функции h(t) для аналоговых и h(kDt) для цифровых систем, которые является реакцией (откликом) систем на единичные входные сигналы: дельта-функцию d(t) для аналоговых и импульс Кронекера d(kDt) для цифровых систем, поступающие на вход систем соответственно при t=0 и k=0. Эта реакция однозначно определяется оператором преобразования:

Импульсный отклик аналоговой системы, как результат операции над дельта-функцией, в определенной степени представляет собой математическую абстракцию идеального преобразования. С практической точки зрения под импульсным откликом аналоговой системы можно понимать математическое отображение реакции системы на импульсный входной сигнал произвольной формы с площадью, равной 1, если длительность сигнала пренебрежимо мала по сравнению со временем реакции системы или с периодом ее собственных колебаний. Под временем (длиной) реакции системы обычно понимают интервал, на котором значения функции h(t) существенно отличаются от нуля после прекращения действия единичного сигнала на ее входе.

Для цифровых систем импульсный отклик однозначно определяется реакцией системы на импульс Кронекера d(kDt)=1 при k=0.

Функцию импульсного отклика называют также весовой функцией системы.

На рисунке 11.2.1 приведен пример импульсного отклика h(t) интегрирующей RC-цепи. При подаче на вход RC-цепи импульса заряда Dq емкость С заряжается до напряжения V_о = Dq/C и начинает разряжаться через сопротивление R, при этом напряжение на ней изменяется по закону v(t) = V_oЧ e^-t/RC = (Dq/C)Ч e^-t/RC. Отсюда, отклик RC-цепи по выходному напряжению на входной сигнал Dq = 1: h(t) = (1/C)Ч e^-t/RC. По существу, импульсным откликом системы h(t) определяется доля входного сигнала, которая действует на выходе системы по истечении времени t после поступления сигнала на вход (запаздывающая реакция системы).

Реакция системы на произвольный сигнал. Если функция импульсного отклика системы известна, то, с учетом принципа суперпозиции сигналов в линейной системе, можно выполнить расчет реакции системы в любой произвольный момент времени на любое количество входных сигналов с любыми моментами времени их прихода путем суммирования запаздывающих реакций системы на эти входные сигналы, как это показано на рис. 11.2.1 для трех входных импульсов. В общем случае произвольный сигнал на входе системы может быть разложен в линейную последовательность взвешенных единичных импульсов:

y(t) = T[s(t)] є T[s(t)d(t-t) dt]. (11.2.2)

На основании принципа суперпозиции линейный оператор Т может быть внесен под знак интеграла, т.к. последний представляет собой предельное значение суммы. При этом операция преобразования действует только по переменной t:

y(t) =s(t) Т[d(t-t)] dt =s(t) h(t-t) dt. (11.2.3)

Это выражение представляет собой интеграл Дюамеля или свертку (конволюцию) входного сигнала с импульсной характеристикой системы. Заменой переменных t-t = t можно убедиться в том, что эта операция, как и положено свертке, коммутативна:

s(t) h(t-t) dt є h(t) s(t-t) dt.

Аналогично, для дискретных сигналов:

y(kDt) =s(nDt) h(kDt-nDt) є h(nDt) s(kDt-nDt). (11.2.3')

В символической форме математического представления:

y(t) = s(t ) _* h(t-t) є s(t-t) _* h(t ) є s(t) _* h(t).

В реальных физических системах импульсный отклик h(t) равен нулю при t<0 (реакция на выходе системы не может опережать входной сигнал) и, как правило, отличен от нуля только на определенном интервале r, по которому и ведется интегрирование или суммирование в выражениях свертки. При обработке данных на ЭВМ требований по односторонности импульсного отклика не предъявляется, равно как и по его размерам вперед и назад от нуля по координатам.

Усиление постоянной составляющей сигнала. Подадим на вход системы постоянный сигнал s(t) = A. При этом сигнал на выходе системы:

y(t) = h(t) s(t-t) dt = Аh(t) dt = АЧ К_пс, (11.2.4)

т.е. площадь импульсного отклика (для цифровой системы соответственно сумма коэффициентов импульсного отклика) является коэффициентом К_пс усиления постоянной составляющей входного сигнала. Если при обработке сигналов должны изменяться только динамические характеристики их формы без изменения постоянной составляющей, а равно и различных постоянных уровней (фона, пьедесталов, региональных трендов и т.п.), то площадь импульсного отклика (сумма коэффициентов) должна нормироваться к единице.

На рис. 11.2.2 приведен пример выполнения свертки рассмотренной нами выше RC-цепью при нормированной к 1 площади импульсного отклика h(t). Входной сигнал s(t) находится на постоянном фоновом значении, в данном случае - нулевом, при этом, как и следовало ожидать, площадь выходного сигнала y(t) равна площади входного сигнала.

Усиление шумов. Критерием качества системы при использовании любого метода обработки информации можно считать выполнение целевого назначения с минимальным усилением шумов (максимальным их подавлением). Допустим, что система имеет нормированный к 1 импульсный отклик h(k). Обозначим через e(k) аддитивный шум с математическим ожиданием M{e(k)}= = 0 и дисперсией D_e = s², который в сумме с сигналом поступает на вход системы. Значения e(k) статистически независимы и некоррелированы с сигналом. С учетом помехи во входном сигнале значение сигнала на выходе системы:

y(k) = h(n)[x(k-n)+e(k-n)].

Математическое ожидание значений выходного сигнала:

M{y(k)} = h(n)[x(k-n)+M{e(k-n)}] = h(n) x(k-n).

Вычислим дисперсию распределения отсчетов выходного сигнала:

D{y(k)} = M{[h(n)[x(k-n)+e(k-n)]-M{y(k)}]²} = M{[h(n) e(k-n)]²}.

Если правую часть последнего выражения представить в виде

M{[h(n) e(k-n)]Ч [ h(m) e(k-m)]},

то в этом выражении математические ожидания всех членов произведения с сомножителями e(n)e(m) при n № m равны 0 в силу статистической независимости значений шума. Остаются только члены с n = m, т.е.:

M{h²(n)e²(n)} = h²(n) M{e²(n)} = D_e h²(n) = s² h²(n). (11.2.5)

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

Пример. Сглаживающий фильтр: y(k) = 0.2x(k-n).

Коэффициент усиления шумов: 5 (0.2²) = 0.2. Дисперсия шумов уменьшается в 1/0.2 = 5 раз.

Для систем с m входами и n выходами аналогично определяются парциальные импульсные отклики h_ij(t), i = {1,2, ... ,n}, j = {1,2, ... ,m}, каждым из которых отображается сигнал на i-м выходе при поступлении сигнала d(t) на j-й вход. Полная совокупность импульсных откликов образует матрицу:

,

а выражение свертки приобретает вид:

=dt.

Определение импульсной реакции требуется, как правило, для рекурсивных систем, так как импульсная реакция для НЦС специального определения не требует:

h(k) =b(n)d(k-n) є b(k).

Если выражение для системы известно в общей форме (11.1.2), определение импульсной реакции производится подстановкой в уравнение системы импульса Кронекера с координатой k = 0 при нулевых начальных условиях, при этом сигнал на выходе системы будет представлять собой импульсную реакцию системы: y(k) є h(k).

Пример. Уравнение РЦС: y_k = x_k + 0.5y_k-1. Входной сигнал: x_k= d_o= {1,0,0,0,...}.

Расчет выходного сигнала при нулевых начальных условиях:

y_o = x_o+0.5 y_-1 = 1+0 = 1 = h_o. y₁ = x₁+0.5 y_o = 0+0.5 = 0.5 = h₁. y₂ = x₂+0.5 y₁ = 0+0.25 = 0.25= h₂

y₃ = x₃+0.5 y₂ = 0.125 = h₃. y₄ = x₄+0.5 y₃ = 0.0625 = h₄, т.д.

Импульсный отклик системы: h_k = {1, 0.5, 0.25, 0.125, ... } є (O.5)^k, k = 0,1,2....

Определение импульсной реакции физических систем обычно производится подачей на вход систем ступенчатой функции u_o(k) = 1 при k і 0, и u_o(k) = 0 при k < 0:

g(k) =h(n) u_o(k-n) =h(n).

h(k) = g(k) - g(k-1), k=0,1,2,...

Функция g(k) получила название переходной характеристики системы (перехода из одного статического состояния в другое).

11.3. Передаточные функции цифровых систем [5,18].

Z-преобразование является удобным методом решения разностных уравнений линейных систем. Применяя z-преобразование к обеим частям равенства (11.1.2), получаем:

Y(z)a_mz^m = X(z)b_nzⁿ, (11.3.1)

где X(z),Y(z)- соответствующие z-образы входного и выходного сигнала. Из этого выражения, полагая a_o = 1, получаем в общей форме функцию связи входа и выхода системы - уравнение передаточной функции системы (или системной функции) в z-области:

H(z) = Y(z)/X(z) =b_nzⁿ(1+a_mz^m). (11.3.2)

Для нерекурсивных систем при a_m = 0:

H(z) =b_nzⁿ. (11.3.3)

При подаче на вход системы единичного импульса Кронекера d_о, имеющего z-образ d(z) = zⁿ = 1, сигнал на выходе системы будет представлять собой импульсную реакцию системы y(k) = h(k), при этом:

H(z) є Y(z) = Y(z)/d(z) = TZ[y(k)] =h(k)z^k, (11.3.4)

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

h(k) Ы H(z). (11.3.5)

Если функция H(z) представлена конечным степенным полиномом, как это обычно имеет место для НЦС, то обратное z-преобразование элементарно. Передаточная функция РЦС также может быть представлена степенным полиномом прямым делением числителя на знаменатель правой части выражения (11.3.2), однако результат при этом может оказаться как конечным, так и бесконечным, т.е. система может иметь либо конечную, либо бесконечную импульсную характеристику. Системы с бесконечной импульсной характеристикой получили название БИХ-систем, с конечной импульсной характеристикой соответственно КИХ-систем. Нерекурсивные системы всегда имеют конечную импульсную характеристику, т.к. длительность импульсной реакции НЦС определяется окном фильтра.

Пример. Передаточная функция рекурсивной системы: H(z) = (1-z⁵)/(1-z).

Делением числителя на знаменатель получаем: H(z) = 1+z+z²+z³+z⁴. H(z) Ы h(n) = {1,1,1,1,1}.

Система имеет конечную импульсную характеристику.

Пример. Передаточная функция: H(z)=1/(1-2z). Методом обратного преобразования: h(n) = 2ⁿ.

Система имеет бесконечную импульсную характеристику.

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

Для конечного по энергии входного сигнала, можно записать:

|y(t)| Ј |h(t)||x(t-t)| dt = A|h(t)| dt.

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

|h(t)| dt < Ґ , (11.3.6)

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

|h(n)| < Ґ . (11.3.6')

Анализ устойчивости может быть проведен по передаточной функции. В устойчивой системе значение H(z) должно быть конечным во всех точках z-плоскости, где |z| Ј 1, а, следовательно, передаточная функция не должна иметь особых точек (полюсов) при z Ј 1 (на и внутри единичного круга на z-плоскости). Полюсы определяются корнями многочлена знаменателя передаточной функции H(z).

Пример.

Передаточная функция РЦФ: H(z) = 1/(1+1.25z). Корни знаменателя: z = -0,8. |z| < 1.

Фильтр неустойчив.

Передаточная функция: H(z) = (1-z³)/(1-0.6z+0.25z²). Корни: z_1,2 = 1.2+1.6j. |z| = 2 > 1.

Фильтр устойчив.

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

Проверка на устойчивость требуется только для рекурсивных цифровых фильтров (систем с обратной связью), нерекурсивные системы всегда устойчивы.

11.4. Частотные характеристики систем [2,5].

Для линейных систем, принимая в качестве сигнала на входе системы собственную функцию x(kDt) = B(w)exp(jwkDt), мы вправе ожидать на выходе системы сигнал y(kDt) = A(w)exp(jwkDt). Подставляя эти выражения в разностное уравнение системы (11.1.2), получаем:

a_m A(w)exp(jwkDt-jwmDt) =b_n B(w)exp(jwkDt-jwnDt).

A(w)exp(jwkDt)a_m exp(-jwmDt) = B(w)exp(jwkDt)b_n exp(-jwnDt).

A(w)a_m exp(-jwmDt) = B(w)b_n exp(-jwnDt). (11.4.1)

Отсюда, частотная передаточная функция системы (частотная характеристика при нормировке к а_о=1):

H(w) = A(w)/B(w) =b_n exp(-jwnDt)[1+a_m exp(-jwmDt)]. (11.4.2)

Нетрудно убедиться, что подстановкой z = exp(-jwDt) в выражение передаточной функции H(z) (11.3.2) может быть получено абсолютно такое же выражение для частотной характеристики, т.е.:

H(w) є H(z) при z = exp(-jwDt).

При обратном преобразовании H(z) во временную область с использованием выражений (11.3.4-5) отсюда следует также, что частотная характеристика системы представляет собой Фурье-образ ее импульсной реакции, и наоборот. При Dt = 1:

H(w) =h(n) exp(-jwn), (11.4.3)

h(n) = (1/2p)H(w) exp(jwn) dw. (11.4.4)

В общем случае H(w) является комплексной функцией, модуль которой R(w) называется амплитудно-частотной характеристикой системы (АЧХ), а аргумент j(w) - фазочастотной характеристикой (ФЧХ).

A(w) = |H(w)| = .

j(w) = arctg(Im H(w)/Re H(w)).

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

Основные свойства частотных характеристик систем:

1. Частотные характеристики являются непрерывными функциями частоты.

2. При дискретизации данных по интервалам Dt функция H(w) является периодической. Период функции H(w) равен частоте дискретизации входных данных F = 1/Dt. Первый низкочастотный период (по аргументу w от -p/Dt до p/Dt, по f от -1/2Dt до 1/2Dt) называется главным частотным диапазоном передачи сигнала. Граничные частоты главного частотного диапазона соответствуют частоте Найквиста ± w_N, w_N = p/Dt. Частота Найквиста определяет предел частотной разрешающей способности системы по обработке данных.

3. Для систем с вещественными коэффициентами импульсной реакции h(nDt) функция АЧХ является четной, а функция ФЧХ - нечетной. С учетом этого частотные характеристики систем обычно задаются только на интервале положительных частот 0-w_N главного частотного диапазона. Значения функций на интервале отрицательных частот являются комплексно сопряженными со значениями на интервале положительных частот.

11.5. Реакция систем на случайные сигналы [2,14].

Если сигнал на входе линейной системы является детерминированным, то, при известных параметрах системы, его соотношение с выходным сигналом является однозначным. Таким же однозначным является соотношение процессов на входе и выходе и для случайных сигналов, однако в силу природы последних явное представление, как входного сигнала, так и отклика системы, не представляется возможным. Для описания отклика системы необходимо использовать статистический подход. При рассмотрении данной темы ограничимся только физически реализуемыми системами с односторонним импульсным откликом h(t) (h(t)=0 при t<0) и соответствующей частотной характеристикой H(f). Если параметры входного сигнала специально не оговариваются, то по умолчанию принимается, что на вход системы поступает реализация случайного стационарного процесса x(t) с нулевым средним и вызывает сигнал y(t) на выходе системы.

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

Допустим, что система имеет импульсный отклик h(t) = exp(-at), t і 0. Квазидетерминированный случайный сигнал стационарен, не обладает свойством эргодичности, но может быть описан в явной математической форме. Зададим сигнал на входе системы следующего вида:

x(t) = A + cos(2t+j),

где A и j - взаимно независимые случайные величины, причем j равномерно распределена в интервале [0, 2p]. Выходной сигнал определится выражением:

y(t) = h(t) _* x(t-t) є h(t)Ч x(t-t) dt = A/3 + [3 cos(2t+j) + 2 sin(2t+j)]/13.

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

Математическое ожидание произвольного случайного стационарного сигнала x(t) на выходе линейной системы определится выражением:

= М{y(t)}= M{x(t-t)}Ч h(t) dt =h(t) dt = К_пс. (11.5.1)

Отсюда следует, что математическое ожидание выходных сигналов системы равно математическому ожиданию входных сигналов, умноженному на коэффициент усиления системой постоянной составляющей. Если система не пропускает постоянную составляющую сигналов (площадь или сумма коэффициентов импульсного отклика системы равна нулю), то случайный выходной сигнал всегда будет иметь нулевое математическое ожидание.

Корреляционные соотношения. Для произведения выходных сигналов y(t) и y(t+t) линейной системы можно записать:

y(t)Ч y(t+t) =h(a)h(b) x(t-a)x(t+t-b) da db.

Если взять математические ожидания от обеих частей этого равенства, то, с учетом соотношения в подынтегральном выражении

M{x(t-a) x(t+t-b)} = -R_x(t-a-t-t+b) = R_x(t+a-b),

получим:

R_y(t) =h(a)h(b) R_x(t+a-b) da db є R_x(t) _* h(t+a) _* h(t-b). (11.5.2)

Таким образом, функция автокорреляции (АКФ) выходного сигнала равна АКФ входного сигнала, свернутой дважды, в прямом и обратном направлении, с импульсным откликом системы, что сохраняет четность АКФ выходного сигнала. Для нецентрированных процессов аналогичное заключение действительно и для ковариационных функций.

Заметим, что для свертки импульсных откликов, производя замену t-b = t, мы имеем равенство:

h(t+a) _* h(t-b) = h(t+a+b) _* h(t) = h(t) _* h(t+g) = K_h(t),

где K_h(t) - функция ковариации импульсного отклика системы. Отсюда:

R_y(t) = R_x(t) _* K_h(t). (11.5.2')

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

Для взаимной корреляционной функции (ВКФ) R_xy входного и выходного сигналов соответственно имеем:

x(t)Ч y(t+t) =h(a) x(t)y(t+t-a) da.

R_xy(t) =h(a) R_x(t-a) da є R_x(t) _* h(t-a). (11.5.3)

т.е. функция взаимной корреляции входного и выходного сигналов равна свертке АКФ входного сигнала с функцией импульсного отклика системы - фильтрации АКФ сигнала этим же фильтром. Заключение действительно и для функций ковариации.

Другая взаимно корреляционная функция R_yx может быть получена из соотношения:

R_yx(t) = R_xy(-t) є R_x(t) _* h(t+a). (11.5.4)

Отметим, что для статистически независимых случайных величин при одностороннем импульсном отклике h(t)=0 при t<0 функция R_xy(t) также является односторонней и равна 0 при t<0, а функция R_yx соответственно равна 0 при t>0.

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

Применяя преобразование Фурье к выражениям (11.5.2), для спектра мощности выходного сигнала получаем:

W_y(f) = W_x(f) |H(f)|². (11.5.5)

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

Аналогично, для взаимного спектра мощности на основе выражений (11.5.3-4):

W_xy(f) = W_x(f) H(f). (11.5.6)

W_yx(f) = W_x(f) H(-f).

Взаимный спектр мощности при одностороннем импульсном отклике является комплексным и содержит как амплитудную, так и фазовую характеристику процесса.

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

H(f) = W_xy/W_x Ы h(t).

Дисперсия выходного сигнала может быть определена с использованием формул (11.5.2,5):

s_y² = R_y(0) =W_x(f) |H(f)|² df є R_x(0)h²(t) dt = s_x²h²(t) dt, (11.5.7)

что полностью соответствует полученной ранее формуле (11.2.5) для цифровой системы.

Если сигнал нецентрированный и значение дисперсии входного сигнала неизвестно, то по аналогичным формулам вычисляется сначала средний квадрат выходного сигнала или так называемая средняя мощность сигнала:

== R_y(0) є h²(t) dt є W_x(f) |H(f)|² df, (11.5.8)

Вывод: средняя мощность выходного сигнала равна средней мощности входного сигнала, умноженной на квадрат площади импульсной реакции системы (для цифровых систем - сумму квадратов коэффициентов импульсного отклика). Для центрированных случайных сигналов средняя мощность равна дисперсии сигналов. Для нецентрированных выходных сигналов:

s_y² = - ² є (-²)h²(t) dt. (11.5.9)

Функция когерентности дает оценку точности принятой линейной модели системы. Когерентность входного и выходного сигналов системы оценивается по формуле:

g_xy²(f) = |W_xy(f)|²/[W_x(f)Ч W_y(f)]. (11.5.10)

Если функции W_x(f) и W_y(f) отличны от нуля и не содержат дельта-функций, то для всех частот f значения функции когерентности заключены в интервале:

0 Ј g_xy²(f) Ј 1.

Для исключения дельта-функций на нулевой частоте (постоянная составляющая сигнала) определение функции когерентности производится по центрированным сигналам. Для линейных систем с постоянными параметрами функция когерентности равна 1, в чем нетрудно убедиться, если в формулу (11.5.10) подставить выражения W_xy и W_y, определенные через W_x в формулах (11.5.5-6). Для совершенно не связанных сигналов функция когерентности равна нулю. Промежуточные между 0 и 1 значения могут соответствовать трем ситуациям:

1. Система осуществляет преобразование x(t) Ю y(t), но в измерениях этих сигналов или одного из них присутствует внешний шум. Так, например, в сигналах, зарегистрированных с ограничением по разрядности, появляется шум квантования (округления значений).

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

3. Выходной сигнал y(t) помимо x(t) зависит еще от каких-то входных или внутренних системных процессов.

Величина 1-g_xy²(f) задает долю среднего квадрата сигнала y(t) на частоте f, не связанную с сигналом x(t).

Использование функций когерентности в практических методах анализа случайных данных подробно рассмотрено в работе /2/.

11.6. Структурные схемы систем [5,20].

Структурные схемы. Алгоритмы обработки сигналов в системах наглядно отображаются в виде структурных схем. Базовые элементы схем и примеры построения структурных схем приводились ранее на рис. 11.1.1-3. Как правило, структурные схемы соответствуют программной реализации систем и будут рассматриваться ниже применительно к цифровым системам, но не определяет аппаратной реализации в специальных радиотехнических устройствах, которая может существенно отличаться от программной реализации, особенно для аналоговых систем.

Графы систем. Наряду со структурной схемой система может быть представлена в виде графа, который отображает диаграмму прохождения сигналов и состоит из направленных ветвей и узлов.

Пример структурной схемы системы с передаточной функцией H(z) = (1+b₁z)/(1+a₁z) и графа, ей соответствующего, приведен на рисунке 11.6.1. С каждым i узлом графа связано значение сигнала x_i(k) или его образа X_i(z), которые определяются суммой всех сигналов или их z-образов входящих в узел ветвей. В каждой ij - ветви (из узла i в узел j) происходит преобразование сигнала в соответствии с передаточной функцией ветви, например, задержка сигнала или умножение на коэффициент.

Соединения систем. Различают следующие соединения систем (рис. 11.6.2-4).

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

H(z) = H₁(z)Ч H₂(z)Ч ...Ч H_N(z).

2. Параллельное соединение. Сигнал подается на входы всех параллельно соединенных систем одновременно, выходные сигналы систем суммируются. Эквивалентная передаточная функция общей системы равна сумме передаточных функций систем, в нее входящих:

H(z) = H₁(z)+H₂(z)+...+H_N(z).

3. Соединение обратной связи. Выходной сигнал первой системы подается на выход системы и одновременно на вход системы обратной связи, выходной сигнал которой суммируется, со знаком плюс или минус в зависимости от вида связи (отрицательной или положительной), с входным сигналом первой системы. Эквивалентная передаточная функция общей системы: H(z) = H₁(z)/(1± H₁(z)H₂(z)).

Схемы реализации систем. По принципам структурной реализации систем различают следующие схемы:

1. Прямая форма. Реализуется непосредственно по разностному уравнению

y_k =b_nx_k-n -a_my_k-m,

или по передаточной функции

H(z) =b_nzⁿ /(1+a_mz^m).

Пример прямой системы приведен на рис. 11.6.5.

2. Прямая каноническая форма содержит минимальное число элементов задержки. Передаточную функцию РЦС можно представить в виде:

H(z) = Y(z)/X(z) = H₁(z)H₂(z),

H₁(z) = V(z)/X(z) = 1/(1+a_mz^m), H₂(z) = Y(z)/V(z) =b_nzⁿ.

Отсюда:

v(k) = x(k) -a_mv(k-m), (11.6.1)

y(k) =b_nv(k-n). (11.6.2)

В разностных уравнениях (11.6.1-2) осуществляется только задержка сигналов v(k). Граф реализации РЦС приведен на рисунке 11.6.6.

3. Каскадная (последовательная) форма. Соответствует представлению передаточной функции в виде произведения:

H(z) =H_i(z).

H_i(z) - составляющие функции типа (1-r_iz)/(1-p_iz) при представлении H(z) в факторизованной форме, где r_i и p_i- нули и полюсы функции H(z). В качестве функций H_i(z) обычно используются передаточные функции биквадратных блоков - фильтров второго порядка:

H_i(z) = (b_0i+b_1iЧ z+b_2iЧ z²)/(1+a_1iЧ z+a_2iЧ z²).

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

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

литература

1. Баскаков С.И. Радиотехнические цепи и сигналы: Учебник для вузов. - М.: Высшая школа, 1988.

2. Бендат Дж., Пирсол А. Прикладной анализ случайных данных. – М.: Мир, 1989. – 540 с.

5. Гольденберг Л.М. и др. Цифровая обработка сигналов: Учебное пособие для вузов. - М.: Радио и связь, 1990.- 256 с.

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

14. Купер Дж., Макгиллем А. Вероятностные методы анализа сигналов и систем. – М.: Мир, 1989.

18. Оппенгейм А.В., Шафер Р.В. Цифровая обработка сигналов. – М.: Связь, 1979. – 416 с.

20. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. – М.: Мир, 1978.

Copyright ©2005 Davydov А.V.

к оглавлению

Знаете ли Вы, как разрешается парадокс Ольберса?
(Фотометрический парадокс, парадокс Ольберса - это один из парадоксов космологии, заключающийся в том, что во Вселенной, равномерно заполненной звёздами, яркость неба (в том числе ночного) должна быть примерно равна яркости солнечного диска. Это должно иметь место потому, что по любому направлению неба луч зрения рано или поздно упрется в поверхность звезды.
Иными словами парадос Ольберса заключается в том, что если Вселенная бесконечна, то черного неба мы не увидим, так как излучение дальних звезд будет суммироваться с излучением ближних, и небо должно иметь среднюю температуру фотосфер звезд. При поглощении света межзвездным веществом, оно будет разогреваться до температуры звездных фотосфер и излучать также ярко, как звезды. Однако в дело вступает явление "усталости света", открытое Эдвином Хабблом, который показал, что чем дальше от нас расположена галактика, тем больше становится красным свет ее излучения, то есть фотоны как бы "устают", отдают свою энергию межзвездной среде. На очень больших расстояниях галактики видны только в радиодиапазоне, так как их свет вовсе потерял энергию идя через бескрайние просторы Вселенной. Подробнее читайте в 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 -> Просвещение от Вильгельма Варкентина - Карим_Хайдаров.