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

ЦСУ_Егоров

.pdf
Скачиваний:
16
Добавлен:
12.02.2016
Размер:
2.4 Mб
Скачать

х (t)

4

 

 

3

 

 

2

 

 

1

1 2 3 4 5 6 7 8 9 10 11

t

 

 

Рисунок 3.1 Непрерывный сигнал

Дискретный по времени непрерывный сигнал (рис. 3.2). Значе- ния сигнала определены только в определенный момент времени и нахо- дятся во всем непрерывно допустимом числовом диапазоне.

х [nT]

5

4

3

2

1 2 3 4 5 6 7 8 9 10 11 nT

Рисунок 3.2 Непрерывный сигнал дискретизированный по времени.

T - период отсчета сигнала (период дискретизации) n=0… .

Дискретный сигнал непрерывного времени (квантование по уровню) (рис. 3.3). Вся шкала возможных значений сигнала соответствует определенным уровням. Отсчеты сигнала производятся по времени непре- рывно.

х (t)

1

2

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 3 4 5 6 7 8 9 10 11 12

Рисунок 3.3 Сигнал, квантованный по уровню. 1-исходный сигнал, 2-квантовый сигнал

Дискретный сигнал по уровню и по времени (рис. 3.4.)

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

22

х[nT]

7

 

6

 

5

1

4

 

3

2

2

 

1

 

1 2 3 4 5 6 7 8 9 10 11 12 nT

Рисунок 3.4 Дискретный сигнал по уровню и по времени 1-исходный непрерывный сигнал, 2-дискретный сигнал по уровню и по времени

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

Цифровым называется дискретный сигнал, закодированный с по- мощью двоичной или двоично-десятичной системой.

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

Под дискретизацией сигнала по времени понимают преобразование модели сигнала функции непрерывного времени в функцию дискретного времени.

Под квантованием сигнала по уровню подразумевают преобразова- ние некоторой величины с непрерывной шкалой в величину, имеющую дискретную шкалу значений.

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

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

23

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

На рисунке 3.5 показан пример преобразования непрерывного сиг- нала в цифровой код. Задано 8 уровней квантования в двоичной системе счисления с дискретностью отсчетов 1 с. Значения сигнала считываются с интервалом в 1 с. и полученным значениям присваивается двоичный код.

Рассмотрим более подробно операции дискретизации сигнала по времени и квантования. Для получения цифрового сигнала из непрерывно- го необходимо взять его отсчет в дискретные моменты времени t=nT, где t - текущее время, Т - интервалы между отсчетами, n - порядковый номер отсчета. Если полученный дискретный сигнал представить в функции дис- кретного времени, то полученная функция называется решетчатой (рис.3.6). В интервалах между отсчетами значение решетчатой функции равно нулю.

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

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

X(t)

111

110

101

100

011

010

001

000

1

2

3

4

5

6

7

t

X(t)

 

Цифровой код

 

 

Уровень “1”

 

 

 

 

 

 

 

Уровень “0”

2

3

4

5

6

7

t

1

Рисунок 3.5 Преобразование непрерывного сигнала в цифровой код

24

Рисунок 3.6 Решетчатая функция

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

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

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

Теорема устанавливает принципиальную возможность полного вос- становления непрерывной функции с ограниченным спектром по ее отсче- там и указывает предельное значение интервала времени между отсчетами, при котором такое восстановление еще возможно. Формулируется она сле- дующим образом: Функция f(t), допускающая преобразование Фурье и имеющая непрерывный спектр, ограниченный полосой частот от f=0 до f= fmax где fmaxmax/2π - максимальная частотная составляющая спектра, пол- ностью определенная дискретным рядом своих мгновенных значений, вы- численных через интервал времени

T =

1

.

(3.1)

 

2 f max

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

25

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

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

Если воспользоваться теоремой В.А. Котельникова при выборе ша- га дискретизации, то вероятность случайного выброса функции на интер- вале времени Т является малой величиной.

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

Предположим, что для некоторого процесса возможное изменение сигнала в частотном спектре составит 0≤ f ≤1000 Гц. В соответствии с тео- ремой В.А.Котельникова интервал дискретизации следует взять ∆Т ≤ 0,5 мс.

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

Теперь если выбрать в fmax=fгр то в соответствии с теоремой Котель- никова частота дискретизации будет fд=2÷4 f гр.

Рассмотрим операцию квантования непрерывной функции по уров- ню. Поскольку числовое значение каждого отсчета должно быть выражено некоторым числом конечных цифр, то приближенно описать бесконечную последовательность возможных значений непрерывного процесса можно с помощью конечного числа уровней квантования. Как бы ни была точна шкала, необходимо сделать выбор между двумя соседними значениями. На рис.3.7 изображен процесс квантования непрерывного сигнала в соответ- ствии с разрешенными уровнями, происходящий в аналого-цифровом пре- образователе АЦП.

26

Рисунок 3.7 Квантование непрерывного сигнала:

а - непрерывный входной сигнал АЦП; б - выходной сигнал АЦП

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

Таким образом, если квантование выполнено верно, то истинным значениям исходного непрерывного сигнала будут соответствовать наибо- лее близкие к ним уровни квантования. Точность приближения к непре- рывному процессу зависит от числа уровней квантования. В современных цифро-аналоговых преобразователях сигнал на выходе представляется в двоичном коде и имеет 1024 уровня квантования (десятиразрядное АЦП) и более. При идеальном преобразовании ошибка преобразования распре- делена равномерно с отклонением σх 0,29 х, где х шаг кантова- ния. В этом легко убедиться. Пусть р(х) - плотность вероятности ошибки равномерного квантования

 

1

 

0 ,5 x x 0 ,5

 

 

 

p(x) =

 

при

x

,

(3.2)

 

 

x

 

 

 

 

0 при других значениях х

 

 

 

где х текущее значение сигнала, соответствующее определенному уровню квантования (см. рис. 3.7). Очевидно, что среднее значение μх равно нулю, поскольку плотность р (х) симметрична относительно значения уровня. Тогда дисперсия ошибки

x 2 p(x)dx =

1

0 ,5

x

 

0 ,25

x

3

 

σ x2 =

 

x 2dx =

 

(3.3)

 

 

 

 

 

−∞

 

 

 

x

0 ,5

x

3

x

 

 

и отклонение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.4)

σ x =

x

2

0.29 x .

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

27

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

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

Как было показано выше при цифровой обработке информации не- прерывный сигнал в АЦП, преобразуется в последовательность чисел, ко- торую можно рассматривать как решетчатую функцию, аргументом кото- рой является дискретное время [nT], где n = 0,1,2,3..., Т период дискре- тизации.

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

3.2. Частотная форма представления детерминированных сигналов

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

ций ψk (t) [1]

f ( x ) = Ckψ k (t ),

(3.5)

k

где Сk- базисные коэффициенты; [t1,t2]-интервал существования сигнала.

Совокупность базисных коэффициентов Сk называется дискретным спектром сигнала. Вычисление спектральных составляющих сигнала су- щественно облегчается при выборе в качестве базиса системы ортогональ- ных функций. Систему функций ψk, ψm называют ортогональной на отрез- ке [t1, t2], если для всех k=0…n, m=0…n за исключением k=m удовлетворя- ется условие:

t2

 

ψ k ( t )ψ m ( t ) = 0 .

(3.6)

t1

 

Эта система функций будет ортонормированной, если для всех k=m справедливо соотношение:

28

t2

 

ψ m2 ( t )dt = 1 .

(3.7)

t1

При анализе систем автоматического управления широко использу- ется представление детерминированных сигналов с применением базисных функций ept. При p=±jω-преобразование Фурье, а при p=s±jω- преобразование Лапласа.

Использование экспоненциальных базисных функций в преобразо- вании Фурье комплексно сопряженными парами (с положительным и от- рицательным параметром ω) позволяет в соответствии со следствием фор- мулы Эйлера:

e jω + e

jω

 

eiω e

iω

 

 

= cosωt;

j

 

 

= sin ωt (3.8)

2

 

2

 

 

 

 

 

 

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

Рассмотрим разложение в ряд Фурье периодических функций. Пусть функция х(t) , заданная на интервале времени t1t t2 и удо-

влетворяющая условиям Дирихле повторяется с периодом T = = t2-t1,

ω

на протяжении времени от − ∞ до +∞, где ω - угловая частота изменения этого сигнала. Эта функция может быть представлена в виде ряда

 

x(t) = Ck e jkωt .

(3.9)

k =−∞

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

Обозначим

Ck

=

ak + jbk

;Ck

=

ak jbk

,

 

 

 

(3.10)

 

 

 

 

 

 

 

 

 

 

тогда

2

 

 

2

 

 

 

 

 

 

 

 

 

e jkωt + ejkωt

 

e jkωt ejkωt

 

 

 

 

 

 

 

x( t ) = C0 +

Cke jkωt

+ Ckejkωt

= C0 + ak

 

 

+ bk j

 

.

(3.11)

2

2

 

k=1

 

 

 

k=1

 

 

 

В соответствии с формулой Эйлера получим окончательное выра- жение для ряда Фурье в тригонометрическом виде:

29

 

x(t) = a0 + ak cos(kωt) + bk sin (kωt)

(3.12)

k=1

Вразложении периодического сигнала в ряд участвуют гармоники с начальной частотой ω, дискретно изменяющейся при изменении k от 0 до ∞.

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

выберем ψ m (t ).

ψ m

(t ) = e

jm ω t

.

(3.13)

 

 

 

Действительно, при k m интеграл произведения базисных функ-

ций

t2

e jkωt e jmωt dt = 0 .

t1

А при значениях k=m под интегральная функция равна 1. Опреде- ленный интеграл на интервале интегрирования, равным периоду функции равен T.

Умножим правую и левую части ряда на вторую базисную функ- цию еjkωt и осуществим интегрирование на интервале t1, t2

t2

 

 

 

x(t )e

jmωt dt =

C k e jkωt e

jmωt dt ,

(3.14)

t1

 

k =−∞

 

 

при m=k получим:

t

2 x (t )e jk ωt dt = TC k ,

t1

C k

=

1

 

t 2

x (t )e jk ω t dt .

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

t

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем выражение (3.16) в виде:

 

 

 

 

1 t2

 

 

e jkωt + e jkωt

 

e jkωt e

jkωt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C k

=

 

 

x(t )

2

j

 

dt =

 

T

t1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 t2 [x(t )cos(kωt )dt jx(t )sin(kωt )]dt .

T t1

В соответствии с обозначением

C k = α k + jbk ,

2

определим aк, bк:

(3.15)

(3.16)

(3.17)

(3.18)

30

t2

αk = T2 x(t )cos(kωt )dt , (3.19)

t1

t2

bk = T2 x(t )sin(kωt )dt . (3.20)

t1

При k=0 получим, в соответствии с (3.16), составляющую

t2

C0 = 1 x(t )dt (3.21)

Т t1

и через принятое обозначение:

α =

2

 

t2

x(t )dt ,

a = 2C .

 

 

(3.22)

 

0

T

 

0

0

 

 

 

 

 

 

 

 

 

t1

 

 

 

 

Окончательно в тригонометрическом виде ряд периодической

функции запишется так:

 

 

 

 

 

 

 

 

 

 

 

x(t ) = α 0 + [α k cos(kωt )+ bk sin(kωt )].

(3.23)

 

2

 

k =1

 

 

 

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

При записи ряда Фурье в комплексной форме функцию

C k

=

1

t 2

x (t )e jk ω t dt

(3.24)

T

 

 

 

t

1

 

 

 

 

 

 

 

 

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

Обозначим в выражении (3.17)

 

1

t2

 

 

 

 

 

 

1

t2

 

 

× x(t) × cos(kωt)dt

через Аk(kω) , а

j ×

× x(t) ×sin(kωt)dt

че-

 

 

T

T

 

 

t

 

 

 

 

 

 

t

 

 

 

 

1

 

 

 

 

 

 

 

1

 

рез

Вk(kω),

тогда коэффициент Сk

в

комплексном виде

равен,

Ck

= Ak (kω) + jBk (kω)

и модуль

этой

комплексной функции

A(kω ) определяет значение частотного спектра

 

 

 

 

A(kω) =

 

 

.

 

 

 

 

(3.25)

 

 

 

Ak2 (kω) + Bk2 (kω)

 

 

 

 

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

31