Файл: Кондратьев, С. Л. Применение метода функционального моделирования для оценки помехоустойчивости систем связи.pdf

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

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

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

Добавлен: 21.10.2024

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

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

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

при каждом данном значении отношения энергии сигнала к спектральной плотности шума

^ = Л8

(2.2.5)

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

Для решения этой задачи необходимо обработать напряже­ ние на входе u(t):

и (0 = с (<, 7) + п (*),

(2.2.6)

где т — вектор неизвестных параметров, т. е. амплитуды, фазы, частоты.

Упростим задачу, полагая, что частота и фаза нам известны или безразличны, т. е. будем интересоваться только амплитудой. Очевидно, что можно предложить несколько методов gi решения задачи, но нам требуется такое gu чтобы оценка амплитуды а была оптимальной, т. е. обеспечивалось р0 ш=Ропт-

Пример 2. В соответствии с теорией Котельникова В. А. потенциальная помехоустойчивость при гауссовом шуме дости­ гается использованием фазовой модуляции ФМ, если в качестве схемы обработки сигнала с шумом применяются согласованные фильтры или подобные им схемы корреляционной обработки.

Полученное выражение

/>ош= 1 - / 4 1 ^ ) ,

«о„т = 2,

(2.2.7)

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

t0=iH+kT.

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

50


пользуя, например, минимальные свойства совокупности о т к . ^ нений Дг при выбранном алгоритме обработки сигнала gi.

Обратим внимание на то, что если выбранный нами критерий качества оценки не учитывает истинного характера процесса, то полученная «оптимальная» процедура оценки gi на самом деле может оказаться хуже «.неоптимальной».

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

Итак,

какую же функцию G отклонений

Д,- следует при­

знать наилучшей?

 

 

 

 

 

 

 

 

Поскольку величина Дг зависит

от

самих

параметров

T1.T2.---.Tfn и результатов

измерений

хх, x2,...,xN,

gj

то

полу­

ченные с помощью алгоритма обработки сигналов

величины

Yi*. Y2*, - • •. Уп*, обеспечивающие

минимизацию

G (Д|, Дэ,... ,ДП ),

будут зависеть как от вида функции

G, так и от

результатов

измерений, т. е.

 

 

 

 

 

 

 

 

 

Ъ* = Т * (О,

. . . , * * ) .

 

 

 

(2.2.8)

Поэтому

алгоритм gj получения

оценки

тг* (алгоритм

оцени­

вания) в общем случае сам является функцией от G, т. е.

 

gj—fiG, xt).

 

 

 

 

(2.2.9)

Поясним это обстоятельство.

Произведем

прямое

измере­

ние и выберем оценку в виде (*):

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

T i * = ^ 2 ^ .

 

 

 

 

(2-2.10)

Тогда, как будет показано

ниже, в минимум обращается

сумма

квадратов

ошибок:

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G, (Al f Д „ . . . ; Д„) = 2 V

= min.

 

 

(2.2.11)

 

 

 

j=i

 

 

 

 

 

 

Если в качестве оценки взять медианное значение (***):

T a * = { * « - i < * i < - * ' » + i } .

(2.2.12)

то в минимум обратится сумма абсолютных

значений оши­

бок:

 

ЛГ

 

О з ^ , Д 2 ) . . . , Дл,) = 2 | А ' И т 1 п -

(2 - 2 - 1 3 )

i-i

 


Полагая, что можно воспользоваться оценкой (**):

 

 

 

 

1

 

 

 

 

 

 

 

 

(2.2.14)

 

 

 

 

 

 

 

 

 

 

 

 

 

убеждаемся в том, что в минимум

обращается

максимальное

из всех значений

Д ь т. е.

 

 

 

 

 

 

 

 

 

 

 

 

G2 (Д„ Д 2 , . . . , Д Л ) = т а х (Д„ Д 2 , . . . , Д„)=т1п .

 

(2.2.15)

Таким образом, рассуждая в обратном порядке, можем

поставить

вопрос,

какую же

процедуру

gt. оценивания следует

 

 

 

 

 

 

использовать,

если

руко

а)

 

 

 

 

 

водствоваться

избранным

 

 

 

 

 

критерием

качества

оцен­

 

j дцелки\

 

 

 

ки

G.

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

 

 

тре­

 

 

 

 

 

 

 

 

 

\Слещенныё~\

\Неслещенные\

 

буется

определить,

какие

 

же

оценки

являются

«хо­

 

 

 

 

 

 

рошими». Теория

оценок

 

 

 

 

 

 

вводит

для

этого

опреде­

 

 

 

 

 

 

ления,

характеризующие

 

 

 

 

 

 

оценку

с

разных

сторон.

 

 

 

 

 

 

Здесь ограничимся только

 

 

 

 

 

 

основными,

замечая,

что

 

 

 

 

 

 

в

основе

 

разбиения

на

 

 

 

 

 

 

классы

опять

положена

 

 

 

 

 

 

общность

 

 

(отношение)

 

 

 

 

 

 

между

оценками

 

(рис.

 

\1остаятельныё\

 

2.1а, б).

Таким

образом,

\асанптотичеекои.

 

 

разбиваем

все оценки

на

\Acunn7oTayecxv-

два

класса: смещенные и

 

оцет

 

 

 

 

 

 

 

 

 

 

 

 

УаффективиостьнЛ е* limen - /

несмещенные.

 

 

 

 

 

 

 

 

 

 

 

 

1.

Оценка

т*

 

назы­

 

 

 

 

 

 

вается

несмещенной, если

 

РИС. 2.1.

 

 

ее

математическое

ожи­

 

 

 

дание

равно

самой

вели­

 

 

 

 

 

 

 

 

 

 

 

 

чине if:

 

 

 

 

 

 

 

 

 

 

 

 

 

Т + 6 = Т .

 

 

 

 

(2.2.16)

т. е. если

среднее

значение

смещения

b

равно

нулю.

При

Ь Ф 0 оценка будет

смещенной.

 

 

 

 

 

 

 

 

 

 

2. Несмещенная оценка называется эффективной, если она обладает минимальной дисперсией, т; е. если:

1) М (-г*)=7*=т;

2) о?,[ т * - 7 * ] 2 = т * 2 - 2 т * Г + Т 3 = т 1 п .

(2.2.17)

62


Заметим, что если определена условная плотность распре­ деления W (.и/ч), то в общем случае эффективность оценки опре­ деляется неравенством Крамера—Рао, т. е. оценка эффективна в общем смысле, если неравенство Крамера—Рао при выбран­ ной оценке т* обращается в равенство:

 

 

(

М |(Г*-Т)2 1 >

6=т—t*.

(2.2,18)

 

м

 

Если смещение Ь не зависит от у или равно нулю, то полу­ чаем

1

(2.2.18а)

 

Этим устанавливается минимально достижимое значение диспер­ сии оценки (т. е. нижняя граница).

Можно показать, что дисперсия минимальна и равна правой части неравенства Крамера — Рао тогда и только тогда, когда функция правдоподобия параметров удовлетворяет следующему условию:

- ^ - Ш

^ ( и / т ) + ( т * - Т ) А ( т ) ,' .

(2.2.19)

где k(y)—произвольная

функция параметра

у, в частности,

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

Если оценка не обладает свойством эффективности, то она называется неэффективной и для характеристики ее качества (близости) .вводится мера эффективности:

Л2

 

'эф < 1 ,

(2;2.20)

которая, равна единице только для эффективной

оценки.

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

ляется как

предел

частости n/N (несмещенной и эффективной

оценки):

"

 

 

Л=Нт -77,

53