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

книги из ГПНТБ / Сакрисон, Д. Лекции об аналоговой связи

.pdf
Скачиваний:
22
Добавлен:
19.10.2023
Размер:
4.63 Mб
Скачать

6.3. Теорема кодирования источников

141

усредненное по этой совокупности качество кодирования можно сделать сколь угодно близким к величине, опре­ деляемой функцией Ra(d). Отсюда будет следовать, что

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

Для определения совокупности кодов введем распре­

деление в пространстве L -векторов. Положим

 

 

 

Ра(и) =

П Р а М ,

 

(6-65)

 

 

 

;=і

 

 

 

Рѵ(и|и) =

П Р ѵ(й/І«/).

(6-

 

раѵ (и, и) =

Рѵ (и I и) Ра (и),

(6.67)

 

Раѵ(Й) =

J PY(« \u)dPa{u),

(6.68)

 

 

 

41

L

 

Раѵ(5)=

J

 

 

 

Рѵ(5|и)гіРа(и) = Д Р аѵ(й/).

(6.6Э)

41X---X<u

 

'='

 

Рассмотрим множество кодов, т. е. наборов, состоя­ щих из М кодовых слов, или L -векторов, порождаемое

независимым выбором каждого кодового вектора со­ гласно распределению PaY(u). Тогда распределение слу­ чайной величины

L

 

D =

d (U, um) =

4- S

d (u i> “ ” )•

(6-70)

 

 

 

i= \

 

 

где

m e {1,

. . . , M} — индекс,

минимизирующий

d (U, u/i), зависит как от

распределения источника а,

так и от распределения PaY(u), порождающего рассмат­ риваемое множество кодов. Заметим, что распределение именно этой величины' (обозначим его Рс) определяет величину ошибки кодирования.

Найдем теперь связь между распределением вероят­ ностей Рс и вероятностной мерой РаѴ, согласно которой

.определяется наше множество кодон.

142 Гл. 6. Кодирование случайных источников сообщений

Л е м м а 6.1 (Галлагер [3, лемма 9.3.1]). Пусть R* и d* произвольные положительные числа, а csy — произ­ вольное распределение. Зададим множество

Л = {и, и: /аѵ(и,

u ) > L R *

или

d ( u , u ) > d * } .

(6.71)

Тогда

 

 

 

 

Рс ф >

сГ} < Раѵ (Л) +

exp { - Me-™*}-

(6-72)

Д о к а з а т е л ь с т в о .

Мы

следуем доказательству

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

Ли как множество векторов и, для которых и, и е Л :

Ли ={и: /ау(и, и) > L R * и л и cf (и, и) > cf*} (6.73)

(заметим, что множество Л„ измеримо). Рассмотрим

условную вероятность*)

 

 

 

Рс{£> >

d*[u} =

P{cf(u, um) >

d\ m== 1, 2, . . . ,

M} =

 

M

 

 

 

1- J dpav(ü)

M

 

=

n

p (d

u^) >

d*}=

.

(6.74)

Очевидно,

что для любых u, u e A cu

 

 

 

 

 

A av(u, u K e iR’,

 

(6.75)

так что из неравенства

(6.74)

следует, что

 

 

Р Д Я > с Г | и } <

1— е

LR' J Л аѵ (u, и) с?Раѵ (и)]Ж.

(6.76)

Воспользуемся неравенством

(1 — ß * ) * < l — x + e -w

•) Здесь и

ниже через

А с

обозначено дополнение к множе­

ству 4 . —

перев,

 

 

Прим,

 

 

6.3. Теорема кодирования источников

143

(см. Галлагер [3, (9.3.22), (9.3.23)]). Обозначив e~LR*

через ß, а через х — интеграл в (6.76), из последнего не­

равенства получим

РСФ > (u}< J А аѵ (и, и) dPay (ü) +exp {—M e -1**}. Au

Интегрируя обе части этого неравенства по мере Ра (и)

на множестве °U,

находим

 

РСФ >

I dPay(u, u) + exp{— M e-LR'} =

 

 

А

 

=

Pav {A} + exp { - M e-W }.

(6.72)

Покажем теперь, что величину U на выходе источ­

ника можно представить с точностью, произвольно мало отличающейся от d, так, что скорость создания цифровой информации будет сколь угодно близка к R(d). Приво­

димая ниже теорема, а также ее доказательство воспро­ изводятся из книги Галлагера [3].

Т е о р е м а 6.5

(Галлагер [3, теорема 9.6.2]). Пусть

Ra Ф) скорость

при заданном искажении, определяе­

мая распределением источника а. Тогда для любых d >

> 0 и

б > 0 существуют такое

(достаточно

большое)

число L

и такой (М + 1, Ь)-код,

что

 

 

M < e x p {L [tfa(d) + 6]}

(6.77)

 

Ee{ 5 } < d +

Ö.

(6.78)

Д о к а з а т е л ь с т в о (Галлагер). Пусть уо — распре­ деление, принадлежащее Ta (d) и такое, что

144 Гл. 6. Кодирование случайных источников сообщений

Рассмотрим множество (L, М)-кодов, порождаемое рас­ пределением Раѵо(й). и применим доказанную лемму, полагая

d' — d + у ,

R' = R a ( d ) + - J ,

 

M = exp{L[/?a(r f)+ |-5 ]}.

 

Тогда, согласно лемме,

 

 

 

 

•Рc { ö > d +

РаѴо(Л) +

ехр{— еМ/4},

(6.80)

где

 

 

 

 

 

PoYo C^) =

Р аѴо 1 1аУо (U! tl) > L

(Jd ) +

-g-j

 

 

или

 

 

 

 

^

^ аѵ» { ~L ^аъ (U’

U) ^

[^“Yo

 

 

 

 

+ p „Vo{ d > ( d + A ) } .

(6.8 I)

Используя независимость nap

(Ut, Oi), 1 = 1 , . . . ,

L, при

совместном распределении, описываемом соотношениями

(6.65) — (6.69), получаем

L

(^i

u) =

^

Iay„ (ul>

dl)

(6.82)

И

 

/=1

 

 

 

 

 

 

 

 

L

 

L

 

 

 

Д = т 2

0 ' =

т ^

^ &

(),

(6.83)

i=1

/=і

 

 

 

причем правые части этих равенств являются суммами независимых одинаково распределенных случайных ве­ личин. Кроме того, в силу задачи 6.2

Е а ѵ о ( I A r y , ,

І } ^ ~ {U',+

Ü )Л<г у „ ° ° »

E«vo {I Dt 1} =

E„y„ {Di) ^ d

< oo.

6.3. Теорема кодирования источников

145

Таким образом, применяя слабый закон больших чи­ сел [10], убеждаемся, что величину Ра% {Л} можно сде­ лать сколь угодно малой за счет выбора достаточно большого L. Учитывая неравенства (6.80) и (6.81), на­

ходим

Pc{ d > rf + | } < ß ( L ) ,

 

(6.84)

где ß(L) стремится к нулю при L - * оо.

 

 

З а м е ч а н и е . Мы уже на полпути

к завершению

доказательства теоремы. Действительно,

пусть

В — со­

бытие, состоящее в том, что D > d - f 6/2.

Тогда,

как мы

показали, вероятность Р (В) можно сделать сколь угодно

малой. Однако это не означает еще, что E { D } < d + 6, поскольку влияние «хвостов» распределения может ока­ заться существенным. Для завершения доказательства осталось учесть влияние этих «хвостов».

Из неравенства (6.84) вытекает, что

 

Еаѵ. 0 } < [d + I) + РаѴо(В) Еаѵ, ф I Я}.

(6.85)

Эта оценка справедлива для кода, состоящего из М слу­

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

точника поставить в соответствие добавленный нулевой вектор, то получим

 

L

Р|СГѴо (^) Еаѵ»

EaVo (В) EaVo {d(Ut, 0) Iß}. (6.86)

 

/ = 1

Рассмотрим произвольную положительную случай­ ную величину V и произвольное множество В\ выберем d0 так, чтобы

{ V > r f0}=o{ß}.

Тогда

P ( ß ) E { P |ß }= f o d P ( y ) < f vdP(v). (6.87)

В' v ^ d ,

146

Гл.

6.

Кодирование случайных источников сообщений

 

Объединяя

неравенства (6.86) и (6.87), находим, что

при V =

d(U, 0)

 

 

 

 

V d Р (V).

(6.88)

Так как E{V }s^/<g< со, то

(6.89)

Следовательно, правую часть неравенства (6.88) можно сделать меньше 6/2, если d0 взять достаточно большим. Выполнение этого условия можно обеспечить, положив L

достаточно большим, поскольку ß(L) —>-0 при L —*■ оо. Это замечание вместе с неравенством (6.85) завершает доказательство теоремы.

З а д а ч а 6.8. Предположим, что мы определили взаимную информацию по-иному, нежели в (6.26) или (6.27), используя вместо In (Л) другую функцию g( Л). Если бы при этом новом определении оба утверждения

(и позитивное, и негативное) теоремы кодирования вы­ полнялись с иной функцией Ra (d), то мы получили бы

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

6.4. С К О РО СТ Ь П РИ З А Д А Н Н О М И СК А Ж Е Н И И Д Л Я ГА У С С О В С К О ГО С Л У Ч А Й Н О Г О П Р О Ц Е С С А

СО В ЗВ Е Ш Е Н Н О Й К В А Д Р А Т И Ч Е С К О Й М Е РО Й И СК АЖ ЕН И Я

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

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

Другими словами, мы предполагаем, что U — выбо­ рочная функция гауссовского случайного процесса U{t) продолжительностью Т с, 0 ^ t Т. Будем также счи­ тать, что процесс V (t) стационарен, имеет нулевое ма­

6.4. Скорость при заданном искажении

147

тематическое ожидание и известную ковариационную функцию Ru{т) = E{Ut+xU t}.

Определим меру искажения, которую мы будем в дальнейшем рассматривать. Обозначим через А однород­ ную линейную систему, действующую в Ь2[0,Т], с им­ пульсной переходной функцией a{t) и частотной харак­ теристикой /4(f). Предположим, что a(t) Ф 0 лишь на

некотором конечном отрезке [0, Га]. Определим взвешен­

ную ошибку между выборочной

функцией U(t) на вы­

ходе источника сообщений и ее

аппроксимацией

O(t):

г

 

 

Е ( 0 = J a ( t - s ) [ U ( s ) - U ( s ) ] d s ,

0 < t ^ T + Ta.

(6.90)

о

 

 

Поскольку на отрезке Та ^ t ^ Т невозможны ложные

переходные эффекты, определим среднее расхождение как

(6.91)

где Т' = Т - Т а.

Весовая функция A(f) была введена для того, чтобы

приспособить меру искажения к учету субъективных оценок ошибки сигнала O(t). Например, если u(t)

сигнал, воспроизводящий музыкальное звучание, то естественно потребовать, чтобы величина |/4(f)|2 была, грубо говоря, обратно пропорциональна S u(f), т. е.

чтобы высокие звуки, имеющие обычно слабую мощ­ ность, звучали более отчетливо. С другой стороны, вне диапазона слышимости |/4(f) |2 может быстро спадать до нуля.

Если же u(t)— развертка телевизионного изображе­ ния, то, вероятно, следует предположить, что выбор a{t)

должен способствовать аппроксимации как самого им­ пульса, так и его производной, чтобы,, во-первых, улав­ ливать изменение интенсивности и, во-вторых, сохранить четкость изображения предметов. Функция |/4(f) |2 мо­ жет, как и раньше, быстро спадать до нуля вне области разрешения применяемой оптической системы.

Опишем подход к нахожденидо R(d) в рассматривае­

мом случае. Он состоит в таком представлении процесса

148 Гл. 6. Кодирование случайных источников сообщений

U(t), O ^ t ^ T , счетным набором независимых гаус­

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

Для этого определим линейные интегральные опера­ торы А и А* равенствами

 

I a(t — s)u (s) ds,

Ta^ . t ^ . T ,

 

 

 

 

(6.92)

 

в остальных случаях,

 

т

 

 

 

[Л*ы] (s) =

J a{t — s)u{t)dt,

0 < s < 7 \

а

 

 

(6.93)

 

 

 

 

)

в остальных

случаях.

Зададим случайный процесс

V{i),

Та ^

t ^ Т, положив

 

V(t) = [AU](t).

(6.94)

В введенных операторных обозначениях расхождение

(6.91) принимает вид

 

гі = Е { ^ | | Л ( Д - Д ) | | [ г а.г ]} .

(6.91)

Обозначим через R v оператор A o R u ° A *, т.

е. инте­

гральный оператор в L2 [Та, Т] с ядром

 

=

Е іУ Л ,

Ta^ t {, t 2^ T .

(6.95)

Поскольку процесс U (t)

стационарен в квадрате

Та ^

^-Т,

функция Rv

зависит лишь от разности

U t2.

Рассмотрим базис разложения Карунена — Лоэва для V(t) на [Та, Т], т. е. введем систему собственных

функций и собственных значений уравнения

6.4. Скорость при заданном искажении

149

С помощью {ср/i} построим также систему функций для разложения U(t). Рассмотрим оператор А* о А, действую­ щий из Ь2[О, Т] в La [О, Т]. Этот оператор симметричен, и

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

Л’фй = -ц - A*RVфй = А*ARuA*yk .

Таким образом, функция А*щ принадлежит множеству

значений оператора

А*о А,

и,

следовательно, функции

0fc (s) = [А* оЛ]“ 1А ’щ = ±

R„A*Фь

 

(6.97)

k

1, 2, . . . , s g [0,

Т],

 

 

принадлежат L 2[0,

Г]. Заметим, что

 

 

 

[ЛѲЛ] (0 =

 

(0 =

Ф& (О,

^ [ Г 0,

Т).

(6.98)

Возьмем теперь произвольную функцию и в области

определения оператора А и положим

 

 

 

v{t) = [Au](t),

 

Ta< t ^ T ,

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

vN( t ) = y i

vkqk 00,

 

 

 

 

 

ft=l

 

 

 

 

 

 

Vk =

{v,

<Pfc},

 

 

 

 

uN (s )=

N

vkBk (s),

О < s < Т.

 

 

s

 

 

 

ft—1

 

 

 

 

 

В силу (6.98)

 

 

 

 

 

 

 

 

[Лил,](0 =

0Лг(О,

 

 

 

а в силу полноты системы {фй} в L2[Tay Т]

 

 

О = lim II V vN ||rr , г] =

Пт ИА (и uN) ||г

ту

(6.99)

Итак, при взвешенном квадратическом критерии оши­ бок любую функцию в области определения оператора А

можно представить с помощью набора {Ѳ/J. Более того, если оператор А ограничен снизу (по норме) числом k, то в силу (6.99) un сходится к и в смысле среднеквадра­

150 Гл. 6. Кодирование случайных источников сообщений

тической интегральной ошибки, и в этом случае функ­ ции ѲЛ порождают область определения оператора А.

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

lim U N (s) =

 

N

 

0 < s < T, (6.100)

lim

2

(s),

0

N->oo

k—\

 

 

где коэффициенты Vk определяются

равенством

 

 

 

т

 

Vk = ( V ,

<pft) =

f dtV(t)Vk(t).

T'a

Полагая

Ta^ t ^ T ,

k=i

в силу разложения Карунена — Лоэва получаем

У(/) = 1.і.т. Ѵлг (/) ‘)

равномерно для всех t ^ [ T a, Т]. Таким образом,

О = lim Е {II V — VN lljr , ri} =

= lim E { ИЛ (£/ Ufj) \\[т, ri}-

(6.101)

Если оператор А ограничен снизу, то из (6.101) сле­ дует, что UN(s) сходится к U{s) и в. смысле среднеквад­

ратической интегральной ошибки.

Выразим величину искажения через координаты Vh-

Для этого представим процесс источника и закодиро­ ванный процесс в виде

U(s) = I t VkQk(s),

t / ( s ) = S ^ 0 * ( s ) ,

0 < в < Г .

ft=i

ft=i

 

>) Общепринятое обозначение для сходимости в средиеквадра-

тическом. — Прим, перев.

Соседние файлы в папке книги из ГПНТБ