
книги из ГПНТБ / Васильев Г.А. Повышение эффективности комплексной автоматизации
.pdfтипов разложений по сравнению с нулевым. Для этого вычислим отношения
|
min t r i i ( z ) / min m 0( z ) , min (Ji(z)/min a o (z), |
где t = |
1, 2, 3, кдр. При этом предполагается, что коли |
чество |
элементов разложения и коэффициент шумового |
воздействия неизменны.
о |
z |
¥ |
6 |
£ 8 |
Рис. 5.24.
Сравнение эмпирических формул для минимальных значений математического ожидания линейной ошибки различных типов разложений (5.7), (5.46), (5.49), (5.52), (5.19) показывает, что все они могут быть представлены в виде
min т . { ( г ) = р min (а) (1 + а£). |
(5.55) |
150
Обобщенная эмпирическая формула для минимальных значений среднеквадратического отклонения линейной ошибки может быть получена в результате сравнения формул (5.8), (5.47), (5.50), (5.53), (5.20):
min а* (г) = р min з,- (а)(1 |
сЩ. |
(5.56) |
||
Эмпирические |
формулы (5.55) и (5.56) позволяют |
|||
легко оценить искомый выигрыш. Он равен |
|
|||
m inm 4 (z) |
m lnm 4 (a) |
minu4 ( z ) __гшпа4(а) |
|
|
min m 0 (z) |
minm0 (a) |
’ m in s0 (z) |
min o0 (a) |
|
вследствие независимости числовых коэффициентов a, b и с от типа разложения.
Полученный результат позволяет сформулировать важное правило. Относительное уменьшение числовых характеристик линейных ошибок сложных модификаций разложений по отношению к разложению нулевого типа определяется качеством выравнивания размера элемен тов разложения и не зависит от величины шумового воздействия.
5.5.3. Оптимальное количество элементов разложения при наличии шума
Результаты, рассмотренные в предыдущих парагра фах, показывают, что минимумы математического ожи дания линейной ошибки и ее среднеквадратического от клонения достигаются практически при одних и тех же значениях аргумента а. Следовательно, эти критерии мо гут быть обобщены в критерий потенциальной точности. Поэтому в дальнейшем будем рассматривать критерий минимума среднеквадратического отклонения линейной ошибки, отождествляя его с критерием потенциальной
точности.
Точность зависит от общего количества элементов разложения, поскольку в формулах (5.55) и (5.56) N входит в выражение для масштабного коэффициента р и для коэффициента шумового воздействия g:
р = ] / 2’kD2/N, £ = 2 yii3mN/Tобз. |
(5.5а) |
Увеличение общего количества элементов разложения приводит к уменьшению р и, следовательно, к умень шению величины линейной ошибки. В то же время с ро-
151
стом N увеличивается коэффициент шумового воздей ствия, и, следовательно, имеет место тенденция увеличе ния числовых характеристик. Очевидно, существует оптимальное количество элементов разложения N0UT, при котором минимальные значения числовых характе ристик линейной ошибки (5.55) и (5.56) достигают ми- нимум-миниморума.
Для определения iVonT подставим в формулу (5.56) значения р и £ из (5.5а):
minзг- (г) = у |
mmз,- (а) ( 1+ by |
— T~^f— h |
|
|
+ C |
7-овз )• |
|
Продифференцируем это выражение по N и приравняем |
|||
полученную производную нулю. Получим |
|
||
д rain з,- (z)fdN = |
1/vV -f- c2]/3 зш[ТоСа= |
0, с = 0,1. |
Отсюда оптимальное количество элементов разложения выражается следующим образом:
|
Л^опт= ^обз/0,35од1. |
(5.57) |
При этом |
коэффициент шумового воздействия £ОПт= 1 0 . |
|
Заметим, что оба параметра ЛГоит и |
| 0пт не зависят от |
|
типа разложения. |
потенциальной точ |
|
Таким |
образом, для достижения |
ности при создании системы должны быть выбраны оп тимальными: общее количество элементов разложения N-, коэффициент шумового воздействия аргумент а, соответствующий минимальным значениям функций чис ловых характеристик выбранного типа разложения.
5.5.4. Коррекция ошибок передачи
Изучение линейных ошибок передачи при наличи шу ма показывает, что математическое ожидание модуля ошибки отлично от нуля и возрастает при увеличении коэффициента шумового воздействия. Поэтому необходи мо рассмотреть вопрос о математическом ожидании фа зового угла линейной ошибки передачи и на основании его исследования сформулировать правила оптимальной коррекции, подобно тому, как это было сделано для кор-
152
рекции ошибок преобразования в гл. 4. Для сокращения объема рассмотрим только метод интегральной коррек ции.
Запишем выражение для математического ожидания фазового угла линейной ошибки разложения третьего типа при воздействии равномерно распределенного шу мового воздействия:
|
|
|
|
м |
|
|
|
Ad3^ ziD__2ш |
|
||||
|
|
|
|
|
2%, |
-arctg-^ ^ j - dx1dxidxidx4-\- |
■ |
|
Ш |
|
I |
ДЦ3ДР30 2 (1 — k 2) A d t |
|
|
|
|
||||
0 |
|
|
|
|||
|
0 |
kD —Ad„ |
|
|
||
|
|
|
|
2 |
|
|
|
|
|
|
|
Ad_ |
|
+ 2 |
4 i , |
2 |
D |
2 |
|
|
|
|
|
2xa |
arct« |
||
f |
f |
f |
I |
|||
|
|
0 |
0 |
ID - A J t |
|
|
|
|
|
|
|
|
(5.58) |
Исходное выражение для математического ожидания фа зового угла линейной ошибки при нормально распреде ленном шумовом воздействии имеет вид:
А4а др3 ID оо |
2 / п 2 \ |
2 х 3 ехр (- |
X 4/ 2 а ^) |
т .. |
= r X |
фН < - И П |
[ A d 3A$3D 2 ( \ — k 2) omV2n |
Об kD — 00
X |
arctg |
} Xs* 3 |
dx,dx2dx3dx4 -f |
|||
|
A§3_ |
*1 |
Г *4 |
|
|
|
|
|
_ |
. о / л |
2 \ |
||
|
2 |
D |
со |
|||
|
|
|
|
2 x 3 e x p ( — х ф 2 j ^ ) |
||
+ 2 I |
f |
,f |
I |
|
|
^ - x |
Ad3A%D2 {\ —k2) s mV2K |
||||||
0 |
0 |
ID |
—oo |
|
|
|
arctg— |
— dxidx2dx3dxi |
(5.59) |
||||
|
|
#1 “ Г #4 |
|
|
Интегрирование выражения (5.58) может быть выпол нено, хотя и сопряжено с чрезвычайно громоздкими вы числениями; интеграл (5.59) не выражается через эле ментарные функции. Поэтому целесообразно исследо вать математическое ожидание фазового угла линейной ошибки при равномерном и нормальном шумовом воз действиях методами математического моделирования.
153
В этом случае может быть использован алгоритм, описанный в и. 5.5.1. Вычисление значения фазового угла случайной линейной ошибки производится по формуле
<Pi= arctg[x2x3/ (Xi + * 4) ],
где j — номер реализации ошибки в серии испытаний; Xi — ошибка, возникающая под влиянием шумового воз действия канала связи, распределенного по равномерно му или нормальному закону.
Результаты анализа функций т (а, S) и т (а, S,),
полученных методами математического моделирования, показывают, что закон распределения шумового воздей ствия не имеет существенного влияния на величину ма тематического ожидания фазового угла линейных оши бок. Поэтому с достаточной для практических нужд сте пенью точности можно считать, что
«У (*. *) = " У К 9 = (a, S),
если только в обоих случаях величина коэффициента шу мового воздействия | неизменна.
Исследование аналитических выражений математи ческого ожидания фазового угла линейной ошибки при равномерном распределении шумового воздействия и со ответствующих им значений, полученных методами мате матического моделирования, при нормальном распреде лении шумового воздействия позволяет указать основные свойства функции:
5) = |
(at; |
limm .(а, 5)=.0; |
(5.60) |
£-►00 v |
|
дяу (а, £)/<?£< AS, |
|
где A — постоянное число.
Изменение параметров секторных разложений в пре
делах: 0 < А < 0 ,3 , 0 < / < 1 |
оказывает меньшее |
влияние |
на изменение функции |
(a, S), чем изменение |
коэффи |
циента шумового воздействия. Поэтому существует воз можность построить обобщенную функцию математичес кого ожидания фазового угла линейной ошибки в ре
зультате усреднения функций я у (а, |) |
для различных |
типов разложений и их параметров при |
const. |
154
На рис. 5.25 изображено семейство обобщенных функций т,у (я, £). Отсутствие индекса i в обозначении
функции означает ее независимость от типа разложения. В качестве параметра семейства взято значение коэффи
циента шумового воздействия. |
Средняя часть |
семейства |
функций т ч (а, |) при 2 < а < |
8 представляет |
основной |
интерес с точки зрения проектирования, так как она со
ответствует области минимальных значений числовых характеристик линейных ошибок и может быть с доста точной для практики точностью аппроксимирована отрез ками прямых:
M ,(*,5) = [0,4 + 0 ,1 3 (E -2 )]lg« +
+ [0 ,3 3 -0 ,1 7 (1 -2 )], |
\ > 0 . |
Поскольку функция т ф(я, 5) при |
&= const монотонна, |
то правила оптимальной интегральной коррекции, изло женные в гл. 4, могут быть использованы и в данном случае. Отличие обобщенной функции т (а, £) от функ
ции строго соответствующей выбранному типу разложе ния приводит к коррекции систематических ошибок, не сколько отличающейся от оптимальной. Неточность кор рекции определяется отношением
Дт1р;(я + )/т ф(яД) =
= И , (а, £) — т фг (я, Щ /тч(я, +
155
Максимальная величина А т . (а, S) примерно на порядок меньше, чем т (а, 5). Поэтому коррекция линейных оши бок передачи с использованием обобщенной функции т 9 (а. S) приводит на практике к удовлетворительным результатам.
5.6. Законы распределения линейных ошибок
В большинстве случаев для проектирования систем передачи радиолокационной информации достаточно рас полагать сведениями о математическом ожидании и среднеквадратическом отклонении линейной ошибки. Используя эти числовые характеристики, можно выбрать оптимальные параметры ситемы, при которых достижи ма ее потенциальная точность, выяснить механизм влия ния шума на точность передачи информации, выполнить коррекцию систематических ошибок.
Однако для полного представления о характере ли нейной ошибки необходимо знать закон ее распределе ния. Начнем рассмотрение с наиболее простого случая: распределения линейной ошибки преобразования при кадровом методе [65]. Так как все элементы разложения в этом случае равны, то закон распределения линейной ошибки в пределах одного элемента разложения дает исчерпывающую картину явления.
156
Элемент разложения кадрового типа изображен на рис. 5.26. В соответствии с обозначениями, принятыми на рисунке, будем считать в дальнейшем, что А^Кдр1< <Д&кдр2, так как это условие не ограничивает общности рассуждений. Минимальная линейная ошибка преобра зования равна нулю, если отметка цели совпадает с вер шиной элемента разложения, в которую переносятся все отметки, попавшие внутрь элемента. В данном случае эта вершина совпадает с началом координат. Макси мальная ошибка равна диагонали элемента разложе ния.
Функция плотности вероятности линейной ошибки
у = |
У |
+ |
х 2, |
разложения кадрового |
типа имеет |
вид |
|||||
|
|
|
У + 6 |
|
|
------------- ----- |
|
|
|||
®кдр(у)— ^ J |
дакдр (•*,, х 2)Ь(у — у |
х, |
х 2 )dx| dx2 , |
||||||||
|
|
О у—г |
|
|
|
|
|
(5.61) |
|||
|
|
|
|
|
|
|
|
|
|
||
где 8 (и) — импульсивная |
функция |
Дирака; |
G— область |
||||||||
изменения случайной величины х 2, такая что ^ х \ -\-х2 |
на |
||||||||||
ходится в e-окрестности |
точки у; |
o^pix,, х а) — функция |
|||||||||
плотности |
вероятности |
системы |
случайных, |
аргументов |
|||||||
х,, |
х 2. |
|
|
|
|
|
|
|
|
|
|
|
Произведя |
замену |
переменных |
t = |
"j/"х] -(- х 2 >вы_ |
||||||
числим внутренний интеграл в формуле (5.61): |
|
||||||||||
|
|
|
|
»+« |
1 |
|
|
|
tdt |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
^'кдр |
|
J А^кдр1^^кДр2 b{y — t) V П —*; |
|
|||||||
|
|
|
|
у—* |
|
У |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ас■*^КкдяДР1а^акдргКД 2 'VУ2 х2 |
|
|
|
|
|||
|
Полученное выражение проинтегрируем по переменной |
||||||||||
х,: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
КДР2 |
y d x 2 |
|
|
|
||
|
|
|
W. |
|
j |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
К Д Р ( й = |
|
j/" |
У2 |
|
|
|||
|
|
|
|
|
|
|
|
|
Этот интеграл вычисляем по отдельным участкам, в пре делах которых у меняется по одному и тому же закону. Таких участков три (рис. 5.26):
157
|
I. 0 ^ г /^ Д $ Кдр1, |
^:У\ |
|||
2. |
Д<7кдр1 ^ |
у ^^Ай?кдр2, |
Ой-^Хг^А^кдрЬ |
||
3. |
Дс/Кдр2 ^ |
< |
^ |
2 |
+ Д< |
|
|
|
Д^кдр1 |
кдр2 |
Y 7 + Дй*?КДркдп1 < x ^2 < A~ “ dи'КДР] 2-
Вычисление функции плотности вероятности на каж дом из этих участков дает:
ку
2А£/кдР1^^кдР2
|
У |
Ad.КДР1 |
(5.62) |
||
®кдр( у ) z |
^^КДР1^^КДР2 arcsin ■ |
У |
|||
у |
Ad,кдР2 |
. |
V y* — Ad2 |
x |
|
А^кдЯ^кдРг |
arcsin- |
a rc s in - |
у |
р |
Семейство нормированных функций плотности вероятно сти изображено на рис. 5.27 при параметре
!Окдр= Ай?кдрг/А<^кдр1-
При оптимальном разложении (аКдр=1) плотность ве роятности линейной ошибки близка к распределению Симпсона. При « Кдр— >оо плотность вероятности линей
ной ошибки приближается к равномерной вследствие преобладающего влияния одной из компонент, распре деленной согласно исходным условиям равномерно.
Для сокращения громоздких вычислений вопрос о распределении плотности вероятности линейных оши бок секторного разложения рассмотрим качественно. Ли нейная ошибка может изменяться от нуля до величины, равной диагонали наибольшего элемента разложения в случае, если шумовое воздействие отсутствует. Распре-
158
деление ошибки внутри t'-ro элемента разложения описы вается выражением (5.62). Параметром функции плот ности вероятности в этом случае является число щ, рав ное отношению сторон г-го элемента разложения. Для разложения секторного типа с определенным значением нормированного аргумента а величина а; является пере менной, так как соотношение между сторонами элемента разложения меняется с изменением его расстояния от начала координат. Поэтому изменяется и вид функции плотности вероятности ошибки внутри элемента разло жения с изменением его номера, как это показано на рис. 5.27, где параметру щ соответствует параметр аКдр.
Следует обратить внимание на то, что при cti<l боль шая сторона элемента разложения ориентирована вдоль радиуса зоны обзора для всех элементов разложения. При cti> 1 ориентация в центральной части зоны обзора сохраняется для всех элементов разложения вплоть до номера /, причем номером / обозначен элемент, имею щий форму квадрата. Для всех последующих элементов большая сторона ориентирована перпендикулярно радиу су зоны обзора.
Условная функция распределения плотности вероят ности линейной ошибки внутри г-го элемента разложе ния должна удовлетворять тождеству
У—УтахУ)
( Wi (yfi) dy— \,
6
где •wi (y/i) определяется формулой (5.62), в которой соотношение между сторонами элемента разложения рав но щ. Такому же тождеству должна удовлетворять иско мая функция распределения линейных ошибок в преде лах всего сектора разложения. Поэтому интегрирование условных функций распределения должно выполняться с весом, пропорциональным площади элементов разло жения, что является следствием принятой гипотезы о равномерном распределении целей в области зоны об зора:
D |
У— г г т а х б ) |
|
|
w (у) = j* |
f Wi (yfi) dydx3, |
|
|
6 |
О |
|
|
где утах(г') — область |
существования |
линейной |
ошибки |
в пределах г'-ro элемента разложения. |
На рис. |
5.28 по- |
159