Анализ результатов и оценка показателей эффективности АСНИ

СОДЕРЖАНИЕ

1. Задание………………………………………………………………..3

2. Введение………………………………………………………………5

3. Информационный расчет ……………………………………………6

4. Нагрузочный расчет…………………………………………….…...19

5. Анализ результатов………………………………………………… 52

6. Литература……………………………………………………………58

Задание на курсовое проектирование

По дисциплине “Управление в реальном времени” на тему: Расчет и проектирование автоматизированной системы научных исследований и обучения (АСНИО) «Оптимизация характеристик средств сбора и первичной обработки измерительной информации»

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

Исходные данные к курсовой работе: Вариант № 101.

Таблица 1

ПАРАМЕТРЫ ОБЪЕКТА ИССЛЕДОВАНИЯ (ОИ)
№ группы датчиков	Количество датчиков в группе	Вид корреляционной функции сигналов на выходе датчиков	Параметр корреляционной функции
j	m­j	Rj	a
			2.5
		sin(at)/(at)
		(1+at)/exp(at)

Таблица 2

МЕТРИЧЕСКИЕ ХАРАКТЕРИСТИКИ АСНИ
Допустимая ошибка восстановления сигнала в % от его шкалы, по критерию оценки погрешности восстановления сигнала: наибольшего отклонения – (м) или среднеквадратического отклонения – (ск)	Время сбора данных ( сек. )	Допустимый резерв загруз-ки ЭВМ в рабочей точке
e0	Tсб	R0
5,0 (ск)		0,10

Основные требования к АСНИ “Оптимизация характеристик средств сбора и первичной обработки измерительной информации”

Система должна обеспечивать одновременный сбор и первичную обработку информации с заданного количества датчиков в течении заданного интервала времени.

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

Резерв по загрузке ЭВМ в рабочей точке на этапе сбора данных должен быть не менее заданной величины.

Аннотация

В курсовой работе произведен расчет и приведено обоснование автоматизированной системы, предназначенной для научных исследований и производства испытаний газотурбинных двигателей (АСНИ ГТД).

Исходные данные приведены в начале пояснительной записки к курсовой работе.

Проект включает следующие этапы:

- Информационный расчет

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

- Анализ результатов

На первом этапе получены реальные частоты опроса датчиков и информационная производительность объекта исследований. Информационный расчет проводился в два этапа. Так как при использовании для восстановления сигналов метода ступенчатой интерполяции частота следования выборок сигнала не позволила сформировать необходимый КТС АСНИ, пришлось провести вторично информационный расчет для метода линейной интерполяции. В результате получена следующая суммарная частота следования выборок с датчиков: C₀=9×12+3×28+4×224+10×386=4948 (гц). По минимуму суммарной производительности системы сбора определена наиболее подходящая разрядность АЦП – 6.

На втором этапе построена адаптивная равномерная циклограмма опроса системы датчиков АСНИ и получена ее рабочая частота Срт=6176 гц. Произведен расчет рабочей нагрузки АСНИ, построена ее потенциальная нагрузочная характеристика. Производительность системы в рабочей точке составила С_s = 7016 гц, а резерв по нагрузке равен R = 0,135. По критерию минимальной стоимости сформирован оптимальный состав аппаратно-программного комплекса АСНИ.

На четвертом этапе произведена оценка показателей эффективности окончательного варианта АСНИ. Информационная избыточность АСНИ в целом равна: h_s = 0,986 Стоимостной дисбаланс АСНИ, характеризующий асимметрию загрузки системы, равен:

Информационный расчет

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

* выбрать способ восстановления сигналов по дискретным отсчетам;

* рассчитать параметры квантования сигналов с датчиков по критерию минимума информационной производительности ОИ;

* сделать предварительный выбор устройств сбора данных.

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

Расчет параметров квантования сигналов осуществляется по следующим формулам:

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

При среднеквадратической ошибке – (ск):

e²(n,Dt)£ D² (1)

При максимальной ошибке – (м):

e(n,Dt)£ D (2)

Здесь D - граница для допустимой ошибки восстановления сигнала, задаваемая в процентах от диапазона сигнала (шкала сигнала) и зависит от дисперсии сигнала - :

· при равномерном распределении сигнала

· при нормальном распределении сигнала

По заданию сигнал распределен по нормальному закону. По закону 3sза достоверные значения с вероятностью большей 90% принимаются только те, которые лежат на далее 3sвлево и вправо от точки математического ожидания. Ширина этого диапазона D именуется шкалой.

где - плотность распределения

амплитуды сигнала

3

- дисперсия сигнала

Рис. 1

D

U -средняя составляющая

напряжения сигнала

U - напряжение сигнала

Вероятность превышения напряжения сигнала некоторого уровня равна

Тогда: D=e₀*D, D=6*s_s, D²=36*s_s² => s_s²=D²/36

D²=36*s_s²*e₀² (3)

Если ошибка составляет 1% от шкалы сигнала с нормальным распределением, то это значит при критерии максимальной ошибки:

и соответственно при критерии СКО:

Среднеквадратическая ошибка при ступенчатой интерполяции имеет вид

(4)

Независимо от выбранного критерия оценки погрешности ошибка восстановления сигнала включает две аддитивные составляющие:

- погрешность квантования сигнала по уровню

(5)

- погрешность дискретизации сигнала по времени

(6)

Здесь n - число двоичных символов отводимых на кодирование одного отсчета, s_s² - дисперсия сигнала, R(t) - корреляционная функция сигнала.

Подставив выражения (3), (4), (5), (6) в (1) получим:

Упростим это выражение:

Разрешив это неравенство, как:

f_i£j(n, e),

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

Для удобства интегрирования корреляционную функцию R(t) целесообразно разложить в ряд Маклорена с точностью до двух первых ненулевых членов ряда.

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

B_i=n_i×f_i

Здесь B_i - информационная производительность i-го датчика, n_i - число двоичных символов отводимых на кодирование одного отсчета с i-го датчика (разрядность АЦП), f_i - частота опроса i-го датчика.

Рис. 2

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

f^opt=B^opt/n^opt

Поскольку датчики внутри групп эквивалентны по частоте, то указанный расчет достаточно произвести только для одного датчика из каждой группы.

Ая группа датчиков

a=1[1/сек].

Вид модели сигнала R(t)=

Таблица 4

n
fi
Bi

n^opt=4

f^opt=13

B^opt=52

Ая группа датчиков

a=2.5[1/сек].

Вид модели сигнала R(t)=

Таблица 5

n
fi
Bi

n^opt=4

f^opt=32

B^opt=128

Я группа датчиков.

a=200[1/сек].

Вид модели сигнала R(t)=sin(at)/(at)

Таблица 6

n
fi
Bi

n^opt=4

f^opt=239

B^opt=954

Ая группа датчиков

a=200[1/сек].

Вид модели сигнала R(t)=(1+at)/exp(at)

Ряд Маклорена

R(0)=1/e⁰=1,

,

e²(n,Dt)£ D²

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

Таблица 7

n
fi
Bi

n^opt=4

f^opt=413

B^opt=1651

После расчета оптимальной производительности каждого датчика рассчитывается оптимальная производительность ОИ (объекта исследования):

Здесь m - число датчиков в системе.

или

Здесь K - число групп однотипных датчиков; m_j - число датчиков в j-ой группе; B_j^opt - оптимальная производительность одного из датчиков в j-ой группе.

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

=9×52+3×128+4×954+10×1651=21178

По полученным данным производится предварительный выбор вариантов и требуемое количество УСД. Также необходимо выполнить технологическое условие, состоящее в том, что разрядность АЦП должна быть одной и той же для всех УСД используемых в АСНИ (в данном случае n=6 или n=10). Отход от некоторых оптимальных параметров квантования приведет к некоторому увеличению информационной производительности ОИ (избыточности) на величину:

Здесь m - число датчиков в системе; B_iⁿ - информационная производительность i-го датчика при n удовлетворяющем технологическому условию; B_i^opt - оптимальная производительность i-го датчика или

Здесь K - число групп датчиков; m_j - число датчиков в системе; B_jⁿ - информационная производительность любого датчика j-ой группы при n, удовлетворяющем технологическому условию; B_j^opt - оптимальная производительность любого датчика j-ой группы.

Для перехода к 6-ти или 10-ти разрядной сетки рассчитываются величины:

Переход осуществляется к той разрядности АЦП которая обладает меньшей разностью DB_n. Таким образом выбирается разрядность АЦП n_АЦП в АСНИ.

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

Результаты расчетов параметров квантования сводятся в таблицу (см. таблицу 4).

Информационная избыточность вносимая УСД характеризуется следующим коэффициентом:

DB₆=9×(68-52)+3×(168-128)+4×(1339-954)+10×(2319-1651)=8484

DB₁₀=9×(112-52)+3×(277-128)+4×(2223-954)+10×(3849-1651)=28043

Так как DB₆ <DB₁₀, то выбираем n_ацп=6.

Результаты расчетов параметров квантования сведем в таблицу 4:

Таблица 8

№ группы датчиков	Количество датчиков в группе	Параметры квантования
J	mj	foi, Гц	nацп

Информационная избыточность вносимая УСД:

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

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

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

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

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

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

Снятие показания с датчика и его первичную обработку будем считать прикладной задачей. Для опроса 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₀=9×12+3×28+4×224+10×386=4948

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

Этап 2

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

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

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

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

C₁ = 12×2⁹ =6144

C₂ = 28×2⁸ =7168

C₃ = 224×2⁵=7168

C₄ = 386×2⁴=6176

Этап 3

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

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

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

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

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

С₄=6176

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

f₀= Срт=С₄=6176 гц.

Этап 4

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

r_i = [f₀ / f_oi]₂

f_i = f₀ / r _i

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

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

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

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

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

Таблица 6

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

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

N_ц = max{r₁,r₂,...,r_m}= 512 тактов.

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

Таблица 7