Файл: Силяков В.А., Красюк В.Н. - Системы авиационной радиосвязи - 2004.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 20.10.2020

Просмотров: 5041

Скачиваний: 108

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
background image

61

фильтрации, так как 

ˆ

λ ( )

t

и 

  R

λ

(

t

) исчерпывающе задают одномерную

нормальную АПВ.

Асимптотические свойства рассмотренного алгоритма фильтрации

отмечены в Прил. 2.

Упрощенные алгоритмы фильтрации

. Алгоритмы (3.52) и (3.53)

можно упростить, если учесть, что гауссовская аппроксимация АПВ
справедлива при большом отношении сигнал/шум в канале [14], кото-

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

 R

λ

(

t

) оценки 

ˆ

λ ( )

t

. Разло-

жим функцию 

F

(

t

λ

) в малой окрестности точки 

ˆ

λ ( )

t

в ряд Тейлора и

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

2

2

2

ˆ

ˆ

( ,

λ)

( , λ)

1

ˆ

ˆ

ˆ

( ,

λ)

( , λ)

λ ( ) λ( )

λ( ) λ( ) .

ˆ

ˆ

2

λ

λ

F t

F t

F t

F t

t

t

t

t

=

+

+

(3.54)

Подставим данное разложение  в (3.52) и (3.53). Тогда

λ

ˆ

ˆ

λ

( , λ)

( ,

λ) ( , λ) λ

ˆ

λ

d

F t

g t

p t

d

R

dt

Λ

=

+

;

(3.55)

2

2

λ

λ

λ

λ

2

ˆ

( ,

λ)

( , λ)

2

( ,

λ) λ

.

ˆ

2

λ

λ

d R

N

g t

F

t

R

p t

d

R

dt

Λ

=

+

+

(3.56)

При рассмотрении стационарного режима работы приемного устрой-

ства в дифференциальных уравнениях (3.55) и (3.56) следует положить

значение 

λ

0

d R

d t

=

. Тогда выражение для дисперсии фильтрации оп-

ределяется алгебраическим уравнением. Если уравнение сообщения яв-
ляется линейным

λ

λ

α λ

( ),

d

n t

d t

= −

+

(3.57)

то выражения (3.55) и (3.56) приводятся к виду

λ

ˆ

ˆ

λ

( , λ)

ˆ

α λ

;

ˆ

λ

d

F t

R

dt

= −

+

(3.58)


background image

62

2

2

λ

λ

λ

2

ˆ

( ,

λ)

2

λ

0.

ˆ

2

λ

N

F t

R

R

+

=

(3.59)

Для отыскания алгоритмов фильтрации при некоторых законах мо-

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

3.5.2. Многомерные уравнения наблюдения и сообщения

Обобщение уравнений (3.52), (3.53) и (3.55), (3.56) на случай много-

мерных наблюдений и сообщений, описываемых соотношениями (3.17)
и (3.18), можно выполнить по методике,  аналогичной использованной
в Прил. 2.  Применение такого подхода при аппроксимации АПВ выра-
жением (3.50) приводит к алгоритму фильтрации, определяемому сис-
темой уравнений  [14]

  

( )

{

}

( )

1

p

1

1

λ

,

,

λ

λ

P

p

p q

q

q

d

F t

M

a

t

R

M

dt

=

=

+

λ

,

( )

( )

( )

( )

p

1

1

µ

1

µ

µ

µ

µ 1

2

µ

ν

1

µ

ν

µ 1

1

,

,

,

2

λ

λ

,

;

,

1... ,

λ

λ

P

p q

pq

q

p

q

P

P

p

q

dR

N

t

a t

a

t

M

M

R

M

R

dt

F t

R R M

p q

P

=

=

=

=

+

+

+

+

=

∑ ∑

λ

λ

λ

λ

ν

где 

{ }

1

M

 – оператор осреднения с апостериорной плотностью вероят-

ностей (3.50).

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

денных соотношений частный случай, а именно, когда принимаемая
смесь 

ξ

(

t

) описывается скалярным уравнением, а сообщение 

λλλλλ

(

t

) явля-

ется векторным многокомпонентным 

[

]

Т

1

2

( )

( ),    

( ),

,

( )

P

t

t

t

t

=

λ

λ

λ

λ

и

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


background image

63

λ

( )

( ) ( )

( ),

d

t

t

t

t

dt

=

+

λ

A

λ

n

где 

1

T

λ

λ

λ

( )

( ),  ...,   

( )

P

t

n

t

n

t

= 

n

– вектор-столбец независимых БГШ, ха-

рактеризуемых  матрицей односторонних спектральных плотностей

{

}

(

)

T

λ

1

1

λ

2

1

2

λ

( )  ( )

δ

2

N

M

t

t

t

t

=

n

n

 A

(

t

) – матрица размерности (

р

×

р

),

имеющая вид

11

1

1

α    α

( )

.

α    α

P

P

PP

t

= 

A

В рассматриваемом случае квазиоптимальный алгоритм фильтрации

определяется системой уравнений (частным случаем которой  является
уравнение (3.55) для оценок компонент (

p

 = 1, ..., 

P

)

1

1

ˆ

ˆ

λ

( ,

)

ˆ

( )

λ

( )

,

ˆ

λ

i

P

P

p

i

pq

q

pq

q

q

q

d

F t

t

R

t

dt

=

=

=

+

λ

α

(3.60)

и системой уравнений (частным случаем которой является уравнение
(3.56)) для корреляций ошибок оценок

µ

µ

µ 1

1

( )

δ

α ( )

( )

α ( )

( )

2

P

P

pq

p

pq

p

q

q

p

R

t

N

t R

t

t R

t

t

=

=

=

+

+

+

ν

ν

ν

2

µ

ν

µ

ν

µ 1

1

ˆ

( , )

( )

( )

.

ˆ

ˆ

λ

λ

P

P

p

q

F t

R

t R

t

= =

+

∂ ∂

∑∑

λ

ν

(3.61)

Здесь 

{ }

λ λ

( )

,   ( ,   

1, ..., )

pq

R

t

R

p q

P

=

=

– матрица взаимных корреля-

ций ошибок оценок компонент 

ˆ

ˆ

и

λ ; δ

1, δ

0,

p

q

pp

pq

λ

=

=

при 

p

q

.

Упрощенные алгоритмы фильтрации

. Моделирование уравнения

(3.61) для большинства используемых в системах связи типов радио-
сигналов оказывается затруднительным и схема устройства, реализую-
щего уравнения (3.60) и (3.61), оказывается весьма сложной. Поэтому


background image

64

для построения подоптимальных алгоритмов используют ряд упроще-
ний [14]. Рассматривают стационарный режим работы, полагая

0

pq

R

t

∂ =

. Вместо вторых производных 

2

ˆ

ˆ

ˆ

( ,

λ) λ

λ

i

j

F t

, входящих

в (3.61), используют их значения, осредненные по времени. Отбрасыва-
ют члены с удвоенными частотами, так как последние сильно ослабля-
ются сглаживающими цепями. При указанных упрощениях получим
систему уравнений фильтрации [6,14]

      

p

1

1

ˆ

ˆ

λ

( ,

)

ˆ

α ( )λ

( )

,

ˆ

λ

P

P

i

pq

q

pq

q

q

q

d

F t

t

R

t

dt

=

=

=

+

λ

(3.62)

( )

µ

µ

µ

ν

ν

1

ν 1

2

µ

ν

µ

ν

µ 1 ν 1

δ

( )

( )

α ( )

( )

2

ˆ

,

( )

( )

0,

ˆ

ˆ

p

P

P

pq

p

q

p

p

i

P

P

p

q

N

t R

t

t R

t

F t

R

t R

t

=

=

=

=

+

+

+

+

=

∂ ∂

∑ ∑

λ

λ

λ

α

(3.63)

где  

( )

2

µ

ν

ˆ

ˆ

ˆ

,

F t

λ

λ

λ

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

Приведенные соотношения (3.62) и (3.63)  являются исходными для

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

3.6.  Фазовая модуляция

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

уравнения наблюдения и сообщения можно записать в следующем виде:

( )

(

)

( )

0

0

ф

ξ

cos ω

λ

,

t

A

t

M

n t

=

+

+

(3.64)

( )

( )

( )

λ

λ

αλ

d

t

t

n

t

dt

= −

+

,

(3.65)


background image

65

где 

А

0

ω

0

 – значения амплитуды и частоты радиосигнала, которые счи-

таются известными; 

λ

(

t

) – передаваемое аналоговое сообщение; 

M

ф

  =

 = 

σ

ф

/

σ

λ

 – известная крутизна характеристики фазового модулятора;

σ

ф

 = (

D

ф

)

0,5

 – среднее квадратическое отклонение фазы, обусловленное

модуляцией; 

σ

λ

 = (

D

λ

)0,5 – среднее квадратическое значение сообще-

ния 

λ

(t).

Уравнения нелинейной фильтрации (3.62) и (3.63) с учетом приня-

тых в них обозначений для модели полезного радиосигнала, входящего
в (3.64), преобразуются к виду [6]

λ λ λ

ˆ

λ( )

ˆ

α λ

,

d

t

R F

dt

= −

+

(3.66)

2

λ

λ λ

λ λ λ λ

4

α

2

0,

N

R

R

F

+

=

(3.67)

где для упрощения записи обозначено: 

( ) ( )

11

11

ˆ

ˆ

ˆ

α

–α, λ

λ

λ.

t

t

=

=

=

Учитывая, что информационный параметр неэнергетический, имеем

(

)

1

λ

0

ф

0

ф

ˆ

( ,

λ)

ˆ

2

ξ( )sin ω

λ

ˆ

λ

F t

F

N

A M

t

t

M

=

= −

+

(3.68)

2

1

2

2

λ λ

0

ф

2

ˆ

( ,

λ)

.

ˆ

λ

F t

F

N

A M

=

=

(3.69)

Подставив выражение для 

λ

F

 из (3.68)  в  (3.66), получим

1

λ λ

ф

0

0

ф

ˆ

λ

ˆ

ˆ

α λ 2

ξ( ) sin(ω

λ),

d

R M N

t A

t

M

dt

= −

+

(3.70)

или в операторной форме

0

0

ф

λ

ˆ

ˆ

λ

ξ( ) sin(ω

λ),

1

K

t A

t

M

T D

= −

+

+

(3.71)

где 

/

D

t

= ∂ ∂

 – оператор дифференцирования; 

λ

1

α

T

=

; произведение

α

N

 определяет мощность аддитивного белого шума 

n

(

t

) в полосе 

α

f

∆ =

частот сообщения

 λ

(

t

) [6,14]. Величину