книги из ГПНТБ / Сакрисон, Д. Лекции об аналоговой связи
.pdf6.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 |
|
>) Общепринятое обозначение для сходимости в средиеквадра-
тическом. — Прим, перев.
