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

книги из ГПНТБ / Шаталов, В. А. Применение ЭВМ в системе управления космическим аппаратом

.pdf
Скачиваний:
34
Добавлен:
20.10.2023
Размер:
8.16 Mб
Скачать

Продолжение

2

 

 

 

 

 

 

т

о с н о в н ы е ф азы

п о л ета К К

О сновны е за д а ч и , р е ш а е м ы е

н а б о р ту

К К

не­

 

„А поллон -

при п о м о щ и Б Ц В М

 

2

 

 

 

 

 

 

20

Отделение командного

отсе-

 

 

ка от служебного отсека

отсе­ Определение положения и скорости в про­

21

Вход

командного

ка

с экипажем

в атмосферу странстве

вектора

скоро­

Земли

 

 

Управление ориентацией

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

К первой группе следует отнести операции сложения, вычита­ ния, сравнения, различных передач. Общее число операций типа сложения определяется соотношением

 

 

N сл = «сл + kxnB-j- k2ncр -f-ksnneV,

(7. 3)

где Пси,

«в, Пер, «пер —число операций

сложения, вычитания

 

 

сравнения, передачи команд;

 

 

 

ku k2, k3— коэффициенты пропорциональности, учи­

 

 

тывающие соотношение между временем

 

 

выполнения стандартной операции сло­

 

 

жения и операциями другого типа.

 

Ко второй группе можно отнести операции типа умножения,

деления,

сдвигов. В этом случае общее

число операций

типа

умножения запишется в виде

 

 

 

 

Му—% + &4«д~Ь ^5«СД)

 

(7-4)

где яу, Яд,

яСд — число операций умножения, деления, сдвига;

 

ki,

k5— коэффициенты пропорциональности.

 

Таким образом, можно считать, что ЦВМ способна выполнить Nсл сложений или Му умножений в единицу времени. Номиналь­ ное быстродействие позволяет с определенной степенью достовер­ ности осуществить сравнение ЦВМ по скорости выполнения вычислительных работ.

Основным условием, из которого следует исходить при опре­ делении номинального быстродействия БЦВМ, является выпол­ нение неравенства

 

 

 

 

 

(7.5)

где

At — время, отводимое на решение требуемой задачи;

 

 

t(p ■— фактическое время решения задачи.

 

 

Величина определяется выражением

 

 

 

 

Ц = ^ в(1 ^l) ! ^гр I ^к.ф

1 ^в

! ^выв ( ^ — 8г).

(7- 6)

где

tBB, Цгв — время, отводимое

на

ввод

исходной инфор­

 

мации и вывод результатов

решения задачи;

170

si, 62 — коэффициенты совмещения

ввода

и

вывода

информации с процессом решения;

считывания

trp— время

групповых обменов

для

программы, контрольных исходных

данных

и программ математического обеспечения из

ВЗУ в ОЗУ;

решения

задачи

tB— время

непосредственного

вычислителем;

 

 

 

/к. ф— время,

отводимое для контроля функциониро­

вания ЦВМ, например при помощи контроль­ ных тестов.

Рассмотрим кратко возможность определения составляющих выражений (7.6). Значение /вв определяется количеством сооб­ щений, поступаемых на вход ЦВМ, и временем передачи одного /-го сообщения. Величина /вв может быть найдена из соотно­ шения

4в= 2 т ^г’

(7-7)

г-1

 

где ti — njqf— время передачи одного /-го сообщения; п — число разрядов в коде сообщения; f — частота передачи разрядов кода; q — число, одновременно передаваемых разрядов; /пг — количе­ ство однотипных сообщений.

Выбор коэффициента si может быть осуществлен следующим образом [5].

Если для работы устройства ввода требуется более 0,15— 0,20 времени решения задачи в вычислителе, то следует выби­ рать структуру ЦВМ с совмещением времени работы устройства ввода и вычислителя. Коэффициент совмещения может изменять­

ся от нуля (при /с = 0 ) до максимального значения,

равного

п гаах= — ,

(7.8)

когда t c = tBв.

Аналогичные рассуждения могут быть проведены при выборе составляющих /выв и s2.

Для оценки времени групповых обменов /гр необходимо задаться структурой и типом внешних запоминающих устройств (ВЗУ). Общее время для выполнения групповых обменов инфор­ мацией между оперативным запоминающим устройством (ОЗУ)

и ВЗУ может быть определено соотношением

 

ДР= V ((3 4 1;. + Af2/)T;.)(1- P ;.),

(7.9)

j

 

где Mij — общее количество чисел и команд,

пересылаемых

из /-го ВЗУ в ОЗУ; М2;- — количество чисел и команд, пересылае-

1/1

— -\

Рис. 7. 1. Временная диаграмма выполнения

команд

 

jon

З У

 

'К.ЗУ

'Ч. ЗУ

мых из ОЗУ в /-ое ВЗУ; Xj — время, приходящееся на пересылку

одного числа из ОЗУ в ВЗУ или наоборот;

■— коэффициент сов­

мещения встречных групповых передач

информации из ОЗУ

в /-ое ВЗУ и из /-го ВЗУ в ОЗУ.

 

Часто можно вычислительный процесс строить так, что выпол­

няется условие

 

М ц = М х = М },

(7. Ю)

тогда время для групповых передач

 

^ р = 2 уИЛ '2 (1 - р >)-

(7Л1)

)

 

Время, отводимое для контроля функционирования, tK.ф существенным образом зависит от применяемого метода контроля правильности функционирования вычислителя.

Время непосредственного решения задачи tB вычислителем определяется последовательностью выполнения команд в вычис­ лителе. Для его оценки необходим анализ временных диаграмм выполнения команд. Рассмотрим принцип разделения времени работы ЗУ и АУ.

Время выполнения команд на ЦВМ с программным управле­ нием слагается из составляющих:

времени выбора из ЗУ команды /к.зу;

времени выбора из ЗУ числа или чисел t4- зу в зависимости от количества адресов в команде k\

времени выполнения операции ton.

На рис. 7.1 для примера приведена временная диаграмма вы­ полнения команды для ЦВМ с одноадресной системой кодиро­ вания. Временная диаграмма выполнения команды состоит из трех частей: времени получения из ЗУ команды — /к. зу, времени считывания или записи числа — 14.зу и времени выполнения опе­ рации.

Таким образом, для определения величины tcР могут быть оценены все составляющие.

При невыполнении условия (7. 5) необходимо провести тща­ тельный анализ составляющих соотношения (7. 6) и внести изме­ нения в структурную схему ЦВМ, которые установили бы требуе­ мый баланс времени.

172

Требуемое быстродействие ЦВМ или непосредственно вычис­ лителя определится соотношением

V Ц В М

v

Ц В М

(7. 12)

СЛ

____ м

у

У

At

 

At

 

 

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

После этого составляется сводная диаграмма фаз и потреб­ ностей каждой из фаз по памяти и быстродействию. Наиболее нагруженные фазы и определят (с учетом необходимого резерва) требования к БЦВМ. В табл. 7.2 представлены суммарные потребности по памяти и быстродействию для фаз полета к Марсу [15].

Т а б л и ц а 7. 2

азыф

Н мо ре

 

П о к а з а т е л ь П о т р е б н а я

З а д а ч а ф азы

п р о д о л ж и ­

п а м я т ь ,

те л ьн о сти

т ы с я ч и

 

 

ф а зы , ч

сл о в

С к о р о с ть в ы ч и сл и ­ т е л я , ты с .

о п ер /с

 

!

 

 

 

 

 

 

 

1

Вывод на орбиту Земли

 

 

0 , 2

2,5

1 0

2

Полет по орбите вокруг Земли

80

8 , 6

п о

3

Вывод

на

траекторию

перелета

0 , 2

4,5

45

 

к Марсу

 

 

 

 

 

 

60

4

Пассивный полет к Марсу

 

119,8

8 , 8

5

Коррекция траектории

 

 

0 , 2

9

140

6

Управление

ориентацией

на марше­

2759,8

9,2

60

 

вом полете

 

 

 

 

 

 

140

7

Коррекция при достижении Марса

0 , 2

9,5

8

Переход на орбиту Марса

 

 

0 , 6

7,6

55

9

Полет по орбите Марса

 

 

974,3

22,3

365

1 0

Вывод на орбиту к Земле

 

 

0 , 2

8,5

70

11

Пассивный полет к Земле

 

 

119,8

8 , 8

60

1 2

Коррекция траектории

 

 

0 , 2

9

140

13

Управление

ориентацией

на

маршевом

6119,8

9

60

 

полете

 

 

 

 

 

 

 

14

Коррекция

при достижении

Земли

0 , 2

6 , 2

2 0 0

15

Управляемый вход в плотные слои

0 , 2

9

60

 

атмосферы

 

 

 

 

 

 

 

173

7,2. ОЦЕНКА ТРЕБУЕМОЙ ТОЧНОСТИ ВЫЧИСЛЕНИЙ И ВЫБОР РАЗРЯДНОСТИ ЦВМ

Одной из важнейших характеристик БЦВМ является длина кода и разрядность машинного слова.

Разрядность определяется следующими факторами: •— точностью решаемых на борту задач;

— характеристиками входной информации, поступающей для хранения и обработки;

требуемой разрядностью команд;

требованиями по быстродействию и массе оборудования.

Рассмотрим задачу выбора разрядности. Предположим, что при помощи БЦВМ реализуется управляющая функция

ll(t)=f(Xi, х2,..., хп)

(7. 131

от поступающих на ее вход аргументов хи х2, .. ., хп. Общая ошиб­ ка на выходе БЦВМ определяется неточностью математического описания (например, линейной аппроксимацией нелинейной за­ висимости) физического процесса, реализуемого при помощи алгоритма БЦВМ, трансформированной ошибкой, зависящей от вида функции f(xit х2, .. хп) и ошибок аргументов на входе БЦВМ и инструментальной ошибки цифровой машины [3].

Перечисленные ошибки независимы, следовательно, для БЦВМ можно составить баланс ошибок, описываемый выраже­ нием

(7. 14)

где оу — среднее квадратичное значение ошибки на выходе БЦВМ; сгы — средние квадратичные ошибки математических аппроксимаций; от — среднее квадратичное значение трансфор­ мированной ошибки; (ти — среднее квадратичное значение инст­ рументальной ошибки ЦВМ.

Проанализируем составляющие соотношения (7. 14). Ошибка математических аппроксимаций обусловлена тем, что управляю­

щие сигналы в ЦВМ формируются не в

точном соответствии

с желаемой управляющей функцией f(xu

х2, . .., хп), а в виде

близкого к этой функции аналога f*(xu х2, . . ., хп), алгоритми­ чески реализуемого в вычислителе. Трансформированная ошибка зависит от неточности входных величин и вида преобразующей функции f(x 1, х2, . , ., х„). Значение трансформированной входной ошибки определяется соотношениями

2

д/ ( х и

х 2,..., х п) d f

( х \ ,

х ъ.,, .х п)

1/2

 

д

х }

д х л

д х , г

(7. 15)

 

U ь

 

 

 

 

174

если величины X j и x h коррелированы, и

 

d f ( х и х ъ . . х п)

а 2

1/2

а

(7. 16)

т

д х j

ВХ

 

 

 

 

в случае независимости этих величин.

В выражениях (7.15, 7. 16) — j — среднее квадратичное значение ошибок аргумента х, на входе БЦВМ, rik — коэффи­ циент корреляции величин Xj, Хъ.. Величина aBxj в свою очередь определяется соотношением

а) + а1г

(7- 17)

где в, — средняя квадратичная ошибка

аргумента на входе

/-го датчика; Ohj — средняя квадратичная

ошибка квантования

преобразователя. Ошибки квантования преобразователя связаны с его разрядностью. Пусть каждой величине Xj соответствует свой преобразователь непрерывной величины в дискретную и уровень квантования /-го преобразователя равен Дх^. Ошибки

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

равновероятному закону, а их предель-

ные значения равны ±

. Следовательно,

средние квадра­

тичные ошибки квантования

преобразователя

будут

 

3.. =

- ^ % - .

(7.18)

Разрядность же Rj преобразователя можно определить из соот­ ношения

• * 7 m a x х 3m in *2,^3 j

Точность работы непосредственно ЦВМ (инструментальная ошибка) сги определяется совокупностью таких характеристик, как разрядность преобразователей, разрядность запоминающих устройств и операционного устройства и количеством разрядов, необходимых для представления выходных величин с требуемой точностью [3].

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

Рзу

(7.20)

Разрядность операционного (арифметического) устройства может быть записана как

Ron= R 4 + S,

(7.21)

где R4 — разрядность чисел, вводимых в ЦВМ;

175

S — количество дополнительных разрядов для компенсации ошибок округления.

В свою очередь разрядность чисел, вводимых в машину, опре­ делим разрядностью преобразователя. Значение 5 можно подсчи­

тать, исходя из следующих рассуждений.

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

~mfx — ед. младшего разряда,

юкр

а ее среднее квадратичное значение равно

окр"

ед. младшего разряда.

 

2 / 3

При независимых однократных округлениях средняя квадра­ тичная инструментальная ошибка сти для последовательности «-вычислений равна

+

32 _I

, . _ j —02 = 3

'■'окр У п

(7. 22)

■V'окр! I

окрokd2 I

1 окр л

 

 

или

ед. младшего разряда.

Для компенсации накопления ошибок округления необходимо ввести дополнительное число разрядов

•S = log2(oH+ 1)

(7.

23)

“/Ц .2 +’ ) '

(7.

24)

 

Таким образом, может быть определена общая разрядность операционного устройства. Разрядность выходных устройств RBЫх обычно выбирают равной разрядности входных преобразователей

изапоминающих устройств. В результате, используя уравнение баланса ошибок и зная общие требования, предъявляемые к ау,

ивычислив значения ом и сгт, можно определить аи, а следова­ тельно, и минимально требуемую разрядность.

Анализ точности решения различных задач на борту показал [1], что можно выделить две существенно различающиеся группы задач, требующие различной точности и соответственно разряд­ ности ЦВМ. При решении задач типа навигации требуемая точ­ ность колеблется в пределах от 27 до 32 разрядов. При решении

176

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

еm a x

(7.25)

где п — число разрядов, отведенное под мантиссу числа.

Если выбрать разрядность БЦВМ, учитывая самые жесткие требования, возникающие при решении задач навигации, то необ­ ходимо выбрать машинное слово, имеющее от 27 до 32 разрядов. Однако в этом случае увеличивается масса оборудования и уменьшается быстродействие за счет увеличения времени про­ бега сигнала переноса в арифметическом устройстве БЦВМ. Одним из путей разрешения данного противоречия является выбор машинного слова, удовлетворяющего основным требова­ ниям, и использование для решения задач, требующих повышен­ ную точность, специальных программ вычислений с повышенной точностью, отводя при этом для хранения каждого числа не одну, а несколько ячеек.

Рассмотрим принцип построения алгоритма сложения чисел с повышенной точностью в ЦВМ с фиксированной запятой. Пусть мы имеем два числа А и В, размещенные каждое в двух ячейках памяти, и А с и Вс и Вы— части чисел А н В, содержащие старшие и младшие разряды соответственно. Пусть результат А' также должен быть получен с двойной точностью и состоит из Хс и Ам, тогда для сложения или вычитания

(7. 26)

где р — единица переноса при вычислении Хж.

При выполнении сложения и вычитания могут быть исполь­ зованы два способа представления результата:

знак числа единый для частей числа, хранящихся в раз­ ных ячейках;

каждая часть числа может иметь свои знаки.

Пример. Рассмотрим двузначные десятичные числа. Пусть А

= 4-43; В —

= +28; А с = 4; Л м= 3 ; В с — 2; В „ = 8. Тогда Хм = 1; р — 1; Х с =

7.

При сложении и вычитании чисел с повышенной точностью каждая со­ ставляющая может получить свой знак. Будем обозначать числа следующей записью: А = + 4 3 = [ + 4 , +3]. Тогда при Л = +43 и В = —28

[+4, +3]+[—2, —8 ]= [+ 2 , —5 ]= [+ 1, +5].

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

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

7

994

177

[+4, + 3 ]= [+ 5 , -7 ];

[—2, —8]== [—3, +2].

Однако операции, проводимые над числами, дают в этом случае одинаковый

результат:

[+ 4, +3]+[—2, —8] = [+2, -5 ]; [+5, —7]+[—2, —8]= [ + 2, -5 ]; [+4, +3]+[—3, +2] = [+1, +5]; [+5, —7]+[—3, +2] = [+ 2, -5 ].

. Рассмотренные примеры показывают, что когда появляется перенос, то он вместе со своим знаком переходит в сумму более старшей части. Часть, в которой сформировано переполнение, сохраняет при этом свой знак. Таким образом

р = 0 или р = ± 2~п.

Примером использования алгоритма вычислений с повышен­ ной точностью может служить ЦВМ, размещаемая на КА «Апол­ лон», в которой применяется алгоритм вычислений с двойной точ­ ностью (сложение, вычитание, умножение и деление) и алгоритм сложения и вычитания с утроенной точностью. При этом данные достигнутого времени выполнения операций с повышенной точ­ ностью приведены в табл. 7. 3.

Т а б л и ц а

7.3

 

 

 

 

 

 

 

У р о в ен ь т о ч н о с т и

В р е м я

вы п олн ения

 

 

 

о п ер а ц и и

 

О б ы ч н ая т о ч н о с т ь

У д в о е н н а я то ч н о с т ь

 

 

 

 

Время

сложения,

мкс

23,4

35,1

Время

умножения,

мкс

46,8

57,5

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

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

Так, в ЦВМ «Аполлон» выбрано 15-разрядное машинное слово. Такая разрядность не обеспечивает решение поставлен­ ной задачи, без специальных мер, так как позволяет кодировать лишь 8 команд и производить адресацию к 4096 ячейкам ЗУ. Число команд и адресов в этом случае расширяется при исполь­ зовании специальных регистров.

178

Таким образом, машинное слово в 18—20 разрядов и способ­ ность машины обрабатывать данные с повышенной точностью, уже позволяют удовлетворить требования, предъявляемые кБЦВМ.

7.3. АНАЛИЗ ХАРАКТЕРИСТИК БЦВМ ПРИ ПОМОЩИ МОДЕЛИ

Определение отдельных требований к БЦВМ при помощи ана­ литических выражений может быть использовано лишь в про­ стейших случаях, • например, если вычислительные бортовые средства состоят из одного процессора. Кроме того, при помощи аналитических выражений крайне трудно с достаточной степенью точности учесть аппаратурные системы вычислительных средств, структуру алгоритмической системы, надежность средств обра­ ботки и т. д.

Указанные недостатки устраняются при использовании для оценки параметров вычислительной системы метода моделиро­

вания [2, 9].

Модель, состав которой приведен на рис. 7. 2, позволяет уста­ новить степень влияния различных характеристик на показатели качества системы и определить, как справляется система обра­ ботки данных с обслуживанием входящего потока информации при различном составе средств, различных схемах обработки и воз­ можных критических и аварийных ситуациях.

Общая модель состоит из частных моделей. Рассмотрим осо­ бенности их построения и основные исходные данные, которые необходимы для создания модели.

1. Модель источников позволяет сформировать входной поток данных, поступающих в вычислительную систему. Под источни­ ками применительно к БЦВМ понимаются различного типа дат­ чики, пульты и другие устройства, с которых поступают данные, а также сигналы, вырабатываемые в самой системе, на решение отдельных задач в соответствии с программой полета. Модель источников, формируя время выдачи сообщений, должна указы­ вать номер задачи, которая обслуживает данный тип сообщений. Необходимо отметить, что точность знания характеристик вход­ ного потока данных в основном определяет степень соответствия модели реальной системе обработки данных на борту КА.

2.Модель программной сети является основой общей модели

иее главной частью.

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

7

179

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