Гл а в а VIII. ОЦЕНИВАНИЕ ПАРАМЕТРОВ ЛИНЕЙНЫХ МОДЕЛЕЙ ДВИЖЕНИЯ КОСМИЧЕСКИХ ОБЪЕКТОВ ПРИ НАЛИЧИИ ОШИБОК ИЗМЕРЕНИЙ В ОПЫТНЫХ ЗНАЧЕНИЯХ КООРДИНАТНЫХ ФУНКЦИЙ
1
В предыдущей главе было получено общее выражение для функции правдоподобия в задачах оценки параметров моделей движения космических объектов при наличии ошибок измерений в опытных значениях координатных функций. Функция правдо подобия была записана при достаточно общих предположениях о виде модели движения космического объекта и о законе рас пределения случайного вектора ошибок измерений. По этой при чине результаты предыдущей главы могут рассматриваться как основа для получения конкретных алгоритмов решения задач оценки параметров, встречающихся в практике эксперименталь ной космической баллистики.
Большой класс задач оценки параметров составляют такие, в которых вектор оцениваемых параметров входит в уравнение модели движения космического объекта линейно. Такие модели движения являются в экспериментальной космической баллис тике основными. В большинстве случаев к ним путем линеари зации могут быть приведены и нелинейные модели. Важное преимущество линейных моделей движения заключается в том, что для оценки их параметров могут быть получены сравнитель но простые алгоритмы, очень часто конечные.
Рассмотрению алгоритмов оценивания параметров такого класса моделей движения космических объектов посвящена на стоящая глава. Здесь, как и в предыдущей главе, полагаем, что опытные значения части координатных функций содержат ошиб ки измерений, поэтому основой для получения этих алгоритмов является общее выражение для функции ,правдоподобия, полу ченное в этой главе.
Вектор ошибок измерений принят распределенным по нор мальному закону с нулевым вектором средних значений и изве
стной корреляционной матрицей. Такое допущение об условиях эксперимента является в экспериментальной космической бал листике общепринятым и здесь принимается без дополнительных
обоснований.
Кроме алгоритмов оценивания параметров модели движения космического объекта, основанных на применении к решению задач оценки параметров при наличии ошибок измерений в опытных значениях координатных функций метода максималь ного правдоподобия, в настоящей главе обсуждаются и некото рые другие алгоритмы. Эти алгоритмы основаны на модифика ции алгоритмов оценивания параметров при точно заданных значениях координатных функций. Рассматривается также воз можность применения для оценки параметров моделей движения космических объектов способа статистик Вальда.
§8.1. МОДЕЛЬ ДВИЖЕНИЯ, СХЕМА ЭКСПЕРИМЕНТА
ИОСНОВНЫЕ ДОПУЩЕНИЯ
Уравнение модели движения космического объекта, линейное относительно вектора оцениваемых параметров, может быть за писано в следующем виде:
|
|
Р |
|
|
|
|
х (Л= 2 ?А ( 0 + 2 |
Ь + Л (0, ' |
|
(8.1.1) |
где x(t) — ордината |
модели движения; %j(t) и |
— ее |
коор |
динатные |
функции; |
<7j — постоянные |
параметры |
модели |
дви |
жения. |
|
|
|
|
|
Значения ординаты модели движения и координатных функ |
ций lj(t) |
(/= 1, 2, ..., |
р) в дискретные моменты времени U (г'=1, |
2, ..., N) доставляются проводимыми в процессе эксперимента измерениями/а значения в те же моменты времени координат ных функций (/= 1, 2, ..., s) известны точно.
Измеренные значения ординаты модели движения и коорди натных функций tj(t) содержат аддитивные ошибки измерений. Вектор ошибок измерений h* удобно представить в виде матри цы-столбца ciV (p+l) элементами:
, ( 8. 1.2)
где е* — подвектор вектора h*, содержащий своими элементами случайные ошибки измерений, имеющие место при измерении в моменты времени ti значений х* ординаты модели движения x{t):
# II * * |
* 1|т |
1 |
(8.1.3) |
®—I] е1$2. • • £лч| |
а б* — подвектор вектора Н*, содержащий своими элементами случайные ошибки измерений, имеющие место при измерении в те же моменты времени значений координатных функций мо дели движения £ j(0 :
|| 3iT ЗГ i |
- • - |
(8.1.4) |
$*) = 118Л3/2 • • • |
1 |
(/= 1 , 2, |
. .., р). |
Корреляционная матрица Ви вектора h* содержит своими элементами дисперсии ошибок измерений, а также изохронные изоканальные и перекрестные корреляционные моменты связи между ошибками измерений:
В принятых обозначениях выражение для плотности вероят ностей случайного вектора Л*, распределенного по нормальному закону с нулевым вектором средних значений, имеет вид
П |
1 . |
|
/?(A*) = (2 itf Т | 5 ЛГ У ехр(|— -Хк * тВ н 1к ^ , |
(8.1.6) |
где n — N(p+ 1) — размерность случайного вектора Л*.
Вектор q неизвестных* постоянных параметров |
модели дви |
жения удобно разбить на два подвектора: |
|
q = \\aT\bx\\. |
(8.1.7) |
Компонентами подвектора а вектора q являются постоянные параметры модели движения, служащие в ее уравнении множи телями при координатных функциях \j(t) :
а = \ а га2 . . . ар ||т,
а компонентами подвектора & вектора q — постоянные парамет
ры модели движения, служащие |
в ее уравнении множителями |
при координатных функциях |
|
|
& H IM 2 |
М т- |
Компоненты qj вектора q(j=U |
2, ..., р, р+ 1, ..., p + s) связа |
ны с компонентами Cj и bj |
векторов а и Ь следующими очевид |
ными равенствами: |
|
|
Q\— ai> |
Qp+\~bi > |
q |
qp+2 ~ b 2, |
ЯР*=аг |
qP+s= bs- |
В этих обозначениях компонент вектора оцениваемых пара метров уравнение модели движения (8.1.1) примет вид
х + := 2 |
00+ 2 |
+• |
(8.1.8) |
J-1 |
7 = 1 |
|
|
В практике обработки результатов измерений могут иметь место различные частные случаи линейной относительно вектора оцениваемых параметров модели движения (8.1.8). Так, при s = = 0 получается уравнение модели движения вида
* 0 0 = 2 + ^ ’ |
(8.1.9) |
7=1 |
|
характерной особенностью которой является то, что опытные значения всех входящих в нее координатных функций достав ляются проводимыми в процессе эксперимента измерениями и содержат, следовательно, ошибки измерений. При р = 0 модель движения (8.1.8) принимает вид
* + = 2 V 0 ( + |
(8.1.10) |
7 = 1 |
|
Особенностью этого частного случая является то, что значе ния в точках U всех входящих в нее координатных функций из вестны точно. Задача оценки параметров с уравнением модели движения вида (8.1.10) является обычной линейной задачей регрессионного анализа. Эта задача рассмотрена в гл. VI.
Очевидно, что наиболее общей является линейная ’ модель движения космического объекта вида (8.1.8). Частными случая ми алгоритма решения задачи оценки параметров с моделью движения вида (8.1.8) будут алгоритмы решения задач оценки параметров с моделями движения вида (8.1.9) и (8.1.10).
§8.2. ФУНКЦИИ ПРАВДОПОДОБИЯ
СРАЗЛИЧНЫМИ ВЕСОВЫМИ ФУНКЦИЯМИ
Впредыдущей главе было показано, что для получения функ ции правдоподобия необходимо записать плотность вероят ности вектора измерений г. Эта плотность вероятности может
быть получена из выражения (8.1.6) для плотности вероятно стей случайного вектора ошибок измерений h*. Реализация век тора п В проведенном эксперименте может быть выражена че
рез вектор измерений z и вектор измеряемых величин t |
следую |
щим образом: |
1 |
h = z — ls. |
(8 .2 .1 ) |
Подстановка (8.2.1) в (8.1.6) дает выражение для плотности вероятности реализации z случайного вектора г*, т. е. для плот ности вероятности измерительной информации, полученной в проведенном эксперименте:
p{Z, £) = (2я) В , ехр ~ ( г - ф В й 1{ г - Ъ 5.2.2)
Выражение (8.2.2) определяет значение плотности вероятно сти вектора измерений z при фиксированном значении измеряе мого вектора £. Для записи функции правдоподобия необходимо иметь плотность вероятности ведтора измерений z, выраженную
через вектор оцениваемых параметров q |
и вектор |, компонен |
тами которого являются точные значения |
в моменты време |
ни ti координатных функций %j(t). Этого можно достичь, если записать через векторы q и %вектор измеряемых величин
Для выражения вектора £ через векторы q и \ воспользуемся уравнением модели движения в виде (8.1.8). Это уравнение, за писанное для каждого из моментов времени ti, дает следующую систему равенств:
7 - 1 • 7 - 1
*2= 2 ^ 7 2 + 2 & ;V > |
(8.2.3) |
7 - 1 |
7 = 1 |
|
РS
*N = 2 ar)N+ 2
В левых частях равенств (8.2.3) стоят компоненты подвекто ра х измеряемого вектора £. Правые части этих равенств пред ставляют собой выражения компонент подвектора х измеряемо го вектора £ через компоненты векторов и §, введенных ранее, а также вектора ф, который введем равенством
ф = |
II ф! j фг I ■• • [ф*Г, |
' |
(8.2.4) |
где ф;- = I ф;-{ф;2 . . . фду||т- |
(7 = 1,2, . . . , s). |
|
|
Из (7.3.3) следует, что второй подвектор измеряемого векто ра £ равен вектору £ и не зависит от вектора оцениваемых па раметров q .
Для записи вектора £ через векторы q и %введем следующие
матрицы: |
|
A e = \ A \ E Np\, |
(8.2.5) |
5 0= | | Я | О лг*хлгЛ , |
I8-2-6) |
где А — прямоугольная матрица, имеющая Np строк и N столб цов;
а\e n
^а2En
apEN
В — прямоугольная матрица, имеющая Ns строк и N столб цов;
b,EN
N
ENp — единичная матрица порядка Np; Ошх Np — нулевая мат рица, имеющая Ns строк и Np столбцов. В выражениях для матриц Л и В EN— единичная матрица порядка N.
С помощью введенных матриц АЕ и В0 измеряемый вектор £ может быть записан через векторы q , \ и г): следующим образом:
Вектор | входит в выражение для измеряемого вектора £ в- явном виде, компоненты же вектора q — через элементы матриц Ае и В0. Равенство (8.2.7) позволяет записать плотность вероят
ности вектора измерений 2 при фиксированных значениях векто ров q и \\
__п_ _ 1 |
|
р ( г , Ъ, д) = {2п) 2 \ B h\ |
2 Х |
|
Хехр - ± ( z - A U - В Ц у ВТг1{2- |
A l l - ВЦ) |
(8. 2.8) |
Плотность вероятности векторам в виде (8.2.8), подставлен ная в выражение (7.3.6), дает возможность получить функцию правдоподобия в задачах оценки параметров линейной относи тельно вектора оцениваемых параметров модели движения при
наличии ошибок измерений в опытных значениях координатных функций:
