Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник

.pdf
Скачиваний:
16
Добавлен:
19.10.2023
Размер:
6.54 Mб
Скачать

насьева, которой весьма поучительным образом восполь­ зовался для обработки наблюдении, выполненных на при­ ливо-отливном течении, академик А. Н. Крылов [43]. Поз­ же в печати высказывались предложения об использова­ нии зависимостей между авансом и отношением крутяще­ го момента к квадрату числа оборотов движителей в ми­ нуту [76] или между авансом и отношением мощности на валу к кубу числа оборотов [85]. Но к .простым и удобным способам обработки наблюдений они не при­ вели.

Обобщив формулу В. И. Афанасьева, можно выска­ зать утверждение, что отношение авансов на двух разных режимах при штиле, неизменных осадке и состоянии под­ водной части корпуса корабля есть некоторая функция от­ ношения числа оборотов движителей:

Vj:Nj

/(NJ-.NJ,). (2.88)

Предположим, |Что при тех же заданных числах оборо­ тов движителей в минуту, как на j-м и }'-м режимах, вы­ полнялись наблюдения и на ряде предыдущих испытаний. Пусть VjC и Vj,c — средние значения полученных тогда скоростей хода, приведенных к стандартным условиям, причем именно эти скорости были введены в качестве при­ ближенных величин в формулу (2.63) для вычисления сво­ бодных членов уравнений поправок (2.60). Тогда, как сле­ дует из выражения (2.88), с той же точностью, как это было раньше, должно удовлетворяться равенство

' j'c

 

(2.89)

 

 

что приведет к еще одной системе уравнений

поправок

вида

 

 

aJWJ + a m W J ,

= vJJ

(2.90)

где

 

 

 

 

(2.91)

а7<Д

Vf9

(2.92)

• 80

Уменьшения числа искомых величин и повышения точ­

ности их оценок можно достичь, аппроксимируя

поправку

лага степенной функцией скорости

хода:

 

. AAJ = х + Vj~x2 + . . .

+VJ-1xk.

(2,93)

Это допустимо только в тех случаях, когда есть уверен­

ность, что ошибки

аппроксимации не превысят допустимых

пределов,

в

частности, когда

поправка, вводимая коррек­

тором лага,

с

достаточной

точностью

аппроксимируется

функцией

вида

(2.93), и имеет смысл, если число режи­

мов,

па

которых

проводятся

испытания,

превышает

чис­

ло к

искомых

величин, при котором удовлетворяется

пер­

вое

условие.

 

 

 

 

 

 

Составлению системы нормальных уравнений должно

предшествовать

оценивание корреляционной матрицы оши­

бок свободных членов уравнении поправок. Организацию этих вычислений облегчает упомянутое в § 1.1 общее пра­ вило [49]: если случайные величины X и Y являются сум­ мами независимых случайных величин и имеют общие

слагаемые, то корреляционный момент случайных

вели­

чин X и У равен сумме

дисперсий их

общих слагаемых.

В уравнениях поправок

(2.60) и (2.67),

относящихся

к од­

ному пробегу, общими являются ошибки измерения про­ дольной проекции абсолютной скорости и ошибки аппро­

ксимации течения,

в уравнениях (2.60)

и (2.77)—ошибки

измерения среднего за пробег числа

оборотов

движите­

лей,

в уравнениях

(2.60) и (2.67),

относящихся

к

разным

пробегам,— ошибки

аппроксимации

течения.

Корреляцион­

ная

матрица

ошибок уравнений

поправок

(2.90)

может

оцениваться

апостериорио по результатам предыдущих хо­

довых испытаний с последующей экстраполяцией на вновь проведенные испытания.

При испытании корабля на нескольких режимах общее число неизвестных в системе уравнений поправок (2.60), (2.67), (2.77) и (2.90) может достигать полутора-двух де­ сятков. Общий путь составления и решения системы нор­ мальных уравнении приводит к трудоемким вычислениям, практически осуществимым только на электронной вычис­ лительной машине. (Следует надеяться, что в недалеком будущем машинная обработка результатов наблюдений по­ лучит широкое распространение.) Если такая возможность отсутствует, вычисления хледует вести по упрощенной со­ кращенной схеме. Для этого надо разбить систему урав-

a i

нений поправок на ряд частично независимых подсистем и решать их способом последовательных приближений, на­ пример, так, как показано в § 3.10. Там же приведены рас­ четные формулы. Сделаем несколько дополнительных за­ мечаний.

Ведя обработку наблюдений, приходится вычислять средние за время пробега значения продольной проекции скорости абсолютного перемещения корабля, числа оборо­

тов

гребных валов в минуту и скорости нарастания отсче­

та

лага. Если требуется точность определения I поправки

лага 0,1—0,2%, эти вычисления должны выполняться с точностью до четвертой-пятой значащей цифры. Логариф­ мическая линейка такой точности не обеспечивает, а бо­ лее совершенные средства вычислений имеются не на всех кораблях. Требуемой точности можно достичь, пользуясь несложными вспомогательными таблицами, помещенными в § 3.9. Там же приведен пример пользования ими.

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

лого хода через

поправки ДУ , и ДУ..

для

двух

со­

седних

режимов.

Как следует

из выражения (2.89),

 

 

A

V

Zl!LAV

1IL.AV

 

 

( 2 .94)

 

 

 

J С

J с

 

 

 

Использование этого соотношения позволит на режиме

самого

малого

хода

обойтись

только

одним

пробегом.

Чтобы

при этом

не

приходилось экстраполировать

про­

дольную проекцию скорости течения, наблюдения следует организовать следующим образом. Сначала на j'-м ре­ жиме надо выполнить два пробега, потом на /-м режиме самого малого хода—один пробег, направление которого противоположно лоследнем-y пробегу на предыдущем режи­ ме, затем на /"-м режиме выполнить два пробега. Направ­ ление первого из них должно быть противоположно на­ правлению пробега на режиме самого малого хода.

Если теперь, воспользовавшись формулой (2.94), выра­ зить поправку AVj в уравнении, поправок (2.60), состав­ ленном по результатам наблюдений на единственном про­ беге режима самого малого хода, сначала через поправ­ ку ДУ., для предыдущего режима, а затем через поправ?

82

ку Д1Л„

для последующего режима и присоединить полу­

ченные два уравнения поправок к

уравнениям, составлен­

ным

по' наблюдениям

на /'-м

и

j"-u

режимах,

получим

три

уравнения

с неизвестней

 

'.,

и три уравнения

с

неизвестной

AV.„.

Их решение,

осуществляемое

по той

же формуле

(2.71),

что и при трех пробегах на одном ре­

жиме, позволит

найти

оценки

искомых величин

Д1Л, и

Д1Л„,

а также

продольной проекции

скорости течения

на

каждом

пробеге, -в том

числе

и на единственном

пробеге

/-го режима. Подставив

оценки

Д1Л, и Д1Л, - в

выраже­

ние (2.94), можно получить искомую поправку ДК,- в пер­ вом приближении.

Оценивание точности результатов испытаний. Вначале рассмотрим эту задачу применительно к традиционным Способам обработки наблюдений [42], [18], [59, гл. 1], [1], когда в принятых нами обозначениях оценки искомых ве­ личин вычисляются по формулам (2.70), (2.71). Из выра­ жений (2.70) и (2.59) следует, что при двух пробегах на режиме истинная ошибка определения скорости хода Vj равна

 

 

-{dx +

d2)kw +

 

{v,x-vT,)],

 

 

 

(2.95)

где

A v x ,

ДК 2 —истинные

ошибки

измерения

средних

 

 

за

пробег

 

значений -продольной про­

 

 

екции

абсолютной

скорости на

первом

 

 

и

втором

 

пробегах;

 

 

 

 

A/vi>

Адго — истинные

ошибки

определения

сред­

 

 

них чисел оборотов движителей в ми­

 

 

нуту

на этих пробегах;

 

 

 

 

VTU

VT2 — истинные

значения

средних

величин

 

 

продольной

проекции скорости

тече­

 

 

ния на перв_ом и втором пробегах.

 

Если длина пробега известна с пренебрежимо малой

ошибкой, а секундомеры перед испытаниями

тщательно

проверены (имеющие

большой

 

и непостоянный

ход

отбра­

ковываются),

то при

должной

организации

наблюдений

систематические ошибки измерений величин V,-, Ni и Jli будут невелики. Возможности же уменьшения случайных ошибок измерений ограничены числом наблюдателей и чи­ слом секундомеров, которыми они могут быть снабжены.

83

Поэтому будем считать, что в величинах V,-, Nt и Л ; слу­ чайные ошибки преобладают над систематическими. Что же касается ветра и течения, то их влияние является ис­ точником преимущественно систематических ошибок. Если учесть это обстоятельство, то, как видно из выражения (2.96), при двух пробегах на режиме дисперсию ошибки оценки скорости хода Vj следует считать равной

D <у})

= - 1 - D (V;) + - f fD

 

+ -J- № + d2f

k\,

+

 

 

+ - f c°i ^

 

 

( 2 - 9 6 )

где D{ V[),

D (N,) — дисперсии

ошибок измерения

средних

 

 

за пробег значений продольной проек­

 

 

ции абсолютной скорости

и числа обо­

 

 

ротов

движителей;

 

 

 

 

 

 

—средняя

 

квадратическая

 

величина

 

 

первой

производной

по

времени от

 

 

средней за время пробега величины

 

 

продольной

 

проекции

скорости

тече­

 

 

ния, мили • ч~2 ;

 

 

 

 

 

 

т—-промежуток

времени

между

средними

 

 

моментами

 

пробегов,

ч.

 

 

 

Второй

член

формулы

(2.96)

выражает влияние

оши­

бок определения

среднего за пробег числа оборотов движи­

телей. На

необходимость уделять "определению числа

обо­

ротов не меньше внимания, чем измерению скоростей, не­

однократно указывал академик

А. Н. Крылов [43], [44].

Но со временем его напоминания

оказались забытыми, и

в работах, посвященных анализу точности испытаний хо­ довых качеств кораблей [67], [86], [74], [59, § 32], [1], эти ошибки не рассматриваются.

Четвертый член формулы (2.96) выражает влияние из­ менчивости течения. Если рассматривать величину у^как стационарную случайную функцию времени и обозначить символом D(Vj) ее дисперсию, символом /"(-с)—нормиро­ ванную корреляционную функцию, это влияние может быть выражено формулой, .предложенной В. С. Болдыре­ вым и Б. Н. Беляевым [8]:

D4 (\/y ) = 4 - D ( z / T ) [ l - r ( , ) ] .

(2.97)

84

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

его первым

членом.

 

 

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

дисперсии ошибки в оценке скорости хода V} при трех про­

бегах на режиме:

 

 

 

D<yJ)

=

±D<yi)

+ ±-fD(N,)

+

+

±

(с!л

+ 2d, +

с!лУ k\, + ± (рц

(2.98)

где о- —средняя квадратическая величина второй произ­

водной по времени от средней за

время

пробега

величи­

ны

продольной

проекции

 

скорости течения, милн - ч - 3 .

 

Последний член формулы (2.98), выражающий

влия­

ние

изменчивости течения,

может

быть

представлен

так­

же

в виде,

предложенном

В. С. Болдыревым и Б. Н. Бе­

ляевым:

 

 

 

 

 

°

 

 

 

 

D,(V.)

= 4- [3 -

4г (•:) +

г (2х)1 D ( О .

(2.99)

 

Приняв

приближенно

 

Jli—Vj,

можно

получить

выра­

жение

и для априорного

оценивания средней квадратиче-

CKOJ'

величины

ошибки

определения поправки лага:'

 

 

 

 

с(Алу ) =

у, J )

100%.

-

(2.100)

Обработка наблюдений, основанная на принципах, из­ ложенных в настоящем параграфе, позволяет производить не только априорное, но и апостериорное оценивание точ­ ности оценок искомых величин. Не останавливаясь на ап­ риорном оценивании, поскольку оно существенно не от­ личается от только что рассмотренного для традиционных способов, отметим, что испытания на двух режимах (по два пробега на каждом) позволяют определить и учесть скорость изменения продольной проекции скорости тече­ ния по времени. Если испытания проводятся при полном штиле, то для достижения тех же точностей оценок иско­

мых

величин,

что и традиционными способами,

достаточ­

но на каждом

режиме делать по два пробега вместо трех,

а на

режиме

самого малого хода ограничиться всего од­

ним

пробегом.

При ветре заметного ухудшения

точности

85

определения аванса, как при обработке наблюдений тра­ диционными способами, происходить не будет.

При апостериорном оценивании точности оценок иско­ мых величин следует всегда исходить из некоторой гипо­ тезы о соотношении точностей отдельных измерений. Огра­ ничимся простейшим случаем, когда нет оснований по­ лагать, что некоторые однотипные измерения на какомлибо режиме были относительно менее точными, нежели на другом. Тогда все уравнения поправок вида. (2.60) и

(2.67) можно считать равновесными. Обозначим

симво­

лом а„ среднюю квадратическую величину ошибки, с

какой

из этих уравнений определяется продольная проекция от­ носительной скорости. Пример ее оценивания показан в

§3.10.

Оценивание точностей оценок аванса и поправки лага требует более сложных вычислений. Но приближенно мож­ но считать, что средняя квадратическая величина ошибки определения скорости Vj, соответствующей при штиле за­ данному *числу оборотов движителей, есть величина при­ мерно того же порядка, что и а0. О порядке ошибки опре­ деления поправки лага дает представление формула (2.100).

Определение расхода топлива и дальности плавания корабля составляет одну из важных задач ходовых испы­ таний. Для определения расхода топлива в единицу вре­ мени обычно рекомендуется применять мерные баки спе­ циальной конструкции. Относительная точность этих изме­ рений мала. Если мерные баки отсутствуют и применяются расходомеры, она становится еще более низкой. Результаты, полученные при кратковременных испытаниях, приходится экстраполировать на все дальнейшее плавание корабля, что ведет к дополнительным ошибкам. Соотношения рас­ ходов топлива на разные нужды, характерные для повсе­ дневной эксплуатации корабля, при испытаниях модели­ руются грубо. Не учитывается влияние таких факторов, как обрастание подводной части корабля, дополнительные потери топлива, возникающие вследствие износа механиз­ мов, и т. д. Все это приводит к значительным системати­ ческим отклонениям фактических расходов топлива ' при эксплуатации корабля от предвычисляемых на основании результатов испытаний. Основные статистические харак­ теристики этих отклонений могут быть оценены примене­

нием

метода последовательного уточнения оценок иско­

мых

величин.

86

Заметим, что суммарный расход топлива за большой промежуток времени (например, от одного пополнения за­ пасов топлива до другого) может определяться со срав­ нительно небольшими относительными погрешностями. Нас обычно интересуют удельные расходы топлива как - на милю пройденного расстояния (при длительных плава­ ниях на постоянном режиме работы энергетической уста­

новки),

так и за

единицу времени

(при плавании

малыми

ходами,

на стоянке

корабля

и т. д . ) . Обозначим

симво­

лом

а

значение

параметра

(пройденного

кораблем рас­

стояния

или времени), на единицу

которого определяется

расход

топлива,

символом

— удельный

расход

топлива

(на

единицу этого параметра)

на /-м режиме. Оценки иско­

мых

величин будем

отыскивать

в виде

 

 

 

 

 

 

lj = tjc +

Xj,

 

 

(2.101)

где

SyC —оценка

искомой величины

по результатам пре­

 

 

дыдущих

испытаний;

 

 

 

 

 

Xj

— искомая

поправка.

 

 

 

 

Тогда зависимость между наблюденным за достаточно большой период времени (от одного пополнения запасов топлива до другого) расходом топлива и оценками иско­ мых величин может быть выражена в виде уравнения по­ правок

апхх + ... + auXj + . . .

+ aimxm -1. = vp (2.102)

где

(2.103)

UJ. — суммарное время, в течение которого энергетическая установка эксплуатиро­ валась в /-м режиме (если £j выражает удельный расход топлива на единицу времени), или пройденное на /-м режиме расстояние (если |у — удельный расход топлива на милю пути);

—фактический общий расход топлива по результатам (-го наблюдения;

а?

т

счнслимый расход топлива на т режи­

мах,

для которых оценки Xj включены в

-число

искомых

величин (вычисляется ис­

ходя из найденных на предыдущих ис­

пытаниях оценок удельных расходов то­

плива на этих

режимах);

Rir — счислимый или

непосредственно изме­

ренный расход топлива на всех осталь­

ных

режимах.

 

Составление и решение системы нормальных уравне­ ний, а также апостериорное оценивание точностей оценок искомых величин могут вестись в соответствии с общими основами способа последовательного уточнения оценок ис­ комых величин: Если вычисления ведутся вручную, то, чтобы избежать излишней их трудоемкости, приходится ограничивать число искомых величин тремя-четырьмя. Это понижает точность их оценок. Поэтому, как и при испыта­ ниях на мерной линии, целесообразна организация цен­ трализованной обработки результатов наблюдений в учре­ ждениях, располагающих современной вычислительной тех­ никой.

§ 2.4. О Ц Е Н И В А Н И Е ТОЧНОСТИ СЧИСЛЕНИЯ

Пусть в течение промежутка времени t велось непре­ рывное счисление пути корабля без корректировки обсер­ вациями. Будем рассматривать вектор r(t) истинной ошиб­ ки счислимого места как случайную функцию времени. Истинной ошибкой счисления, накопившейся за время /, назовем вектор

 

с(0 = г ( 0 - г ( 0 ) ,

 

(2.104)

где /'(0)—истинная ошибка счислимого

места в началь­

ный момент

времени.

 

 

Если счисление ведется по показаниям компаса

и от­

носительного

лага, то ошибки счисления

являются

след­

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

88

нне, не прерываемое обсервациями, можно считать прак­

тически стационарным

[7], [69], [57, § 25], [39, § 13].

 

В

§ 2.1 уже упоминалось, что истинная ошибка

изме­

рения

курса корабля

в любой момент времени

может

быть

представлена как. сумма трех составляющих: ошибки

определения поправки компаса, изменения среднего зна­ чения истинной ошибки курсоуказания за промежуток вре­ мени, прошедший после определения поправки компаса, и случайного отклонения мгновенного значения ошибки кур­ соуказания от ее среднего значения. Первые две ошибки в течение всего времени плавания по счислению остаются постоянными и проявляют свойства систематических оши­ бок. Подобным образом и ошибка поправки лага может быть представлена как сумма трех составляющих: ошибки (Определения поправки лага, изменения среднего значения поправки лага, происшедшего после определения, и слу­ чайного отклонения поправки лага в рассматриваемый мо­ мент времени от ее среднего значения.

Рассмотрим влияние постоянных' составляющих этих ошибок. Пусть поправка компаса известна с постоянной систематической ошибкой Дк ь поправка лага — с постоян­ ной ошибкой ДлЬ

Пусть

Ур (0—истинное значение вектора скорости относитель­ ного перемещения корабля в момент времени г;

V (t)—счислимое

значение вектора относительной ско­

рости;

 

у — единичный вектор, направление которого полу­ чено вращением вектора V (t) на 90° по ча­ совой стрелке.

Если отвлечься от остальных ошибок измерений, то

Fp c (0 = vp (0 + т ^ - й р с ( 0 + ^ v(C (о/ (2.Ю5)

Истинная ошибка счисления, накопившаяся за время г, будет равна

СЛ0

=

l[v^t)~v9{t)]dt^

 

 

о

 

 

=

H

^ + W 5 P c ? c . •

(2.106)

89

Соседние файлы в папке книги из ГПНТБ