Нагрузочный расчет

Цель нагрузочного расчета - формирование рабочей нагрузки АСНИ во времени между устройствами системы. При этом необходимо решить следующие задачи:

* выбрать режим сбора и первичной обработки данных;

* рассчитать параметры временной диаграммы работы АСНИ по критерию минимума непроизводительных затрат при ограничении на скорость сбора данных;

* рассчитать требуемый объем оперативной памяти;

* выбрать ЭВМ, состав программного обеспечения и интерфейс (предварительно) из заданного набора альтернативных вариантов.

Снятие показания с датчика и его первичную обработку будем считать прикладной задачей. Для опроса i-го датчика каждый раз запускается i-ая прикладная задача. В начале расчета подсчитывается С₀ - суммарная частота запуска прикладных задач по следующей формуле:

Здесь m - число датчиков в системе, f_oi - частота опроса i-го датчика (получены на этапе информационного расчета).

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

r_n(C)=1-r_д(C)

Здесь r_д(C) - приведенные затраты процессорного времени на диспетчеризацию.

r_д(С)=1-r_п(С)

0£r_n(C)£1

Конкретный вид ПНХ определяется выбранной операционной системой (ОС). Операционная система выбирается в соответствии с запросами системы и с учетом критерия минимизации стоимости системы. Возможные ОС представлены в П1 пункт 6.2.3 таблица 6. Наиболее быстродействующей является 3-я ОС.

 

Рис. 4

Каждому рабочему режиму АСНИ на плоскости нагрузочной характеристики соответствует рабочая точка (РТ) с координатами: C_рт и r_рт.

r_рт=С_рт(t_усд+t_инт+t_по)

 

t_усд - задержка вносимая УСД.

Здесь С_усд - максимальная скорость преобразования, выбирается из [1] пункт 6.2.3 таблица 2.

t_инт - задержка вносимая интерфейсом при обмене данными между ЭВМ и УСД, выбирается из П1 пункт 6.2.4 таблица 4.

t_по - затраты процессорного времени не связанные с обменом (подготовительные операции, первичная обработка), выбирается из П1 пункт 6.2.4 таблица 5.

При выборе ОС значения t₁ и t₂ в П1 таблица 6, а также t_по даны для 3-ей ЭВМ. Для их пересчета на выбранный вариант ЭВМ используются формулы:

Здесь

, , - значения для выбранной ЭВМ;

, , - значения для 3-ей ЭВМ из П1 таблицы 5 и 6;

, , , - индивидуальные коэффициенты производительности для выбранной и 3-ей ЭВМ из П1 пункт 6.2.3 таблица 3.

Для выполнения ограничений на верность восстановления сигнала необходимо, чтобы на этапе сбора данных частоты опроса датчиков были равномерными и лежали в заданных пределах:

f_i ³ f_oi

Здесь f_oi - частота опроса i-го датчика при выбранной разрядности АЦП (n_ацп), полученные на этапе информационного расчета; f_i - реальная частота опроса i-го датчика. f_i может быть больше f_oi для выбранной разрядности АЦП, но не меньше, так как в этом случае не будет обеспечено восстановление сигнала с заданной точностью.

Варьируя значение f_i, мы перемещаем РТ по плоскости ПНХ (по координате С_рт), тем самым выбирая выгодное для нас положение РТ, при котором обеспечивается восстановление сигнала с заданной точностью и минимальная стоимость используемого оборудования. Изменение стоимости происходит за счет выбора различных устройств (движение по координате r_рт).

При выборе f_i и устройств необходимо обеспечить выполнение условия:

r_n(С_рт) ³ r_рт(С_рт)+R₀

или

R₀³ r_n(С_рт)-r_рт(С_рт)

Здесь R₀ - допустимый резерв загрузки ЭВМ в РТ заданный в техническом задании.

Для удобства обозначим:

=f₀

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

Если пронумеровать все такты временной диаграммы элементами натурального ряда чисел, то последовательность номеров тактов, в которых вызывается i-ая задача, можно рассматривать как класс вычетов j _i по модулю r _i. Здесь j_i (начальная фаза) - номер такта временной диаграммы, в котором i-ая задача вызывается в первый раз; r_i (тактовое расстояние) - расстояние между соседними моментами запуска i-ой задачи, выраженное в тактах временной диаграммы. r_i и j_i - целые числа. Тогда частоты запуска задач:

f _i = f₀ / r _i

Рис. 5

Выбор параметров временной диаграммы f₀, j=[j₁,j₂,...,jm], r=[r₁,r₂,...,r_m], где m - число датчиков в системе, следует проводить по критерию минимума суммарной загрузки процессора. Чем меньше загрузка процессора, тем менее производительный процессор можно использовать в АСНИ, а следовательно снизить ее стоимость. Для поиска минимальной загрузки процессора необходимо решить следующую задачу:

При следующих ограничивающих условиях:

r_i £ f₀/f_oi - ограничение на погрешность восстановления, вытекающее из требования f_i ³ f_oi и f_i=f₀/r_i;

j_i(mod r_i)¹j_j(mod r_j) - требование, согласно которому в каждом такте временной диаграммы должно начинаться выполнение не более одной задачи;

н.о.к (r₁,r₂,r₃,...,r_m)£N₀ - ограничение на длину циклограммы, накладываемое оперативной памятью (н.о.к - наименьшее общее кратное).

На этапе предпроектного анализа целесообразно использовать следующий подход к выбору параметров временной диграммы. Примем тактовое расстояние r_i равным ближайшей к f₀/f_oi степени числа 2, меньшей f₀/f_oi, т.е. r_i=[f₀/f_oi]₂=2^a. При этом задача сведется к нахождению величины f₀, минимизирующей суммарную загрузку процессора

r_S(f₀)=r_рт(f₀)+r_д(f₀)

при следующем ограничении:

Здесь К_з.ц. - коэффициент загрузки циклограммы, характеризует долю ненулевых элементов в циклограмме, М - число датчиков в системе.

Функция r_S(f₀) имеет пилообразный характер, причем, локальные минимумы наблюдаются в “особых” точках, имеющих следующие значения:

S^(k,i) = f_oi × 2^k, i=1,2,...,m k=1,2,....

Значение частоты f₀, обращающее r_S в минимум, лежит на интервале [C₀, 2C₀] в одной из особых точек. Напомним, что:

Ограничивающее условие К_з.ц. можно записать в следующем виде:

Здесь K - число групп датчиков, М_j - число датчиков в j-ой группе, (2^к)_j - тактовое расстояние кратное степени числа 2 для j-ой группы.

Алгоритм определения параметров временной диаграммы состоит из следующих этапов:

Этап 1

Вычисление области поиска рабочей частоты циклограммы f₀: [C₀, 2С₀] – значения тактовой частоты циклограммы из указанного интервала должны выбираться по возможности наименьшими, что снизит требования на быстродействие КТС.

C₀=10×164+2×80+9×43+10×46=2647

Область поиска рабочей частоты циклограммы: [2647; 5294]

Этап 2

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

Выбор предпочтительных частот: С₀ £ f_oi×2^a £ 2C₀

Для каждой группы датчиков вычисляется: С_j = f_j × 2^a

Здесь f_j - частота опроса f_oi любого датчика j-ой группы.

C₁ = 164×2⁵ =5248

C₂ = 80×2⁶ =5120

C₃ = 43×2⁶=2752

C₄ = 46×2⁶=2944

Этап 3

Из рассчитанных для каждой группы С_j выбирается наименьшая С_j и проверяется в качестве тактовой частоты циклограммы - f₀. Вычисляются тактовые расстояния для каждой группы датчиков - значения (2^к)_j=[f₀/f_oi]_2.

По условию построения равномерной адаптивной циклограммы тактовые расстояния в нашем случае выбираются кратными степени двойки, что может иметь то преимущество, что в качестве задатчика циклограммы можно использовать двоичный счетчик. Если для выбранной частоты циклограмма может быть построена - удовлетворяется условие К_з.ц. £ 1, то f₀ считается допустимой и квазиоптимальной для всех особых точек r_S(f₀). В противном случае она отбрасывается и этап 3 повторяется для других по порядку возрастания С_j. При таком подходе полученная f₀ является минимальной из возможных, где загрузка процессора также минимальна.

Начинаем подбор с меньшей частоты C₃=2752:

Выбор исследуемой тактовой частоты циклограммы (2752 гц.) позволяет получить по крайней мере для одной дроби искомый результат, совпадающий со значением равным степени двойки, что уменьшает процент пустых тактов в циклограмме и позволяет считать это значение частоты наиболее предпочтительным. Выбор в качестве тактовых расстояний степеней двойки (4,32,16,256), меньших расчетных значений, обуславливается необходимостью обеспечения заданной точности восстановления сигнала.

Вычислим коэффициент заполнения циклограммы:

 

Так как условие К_з.ц. £ 1 не выполняется, то переходим к следующей частоте: С₄=2944

Таким образом в качестве С_рт выбирается С₁:

f₀= С_рт=С₁=2944 гц.

 

Этап 4

Вычисляются тактовые расстояния и частоты запуска для каждой прикладной задачи:

r_i = [f₀ / f_oi]₂

f_i = f₀ / r _i

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

Вычисляется длина циклограммы:

N_ц = max{r₁,r₂,...,r_m}

Номера датчикам назначаются следующим образом: 1-ая цифра - номер группы, которой принадлежит датчик, 2-ая цифра - номер датчика в группе.

Параметры временной диаграммы сведем в таблицу:

 

Таблица 9

№ датчика	Параметры циклограммы	Частота опроса, Гц
i	ri	ji	foi	fi

Длина циклограммы:

N_ц = max{r₁,r₂,...,r_m}= 64 такта.

Построим циклограмму (табл. 7):

Таблица 10

 

Заполнение циклограммы начинают с отведения тактов для датчиков с меньшими тактовыми расстояниями. Датчикам присваиваются двухзначные номера, где первая цифра означает номер группы датчика, а вторая цифра номер датчика в группе.

 

После построения циклограммы опроса рассчитываются r_рт и r_п для выбранного оборудования, и проверяется условие:

r_п ³ r_рт + R₀ (7)

Если необходимо делаются замены в составе аппаратно-программного обеспечения и вновь проверяется условие. Таким образом выбирается необходимое оборудование.

С целью уменьшения трудоемкости расчетов целесообразно на начальном этапе исключить из рассмотрения сочетание программно-аппаратных средств АСНИ, заведомо недопустимых по производительности. Для этих вариантов не выполняется хотя бы одно из следующих неравенств:

r_п(С₀) ³ r_рт(С₀) + R₀ (8)

МV_зсд+V_ос+N_ц+С₀ n₀ T_сб £ maxV_озу (9)

Здесь М - число датчиков в системе; V_зсд и V_ос - объем занимаемой ОП выбранными ЗСД и ОС соответственно; Т_сб - время сбора данных (задается в техническом задании); V_озу - имеющийся объем ОП в выбранной ЭВМ; n₀=[n_ацп/8], где[ ] означает дополнение до большего целого.

Объем занимаемой оперативной памяти ЭВМ в РТ рассчитывается по формуле:

V_озу= МV_зсд+V_ос+N_ц+С_рт n₀ T_сб

Также необходимо чтобы объем занимаемой памяти был меньше максимально допустимого объема ОП выбранной ЭВМ:V_озу£ maxV_озу.

Задачу выбора работоспособного варианта системы сбора и обработки данных по проектированной системе можно было бы решить простым перебором, но такой подход даже в данном случае (ЭВМ-7, интерфейс-3, УСД-4, ОС-3) предполагает проверку 7*3*4*3=252 возможных вариантов решения, что является достаточно трудоемкой задачей.

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

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

Шкала изменения стоимостей элементов системы представлена в таблице 8:

Таблица 8

ЭВМ _(без ОЗУ)

5.6

ОЗУ

5.6.7

1.2.3.4

УСД

ОС (с настройкой)

Интерфейс

Стоимость

10000 С у.е.

 

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

Также, рекомендуется проверять условие (9) на любом этапе нагрузочного расчета.

 

Так как для нашего варианта имеется только один алгоритм ЗСД, то можно рассчитать объем ОП занимаемый ЗСД:

МV_зсд=21×200=4200 байт

Также можно рассчитать объем ОП занимаемый циклограммой и регистрируемой измерительной информацией на этапе сбора данных (так как при выборе оборудования они остаются неизменными):

С_рт×n₀×T_сб+N_ц=2944×1×50+64=147264 байт

Минимальный объем ОП занимаемый ОС:

V_ос=5×1024=5120 байт

Таким образом минимальный требуемый объем ОП:

V_озу=МV_зсд+V_ос+N_ц+C_рт×n₀×T_сб=4200+5120+147264=156584 байт

V_озу > 128 кб=131072 байт

Поэтому выбор ЭВМ возможен только между 5-ой и 7-ой машиной, так как их максимальный объем памяти равен 256 кб. Дальнейший подбор устройств КТС рекомендуется проводить в такой последовательности, когда сначала подбираются более дорогостоящие составляющие КТС и ПО: ЭВМ, интерфейс, УСД и ОС. Так как все параметры задержек задаются в методическом пособии для третьей машины, соответствующие параметры для других ЭВМ определяются с помощью коэффициентов пересчета (см. начало нагрузочного расчета). В нашем случае, так как 5-ая ЭВМ дешевле, то целесообразнее выбрать ее. Проверим выполнение условия (7) для 5-ой ЭВМ:

Пересчитаем t_по на 5-ю машину:

мкс

Выберем самое быстродействующее оборудование и ПО:

* параллельный интерфейс сек

* 3-й вариант УСД сек

* 3-й вариант ОС (пересчитаем на 5-ю ЭВМ)

сек

Так как резерв по нагрузке для выбранных устройств значительно превышает допустимый, то выберем более дешевые и менее быстродействующие устройства. Так как львиная доля стоимости приходится на сеть связи, то желательно значительно снизить стоимость интерфейса. Возьмем самый дешевый 3-й вариант интерфейса, а также более дешевое УСД - 1-й вариант. Проверим выполнение условия (7):

t_усд=100 мкс

t_инт=100 мкс

Так как резерв больше допустимого, попробуем изменить вариант ОС. Возьмем 2-ую ОС:

Пересчитаем параметры ОС на 5-ю ЭВМ:

сек

сек

Так как резерв превышает допустимый, попробуем взять 1-й вариант ОС:

Пересчитаем параметры ОС на 5-ю ЭВМ:

сек

Резерв по нагрузке ЭВМ в рабочей точке обеспечивается, поэтому остановимся на этом варианте ОС, где резерв R=0,1491, что на 32% больше допустимого.

Рассчитаем объем ОЗУ необходимый для нашей АСНИ и число блоков памяти по 16 кб=16384 байт:

Как указывалось выше

МV_зсд=4200 байт

N_ц+C_рт×n₀×T_сб=147264 байт

Для выбранного варианта ОС: V_ос=10 кб=10240 байт

Таким образом требуемый объем памяти:

V_озу=МV_зсд+V_ос+N_ц+C_рт×n₀×T_сб=161704 байт

Количество блоков памяти: h= [V_озу / 16 кб]=[161704 / 16384]=[9,87]=10

Таким образом получим следующий состав аппаратно-программного обеспечения:

* 1-й вариант УСД t_усд=100 мкс

* 5-й вариант ЭВМ

* 3-й вариант интерфейса (последовательный) t_инт=100 мкс

* 1-й вариант алгоритма ЗСД t_по=21,8 мкс

* 1-ый вариант ОС r_п(С)=1-67,2×10^-6×С

* 10 блоков памяти по 16 кб h=10

Основные величины:

С_рт=2944

r_рт(С_рт)=0,6530

r_п(С_рт)=0,8021

Когда окончательно выбрано оборудование, рассчитываются следующие величины:

t - величина средних затрат процессорного времени на однократное выполнение одной задачи. Численно равна тангенсу угла наклона отрезка прямой соединяющей РТ с началом координат.

t = r_рт / С_рт

С_s – производительность системы в РТ, является проекцией точки пересечения ПНХ с прямой проходящей через начало координат и РТ. Определяется как корень следующего уравнения:

tC_s = r_п(C_s)

C_max - теоретический предел производительности системы. С_max = C_s при t®0. Определяется из следующего уравнения:

r_п(C_max)=0

Приведенные затраты процессорного времени на диспетчеризацию в РТ:

r_д(C_рт)=1-r_п(C_рт)

По полученным данным вычерчивается график ПНХ.

При t®0 С_s®С_max:

1-67,2×10^-6×С

Решив уравнение получим С_max=14880

r_рт(С_рт)=1-0,8021=0,1979

Построим ПНХ (рис.6):

Рис. 6

Основные величины и выражения:

* Вид ПНХ: r_п(С)= 1-67,2×10^-6×С

* Суммарная частота запуска прикладных задач в РТ: С_рт=2944

* Производительность системы: С_s=3747

* Теоритический предел производительности системы: С_max=20726

* Резерв загрузки ЭВМ в РТ: R_рт=0,1841

* Загрузка процессора в РТ: r_рт(С_рт)=0,6530

* Максимальная возможная загрузка процессора в РТ: r_п(С_рт)=0,8021