
книги из ГПНТБ / Скворцов, М. И. Теория и практика решения задач кораблевождения с учетом влияния систематических ошибок учебник
.pdfнасьева, которой весьма поучительным образом восполь зовался для обработки наблюдении, выполненных на при ливо-отливном течении, академик А. Н. Крылов [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