Файл: Брандин, В. Н. Основы экспериментальной космической баллистики.pdf

тельную информацию. Так v — подвектор вектора г , имеющий своими компонентами измеренные значения щ в точках ^ орди­ наты модели движения космического объекта x(t), a w — под­ вектор вектора измерений z, компонентами которого являются измеренные значения в точках ti координатных функций

модели движения космического объекта l\(t), £г(0> •••> £р (0- Это может быть представлено следующей записью:

Вектор z представляет'собой измеренное значение вектора £, компонентами которого являются точные значения в точках U ординаты модели движения космического объекта x(t) и ее ко­ ординатных функций h(t), l2(t), ..., l P(t):

С=11* т ГГ;

X = I х {х 2.

(7.3.3)

Поскольку все компоненты вектора £ измеряются с ошибка­ ми, полная плотность вероятности вектора г, как показано выше (см.. § 7.2), должна быть записана в виде

Здесь p(z, £) — плотность вероятности вектора г при фикси­ рованном положении измеряемого вектора £. Поскольку вектора содержит всю полученную в опыте измерительную информацию, эта плотность вероятности позволяет учесть все возможные кор­ реляционные связи между ошибками измерений. Функция ср(£) в выражении (7.3.4) — по-прежнему весовая функция, приписы­ вающая различные в общем случае веса различным гипотезам о положении измеряемого вектора £ в области его возможных зна­ чений. Область Z возможных значений измеряемого вектора £ образуется непрерывным множеством векторов, удовлетворяю­ щих своими одноименными по второму индексу компонентами модели движения (7.1.1). Такое определение области Z возмож­ ных значений измеряемого вектора £ позволяет записать следую­ щее соотношение между его подвекторами х и %:

■ * = / ( 5 , Я \

202

отк уд а следует, что в области Z

C=C(S, я).

Последнее равенство позволяет записать выражение для пол­ ной плотности вероятности вектора измерений г в следующем виде:

Функция правдоподобия, соответствующая плотности веро­ ятности вектора измерений z в форме (7.3.5), имеет вид

бия; k — положительный нормирующий коэффициент, не завися­ щий от вектора варьируемых параметров функции правдоподо­ бия i' .

Функция правдоподобия (7.3.6) может служить основой для получения алгоритмов обработки результатов измерений в зада­ чах оценки параметров моделей движения космических объектов при наличии ошибок измерений в опытных значениях координат­ ных функций. Методика записи этой функции правдоподобия является достаточно общей. С ее помощью могут быть учтены все возможные корреляционные связи между ошибками измере­ ний, т. е. изохронные, изоканальные и перекрестные.

При записи функции правдоподобия (7.3.6) не накладывалось никаких ограничений ни на вид модели движения космического объекта, ни на закон распределения ошибок измерений. Поэтому изложенную выше методику можно считать достаточно универ­ сальной. Следует отметить, однако, что функция правдоподобия (7.3.6) записана еще не в окончательном виде. Для приведения ее к виду, удобному для практического использования, необхо­ димо вычислить интеграл по области Z. Успех вычисления этого многомерного интеграла зависит от вида плотности вероятности вектора измерений z, вида модели движения космического объ­ екта, а также вида весовой функции ср(|, q).

§ 7.4. ВЕСОВЫЕ ФУНКЦИИ

Весовая функция введена 'в выражение для плотности веро­ ятности вектора измерений .г для того, чтобы можно было учиты­ вать различную в общем случае возможность получения вектора измерений z при измерении различных векторов £ из области Z их возможных значений. При известной модели движения косми­ ческого объекта измеряемый вектор £ полностью определяется вектором |, компонентами которого являются измеряемые значе-

203

ния координатных функций Поэтому весовая функция ф(^, q) определяет значимость гипотезы \ о положении в области Z измеряемого вектора £. При этом весовую функцию можно запи­ сать в виде, содержащем явно вектор | и вектор оцениваемых параметров q, что и сделано в § 7.3. Весовые функции должны выбираться, как правило, с учетом конкретного вида модели дви­ жения космического объекта, закона распределения ошибок из­ мерений и особенностей проводимого эксперимента. Вместе с тем по поводу вида весовых функций можно высказать ряд предва­ рительных общих соображений.

Поскольку весовая функция в виде <р(|, q) приписывает веса различным гипотезам § о положении в области Z измеряемого вектора £, прежде всего можно ввести такую весовую функцию, которая всем гипотезам о положении измеряемого вектора при­ писывает одинаковые веса. Талую весовую функцию можно наз­ вать равномерной и ввести следующим образом:

Равномерная весовая функция задает одинаковые веса раз­ личным гипотезам о положении измеряемого вектора либо во всей области Z, либо только в части ее. Определение области задания весов этой функцией производится в соответствии с конкретными условиями эксперимента. Учитывается область за­ дания весов при вычислении многомерного интеграла путем на­ значения соответствующих пределов интегрирования.

Очевидно, что равномерная весовая функция является край­ ним случаем возможных весовых функций, поскольку с ее по^ мощью неделается практически никаких различий между раз­ ными гипотезами о положении измеряемого вектора из области его возможных значений. Другим крайним случаем возможных весовых функций является, по-видимому, такая, которая припи­ сывает бесконечный вес действительному значению измеряемого вектора в проведенном эксперименте и веса, равные нулю, — всем остальным значениям измеряемого вектора из области его воз­ можных значений. Такой весовой функцией может служить дель­ та-функция с векторным аргументом \ —

Сравнивая введенный здесь вектор' %й с вектором |, введен­ ным ранее равенством (7.3.3), можно сделать вывод, что эти век-

204

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

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

Вектор §6 содержит своими компонентами зафиксированные дельта-функцией точные значения в моменты времени ti коорди­ натных функций £3-(/). Эти значения в области возможных зна­ чений являются единственными. В связи с этим при использова­ нии дельта-функции в качестве весовой возникают дополнитель­ ные трудности, связанные с увеличением вектора оцениваемых параметров. Это увеличение математически является следствием фильтрующего свойства дельта-функции, заключающегося в том, что

оо

— оо

где /( |) — произвольная скалярная функция вектора |.

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

Между рассмотренными выше крайними .случаями весовых функций находится класс таких весовых функций, которые раз­ личным гипотезам о положении в области возможных значений измеряемого вектора приписывают различные ограниченные и не равные нулю веса. Из этого класса весовых функций представ­ ляет интерес функция, определяемая равенством

Эта весовая функция гипотезе \ о положении в области Z из­ меряемого вектора £ приписывает вес, равный плотности веро­ ятностей случайного вектора .г* в точке z при расположении центра рассеивания вектора z* в точке £, определяемой векто­ ром \ и моделью движения космического объекта. Исходя из физического смысла весовой функции определяемая равенством (7.4.3) функция может быть названа вероятностной.

По смыслу задачи оценки параметров модели движения кос­ мического объекта вероятностная весовая функция является наиболее корректной из указанного промежуточного класса функций, поскольку она наиболее полно учитывает конкретные условия эксперимента, информация о которых содержится в вы­ ражении для плотности вероятностей ошибок измерений.

205

§ 7.5. СМЕЩЕНИЕ ОЦЕНОК НЕИЗВЕСТНЫХ ПАРАМЕТРОВ

МОДЕЛИ ДВИЖЕНИЯ КОСМИЧЕСКОЕО ОБЪЕКТА

Одним из основных статистических свойств оценок неизвест­

ных параметров является смещение. Оценку q вектора оцени­ ваемых параметров q можно рассматривать как реализацию

случайного вектора q *, случайность которого определяется слу­ чайностью вектора ошибок измерений h*. Смещение вектора оце­ ниваемых параметров определяется как разность между матема­

тическим ожиданием случайного вектора q* и точным значением <7 вектора оцениваемых параметров. При таком определении сме­ щение является вектором той же размерности, что и вектор оце­ ниваемых параметров. Компонентами вектора смещений будут смещения оценок отдельных параметров модели движения кос­ мического объекта.

Выше было указано, что при подстановке в функцию правдо­ подобия вместо вектора варьируемых параметров вектора оце­ нок неизвестных параметров эта функция принимает наибольшее свое значение. Символически это условие можно записать так:

Точи наибольшего значения функции правдоподобия являет­ ся точкой ее экстремума, поэтому равенству (7.5.1) соответству­ ет уравнение правдоподобия вида

XQ')

л = 0. (7.5.2) dq’ q’=q

В общем случае уравнение правдоподобия (7.5.2) векторное. Неизвестным в этом уравнении является вектор q-

Следует отметить, что если равенство (7.5.1) определяет оцен­ ку вектора q однозначно, то уравнение правдоподобия (7.5.2) этим свойством не обладает. Последнее связано с тем, что ра­ венство (7.5.2) справедливо не только для точки наибольшего значения функции правдоподобия, но и для всех экстремумов этой функции, которых в общем случае может быть несколько. По этой причине отыскание оценки вектора неизвестных пара­ метров является задачей отыскания такого корня уравнения правдоподобия, при котором функция правдоподобия принимает наибольшее свое значение. Иначе эта задача может быть сфор­

мулирована как задача • отыскания абсолютного максимума функции правдоподобия.

Уравнение правдоподобия может быть либо линейно, либо нелинейно относительно искомого вектора оценок. В первом случае уравнение правдоподобия имеет конечное аналитическое ре-