рефераты конспекты курсовые дипломные лекции шпоры

Реферат Курсовая Конспект

Стационарные линейные системы

Стационарные линейные системы - раздел Математика, Сигналы И Линейные Системы ...

СИГНАЛЫ и ЛИНЕЙНЫЕ СИСТЕМЫ

Тема 11: СТАЦИОНАРНЫЕ ЛИНЕЙНЫЕ СИСТЕМЫ

Марк Туллий Цицерон. О гадании. Римский философ и политик, I в.д.н.э. Геология, наряду с медициной и богословием, относится к точным наукам. Есть геологический объект - должен быть сигнал.…

Введение

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

11.1. Линейные системы [1,2,9,14,18,20]

Общие понятия систем.Безотносительно к назначению и исполнению система всегда имеет вход, на который подается входной сигнал или входное воздействие, в общем случае многомерное, и выход, с которого снимается обработанный выходной сигнал. Если устройство системы и внутренние операции преобразований принципиального значения не имеют, то система в целом может восприниматься как “черный ящик”, в формализованном виде. Формализованная система представляет собой определенный системный оператор (алгоритм) преобразования входного сигнала – воздействия s(t), в сигнал на выходе системы y(t) – отклик или выходную реакцию системы. Символическое обозначение операции преобразования (трансформации):

y(t) = T[s(t)].

Системный операторT - это правило (набор правил, алгоритм) преобразования сигнала s(t) в сигнал y(t). Для общеизвестных операций преобразования сигналов применяются также расширенные символы операторов трансформации, где вторым символом и специальными индексами обозначается конкретный вид операции (как, например, TF - преобразование Фурье, TF-1 - обратное преобразование Фурье).

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

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

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

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

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

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

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

T[a(t)+b(t)] = T[a(t)]+T[b(t)].

Принцип однородностиили пропорционального подобия требует сохранения однозначности масштаба преобразования при любой амплитуде входного сигнала:

T[c ´ a(t)]= c ´ T[a(t)].

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

Примеры.

1. Система y(t) = a2t. y(t1) = a2t1, y(t2) = a2t2, y(ct) = a2ct.

y(t1+t2) = a2(t1+t2) = a2t1+a2t2 = y(t1)+y(t2). Система аддитивна.

cy(t) = ca2t = a2ct = y(сt). Система однородна. И в целом линейна.

2. Система y(t) = at2. y(t1) = at12, y(t2) = at22, y(ct) = a(ct)2 = ac2t2.

y(t1+t2) = a(t1+t2)2 ¹ y(t1)+y(t2)= at12+at22. Система не аддитивна.

с y(t) = с at2 ¹ y(сt) = ac2t2. Система неоднородна. И в целом нелинейна.

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

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

Рис. 11.1.1. Графика системных операций

y(t) = b ´ x(t),

y(t) = x(t-Dt),

y(t) = a(t)+b(t).

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

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

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

z(x,y) = c(x,y)×u(x,y),

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

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

s(x,t) = T[a(x,t)], T[a(x-Dx,t-Dt)] = s(x-Dx,t-Dt).

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

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

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

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

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

am= bn. (11.1.1)

При нормировке к ао = 1, отсюда следует:

y(t) =bn-am. (11.1.1')

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

Аналогичная связь в цифровой системы описывается разностными уравнениями:

am y((k-m)Dt) = bn s((k-n)Dt). (11.1.2)

y(kDt) =bn s((k-n)Dt) -am y((k-m)Dt). (11.1.2')

Последнее уравнение можно рассматривать как алгоритм последовательного вычисления значений y(kDt), k = 0,1,2, …, по значениям входного сигнала s(kDt) и предыдущих вычисленных значений y(kDt) при известных значениях коэффициентов am, bn и с учетом задании определенных начальных условий - значений s(kDt) и y(kDt) при k < 0. Интервал дискретизации в цифровых последовательностях отсчетов обычно принимается равным 1, т.к. выполняет только роль масштабного множителя.

Нерекурсивные цифровые системы. При нулевых значениях коэффициентов am уравнение (11.1.2') переходит в уравнение дискретной свертки x(k) с оператором bn:

y(k) =bn x(k-n). (11.1.3)

Рис. 11.1.2. Пример НЦС.

При установленных значениях коэффициентов bn значения выходных отсчетов свертки для любого аргумента k определяются текущим и "прошлыми" значениями входных отсчетов. Такая система называется нерекурсивной цифровой системой (НЦС). Пример простейшей НЦС приведен на рис. 11.1.2. Интервал суммирования по n получил название "окна" системы. Окно системы (11.1.3) составляет N+1 точку, система является односторонней каузальной, причинно обусловленной текущими и "прошлыми" значениями входного сигнала, выходной сигнал не опережает входного. Каузальная система может быть реализована аппаратно в реальном масштабе времени. При k<n проведение обработки входных данных возможно только при задании определенных начальных условий для точек x(-k), k=1,2,..,N. Как правило, в качестве начальных условий принимаются нулевые значения или значения отсчета х(0). Применяется также четное или нечетное продление функции x(k) на интервал отрицательных значений k. Если при обработке данных начальные интервалы массивов x(k) существенного значения не имеют, то обработку можно начинать с отсчета k=N.

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

y(k) =bn x(k-n). (11.1.3')

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

Техника выполнения свертки в координатной области не отличается от техники выполнения обычной дискретной свертки двух массивов данных.

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

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

Пример.Сглаживание данных скользящим симметричным П-окном размером 5 отсчетов.

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

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

Проверьте результат (выполните фильтрацию сигнала).

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

 

Рис. 11.1.3. Пример РЦС.

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

Рис. 11.1.4. Техника вычислений в РЦС.

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

Пример.Уравнение РЦС: yk= boxk+a1yk-1, при bo=a1=0.5, y-1=0. Входной сигнал: xk = {0, 10, 0, 0, 0,....}.

Расчет выходного сигнала:

уo = 0,5xo + 0,5y-1 = 0 y1 = 0,5x1 + 0,5yo = 5 y2 = 0,5x2 + 0,5y1 = 2.5 y3 = 0,5x3 + 0,5y2 = 1.25

y4 = 0,5x4 + 0,5y3 = 0.625 y5 = 0,5x5 + 0,5y4 = 0.3125 и т.д.

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

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

Пример.Уравнение РЦС: yk = boxk+a1yk-1+a2yk-2, при bo=a1=a2=1. xk = {0, 1, 0, 0, 0,....}, y-1 = y-2 = 0.

Расчет выходного сигнала:

yo = xo+y-1+y-2 = 0. y1 = x1+yo+y-1 = 1. y2 = x2+y1+yo = 1. y3 = x3+y2+y1 = 2.

y4 = x4+y3+y2 = 3. y5 = x5+y4+y3 = 5. y6 = x6+y5+y4 = 8. y7 = x7+y6+y5 = 13. и т.д.

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

Пример.Интегрирование данных по формуле трапеций.

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

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

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

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

11.2. Импульсная характеристика системы [1,5]

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

y(t) = T[d(t)] º h(t). (11.2.1)

y(kDt) = T[d(kDt)] º h(kDt). (11.2.1')

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

Рис. 11.2.1. Импульсный отклик системы h(t), входной сигнал s(t) и выходная реакция системы y(t).

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

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

На рисунке 11.2.1 приведен пример импульсного отклика h(t) интегрирующей RC-цепи. При подаче на вход RC-цепи импульса заряда Dq емкость С заряжается до напряжения Vо = Dq/C и начинает разряжаться через сопротивление R, при этом напряжение на ней изменяется по закону v(t) = Vo×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.

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

Усиление шумов.Критерием качества системы при использовании любого метода обработки информации можно считать выполнение целевого назначения с минимальным усилением шумов (максимальным их подавлением). Допустим, что система имеет нормированный к 1 импульсный отклик h(k). Обозначим через e(k) аддитивный шум с математическим ожиданием M{e(k)}= = 0 и дисперсией De = s2, который в сумме с сигналом поступает на вход системы. Значения 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)}]2} = M{[h(n) e(k-n)]2}.

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

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

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

M{h2(n)e2(n)} = h2(n) M{e2(n)} = De h2(n) = s2h2(n). (11.2.5)

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

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

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

Для систем с m входами и n выходами аналогично определяются парциальные импульсные отклики hij(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).

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

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

yo = xo+0.5 y-1 = 1+0 = 1 = ho. y1 = x1+0.5 yo = 0+0.5 = 0.5 = h1. y2 = x2+0.5 y1 = 0+0.25 = 0.25= h2

y3 = x3+0.5 y2 = 0.125 = h3. y4 = x4+0.5 y3 = 0.0625 = h4, т.д.

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

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

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

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

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

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

Y(z)amzm = X(z)bnzn, (11.3.1) где X(z),Y(z)- соответствующие z-образы входного и выходного сигнала. Из этого… H(z) = Y(z)/X(z) =bnzn(1+amzm). (11.3.2)

Пример.

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

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

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

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

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

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

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

am A(w)exp(jwkDt-jwmDt) =bn B(w)exp(jwkDt-jwnDt). A(w)exp(jwkDt)am exp(-jwmDt) = B(w)exp(jwkDt)bn exp(-jwnDt). A(w)am exp(-jwmDt) = B(w)bn exp(-jwnDt). (11.4.1)

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

Квазидетерминированный сигнал в какой-то мере позволяет оценить сохранение однозначности преобразования системой случайных сигналов. Допустим, что система имеет импульсный отклик h(t) = exp(-at), t ³ 0.… x(t) = A + cos(2t+j),

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

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

Wy(f) = Wx(f) |H(f)|2. (11.5.5)

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

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

Wxy(f) = Wx(f) H(f). (11.5.6)

Wyx(f) = Wx(f) H(-f).

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

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

H(f) = Wxy/Wx Û h(t).

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

sy2 = Ry(0) =Wx(f) |H(f)|2 df º Rx(0)h2(t) dt = sx2h2(t) dt, (11.5.7)

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

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

== Ry(0) º h2(t) dt ºWx(f) |H(f)|2 df, (11.5.8)

 

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

sy2 = - 2 º (-2)h2(t) dt. (11.5.9)

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

gxy2(f) = |Wxy(f)|2/[Wx(f)×Wy(f)]. (11.5.10)

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

0 £ gxy2(f) £ 1.

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

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

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

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

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

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

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

Графы систем. Наряду со структурной схемой система может быть представлена в виде графа, который отображает диаграмму прохождения сигналов и состоит… Пример структурной схемы системы с передаточной функцией H(z) =… Соединения систем. Различают следующие соединения систем (рис. 11.6.2-4). Рис. 11.6.2.

Литература

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

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

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

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

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

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

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

Главный сайт автора ~ Лекции по сигналам ~ Практикум

О замеченных опечатках, ошибках и предложениях по дополнению: davpro@yandex.ru. Буду благодарен.

Copyright ©2005 Davydov А.V.

 

– Конец работы –

Используемые теги: стационарные, ные, системы0.059

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

Что будем делать с полученным материалом:

Если этот материал оказался полезным для Вас, Вы можете сохранить его на свою страничку в социальных сетях:

Еще рефераты, курсовые, дипломные работы на эту тему:

Лекция 1. Тема: Операционная система. Определение. Уровни операционной системы. Функции операционных систем. 1. Понятие операционной системы
Понятие операционной системы... Причиной появления операционных систем была необходимость создания удобных в... Операционная система ОС это программное обеспечение которое реализует связь между прикладными программами и...

Микропроцессорные системы: система ДЦ-МПК, система "Юг"
Использован практический опыт внедрения линейных пунктов управления (ЛПУ) на 60 станциях в увязке с ЭЦ-4, ЭЦ-9, МРЦ-12, МРЦ-13. Выполнен переход на… В состав аппаратуры центрального пункта управления (ПУ) входят IBM-совместные… Круглосуточный режим работы аппаратных средств ПУ обеспечивается источниками бесперебойного питания, а также системой…

ТЕЛЕКОММУНИКАЦИОННЫЕ СИСТЕМЫ. СИГНАЛЫ И КАНАЛЫ ЭЛЕКТРИЧЕСКОЙ СВЯЗИ. СИСТЕМЫ СВЯЗИ С ЧАСТОТНЫМ РАЗДЕЛЕНИЕМ КАНАЛОВ. ЦИФРОВЫЕ СИСТЕМЫ ПЕРЕДАЧИ
Лабораторные работы часа... Практические занятия часа... Всего аудиторных занятий часов...

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

Совместные, несовместные линейные системы. Матричная запись линейной системы.
Совместные несовместные линейные системы Матричная запись линейной системы... Невырожденные линейные системы... Решение невырожденных систем по правилу Крамера и с помощью обратной матрицы...

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

Экспертные системы. Классификация экспертных систем. Разработка простейшей экспертной системы
Глава 2. Структура систем, основанных на знаниях. 1. Категории пользователей экспертных систем. 2.2. Подсистема приобретения знаний. 3. База… ЭС выдают советы, проводят анализ, дают консультации, ставят диагноз. Практическое применение ЭС на предприятиях способствует эффективности работы и повышению квалификации специалистов.

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

Нервова система (Нервная система)
У ц лому функц нервово системи полягають в нтеграц ц л сного орган зму включаючикоординац ю та регуляц ю функц й ус х його орган в та доц льн й… До складу входять нерви, нервов вузли та сплет ння,головний та спиннний мозок,… З п двищенням р вня орган зац у ф логенез тварин в дбувалася концентрац я нейрон в у певних частинах т ла, звичайно б…

0.04
Хотите получать на электронную почту самые свежие новости?
Education Insider Sample
Подпишитесь на Нашу рассылку
Наша политика приватности обеспечивает 100% безопасность и анонимность Ваших E-Mail
Реклама
Соответствующий теме материал
  • Похожее
  • По категориям
  • По работам