
книги из ГПНТБ / Свешников А.А. Вероятностные методы в прикладной теории гироскопов
.pdfП Р Е Д И С Л О В И Й |
7 |
собе изложения неизбежно были бы повторения, поскольку раз личные гироскопические устройства часто описываются одинако выми уравнениями движения и для их исследования применяются тождественные математические приемы.
' ЧІІам представилось более рациональным сперва дать общий обзор основных типов ГУ и их уравнений движения, выделив та ким образом типичные математические задачи, которые возникают
вгироскопии, а затем изложить методы решения этих задач, ил люстрируя их применение анализом конкретных ГУ.
Вцелях сокращения объема книги и представления материала
вболее обозримой форме вывод уравнений движения различных гироскопических устройств не приводится, поскольку его можно найти в многочисленных книгах но гироскопии, изданных в по следнее время (в тексте даны необходимые ссылки на соответствую щие источники).
Авторы сочли невозможным приводить в книге изложение об щих вопросов теории вероятностей и теории случайных функций, предполагая, что эти сведения известны читателям. Однако было признано целесообразным дать общую сводку формул теории веро ятностей и теории случайных функций для того, чтобы облегчить чтение дальнейшего текста и избавиться от необходимости услав ливаться в обозначениях при изложении основного материала книги.
Некоторые специальные математические приемы, используе мые в книге, поясняются по мере их применения.
Предлагаемая книга является совместной работой обоих авто ров, и после ее написания оказалось невозможным указать, какие разделы книги написаны каким автором, хотя и можно отметить, что при изложении основных особенностей гироскопических устройств и их уравнений движения основную роль играл G. С. Ривкин, а изложение вероятностных методов исследования в основном было выполнено А. А. Свешниковым. Решение приме ров производилось авторами совместно.
Авторы стремились снабдить книгу полным перечнем известных им источников по применению вероятностных методов в гироскопии и указать общие труды по гироскопии, теории вероятностей и тео рии случайных функций, в которых читатель может найти разъяс нение формул, приводимых в книге без выводов.
В книге принята следующая нумерация глав, параграфов и пунктов: параграфы имеют свою нумерацию в каждой главе, так же как и пункты внутри параграфа. В номере параграфа первая цифра означает номер главы (например, § 3.4 означает 4-й параграф 3-й главы), пункты внутри каждого параграфа отмечаются после довательными номерами (например, пункт 1 § 3.4 означает 1-й пункт 4-го параграфа 3-й главы). При нумерации формул указы вается номер главы и порядковый номер в данной главе (например,
8 П Р Е Д И С Л О В И Е
формула (3.15) означает 15-ю формулу 3-й главы). При ссылках на формулы данной главы первая цифра для простоты опускается.
При написании книги авторы встретились с существенной труд ностью, связанной с выбором единой системы обозначений, так как в книгах по гироскопии и теории вероятностей используются различные обозначения. Более того, даже в различных работах по гироскопии не выдерживается единая система обозначений. По этому, для того чтобы упростить читателю обращение к различным источникам, на которые имеются многочисленные ссылки в тексте, в предлагаемой книге учитывались обозначения, принятые в ра боте, откуда была взята та или другая формула, отступая, однако, от этого принципа в тех случаях, когда по тем или иным причинам это представлялось целесообразным. Это привело к необходимости каждый раз пояснять принятые обозначения. Также не удалось полностью выдержать принятый во многих работах по теории вероятностей принцип, согласно которому большие латинские буквы обозначают случайные величины и функции, а малые ис пользуются для обозначения их реализаций.
Авторы выражают благодарность всем лицам, участвовавшим в обсуждении отдельных вопросов, изложенных в книге, и особенно благодарны Д. М. Климову, прочитавшему всю рукопись и сде лавшему ряд замечаний.
Авторы считают также приятным долгом поблагодарить А. Ю. Ишлинского, интерес которого к применению вероятност ных методов в гироскопии способствовал написанию этой книги.
Г Л А В А 1
СВОДКА ОСНОВНЫХ ФОРМУЛ ТЕОРИИ ВЕРОЯТНОСТЕЙ И ТЕОРИИ СЛУЧАЙНЫХ ФУНКЦИЙ
§ 1.1. Основные формулы теории вероятностей
Общий курс теории вероятностей предполагается известным читателю в объеме [20], [14], [1в], поэтому в данном параграфе при ведены только основные формулы теории вероятностей, которые используются в дальнейшем, и установлены основные обозначения.
Будем, по возможности, обозначать случайные величины боль шими буквами конца латинского алфавита X, Y, Z и т. д., отме чая их в случае необходимости индексами, например Xl7 Х 2, Х 2. Исключение из этого правила будет сделано только для тех детер минированных величин, для которых традиционно в гироскопии приняты большие буквы, а также для случайных величин, обозна чаемых греческими буквами.
Для обозначения возможных значений случайных величин или значений, получаемых в результате опыта, мы будем (по воз можности) использовать малые буквы того же наименования, т. е., например, х, у, z и т. д.
Полной характеристикой случайной величины X является ее функция распределения F(x), равная вероятности того, что слу чайная величина X будет меньше выбранного значения х. Функ
ция распределения определяется равенством |
|
F(x) = P (X < х}. |
( 1. 1) |
Плотностью вероятности /(ж) случайной величины X называется производная от функции распределения F(x) по ее аргументу, т. е.
( 1. 2)
Если случайная величина X непрерывна, то плотность ве роятности / (х) является или непрерывной функцией своего ар гумента или имеет только такие разрывы, что для любого беско нечно малого интервала (х, ж+Дж) интеграл от плотности / (х) имеет порядок величины Дж, т. е.
х-^Ах
^ / (хг) dx1= О (Дж) |
(1.3) |
10 О С Н О В Н Ы Е Ф О Р М У Л Ы [ГЛ. 1
где О (Ах) обозначает величину, удовлетворяющую условию
г |
О (Ах) |
(1.4) |
н т |
—і— - = const. |
|
Дх+0 |
Аж |
|
Если X является дискретной случайной величиной, могущей принимать значения Xj с вероятностями р . (/ = 1,2, . . ., гі), то f(x) распадается на сумму дельта-функций
/( * ) = |
2 |
Pjb(x — Xj), |
(1.5) |
|
|
|
J = 1 |
|
|
где функция 8 (я) определяется равенствами |
|
|||
о |
при |
X |
О, |
( 1. 6) |
8(* )= оо |
при |
X = О, |
Так как дельта-функция будет встречаться в дальнейшем изло жении, напомним ее основные свойства. Для любой непрерывной
в точке х = у |
функции ср (х) |
справедливо равенство |
|
||||
|
|
ъ |
|
|
|
|
|
|
|
I <р(я)8(у — x)dx = |
<?(у), |
(1.7) |
|||
|
|
а |
|
|
|
|
|
если а |
у <Г Ъ. |
|
|
|
|
|
|
Для |
дельта-функции справедливо интегральное представление |
||||||
|
|
|
|
СО |
|
|
|
|
|
8 (ж) = |
^ |
eiwxdiо, |
(1.8) |
||
|
|
|
|
— СО |
|
|
|
где интеграл |
понимается в смысле |
главного значения, т. |
е. как |
||||
|
|
|
|
L |
|
|
|
|
|
lim |
I eiu,xdu>. |
|
|
||
|
|
£-»•со |
". |
|
|
|
|
Наконец, дельта-функцию можно дифференцировать любое |
|||||||
число раз, причем |
|
|
|
|
|
||
|
|
|
|
|
|
СО |
|
|
|
8'", (*) = й Г 8 (* )!= ■ £ |
\ ^ d . , |
(1.9) |
|||
|
|
|
|
|
—СО |
|
|
а для любой функции ср (х), |
|
имеющей в точке х —у непрерывную |
|||||
производную до п-го порядка, справедливо равенство (а < |
у <С Ъ): |
||||||
|
|
ь |
|
|
|
|
|
|
|
^ ср (х) Ь(п) (у — x)dx = |
ср(я) (у). |
(1.10) |
а
§ 1.1] |
Т Е О Р И Я В Е Р О Я Т Н О С Т Е Й |
11 |
Функция распределения F (х) и плотность вероятности / (х) удовлетворяют следующим очевидным соотношениям:
0 < F ( x K l , |
F(— co) = |
ü, F ( + оо )= 1, |
(1.11) |
|
X |
|
|
f {х) > о, / (± ОО) = 0, |
(*) = j |
/ (*i) <telt |
|
|
— СО |
|
|
|
|
со |
|
|
|
^ f(x)dx = 1. |
(1.12) |
Вероятность попадания случайной величины X в интервал (а, Ь) может быть выражена как через функцию распределения, так и через плотность вероятности формулами
|
ь— |
|
р {а |
< 6} = F (b) — F(ä)= j f(x)dx, |
(1.13) |
|
а— |
|
причем для дискретных случайных величин существенно, что правая граница не включена в интервал, а для непрерывных слу чайных величин это не имеет значения, так как в соответствии с (3) в этом случае
Р { Х = 6 } = 0 . |
(1.14) |
В теории вероятностей принято условно считать, что вид функ ций распределения и плотности вероятности определяются аргу
ментами этих |
функций. |
Таким образом, функции F (х), F (у), |
F (z) и т. д. и / |
(х), / (у), / |
(z) и т. д. будут обозначать не одни и те же |
функции, взятые при различных аргументах, а функции распре деления и соответственно плотности вероятности случайных вели чин X, У, Z и т. д., которые могут иметь совершенно различные математические выражения. Мы также будем пользоваться этими обозначениями, отступая от них только тогда, когда это может привести к недоразумениям. В этих случаях функции F (х), / (х) и т. д. мы будем снабжать индексами, указывающими на наиме нование случайной величины.
Вместо полной характеристики случайной величины ее функ цией распределения или плотностью вероятности иногда доста точное представление о свойствах случайной величины можно получить, указав необходимое число первых моментов этой вели чины, т. е. указав математические ожидания целых степеней этой величины (начальные моменты) или математические ожидания целых степеней разности случайной величины и ее первого момента (центральные моменты). Начальный момент порядка у мы будем
12 |
О С Н О В Н Ы Е Ф О Р М У Л Ы |
[ГЛ . 1 |
обозначать m-, центральный момент того же порядка будем обозна чать f i т. е.
СО |
|
СО |
|
nij — I |
xjf(x)dx, |
pj == J (x — m1Yf(x)dx. |
(1.15) |
—CO |
—CO |
|
|
В том случае, |
когда речь идет о нескольких случайных величи |
нах, моменты mj и р . мы будем снабжать сверху индексом, указы вающим на наименование случайной величины, т. е. будем писать, например, т)х), р]ж) и т. д. В приложениях наиболее часто нахо дят применение первый начальный и второй центральный моменты. Первый начальный момент называют просто математическим ожи данием и обозначают большой буквой М перед случайной величи ной или малой буквой, соответствующей наименованию случайной величины с чертой сверху, т. е., например, М [X] или х для слу чайной величины X, М [F] или у для случайной величины Y и т. д. Второй центральный момент называется дисперсией и
обозначается (для случайной величины X) D [X] |
или сг; поло |
|||
жительный корень из дисперсии, |
обозначаемый ах, носит название |
|||
среднего |
квадратического отклонения. |
|
||
Таким |
образом, в соответствии с (15) имеем |
|
||
|
СО |
|
|
|
|
М [X] = х = |
( |
xf(x)dx, |
(1.16) |
|
— СО |
|
|
|
|
|
СО |
|
|
|
D[X] = o2 = |
j |
(.T- x f f { x ) d x . |
(1.17) |
Между дисперсией случайной величины и вторым начальным моментом существует соотношение
а| = тп2 — Ж2. |
(1.18) |
Интегралы (15) для отдельных законов распределения могут рас ходиться, начиная с некоторого значения /. В этом случае говорят, что случайная величина не имеет моментов, начиная с /-го.
Наиболее часто в приложениях встречается нормальный закон распределения, плотность вероятности для которого определяется формулой
|
{х—х)г |
/(*) = — |
(1.19) |
О* V2тс |
|
а функция распределения имеет |
вид |
( 1 . 2 0 )
г М = 4 [ і + ф ( і 7 г ) ;
§ 1.1] Т Е О Р И Я В Е Р О Я Т Н О С Т Е Й 13
где через Ф (х) обозначена интегральная функция |
Лапласа, |
опре |
||
деляемая равенством |
X |
£2 |
|
|
|
|
|
||
Ф (X) = |
j* е |
2 dt, |
(1.21) |
|
|
о |
|
|
|
таблицы которой можно найти в ряде работ (см., |
например, |
I11], |
[83]). Нормальный закон полностью определяется двумя момен тами: математическим ожиданием и дисперсией. Все моменты для этого закона существуют, причем все центральные моменты нечет ного порядка равны нулю, а четного порядка выражаются через дисперсию.
Совокупность случайных величин или, как говорят в теории вероятностей, система случайных величин может быть охарак
теризована |
функцией |
распределения (например, |
F (х , у, z) = |
|
-=Р {X И х, |
У <С у, |
Z <С z} — для системы случайных величин |
||
X, Y, Z) |
или многомерной полотностью вероятности (/ (х, у, z) = |
|||
=<d3F (X, |
у, |
z)/dx ду dz — в данном примере). |
вероятности |
|
Для независимых |
случайных величин плотность |
системы распадается на произведение плотностей вероятностей
отдельных случайных величин. Например, если X, Y, |
Z — не |
зависимы, то |
|
f(x, у, Z) = f(x)f(y)f(z). |
(1.22) |
Если значение одной или нескольких случайных величин, вхо дящих в систему, закрепить, то плотность вероятности оставшейся подсистемы случайных величин (условная плотность вероятности), например для системы двух случайных величин, определится формулой
/ ( * Ы = % Г ’ |
(1-23) |
где / (х I у) — условная плотность вероятности случайной вели чины X при фиксированном У; / (у) — плотность вероятности У, связанная с двумерной плотностью вероятности / (х, у) формулой
СО |
|
f { y ) = \ fH , y)dx. |
(1.24) |
— СО
Понятие моментов случайных величин сохраняется и примени тельно к системе величин. Наиболее важными для приложений являются первые и вторые моменты. Первые начальные моменты (математические ожидания) систем случайных величин Ху, Х 2, ■■■
. . Х п в соответствии с (16) определяются формулой
СО |
СО |
xj = $ xjf (xj) dxj = j |
j Xjf(xv x2, . .., xn) dxxdx2 . . . dxn, (1.25) |
— CO |
— CO |
14 |
|
О С Н О В Н Ы Е Ф О Р М У Л Ы |
[ГЛ. 1 |
а |
вторые |
центральные моменты (дисперсии), в соответствии |
|
с |
(17), — формулой |
|
|
|
СО |
|
|
|
~ 5 {X J |
f (x j) dx j = |
|
|
— СО |
со |
|
|
|
|
|
|
= $ |
• • • \ {Xj — Xj f f { xv х2, . . .,хп) dxxdx2 . . . dxn. |
(1.26) |
— СО
Кроме математических ожиданий и дисперсий для систем слу чайных величин весьма важными характеристиками являются кор реляционные моменты кх .Х[, определяемые формулой
kXjXl=.kJl = М [(Xj — Xj) (Xl — xt)} =
СО
= Jj (Xj — Xj) ( xt — x {) f (Xj, Xl) d X j d x t =
—■•Со
CO
= J • • • $ (Xj — Xj) (X , — Xi)f(xlt x2, . .., xn) dxxdx2 . . . dxn. (1.27)
— CO
При l = j формула (27) совпадает с (26) и, следовательно, kXjXj — = kjj = D [X .|. Совокупность корреляционных моментов называется корреляционной матрицей.
Закон распределения нормальной системы случайных величин полностью определяется их математическими ожиданиями и корре ляционной матрицей. Плотность вероятности нормальной системы
случайных величин Х х, |
Х 2, . |
. ., |
Х п имеет вид |
/ (Х\, х2, . . ., хп) -- |
|
|
|
1 |
1 |
|
A j l (Xj — X j ) ( x l — x l)\, (1.28) |
' (2ті)я/2 ѵ'д exp |
2Д |
^ |
j, i=i
где Д — определитель, элементами которого являются элементы корреляционной матрицы kjlt т. е.
* и |
*12 |
• |
• |
*1я |
|
д = *21 |
*22 |
• |
• |
*2я |
(1.29) |
*ЯІ |
*я2 |
• |
• |
*яя |
|
а A j t •— алгебраическое дополнение элемента kjt этого опре делителя.
Кроме функции распределения или плотности вероятности закон распределения случайных величин полностью определяется их характеристической функцией.
§ 1.1] |
Т Е О Р И Я В Е Р О Я Т Н О С Т Е Й |
15 |
Характеристической функцией Е (и) одной случайной вели чины X называется математическое ожидание е'иХ, т. е.
|
|
|
|
|
СО |
|
|
|
Е (и) = |
М [е,мХ] = |
^ |
eiuxf (X) dx, |
(1.30) |
||
где и — вещественная |
переменная. |
|
|
харак |
|||
Аналогичным образом для систем п случайных величин |
|||||||
теристическая |
функция Е (ult |
щ, . |
. |
.,1и„) будет |
|
||
Е{иѵ н2, |
|
і 2 |
ujXj |
|
|
|
|
е J - 1 |
|
|
|
|
|||
S |
г 2 |
uj xj |
|
|
|
|
|
- e i |
|
/( * i. |
Ха, |
|
xn) dxxdx2 .. . dxn. |
(1.31) |
Между плотностью вероятности и характеристической функцией существует однозначное соответствие. Так, для одной случайной величины имеем
СО |
|
$ e~iuxE(n)du. |
(1-32) |
— СО
Для системы п случайных величин
/ (^ij |
• >Хп)--- |
|
|
СО |
*' 2 uJxj |
|
|
|
|
|
і—1 Е (н 1; и ,, . ,,u^dux . . . dun. (1.33) |
|
— СО |
|
Для нормальных случайных величин характеристические функции имеют вид
Е[(и) = exp I — у а|ц2 + iuz^ , |
|
|
П |
П |
(1.34) |
— у 2 |
kjlujul + i S |
ujxj |
1 j, 1=1 |
i=1 |
|
Моменты случайной величины (если они существуют) могут быть получены путем дифференцирования характеристической функции Е (и) по формулам
1 dJ'E(u) |
Р , = J _ Ü [е~іихЕ (и)] ||(=„. |
(1.35) |
т ; |
||
ди* м-0 |
U дці |
|
16 |
О С Н О В Н Ы Е Ф О Р М У Л Ы |
[ГЛ . 1 |
Аналогичные формулы существуют и для системы случайных величин. Например, корреляционный момент к}^ может быть вы числен по формуле
К г |
д* |
[e-iujXj-iuixiE (р,ѵ щ |
яя)]|Мі=%=...=Ив=о. |
(1.36) |
|
dujdui |
|||||
|
|
|
Математическое ожидание и дисперсия линейной функции Z случайных величин X и Y,
Z = aX + bY + c, |
(1.37) |
|
вычисляются по формулам |
|
|
z — ах + by + |
с, |
(1.38) |
D [Z] = аЮ [X] + |
fe2D [Y] + 2abkxy. |
(1.39) |
Последние формулы имеют очевидное обобщение на линейные функции любого числа случайных величин.
Моменты случайной величины Z, являющейся полиномом более высокой степени, чем первой от заданной системы случайных вели чин, просто вычисляются только для нормальных величин. На пример, если Z = X 1 X 2, то, используя (34) для системы Х г, Х 2, получим
D [ZJ = D т |
D [Х2] + |
к\2+ |
**D [Х2] + |
x*D [XJ + |
2хЛ к12. |
|
(1.40) |
|||
Аналогичным образом для четырех нормальных |
величин |
имеем |
||||||||
М (AjAg-Xg-X^) = k lJc3i -f- к13к^ -)- кгік23 -f- k12x3xi -j- к13х^4 -|- |
|
|
||||||||
-j- k^x^x^ -j- k2 3 x ^ 4 |
к2 іхгх3 -f- к.мхгх 2 -|- x^x2 X^X4 . |
(1.41) |
||||||||
Плотность вероятности |
суммы случайных |
величин Z = Х |
Y |
|||||||
определяется формулой |
|
|
|
|
|
|
|
|
||
|
СО |
|
|
|
00 |
|
|
|
|
|
f A z) = \ |
f(z — y> |
y)dy — |
$ |
f (х, |
z — X) dx |
(1.42) |
||||
или для независимых |
величин |
|
|
|
|
|
|
|||
|
00 |
|
|
|
со |
|
|
|
|
|
/,(*) = |
J f A z — y)fy(y)dy = |
5 |
f A x)fy(z — x) dx- |
(143) |
||||||
|
— со |
|
|
— со |
|
|
|
|
|
Для случайной величины X, закон распределения которой отличается от нормального, часто с успехом может быть исполь зовано разложение плотности вероятности / (х) по производным от плотности вероятности нормального закона распределения