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

книги из ГПНТБ / Филиппов, Л. И. Основы теории радиоприема дискретных сигналов (синтез и анализ) [монография]

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

1. 3. Основные методы отождествлений. Прямой метод

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

При прямом зондировании канал моделируется импульсным откликом h (v, t), где t — момент воздействия 8-импульса или другого испытательного импульса. Для многих практически ис­ пользуемых каналов преобразование Фурье Ft [h (у, t) ] при у= const сосредоточено в ограниченной полосе частот А /Лв окрест­ ности нулевой частоты. При этом h (t, у) можно представить в виде

h(u, t) = ^ lh(jih, v) sin [(<ь/тд) (t jzh)]

(*hh) (t — i*h)

 

 

3

 

 

sin [(tv/ta) (t j4 )]

 

( 1. 21)

(тс/тЛ) [l — jzh)

>

 

3

 

 

где

 

( 1. 22)

h. (v) = h (;ЧЛ, v); тЛ=1/Д /Л.

 

Из (1. 21) следует что h (v, t) полностью описывается своими мгно­ венными значениями h3 (у), отсчитанными через интервал времени

тл= 1 /А /л-

Если входное колебание есть последовательность зондирующих импульсов

«1V) =

2

8 (t — foA),

 

 

S2 (t) = J h(t, y) sx (t v)dv =

 

 

 

 

sin [(njтд) (t jzh)\

— ftt*).

(1.23)

Ыч)

]'ч)

 

 

Пусть продолжительность

отклика не превышает

величины

\Ч1|1 («память канала») и

 

 

 

(1-24)

 

 

 

 

Такой канал будем называть каналом с медленно изменяющимися параметрами. При этом в сумме (1. 23) по / слагаемые можно считать постоянными, равными единице для k = j и равными нулю для kj^j. Таким образом,

=

С125)

3

3

31

а,.,

На выходе канала получается последовательность выборок им­ пульсного отклика.

Практическое использование прямого метода затруднено огра­ ничением максимальной (пиковой) мощности зондирующих пере­ датчиков.

1.4.Метод сравнения с физической моделью

1.В этом методе предполагается наличие четырехполюсника

с откликом йи, зависящем от множества параметров {a f}, так что

/ги = /гм (i, (Х0, ос-р а2, •••)>

причем параметры (af) можно изменять путем внешнего воздей­ ствия.

Отождествляемая система с импульсным откликом h (t) и четырехполюсник hK (t, a0, a1; . . .) подвергаются воздействию одинаковых входных колебаний % (t), а разность выходных коле­ баний

е =

(*)-*,*(£ )

(1.26)

используется для управления параметрами ( а;} до получения минимально возможной величины е. После этого параметры ото­ ждествляемой системы считаются равными параметрам ( а4}.

Может быть использовано большое разнообразие моделей, с которыми производится сравнение. Эйкгофом предложена

обобщенная модель и введена обобщенная ошибка. Обобщенная модель является линейной системой с набором фильтров {hi (t)} и [fj (0} (четырехполюсников с соответствующими импульсными откликами), потенциометров и сумматоров (рис. 1. 1, а).

Ошибка е в обобщенной схеме равна:

е = 2

aflt + 2

Ру*>у = — si (0 +

2

« Л + 2

Ру«у■ (1•27)

1 = 0п

у=О

3 3

<=1

/= о

 

Последнее равенство написано с учетом того, что один из параме­ тров в данном случае а 0, может быть выбран произвольно, например, a0= —1. Параметры {а.} и {(34} регулируются до по­ лучения минимума некоторого критерия Е, являющегося функцио­ налом обобщенной ошибки е: E = F (е). Если положить |30= —1, Рх= р2= . . . = (Зт = 0 , то модель упрощается и обобщенная опшбка переходит в разность выходов отождествляемой системы и модели

(рис. 1. 1, б).

Из (1. 27) видно, что обобщенная ошибка линейна относительно

каждого из параметров {af} и

(^.}, так что

 

де.

дь

(1. 28)

да{

 

 

 

Отсюда следуют две возможности осуществления метода: с исполь­ зованием физической модели и с численным решением уравнений.

32

о

В

В первом случае конструируется схема авторегулирования с надлежащим выбором зависимостей между де/да. и а{, а также между deldfij и (З^.. Во втором случае производится решение системы линейных уравнений де/да.=0, де/д$.=0, в которых неизвест­ ными являются параметры (af), {(ЗД.

2. Для того чтобы минимизировать величину функции Е, необходимо знание производных по регулируемым параметрам. Они могут быть получены способом решения «уравнений чувстви­ тельности», путем использования двух моделей с различными параметрами и одной модели с изменяющимися параметрами.

Пусть система описывается набором уравнений

~df~ = fi(X

*2» • • ч

К1>

Д2’ ‘

• ■» ат)

^ 9 )

и имеет решения

 

 

 

 

 

Xi0 =

<ZV a 2> • '

•> am)-

 

(!• ^ 0)

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

величины

xi по

отношению

к параметру ak называется соотношение dxjdak= u ih. Производные dxi0/dа. называются коэффициентами влияния параметров. Они

могут быть получены одновременно с решением основной системы (1. 29). Действительно, вычисляя д/дак (dxjdt) и учитывая (1.30), получаем

 

dU

I dU „

I

I dfi

(1.31)

dt

дхх

* дх% 2к "т”

* дхп пкш

 

Решение этих уравнений чувствительности определяет коэффи­ циенты влияния параметров. Из вида уравнений (1. 31) может быть построена решающая схема, обеспечивающая получение коэффи­ циентов влияния.

Способ использования двух моделей состоит в построении моделей, имеющих параметры (аД и ( а{-\-Ва.). Частные произ­ водные по параметрам в этом случае аппроксимируются разностью выходов обеих моделей.

3 Л. И. Филиппов

33

Способ использования Модели 6 изменяющимися параметрами следует из того, что

d E

д Е da , d E

 

, л OON

d t ~

да d t " С d t

>

l 1, 6 *>

где а — один из параметров. Если изменять параметр по извест­ ному закону и измерять результат этого изменения dE/dt, то можно получить необходимую информацию о влиянии а на Е, т. е. дЕ/da. Величина dEldt при этом действует как неизбежная помеха, добавляющаяся к помехе от шумов.

Полученные одним из указанных способов частные производ­ ные дЕ/да используются для минимизации ошибки в пространстве ее параметров. Метод сравнения с моделью дает большое разно­ образие вариантов.

1.5. Корреляционно-фильтровые методы

Эти методы состоят в изучении импульсного отклика h (t, v) или коэффициента передачи К (ш, г) с помощью взаимокорреляционной обработки, которая состоит в умножении принятого колеба­ ния s2 (t) на образец излученного, но сдвинутого по частоте и времени, s (t— у Q), с последующим интегрированием для единич­ ных сигналов или низкочастотной фильтрацией для периодиче­ ских сигналов.

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

h (*) = $ й (У) (* — У) dy = y Re J h (у) (t — у) dy. (1. 33)

Взаимокорреляционная обработка узкополосных в радио­ техническом смысле колебаний (Аш ш0) дает взаимокорреляционную функцию, представимую в виде

г (т:,

2) =

J s2 (t) Sj (t — т,

2) dt = y Re J s2 (t) s* (t — x) e-^ -^ d t =

 

 

 

= 4

Re J А(у) Ф (x — y, Q)dy,

(1.34)

где

Ф (x,

£2) — комплексная

обобщенная функция автокорреля­

ции

сигнала sx (t) [19].

 

 

 

Целесообразно ввести также две модификации обобщенной взаимокорреляционной обработки, отличающиеся распределением

запаздывания х и смещения по частоте Q между сигналами:

 

v (х, 2) = ^ s2 (t) Sj (t — x, 2) dt =

 

= — ^ Im ^й(у)Ф (х — у, 2 )dy = —f(t, 2);

(1.35)

34

U(t, 2) == J s2 (t, — 2) Sj_(t z)dt = r (x, 2) cos 2x —

' — y(x, 2) sin 2x = r (x, 2) cos 2x -f-f (x, 2) sin 2x.

(1.36)

Из (1.36) видно, что r(x, 2) равно x/2 действительной части преобразования Фурье для произведения i (t) = tf2 {t ■=)ejSx- Рассмотрим также преобразование Фурье для произведения дей­ ствительных сигналов:

Gz( 2 , х ) = J s 2 ( f ) Sj (t — x) e -^ d t =

 

= j h (y) { J Sj (t — y) s± (t — x) e~J'stdt^ dy.

(1. 37)

Его можно выразить через обобщенную функцию автокорреля­ ции Ф (х, 2):

Сг(2. =

\ h {y )4 (? - y ,

2 )dy +

 

 

+

- 2 ) < Ц .

(1.38)

Вычислим также колебание гф(t, х, Q) на выходе узкополос­ ного фильтра (резонансного контура) при подаче на его вход колебания z (t). Полагая приближенно, что для такого фильтра

7? (2) = 2u 1 [8 (2 — 2 Х) - f 8 (2 - f 2J],

(1. 39)

получим методом преобразований Фурье

 

 

 

Х) * ( В ) eJSidQ =

 

 

=

1 { J h (у) Re 1е-Ув(^)ф (х -

у, 2)] dy +

 

+

J Н(у) Re [е-№ > Ф * (х - у,

- 2 ) ]

(1. 40)

Теперь, пользуясь соотношениями (1. 34)—(1. 40) и задаваясь различными видами функции Ф (х, Q), можно получить ряд ва­

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

 

 

1. Пусть

 

2) =

2it8 (2) е7“»т

(1.41)

ф(х,

(гармоническое колебание

с частотой

ш0). Подставляя

(1. 41)

в (1. 34) и вычисляя результат для х=0

и х = —п/2ш0, получаем

г (0,

2) =

(2) Re Й. (и)0);

(1.42)

г ( ~ 2 ^

S) =

(S) Im Я К )-

(Ф 43)

Этот результат поддается простой интерпретации: из (1. 34) видно, что г (х, й), как функция £2, есть действительная часть спектра произведения сигналов. Из (1. 42) видно, что в рассматри-

3* 35

I

Рис. 1.2. Блок-схема метода параллельного корреляционного анализа

ваемом случае спектр содержит одну составляющую на нулевой частоте. Соотношения (1. 42) позволяют определить действитель­ ную и мнимую части коэффициента передачи канала.

2.Пусть теперь Ф (х, Q) имеет S-образный всплеск в окрест

ности т= т,. Q =Q j (фазоманипулированный сигнал с большим основанием и малой длительностью элемента или сигнал с линей­ ной частотной модуляцией). При этом в (1. 34) и (1. 35) функцию h (у) можно заменить ее значением в точке (хх> Qx) и вынести за знак интеграла. Обозначая

В = \Ф(х, a,) dr,

(1.44)

получим

 

г (г, Q1) = -l-h('c — x ^ R ei?— -|-Я(х— x^Im-B;

(1.45)

v(r, 2 t) = — -i£ (x — x J R e # — \h(x — x^Im R .

(1.46)

Отсюда получаем блок-схему метода отождествления, приведен­ ную на рис. 1.2. Этот метод можно назвать методом параллельного корреляционного анализа. Схема разделения решает систему двух

уравнений (1. 45), (1. 46) с неизвестными h ( x j и %( x j . В частности, если сигналы таковы, что Qx=0, то 1ш 5=0 и схема разделения не требуется.

Рассмотрим далее обработку по (1. 36). С учетом формул (1.45)

и (*, S t) = h (х — Xj) -i- j В |cos (2xx — arg В) -j-

 

+ й(х — i1)^-\B\ sin (QjX — argB).

(1.47)

Операцию умножения на sx (£— x) и интегрирования по £ в (1. 36) можно осуществить с помощью согласованного с sx (£) фильтра (для всех х). Тогда на выходе получим и (х, й х).

Если, в частности, Qx=0, то схема упрощается до одного согласованного фильтра.

36

Рис. 1.3. Блок-схема спектрально-корреляционного метода отождествления

3. Пусть, наконец,

 

W (х, 9 ) = ехр |/о)0т -f- jQ - 1 8Х_3 (ят — Й)}

^

(сигнал с линейной частотной модуляцией в окрестности

й = 0

при большой длительности сигнала). В дополнение к двум преды­

дущим вариантам

здесь

возможен

еще

один. Пусть в (1. 40)

h (у) = 0 при

тх ;> у > т2

( т2 > Tj),

а в образец сигнала введена

задержка х <

тх. Тогда

 

 

 

 

 

2Ф(<> х,

Q) =

^ h {x — -5 -)cos[fi {t — *) +

(% —

(1.49)

Домножив 2ф на cos (Й (t—т ) + [ cu0—(й/2)]

(

й /я )}

и профильтро­

вав результат, получим h [т —(й /я )]. Этот метод можно назвать

спектрально-корреляционным (рис. 1.3).

z

Если же в (1. 49) произвести интегрирование по й, то получаем

(t) (tx, й). Если подобрать начало отсчета так, чтобы

т= т1= 0 , то

 

^2

2(*) = т \ h (у) cos [уК+ “ о) — х]dy-

(1.50)

о

 

 

Если подобрать параметры так,

что яг|/2 2к, то

 

=

о).

(1.51)

 

Это — широко распространенный принцип измерителя частот­ ных характеристик. Он дает модуль коэффициента передачи

вопределенном диапазоне частот.

4.До сих пор предполагалось, что канал не изменяется со вре­ менем. Для канала с переменными параметрами в качестве зон­

дирующих

необходимо использовать периодическое

колебание

с периодом

Та <С 1/Дfa= x a, но Т0 >

т:кап. При этом спектр коле­

бания z (t)=s2 (t) 5Х [if— х, к (2-к/Т0) ]

будет состоять

из непере-

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

37

дельного корреляционного

анализа

останутся

в

силе, если

и 1 тл

В первом случае будем иметь на выходе^ Re it ( <о0, t)

(ш0, t),

а в схеме рис. 1.2 получим h (tk, t) и h (xft, t).

Схема на согласованных фильтрах сохраняется при условии

Тв <^ 1/Д /а. Схема рис.

1.3

останется в силе,

если

полоса пропу­

скания

узкополосного

фильтра

будет

удовлетворять •условию

Ч <

Д/Ф<

(2/Гс) - Д / а.

При

этом

на

выходе

получим

h [ч-к(2к/Т0), *].

1.6. Примеры некоторых простых каналов

Для исследования каналов в ограниченной полосе частот наи­ более удобным является моделирование их функциями h (у, t) [или h (v, £)] и 0J (у, Q) [или б2 (у, Q)]. Целесообразно рассмотреть аналитическую запись этих моделирующих функций для трех, типов простейших каналов, которые могут встретиться на прак­ тике:

1. Канал, вносящий только запаздывание на величину т3. В соответствии с соотношениями (1. 9), (1. 12) и (1. 15) и определе­ ниями получаем:

h(v, t) = 8 (у ч:8),

h(v, i) = 8(y — та) + - Д ;Ч _ _ = 281(у — %),

0, (у, 2) = Ft [h (y, *)] = 2n8 (у - T.) 8 (2),

(1' 52)'

б2(у, 2) = F,[/i(y, £)] = 2u8i (y — xa) 8(2).

2. Канал, вносящий только сдвиг по частоте на величину 2 0:

h (у,

t) =

8 (у) cos Q0t — -i- ^ sin 2 0£,

 

 

h(v,

t) =

28j (y) e-W,

 

 

(у ,

2 ) =

2n [8 ( 2 — 2 0) 8t (y) —f—8 (Q j—2 0) 8X(y)J,

^

>

62(y,

2) =

2и8г (y) 8 (2 — 2 0).

 

 

3, Канал, вносящий запаздывание и сдвиг по частоте:

 

h (у ,

t) =

8 (у — т з) cos 2 0 (t — т а) — 1

sin 2 0 (f xt),

Л (у, t) = 28х (у — хв)

 

б, (у,

Q) =

( 2

- 2 e)81( y _ x , ) e - ^ . +

(1.54)

 

 

 

+ 8 ( 2 + 2 0)[81( у - г в) е- ^ Г ,

б2 (у,

2) =

28х (у — ха) 8(2 — 2„) е -^ з .

 

Из рассмотрения этих соотношений видно, что наиболее удоб­ ной с точки зрения последующих вычислений является аналити­ ческая форма представления функций из-за наличия 8-функций, упрощающих интегрирование иа основании соотношения (В. 17).

38

1.7. Автокорреляционные модели каналов

■Для ряда каналов отклики h (у, t) и преобразования Фурье

£ ( со, t) и 0 (Q, у) являются реализациями случайных процессов. Для их полного описания нужны многомерные законы распреде­ ления. Однако большинство каналов близки к гауссовым. В этом случае знание автокорреляционной функции эквивалентно зна­ нию многомерного распределения [18].

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

автокорреляционные

функции для h (у, t),

К ( ш, t)

и 6а (у, &):

7Д Д Р

^2’

^1»

Д) ---^ (Ур Д) Л- (У2* A)J

 

Kjf(u)j,

COg,

Д,

Д)

К (^1> Д)

(^2*

Д)>

(1*

^ о Д у у2, 2 Х,

2 2) =

М У1> Ql) Оо(г^2,

2 2)-

 

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

Теперь мы можем определить четыре основных типа каналов со случайными свойствами:

1)

стационарный канал коррелированных отражателей,

2)

нестационарный канал некоррелированных отражателей,

3)стационарный канал некоррелированных отражателей,

4)нестационарный канал коррелированных отражателей. Наиболее сложный случай из них четвертый.

Стационарным называется канал, "для которого процессы

являются стационарными во времени t, так что

 

7Д (у-p у21 Д» *

(Дл

Д

ДЬ

 

 

Д)—

(1.56)

Rk{<0v Ш2,

Д, !,s)==Rk(uV у2> Д

Д)-

 

Можно найти., что

 

00

 

 

 

■Rg К , ^2* ^1» Q2) = 2-nrS ( 2

j — S2)

 

 

(1.57)

\ Rh(vv

 

x)e~^xdx,

— СО

где через х обозначена разность Д—Д.

Из (1. 57) видно, что стационарный канал является не корре­ лированным по частотному рассеянию Q (дельта-коррелирован­ ным) .

Каналом некоррелированных отражателей называется такой,

для которого реализации h (у) некоррелированные, т. е.

 

^hiPv Ду

Д) — ДЛДр Д1 Д) ^ Д] ДД

(1.58)

При этом можно найти,

что

 

39

Rk К ,

ш2, tv t,) =

jj ^ Rh (vv y2, tL, U) e

x

p

=

 

 

“—CO

 

 

 

 

0 0

 

 

 

 

=

5 Rh(v!, tv

ts) e x p ^ .K - ^ )]^ = RK(coj — <o2,

tv *2).

(1.59)

 

—CO

 

 

 

 

Отсюда видно, что канал некоррелированных отражателей яв­

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

по ш (но не t)

(Rk зависит

от разности

«1 — ш2).

 

 

 

 

Стационарным каналом некоррелированных отражателей яв­

ляется канал, для которого

 

 

 

R/i (У1> У2>

^г)= Rh ivv

^2) ^ (vi

^2)-

•^0)

При этом, как можно убедиться,

 

 

 

R]c (ш1> ш2’

*£> ^2) == R]c (Ш1

Ш2’ ^1

^)-

(1.61)

Таким образом, этот канал стационарен по обеим переменным.

Для этого канала автокорреляционная функция процесса $2 (у, Q) имеет вид

Rt (vV v2, Q2) = 2nC(vv Ql)b(v1 — v2)b(Q1 — Q2),

где

0 0

 

С (vv 2 Х) = ^ Rh (иг, х) e~->Qxdx

(1.62)

является преобразованием Фурье от В.Лпо переменной x = tt12 и

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

Для эргодических процессов h (у, t), заменяя усреднение по ансамблю усреднением по времени, можно получить

Т/2

RH{v,

х) = -гг lim

J

h (у,

t) h* (у,

t х) dt,

 

 

 

д / 24-со —272

 

 

 

 

где А/ — полоса, в которой изучается процесс h (у, t).

 

Тогда функция рассеяния

 

 

 

 

 

 

С (уг,

2j) = И т 4 - 1022' К ,

2 Х) \\

(1. 63)

 

 

 

Т/2

со х

 

 

 

 

 

 

 

 

 

 

где

^2t (vv

2 г) =

f

h (yL,

t)e~iQlidt.

(1.64)

—272

Тз (1. 63) и (1. 64) видно, что функция рассеяния есть предел при h -> 00 квадрата модуля спектральной функции реализации (у, t). Следовательно, С ( у , Й) является спектральной плотностью средней мощности реализации (как функции текущего времени).

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