Pz_kr_Otsenka_Parametrov_Signala
.pdf
Если дисперсия D Т стремится к нулю с увеличением времени наблюдения Т,
то условие сходимости по вероятности будет выполняться, т.е., если выполняется сходимость в среднеквадратическом, то выполнятся сходимость по вероятности,
и оценка будет состоятельной.
3. Эффективность оценки. Различные функции f y t дают оценки с различными дисперсиями D и смещениями ∆, поэтому желательно выбирать
функцию f y t таким образом, чтобы погрешность ε2 была бы наименьшей.
Шведский математик, актуарий и статистик Х. Крамер и индийский, американский математик и статистик Рао доказали, что дисперсия оценки не может быть меньше некоторой величины (неравенство Рао-Крамера)
|
|
|
|
|
|
|
|
d 2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
D эфф. , |
|
(25) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
||||||
|
|
|
|
|
M |
|
|
|
ln f y |
|
|
y t |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где M - математическое ожидание оценки. Если оценка несмещенная, |
|||||||||||||||||||||
то |
d |
1 и дисперсия оценки примет вид |
|
|
|
|
||||||||||||||||
d |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
1 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y t |
|
|
||||
|
|
|
|
|
|
D M |
|
|
ln f y |
|
|
|
|
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Выражение
J M ln f y
|
2 |
|
2 ln f |
|
|
y t |
|
|||
|
|
|
|
y |
|
|
|
|
|
|
y t |
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
2 |
|||||||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
1 |
(26) |
|
|
||
D эфф. |
||
|
получило название информации по Фишеру. Таким образом, дисперсия любой
несмещенной оценки не меньше величины, обратной информации по Фишеру. |
|||||||
|
|
|
|
|
|
|
|
4. Достаточность оценки. Оценка называется достаточной (достаточная |
|||||||
статистика), если |
совместную |
|
|
можно |
|||
плотность распределения f y t , |
|||||||
представить в виде |
|
|
|
|
|
||
f y t , |
|
f , h y t , |
|
|
|
|
|
|
|
|
|
|
|||
где функция h y t не зависит от , а вся |
|
|
|||||
информация о параметре |
|||||||
содержится в cтатистике. |
|
|
|
|
|||
Необходимым |
условием |
существования |
достаточной статистики |
||||
(достаточной оценки) является возможность представления совместного распределения в виде:
f y t , exp{A B( y(t)) C( y(t)) D( )}.
11
Из сказанного следует, что если оценка эффективная, то она будет достаточной. Обратное не всегда верно. Иногда трудно бывает найти эффективную оценку, хотя она является достаточной.
Среди всех методов оценки параметров метод максимального правдоподобия обладает рядом достоинств:
1.В случае оценки одного параметра оценка наибольшего правдоподобия оказывается всегда состоятельной.
2.При большом времени наблюдения распределение оценки является
приближено нормальным, с центром в точке λ и дисперсией D , определяемой
формулой (25) , т.е. оценка является асимптотически эффективной.
3. Если оценка является эффективной, то решение уравнения
правдоподобия (12) является единственным.
4. Оценка является достаточной статистикой, т.е. метод максимального правдоподобия использует всю информацию о параметре, содержащейся в реализации y(t).
6. Оценка параметров сигнала по методу максимума функционала правдоподобия
6.1. Общая теория метода
Предположим, y(t) – наблюдаемая на интервале 0,T реализация процесса,
представляющая сумму нормального случайного процесса n(t) c нулевым средним и известной корреляционной функцией Rn(t,u) и детерминированного сигнала s(t;λ1,…λm), зависящего от неизвестного векторного параметра λ=(λ1,…λm). Нужно оценить неизвестные параметры детерминированного слагаемого по критерию максимального правдоподобия.
Запишем логарифм функционала отношения правдоподобия:
|
T |
|
1 |
|
|
|
|
|
|||
|
|
|
|||
ln F[ y(t) |
] 2 V t; y t |
|
s t; dt , |
(27) |
|
|
|||||
|
|
|
2 |
|
|
0 |
|
||||
где V(t;λ) – решение неоднородного линейного интегрального уравнения
T |
1 |
|
|
|
|
Rn (t, u)V (u; )du |
s(t; ) , |
t T . |
(28) |
||
2 |
|||||
0 |
|
|
|
||
|
|
|
|
Частные производные по параметрам от логарифма функционала отношения правдоподобия равны
|
ln F y t |
|
|
2T |
|
V u; y t s t; dt , i=1,…,m. |
(29) |
||
|
|||||||||
|
|
||||||||
i |
|
|
|
|
i |
|
|
|
|
|
|
|
0 |
|
|
|
|||
|
|
|
|
|
|
|
|
12 |
|
Из этого непосредственно получаем систему уравнений максимального правдоподобия
Tn |
|
|
V u; y t s t; dt 0 , |
|
|
|
|
|
||||
|
|
|
|
|
i=1,…,m. |
(30) |
||||||
|
|
|
|
|
||||||||
0 |
i |
|
|
|
|
|
|
|||||
Решая систему уравнений относительно неизвестных λ1,…,λm, можно найти |
||||||||||||
|
|
|
|
|
|
|
|
|
|
g |
i y t . |
|
оценки максимального правдоподобия параметров |
i |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Если процесс n(t) представляет белый шум со спектральной плотностью |
||||||||||||
N0/2, то из (28) следует, что |
|
|
|
|
|
|
||||||
|
|
|
|
V t; |
2 |
|
s t; , |
|
|
|
|
(31) |
|
|
|
N 0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||
и система (30) значительно упрощается: |
|
|
|
|
|
|||||||
|
T |
|
s t; y t s t; |
dt 0 , |
i=1,…,m. |
|
(32) |
|||||
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|||
|
0 |
i |
|
|
|
|
|
|
||||
Предположим теперь, что функция s(t; λ1,…,λm) может быть представлена в |
||||||||||||
виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
s(t; 1,..., m ) j s j (t) , |
|
|
|
|
(33) |
|||
j 1
где sj(t), j=1,…,m,–известные функции.
Найдем совместные оценки максимального правдоподобия параметров λ1,…,λm. Подставляя (33) в правую часть (28) , заменим интегральное уравнение (28) системой уравнений
T |
1 |
|
|
|
|
Rn t,u Vi u du |
si t , |
t T , i=1,…,m |
(34) |
||
2 |
|||||
0 |
|
|
|
||
|
|
|
|
причем функция V(t;λ1,…,λm), от которой зависит функционал отношения правдоподобия, равна
|
|
|
m |
|
|
|
V t; 1,...,m jVj t . |
(35) |
|||||
|
|
|
j 1 |
|
|
|
Из (35) находим |
|
|
|
|
|
|
|
|
V t; ,..., |
|
|
V t |
(36) |
|
|
m |
||||
|
i |
1 |
|
i |
|
|
|
|
|
|
|
|
|
и, подставляя (36 2.1.10) в (30 2.1.4) получаем систему уравнений максимального правдоподобия
T |
|
m |
|
|
Vi t y t j s j t dt 0 , i=1,..,m, |
(37) |
|||
0 |
|
j 1 |
|
|
или
13
m |
T |
T |
|
j Vi t s j t dt |
Vi t y t dt , i= 1,…,m. |
(38) |
|
j 1 |
0 |
0 |
|
Введем обозначения |
|
|
|
|
T |
|
|
|
sTij Vi t s j t |
dt , |
(39) |
|
0 |
|
|
|
T |
|
|
|
yTi Vi t y t dt |
(40) |
|
|
0 |
|
|
и представим систему линейных уравнений (38) в виде |
|
||
|
m |
|
|
sTij j yTi , |
i=1,…,m |
(39) |
|
|
j 1 |
|
|
или в матричной форме |
|
|
|
|
STλ=YT, |
(40) |
|
где ST – матрица размером m m, элементы которой равны sTij , а λ и YT – векторыстолбцы, элементы которых λi и yTi соответственно.
T
Полагая, что для всех j s2j t dt и что Rn(t,u) – положительно определенная
0
функция, приходим к заключению, что существует матрица sT1обратная матрице sT. Тогда решение уравнения (40) приводит к следующим оценкам максимального
правдоподобия неизвестных параметров: |
|
|
||||||
|
|
|
|
|
|
|||
|
|
|
sT1 YT . |
|
(41) |
|||
Если случайная составляющая наблюдаемого процесса представляет белый |
||||||||
шум с интенсивностью N0 , то из (34) следует, что |
|
|||||||
N0Vi(t)=si(t), |
|
i=1,…,m, |
|
|
||||
и, следовательно, матрицы ST и XT преобразуются к виду |
|
|||||||
|
|
T |
t s j t dt |
|
|
|
|
|
|
|
|
|
|
|
|||
ST |
|
si |
, |
|
i,j=1,…,m |
(42) |
||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
YT |
|
si t y t dt |
, |
i=1,…,m. |
(43) |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
6.2.Совместные оценки амплитуды и фазы гармонического сигнала
Рассмотрим простой и наиболее распространенный пример сигнала s(t;λ) –
гармонический сигнал с известной частотой ω0 и неизвестными амплитудой A и фазой φ :
s(t; λ1, λ2)=Acos(ω0t – φ)=A cosφ cosω0t + A sinφ sinω0t |
(44) |
|
14 |
и, следовательно |
|
|
|
|
|
|
|
|
|
|
λ1=A cos φ, |
|
|
λ2=A sin φ, |
|
(45) |
|||||
s1(t)=cos ω0t, |
|
s2(t)=sin ω0t. |
|
(46) |
||||||
Выписываем элементы вектора XT и матрицы sT: |
|
|
||||||||
T |
|
|
|
|
|
|
T |
|
|
|
YT 1 2 V1 t y |
t dt , |
|
|
|
YT 2 |
2 V2 t y |
t dt , |
(47) |
||
0 |
|
|
|
|
|
|
0 |
|
|
|
T |
|
|
|
|
|
|
T |
|
|
|
sT 11 2 V1 t cos 0tdt , |
sT 12 |
2 V1 t |
sin 0tdt , |
(48) |
||||||
0 |
|
|
|
|
|
|
0 |
|
|
|
T |
|
|
|
|
|
|
T |
|
|
|
sT 21 2 V2 t cos 0tdt , |
|
sT 22 2 V2 t sin 0tdt , |
(49) |
|||||||
0 |
|
|
|
|
|
|
0 |
|
|
|
где V1(t) и V2(t) представляют решения интегральных уравнений |
|
|||||||||
T |
|
1 |
|
|
|
|
|
|
|
|
Bn t, u V1 |
u du |
cos 0t , |
|
t T, |
|
(50) |
||||
2 |
|
|
||||||||
0 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||
T |
|
|
|
1 |
|
|
|
|
|
|
Rn t, u |
V2 u du |
sin 0t , |
t T |
(51) |
||||||
|
||||||||||
0 |
|
2 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
и Rn(t,u) – корреляционная функция аддитивного, нормального случайного процесса.
Подставляя (47), (48), (49) в (50) и решая систему двух линейных относительно λ1 и λ2 уравнений, получаем оценки максимального правдоподобия этих параметров
TT
V u, v y u sin 0vdudv
1 |
T T |
|
, |
(52) |
||
|
|
|
0 |
0 |
|
|
|
|
V u, v cos 0u sin 0vdudv |
|
|
||
|
0 |
0 |
|
|
|
|
TT
V u, v y u cos 0vdudv
|
|
|
0 0 |
|
, |
(53) |
T T |
|
|||||
2 |
|
u, v cos 0u sin 0vdudv |
|
|||
|
|
V |
|
|||
|
|
0 |
0 |
|
|
|
где V (u,v) V1(u)V2 (v) V1(v)V2 (u) . |
(54) |
|||||
6.3. Оценка амплитуды и фазы гармонического сигнала, при целом отношении T/ T0
Если случайная составляющая наблюдаемого процесса представляет белый шум с интенсивностью N0, то из (50), (51) следует
V t N |
1 cos t, |
V t N 1 sin t |
(55) |
|||
1 |
0 |
0 |
2 |
0 |
0 |
|
V u,v N 1[cos usin |
v cos |
vsin u] |
(56) |
|||
0 |
|
0 |
0 |
0 |
0 |
|
|
|
|
|
|
|
15 |
|
Подставляя (56) в (52) и (53) и полагая ω0Т = |
2 |
T =2πk, где k = |
T |
– целое |
||||||||||||||||
|
T |
T |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
0 |
|
|
||
число, находим после простых вычислений |
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
2 |
T y(u) cos udu, |
|
2 |
|
T |
y(u) sin udu |
|
|
|
(57) |
|||||
|
|
|
|
|
T |
T |
|
|
|
|
|||||||||||
|
|
|
|
|
1 |
|
|
|
0 |
2 |
|
|
0 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
Оценки максимального правдоподобия (52, 53) параметров λ1 |
и λ2 |
могут быть |
||||||||||||||||||
использованы |
для |
получения |
оценок |
|
|
|
амплитуды |
|
и |
фазы: |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
2 |
T |
|
2 |
2 |
T |
|
|
|
2 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
12 22 |
|
y(u) cos 0udu |
|
|
y(u) sin 0udu |
|
|
|
|
|
|
|
(58) |
|||||||
|
|
|
|
|
|
|
|
||||||||||||||
|
|
T |
0 |
|
|
T |
0 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
y(u) sin 0udu |
|
|
|
|
|
|
|||||
|
|
|
arctg |
arctg |
0 |
|
|
|
|
|
|
|
|
|
(59) |
||||||
|
|
|
1 |
T |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
y(u) cos 0udu |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
tg |
2 |
, |
cos |
|
|
|
|
1 |
|
|
, |
sin |
|
|
|
2 |
|
. |
(60) |
||
|
|
|
|
|
|
|
|
||||||||||||||
|
1 |
|
|
|
|
2 |
2 |
|
|
|
2 |
2 |
|
||||||||
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
1 |
|
|
2 |
|
|
|
На рисунке показана структурная схема приемника для оценки амплитуды гармонического сигнала, построенная по формуле (58).
|
|
T |
|
|
|
х |
|
Кв |
1 |
|
|
0 |
|
|
|
Г1 |
cos 1t |
|
|
|
|
|
|
|
|
|
|
|
|
y(t) |
900 |
sin 1t |
|
|
|
|
T |
|
2 |
|
х |
|
Кв |
|
|
|
0 |
|
|
Рис. 1. Структурная схема приемника.
16
|
|
|
Из формул (52) и (53) видно, что оценки 1 |
, 2 |
некоррелированы, а, |
следовательно, и независимы, т.к. распределение каждой из этих случайных величин нормальное. Дисперсии этих оценок равны:
|
D 1 D 2 N0 |
T . |
(61) |
Средние значения этих оценок в силу их несмещенности равны |
|
||
|
|
|
|
|
m 1 Acos , |
m 2 Asin . |
(62) |
|
|
|
|
Оценки A и получены при условии, что k является целым числом. |
|
||
Задание 1. Нарисовать структурную схему приемника для оценки фазы гармонического сигнала. Схему построить по формулам (58) или (59).
7. Потенциальная точность оценки параметров сигнала
Под потенциальной точностью оценок параметров радиосигнала понимают нижнюю границу Рао – Крамера для дисперсии неслучайного параметра, то есть оценки максимального правдоподобия. Потенциальная точность характеризует тот предел точности оценивания, который может быть достигнут только в результате обработки ряда значений Y, то есть без учета априорной информации.
Потенциальная точность оценки параметра радиосигнала обратно пропорциональна информации по Фишеру (26), элементы которой вычисляются как среднее значение вторых производных от функции правдоподобия по оцениваемому параметру. То есть
|
|
2 |
ln f y |
|
y t |
|
|
|
|
|
|||||
D |
M |
|
|
|
|
|
|
|
2 |
||||||
эфф. |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
. |
(63) |
||
|
|
|
|
|
|
|||
|
2 |
ln f y |
|
y t |
|
|||
|
|
|||||||
M |
|
|
|
|
|
|
|
|
|
2 |
|
|
|||||
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
Так, например, потенциальная точность оценки амплитуды (А) радиосигнала определяется соотношением
D |
|
|
|
|
|
1 |
|
|
|
|
1 |
|
N0 |
. |
(64) |
|
|
|
|
2 |
ln f y А y t А |
T |
|
|
|
||||||||
|
|
|
2 |
|
|
T |
|
|||||||||
|
А. |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
2 S (t)dt / N0 |
|
|
|
|
|||
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
А2 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Потенциальная точность оценки фазы ( 0 ) радиосигнала определяется соотношением
D |
|
|
|
1 |
|
|
|
N0 |
|
1 |
, |
(65) |
||
2 ln f |
|
|
y t 0 |
|
|
|||||||||
0 |
|
y |
|
|
2E 2q |
|
|
|||||||
|
|
|
|
|
||||||||||
|
|
M |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
0 |
2 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
17 |
где q = E/N0 – отношение сигнал/помеха.
Потенциальная точность оценки задержки ( з ) радиосигнала по огибающей определяется соотношением
D
з
|
|
где 2 0 |
2 |
|
|
|
1 |
|
|
|
|
N0 |
|
1 |
, |
(66) |
||
2 ln f |
y |
|
y t |
з |
2E 2 |
2q 2 |
||||||||
|
|
|
|
|
||||||||||
|
M |
|
|
з |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
з |
2 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
U ( j ) 2 d / U ( j ) 2 d - нормированный второй момент
энергетического спектра сигнала, используемый в качестве меры ширины спектра сигнала; 0 - несущая частота радиосигнала.
Потенциальная точность |
оценки |
задержки |
( з ) |
синусоидального |
|||||||||||||
радиосигнала по фазе принятого радиосигнала определяется соотношением |
|||||||||||||||||
D |
|
|
|
|
1 |
|
|
|
|
N0 |
|
|
1 |
, |
(67) |
||
|
|
|
|
|
y t |
з |
|
2 |
2 |
||||||||
з |
|
2 ln f |
y |
|
|
2E |
|
|
2q 0 |
|
|||||||
|
|
|
|
|
|
||||||||||||
|
|
M |
|
|
|
з |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
з |
2 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
то есть обратно пропорциональна отношению сигнал/помеха и квадрату |
|||||||||||||||||
несущей частоты радиосигнала. |
Обычно |
0 |
|
, поэтому |
точность оценки |
||||||||||||
задержки ( з ) синусоидального радиосигнала по фазе существенно выше точность оценки задержки ( з ) радиосигнала по огибающей.
Потенциальная точность оценки задержки ( з ) синусоидального радиосигнала по частоте и фазе сигнала дает лучший результат и определяется соотношением
|
|
|
|
D |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
(68) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
2q( 0 |
2 |
|
2 |
) |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
з |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Потенциальная точность оценки частоты ( 0 ) радиосигнала определяется |
|||||||||||||||||||||||||
соотношением |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
N0 |
|
|
1 |
|
, |
(69) |
|||
|
|
|
|
|
2 ln f |
|
|
|
y t 0 |
2E |
2 |
2q |
2 |
||||||||||||
|
|
0 |
|
|
|
y |
|
0 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
0 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
T |
|
|
|
1/ 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где (2 t)S2 (t)dt / S2 (t)dt |
- среднеквадратическая длительность сигнала. |
||||||||||||||||||||||||
0 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Таким |
образом, |
|
потенциальная |
|
|
аддитивного |
|
белого |
гауссов- |
||||||||||||||||
ского шума N0 частоты ( 0 ) радиосигнала обратно пропорциональна отношению сигнал/помеха и квадрату среднеквадратической длительности сигнала.
18
Потенциальная точность совместной оценки частоты ( 0 ) и задержки ( з ) радиосигнала по огибающей определяется соотношением
D |
D |
|
1 |
|
. |
(70) |
|
|
|||||
(2q ) |
2 |
|||||
0 |
з |
|
|
|
|
Отсюда следует, что повышение точности совместных оценок можно достигнуть увеличением отношения сигнал/помеха и (или) произведением эффективной длительности сигнала на эффективную ширину его спектра.
Задание 2. Построить зависимости точности (дисперсии) оценок фазы, задержки, частоты, совместной оценки частоты и задержки радиосигнала при значениях q= E/N0 = 1, 5, 10, 30, 50, 100, 1000 и фиксированных значениях: =1000, 0 =1000000,= 0,01.
Зависимость точности оценки амплитуды радиосигнала построить при длительности радиоимпульса Т= 1,0 мкс; 10,0 мкс; 100,0 мкс; 1,0 мс; 10,0 мс; 20,0 мс и фиксированном значении спектральной плотности мощности помехи 5x10-5 Вт/Гц.
Задание 3. Алгоритм работы приемника системы передачи информации по оценке параметра сигнала
Постановка задачи.
Есть две системы передачи информации:
-аналоговая система передачи информации, осуществляющая передачу информационного параметра λ в аналоговой форме с использованием двух видов модуляции несущего колебания u0 (t): амплитудной АМ, частотной ЧМ;
-цифровая система передачи информации, осуществляющая передачу информационного параметра λ без помехоустойчивого кодирования с использованием двухпозиционной фазовой манипуляции ФМн радиоимпульсами длительностью Т.
По каналу каждой системы передачи информации, имеющему полосу пропускания F, требуется передать информацию от датчика, измеряющего параметр λ удалённого объекта контроля. Параметр λ имеет нормальное распределение, усечённое границами ± 3а, с M(λ)=0 и с дисперсией σ2 = D (λ). В канале действует помеха, представляющая собой аддитивный белый гауссовский шум с равномерной спектральной плотностью мощности N0.
Приемники систем передачи определяют максимально правдоподобную оценку передаваемого параметра . Критерием сравнения качества приема информационного параметра λ по каналам систем передачи информации является
19
:^
точность оценки передаваемого параметра D( ), которая зависит от соотношения средней мощности принимаемого сигнала к мощности аддитивного белого
гауссовского шума (Pс/Рш). |
|
|
|
|
|
|
||
Требуется |
сравнить |
между |
собой |
зависимости |
точности оценки |
|||
|
|
|
^ |
|
|
|
|
|
передаваемого |
параметра |
D( )= |
f(Pс/Рш) |
для двух систем передачи при |
||||
изменении отношения Рс /Рш в диапазоне 2-25 и сделать выводы. |
||||||||
Исходные данные |
|
|
|
|
|
|
||
Поток АТ: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F, кГц |
|
|
T, мс |
|
|
D( ) |
N0, Вт/Гц |
|
3,5 |
|
|
2,5 |
|
|
8 |
5x10-3 |
|
Поток АС: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F, МГц |
|
|
T, мкс |
|
|
D( ) |
N0, Вт/Гц |
|
1,2 |
|
|
1,8 |
|
|
9 |
5x10-6 |
|
Значение параметра a при использовании амплитудной модуляции определяется соотношением
a |
|
1 |
|
. |
|
||||
D( ) |
При использовании частотной модуляции значение параметра a определяется соотношением
a |
|
F |
|||
|
|
|
|
. |
|
|
|
|
|
||
3 |
|
D( ) |
|||
Алгоритм работы приемника
Решающим правилом для приемника является уравнение правдоподобия:
|
|
|
|
df ( y(t) / ) |
|
0 |
|
|
(71) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
d |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Помеха n(t) на входе приемника является аддитивной и гауссовской, |
||||||||||||||
следовательно |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
1 |
T |
|
|
|
||
|
f ( y(t) / ) exp |
|
|
|
|
( y(t) s(t, ))2 dt |
(72) |
|||||||
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
N0 |
0 |
|
|
|
||
Подставляя выражение для f ( y(t) / ) в (71), получим |
||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
d |
T |
|
|
|
|
|
|
|
|||||
|
|
exp |
|
0 |
( y(t) s(t, )) |
dt |
|
0 |
||||||
|
|
|
||||||||||||
|
d |
|
N0 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
20 |
