Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Электроника и схемотехника, учебное пособие, Щ.А.С., О.А.Г

..pdf
Скачиваний:
61
Добавлен:
13.10.2024
Размер:
17.78 Mб
Скачать

431

Рис. 9.12

Работа такого оператора описывается выражением

Н- 1 : x(t) T x(n) (t) , nT t (n 1)T .

(9.32)

Формируемая им последовательность импульсов аппроксимирует x(t)

следующим образом:

x(t) T x(t)

(t nT ) T x(nT ) (t nT )

.

(9.33)

n 0

n 0

 

 

 

 

nT

t (n 1)T

Таким образом, результирующее колебание на каждом интервале

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

Фиксация нулевого порядка (H0 ). Наиболее распространенным восстанавливаю- щим оператором является схема фиксации нулевого порядка (рис. 9.12, б). Ее работа опи-

 

432

 

сывается выражением

 

 

H0 :

x(t) = x(n), пТ <t ≤ (п+1)Т.

(9.34)

Экстраполирующая фиксация первого порядка (H1 ). Работа этой цепи описывает-

ся следующим выражением:

H1 :

x(t) = x(n)+[x(n)-x(n -1)] , пТ <t ≤ (п+1)Т.

(9.36)

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

го интервала использовать разность x(n)-x(n -1) (рис. 9.12, в).

Многоугольная фиксация (Hр ). Этот оператор, представляющий собой линейный интерполятор, выполняет операцию (рис. 9.12, г).

Hр :

x(t) = x(n -1)+[x(n)-x(n -1)] (t-nT)/T, пТ <t ≤ (п+1)Т .

(9.37)

Данный оператор является физически реализуемым вариантом многоугольной фиксации. Нереализуемый вариант описывается выражением x(t) = x(n)+[x(n+1)-

(n)] (t-nT)/T, пТ <t ≤ (п+1)Т .и может использоваться при расчете цифровых фильтров для аппроксимации входного сигнала.

9.5.6. Частотная характеристика фиксирующих операторов

Спектральные характеристики фиксирующих операторов могут быть исследованы с помощью преобразования по Лапласу импульсной характеристики оператора, в кото-

рую вместо s следует подставить j . Нетрудно найти преобразования по Лапласу рассмотренных выше операторов:

 

H 1 (s)

T

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H 0 (s)

1

e sT

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H1

(s)

1

Ts

 

1

e sT

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

T

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H p (s)

 

1 1 e sT

 

 

1

[H

 

(s)]2

 

T

 

 

 

s

 

 

 

 

T

0

 

 

 

 

 

 

 

 

 

 

 

 

.

(9.38, а)

(9.38, б)

(9.38, в)

(9.38, г)

433

Подставляя s j , находим частотные характеристики

 

H 1 (

)

 

 

T

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

T

 

 

 

j T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H0 (

)

 

sin

 

 

e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin( T / 2

2

H

(s)

T

1

 

( T )

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T / 2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H p ( ) T

sin(

T / 2)

 

 

 

 

 

 

 

 

 

 

 

 

T / 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.39, а)

(9.39, б)

(9.39, в)

(9.39, г)

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

приведены на рис. 9.13. Из этих графиков видно, что многоугольный оператор Нр имеет плоскую фазовую характеристику, приводящую к задержке на Т секунд; следовательно, такой оператор-, удобен для использования в разомкнутых восстанавливающих схемах (без обратных связей).

434

435

Рис. 9.13

Фиксирующий оператор нулевого порядка H 0 вносит задержку T/2 секунд. Экстрапо-

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

тя оператор импульсной аппроксимации в диапазоне частот 0 /T имеет идеаль-

ную характеристику, на частотах, превышающих / T , последняя отличается от желаемого нулевого значения.

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

один отсчѐт из Asin(

t) берѐтся в момент времени t

0 , а следующий – в момент време-

ни t

T .

 

 

 

 

 

 

 

 

 

Погрешность

 

x

y равна

Asin( t) ,

0 t

T . Из последнего выражения ясно,

что

MAKC

A T

. Пусть

представляет собой максимальную погрешность, выражен-

 

 

 

 

 

 

 

 

 

 

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

628 fT , %.

 

Если обозначить число отсчѐтов за период через

p ( f T ) 1

 

 

, отсчѐт/период, то полу-

чим соотношение p

628/ p , %, или p 628/

, отсчѐт/период. В табл. 9.2 приведено

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

 

 

 

Таблица 9.4

 

 

 

 

 

Точность, %

Фиксация нулевого

Фиксация первого

 

Фиксация мно-

 

порядка

порядка

 

гоугольника

 

 

 

 

 

0,01

62800

628

 

223

 

 

 

 

 

0,1

6280

199

 

71

 

 

 

 

 

0.,5

1256

89

 

32

 

 

 

 

 

436

1,0

628

63

22

 

 

 

 

В случае экстраполирующей фиксации первого порядка максимальная погрешность

имеет место вблизи пикового значения, поэтому ради удобства положим x

Asin( T ).То-

гда для экстраполирующей функции первого порядка

2 A [1

cos( t)]

A

2t 2

 

 

 

и, сле-

 

MAKC

A

2T 2

100 ( T )2

p 62,8 (

) 1 / 2

 

 

довательно,

 

или

, % или

 

, отсчѐт/период. В

 

 

 

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

9.6. Линейные дискретные системы

Обработка дискретных сигналов линейными фильтрами или весовыми последовательностями производилась еще в 1600 г. Она началась с работ математиков, занимавшихся составлением математических таблиц, и астрономов, которые занимались определением орбит небесных светил. Труды Найпера, Грегори, Ньютона и Бернулли, Эйлера, Тейлора, Лагранжа, Лапласа и Гаусса используются в классических методах численного анализа, применяемых и в настоящее время для численного интегрирования, интерполяции, дифференцирования и т. д. Применения, в которых обычно используют узкополосные сигналы при отсутствии шумов, вытекают из характерных особенностей указанных трудов, обосновывающих полезные способы фильтрации, что становится особенно ясным при их рассмотрении в частотной области. Например, легко показать, что для классиче-

ских симметричных интерполяционных фильтров первые п последовательных производных их амплитудно-частотных характеристик имеют нулевые значения (поведение при нулевой частоте весьма сходно с поведением фильтра нижних частот Баттерворта).

Однако полезность классических методов снижается, если обрабатываемые данные содержат нежелательные шумовые компоненты. Развитие в 1940 и 1950 гг. системы управления огнем и радиолокации дало сильный толчок развитию более совершенных методов обработки сигналов. Появилась богатая литература по дискретным системам управления, начиная с работы Гуревича [3] и др. и кончая потоком литературы по дискретным системам. Эти работы, рассматривая подробно теорию дискретных систем с помощью z- преобразования, в основном ограничиваются анализом систем управления и исследованием методов их оптимизации. Они могут быть использованы для решения многих задач цифровой фильтрации.

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

437

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

9.6.1. Передаточная функция импульсной системы

Реальный импульсный элемент (ИЭ) удобно заменить двумя элементами: идеальным импульсным элементом (ИИЭ) и фиксирующим элементом (ФЭ). Идеальный импульсный элемент замыкает систему на «мгновение» в момент съѐма импульса x(t) , а фиксирующий элемент задерживает (фиксирует) снятое значение импульса x(t) в течение времени T . ФЭ будет далее относится

кнепрерывной части системы. При таком представлении входной сигнал

x(t) ИИЭ оказывается последовательностью дельта-функций, площадь каж-

дой из которых принимается равной значению непрерывного сигнала x(t) на

входе ИИЭ в момент съѐма импульсов, которые можно пронумеровать. Такая последовательность дельта-функций называется решѐтчатой функцией

x (t)

x(t) (t nT )

x(nT ) (t nT ) .

 

n 0

n 0

(9.40)

На рис. 9.14 изображѐн простейший случай, когда x(t) представляет собой функцию Хевисайда вида x(t) 10 (t) , а выходной сигнал x (t) ИИЭ равен

x (t) (t nT ) и является периодической последовательностью дельта-

n 0

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

438

Рис. 9.14

Выведем передаточную функцию импульсной системы (рис. 9.15), непрерывная часть которой имеет функцию веса w(t).

Рис. 9.15

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

фиктивный импульсный элемент, работающий синхронно и синфазно с входным импульсным элементом, его выходной сигнал y(t) также является импульсным и существует только в моменты съѐма импульсов.

Передаточная функция импульсной системы, получаемая путѐм свѐртки и z-преобразования, действительна только для моментов съѐма. Используя принцип суперпозиции, запишем значение выходного сигнала y(t) в произ-

439

вольный момент съѐма t=nT в виде

y(nT ) x(0) w(nT ) x(T ) w(n 1)T x(2T ) w(n 2)T ... x(nT ) w(0) , (9.41)

где w(nT) - импульсная переходная функция или функция веса системы. Об-

щий вид слагаемых реакции y(nT)

x(mT ) w[(n m) T ], m n ,

поэтому

n

y(nT ) x(mT ) w[(n m) T .

m 0

(9.42)

Это – формула свѐртки для импульсных функций, аналогичная формуле свѐртки для непрерывных функций. Решѐтчатая функция имеет вид

 

 

n

y (t)

y(t) (t nT )

x(mT ) w[(n m)T (t nT ) .

 

n 0

n 0 m 0

(9.43)

Применяя преобразование Лапласа к решѐтчатой функции, получим

L[ y (t)]

x(mT ) w[(n m)T ] (t nT ) e st dt

0 n 0

m 0

x(mT ) w[(n m)T ] (t nT ) e st dt.

0 m 0

n 0

В силу фильтрующего свойства дельта-функции получим

L[ y (t)]

x(mT ) w[(n m)T ] e nTs

 

x(mT ) e mTs

m 0

n 0

m 0

n 0

w[(n m)T ] e (n m)Ts .

(9.44)

Введѐм обозначения

z esT , X (z)

x(mT ) e mz , W (z)

w(kT) e kz , k n m .

n 0

 

k 0

(9.45)

Тогда выражение (9.44) примет вид

W (z)

440

L[ y (t)] X (z) W (z) ,

(9.46)

 

 

где W (z)

Y (z)

.

 

 

X (z)

(9.47)

 

 

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

нейной системы. Передаточная функция замкнутой системы (рис. 9.16) получается из еѐ уравнений

X (z) X BX (z) Y (z) , Y (z) W (z) X (z) .

Рис. 9.16

Исключая из уравнений X (z) , получим передаточную функцию замкну-

той системы

W0 (z) 1 W (z) .

(9.48)

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

 

 

1

e

sT

H

0 (s)

 

 

 

 

.

 

 

s