Методы и средства обработки аналоговых сигналов

Курсовая работа по "Системы реального времени"

СОДЕРЖАНИЕ

Задание

Аннотация

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

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

Топологический расчет

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

Литература

ЗАДАНИЕ НА КУРСОВОЕ ПРОЕКТИРОВАНИЕ

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

Исходные данные к проекту: Вариант № 17.

Таблица 1

ПАРАМЕТРЫ ОБЪЕКТА ИССЛЕДОВАНИЯ (ОИ)
№ группы датчиков	Количество датчиков в группе	Вид корреляционной функции сигналов на выходе датчиков	Параметр корреляционной функции
j	mj	Rj
1	9	exp(-22)	20
2	2	exp(-||)	80(0,08)
3	5	sin()/()	30
4	11	exp(-||)(cos()+sin())	0,8

Таблица 2

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

Таблица 3

КООРДИНАТЫ ДАТЧИКОВ
Порядковый номер датчика	№ датчика в циклограмме	Координаты датчиков
		Xi	Yi
1	11	47	77
2	12	25	80
3	13	72	82
4	14	60	75
5	15	85	55
6	16	75	72
7	17	57	87
8	18	87	71
9	19	72	51
10	21	75	35
11	22	71	26
12	31	80	40
13	32	56	26
14	33	82	20
15	34	60	15
16	35	71	14
17	41	36	17
18	42	22	18
19	43	20	25
20	44	11	54
21	45	22	40
22	46	30	22
23	47	43	85
24	48	27	52
25	49	40	25
26	410	15	35
27	411	25	60

Основные требования к АСНИ «Методы и средства обработки аналоговых сигналов»

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

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

Задание получила (число, месяц, год)

Подпись_____________

Задание выдал (преподаватель)

Подпись_____________

АННОТАЦИЯ

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

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

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

Информационный расчет (для интерполяции 1 порядка)

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

Топологический расчет

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

На первом этапе получены реальные частоты опроса датчиков и информационная производительность объекта исследований. Информационный расчет проводился для метода линейной интерполяции. В результате получена следующая суммарная частота следования выборок с датчиков:

C₀=950,51661+2632,708+513,378019+11•3,454699=1825 (Гц).

По минимуму суммарной производительности системы сбора определена наиболее подходящая разрядность АЦП - 6.

На втором этапе построена адаптивная равномерная циклограмма опроса системы датчиков АСНИ и получена ее рабочая частота Срт= 1748,786688 (Гц). Произведен расчет рабочей нагрузки АСНИ, построена ее потенциальная нагрузочная характеристика. Производительность системы в рабочей точке составила С_s =2538,07 ( Гц), резерв по нагрузке равен R = 0,10. По критерию минимальной стоимости сформирован оптимальный состав аппаратно-программного комплекса АСНИ. Окончательный вариант АСНИ состоит из 4-х 8-и канальных 6-и разрядных УСД, ЭВМ №1, ОС №1, последовательного интерфейса и 3-ех блоков памяти по 16 кбайт каждый. Общая стоимость системы Q_s=16834,735 у.е.

На третьем этапе расчета определена топология сети и размещение АСНИ в монтажном пространстве. Уточнено количество УСД:4.

На четвертом этапе произведена оценка показателей эффективности окончательного варианта АСНИ. Информационная избыточность АСНИ в целом равна:

_s = (_усд +1)( _вд +1)( _д +1) - 1 = (0,1748+1)(0,1648+1)(0,33+1) - 1 = 0,8199.

Стоимостной дисбаланс АСНИ, характеризующий асимметрию загрузки системы, равен:

ИНФОРМАЦИОННЫЙ РАСЧЕТ

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

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

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

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

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

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

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

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

При среднеквадратической ошибке - (ск): ²(n,t) ² (1)

При максимальной ошибке - (м): (n,t) (2)

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

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

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

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

Тогда: =₀*D

D=6*_s

D²=36*_s² => _s²=D²/36

²=36*_s²*₀²

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

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

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

(4)

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

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

(5)

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

(6)

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

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

,

где

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

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

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

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

B_i=n_if_i

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



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

f^opt=B^opt/n^opt

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

РАСЧЕТ

1-ая группа датчиков.

=20 [1/сек].

Вид модели сигнала

R()=exp(-²²)

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

²(n,t) ²

Таблица 4

n	fi	bi	ci
5	50,51661	252,5831	454,6495
6	44,94288	269,6573	404,4859
7	43,95093	307,6565	395,5584
8	43,71941	349,7553	393,4747
9	43,66247	392,9623	392,9623
10	43,6483	436,483	392,8347

n^opt=5

f^opt=50,51661

B^opt=252,5831

С₁=454,6495

Рис.4

2-ая группа датчиков.

=80 [1/сек].

m=2

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

Вид модели сигнала:

R()=exp*(-||).

²(n,t) ²

Таблица 5

n	fi	bi	ci
5	632,708	3163,54	1265,416
6	611,4564	3668,739	1222,913
7	606,3648	4244,553	1212,73
8	605,1051	4840,84	1210,21
9	604,7909	5443,119	1209,582
10	604,7125	6047,125	1209,425

n^opt=5

f^opt=632,708

B^opt=3163,54

С₁=1265,416

Рис. 5

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

=30 [1/сек];

m=5;

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

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

²(n,t) ²

Таблица 6

n	fi	bi	ci
5	13,37803	66,89014	66,89014
6	13,26425	79,58549	66,32124
7	13,23655	92,65584	66,18274
8	13,22967	105,8373	66,14834
9	13,22795	119,0516	66,13976
10	13,22752	132,2752	66,13761

n^opt = 5

f^opt = 13,37803

B^opt = 66,89014

С₁= 66,89014

Рис. 6

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

=0,8 [1/сек];

m=11;

Вид модели сигнала

R()=.

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

Таким об разом:

²(n,t) ²

Таблица 7

n	fi	bi	ci
5	3,454699	17,2735	38,00169
6	3,415579	20,49347	37,57137
7	3,406072	23,8425	37,46679
8	3,403712	27,22969	37,44083
9	3,403123	30,6281	37,43435
10	3,402975	34,02975	37,43273

n^opt = 4

f^opt = 3,454699

B^opt = 17,2735

С₁= 38,00169



Рис. 7

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

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

Таблица8

№ группы датчиков	Количество датчиков в группе	Параметры квантования
j	mj	foj, Гц	nацп
1	9	51	6
2	2	633	6
3	5	13	6
4	11	4	6

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

НАГРУЗОЧНЫЙ РАСЧЕТ

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

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

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

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

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

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

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

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

_n(C)=1-_д(C)

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

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

0_n(C)1

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

Рис. 8

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

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

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

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

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

_по - затраты процессорного времени не связанные с обменом (подготовительные операции, первичная обработка).

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

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

, , - значения для 3-ей ЭВМ;

, , , - индивидуальные коэффициенты производительности.

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

f_i f_oi

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

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

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

_n(С_рт) _рт(С_рт)+R₀ или

R_{0 n}(С_рт)-_рт(С_рт)

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

=f₀

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

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

f _i = f₀ / r _i

Рис. 9

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

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

r_i f₀/f_oi - ограничение на погрешность восстановления, вытекающее из требования f_i f_oi и f_i=f₀/r_i; _i(mod r_i)_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. При этом задача сведется к нахождению величины f₀, минимизирующей суммарную загрузку процессора

(f₀)=_рт(f₀)+_д(f₀)

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

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

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

S^(k,i) = f_oi 2^k,

i=1,2,...,m k=1,2,....

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

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

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

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

ЭТАП 1

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

РАСЧЕТ:

C₀ = 9•44.94288+2•611.45+5•13.26425+11•3.415579=1731.291

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

ЭТАП 2

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

Выбор предпочтительных частот:

С₀ f_oi2 2C₀

Для каждой группы датчиков вычисляется:

С_j = f_j 2

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

РАСЧЕТ:

C₁ = 44.942882⁶ = 2876.34432 [1731.291; 3462.583]

C₂ = 611.452² = 2445.8 [1731.291; 3462.583]

C₃ = 13.264252⁸ = 3395.6 [1731.291; 3462.583]

C₄ = 3.4155792⁹ = 1748.786688 [1731.291; 3462.583]

ЭТАП 3

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

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

Расчет:

Начинаем подбор с C₄ =1748,786688:

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

Так как условие К_з.ц. 1 выполняется, то в качестве С_рт выбирается С₄:

f₀ = С_рт = С₄ = 1748,786688 Гц.

ЭТАП 4

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

r_i = [f₀ / f_oi]₂

f_i = f₀ / r _i

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

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

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

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

РАСЧЕТ:

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

Таблица 9

№ Датчика	Параметры Циклограммы	Частота Опроса, Гц
i	ri	i	foi	fi
11	32	3	44,9428	44,9428
12	32	4	44,9428	44,9428
13	32	7	44,9428	44,9428
14	32	8	44,9428	44,9428
15	32	11	44,9428	44,9428
16	32	12	44,9428	44,9428
17	32	15	44,9428	44,9428
18	32	16	44,9428	44,9428
19	32	19	44,9428	44,9428
21	4	1	13,26425	719,08608
22	4	2	13,26425	719,08608
31	128	20	611,45	22,44144
32	128	23	611,45	22,44144
33	128	24	611,45	22,44144
34	128	27	611,45	22,44144
35	128	28	611,45	22,44144
41	512	31	3,415579	5,61786
42	512	32	3,415579	5,61786
43	512	52	3,415579	5,61786
44	512	55	3,415579	5,61786
45	512	56	3,415579	5,61786
46	512	59	3,415579	5,61786
47	512	60	3,415579	5,61786
48	512	63	3,415579	5,61786
49	512	64	3,415579	5,61786
410	512	68	3,415579	5,61786
411	512	71	3,415579	5,61786

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

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

Построим циклограмму:

Таблица 10

21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19	31	21	22	32	33	21	22	34	35	21	22	41	42
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19	43	21	22	44	45	21	22	46	47	21	22	48	49
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19	410	21	22	411		21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19	31	21	22	32	33	21	22	34	35	21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19	31	21	22	32	33	21	22	34	35	21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19	31	21	22	32	33	21	22	34	35	21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22
21	22	11	12	21	22	13	14	21	22	15	16	21	22	17	18
21	22	19		21	22			21	22			21	22

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

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

_{п рт} + R₀ (7)

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

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

_п(С₀) _рт(С₀) + 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], учитывающую характер изменения целевой функции при выборе последовательности шагов нахождения работоспособной системы минимальной стоимости (в нашем случае).

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

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

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

РАСЧЕТ:

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

МV_зсд=27200=5400 байт

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

С_ртn₀T_сб+N_ц=1748,786688 115+512= 26732 байт

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

V_ос=101024=10240 байт

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

V_озу=МV_зсд+V_ос+N_ц+C_ртn₀T_сб = 5400+26732+10240=42581 байт

V_озу < 64 кб=42581байт

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

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

мкс

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

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

сек

3-й вариант УСД

сек

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

сек

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

Проверим выполнение условия (7):

_усд=100 мкс

_инт=100 мкс

Так как нам необходимо добиться максимальной приближенности к резерву по нагрузке, то возьмем более дешевое ОС(2-й вариант):

Так как нам необходимо добиться максимальной приближенности к резерву по нагрузке, то возьмем самое дешевое ОС(1-й вариант):

В курсовом проекте рассмотрена только 1 ЭВМ, остальные шесть ЭВМ также проверены. На основе проверки всех семи ЭВМ выбрана наиболее оптимальное 1 ЭВМ с оптимальными величинами УСД, интерфейса и ОС.

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

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

МV_зсд = 5400 байт

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

Для выбранного варианта ОС:

V_ос = 5 кб = 5210 байт

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

V_озу = МV_зсд+V_ос+N_ц+C_ртn₀T_сб = 26732+5400+5210=37532 байт

Количество блоков памяти:

h= [V_озу / 16 кб]=[ 37532/ 16384]=[5,04]=3

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

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

1-ой вариант ЭВМ

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

1-й вариант алгоритма ЗСД _по=27 мкс

1-ый вариант ОС _п(С)=1-16810^-6С

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

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

С_рт=1748

_рт(С_рт) = 0,396796

_п(С_рт)=0,706336

Рассчитываются следующие величины:

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

= _рт / С_рт

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

C_s = _п(C_s)

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

_п(C_max)=0

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

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

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

РАСЧЕТ:

При 0 С_sС_max:

Построим ПНХ



Рис. 10

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

Вид ПНХ: _п(С) = 1-16810^-6С

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

Производительность системы: С_s = 2538,02

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

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

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

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

ТОПОЛОГИЧЕСКИЙ РАСЧЕТ

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

определить топологию сети связи АСНИ;

выбрать модель трассировки линии связи;

разместить оборудование АСНИ в монтажном пространстве по критерию минимума стоимости сети связи.

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

Рис. 11

В АСНИ датчики могут соединятся как непосредственное с УСД, так и через разъемы (псевдообъекты), тоже относится и к внешним устройствам (МР и ПУ) в смысле их соединения с ЭВМ. Аналогичным образом соединяются между собой УСД и ЭВМ (рис.12):

Рис. 12

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

ЭТАП 1. ОПРЕДЕЛЕНИЕ ТОПОЛОГИИ СЕТИ СВЯЗИ

Выбирается структура связи между УСД и ЭВМ в зависимости от выбранного варианта интерфейса. Различают древовидную и кольцевую схемы соединений (рис.9):

Рис. 13

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

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

Для выбранного типа интерфейса АСНИ необходимо использовать кольцевую схему соединения УСД и ЭВМ.

ЭТАП 2. ВЫБОР МОДЕЛИ ТРАССИРОВКИ ЛИНИЙ СВЯЗИ

Также необходимо выбрать модель трассировки межобъектных линий связи. На практике в качестве модели трассировки наиболее часто используют ортоганальную и евклидову метрики:

- ортоганальная метрика

- евклидова метрика

Здесь l_ij - длина линии связи между i-ым и j-ым объектами.

Так как линии связи в ортоганальной метрике несомненно короче, то за модель трассировки целесообразнее принять именно ее.

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

Приводятся графики дерева модели сети связи и размещения датчиков и устройств в монтажном пространстве.

В качестве модели трассировки линии связи примем евклидову метрику, так как она позволяет проложить более короткую линию связи, чем при евклидовой метрике.

Модель сети связи представим в виде следующей кольцевой иерархической системы (рис.14):

Рис.14

Диаграмма размещения датчиков и устройств в монтажном пространстве (рис.15):



Рис.15

ЭТАП 3. РАЗМЕЩЕНИЕ ОБОРУДОВАНИЯ АСНИ В МОНТАЖНОМ ПРОСТРАНСТВЕ ПО КРИТЕРИЮ МИНИМИЗАЦИИ СТОИМОСТИ СЕТИ СВЯЗИ

В нашей задаче имеется только 4уровня в иерархической модели сети связи. Для решения задачи используем следующий графический метод. Оптимизация проводится дважды: сначала минимизируется линия связи по координате X, а затем таким же образом по координате Y.

Введем понятие удельной стоимости сети связи:

Здесь x - искомая оптимальная координата объекта X, который соединен с объектами X_i координаты x_i которых известны. S_i - удельная стоимость линии связи от объекта X к объекту X_i. m - число объектов связанных с объектом X.

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

Удельная стоимость записывается для каждой координаты, уровня и фрагмента отдельно. При этом верхний индекс при Q указывает на координату, первый нижний индекс на уровень, а второй нижний индекс на фрагмент (или объект).

Вначале, составляется выражение Q для 0-го уровня, что в общем виде можно записать как:

Здесь j - номер фрагмента.

В этом случае Q будет характеризовать стоимость сети связи между объектами 0-го уровня и объектами 1-го уровня. j - номер объекта на первом уровне, который связан с объектами уровня 0.

Рассмотрим координату Х (иерархическая модель сети связи приведена на рисунке 14):

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

Так как в выражении для Q разность (x-x_i) стоит под знаком модуля, то при определении производной необходимо следить чтобы эта разница всегда была больше 0, т.е. если значение производной определяется на промежутке где x>x_i, то эта разность запишется как (x-x_i), а если значение производной определяется на промежутке где x<x_i, то эта разность запишется как (x_i-x).

Разность (x-x_i) назовем узловой точкой с координатой x_i. Для каждой Q координатная ось X разбивается узловыми точками на ряд отрезков в каждом из которых производная имеет постоянное значение. Такая координатная ось изображается для всех Q 0-го уровня. На эту ось наносятся значения координат узловых точек и значения производной в образуемых отрезках. Покажем это на примере, но для начала рассмотрим производную и ее график:

Определим диапазон изменения производной. Примем, что x>x_i, тогда

Пусть теперь x<x_i, тогда

Т.е. диапазон изменения производной [-40; +40].

Пусть теперь x<x_i, тогда

Т.е. диапазон изменения производной [-36; +36].

Пусть теперь x<x_i, тогда

Т.е. диапазон изменения производной [-32; +32].

Пусть теперь x<x_i, тогда

Т.е. диапазон изменения производной [-12; +12].

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

Рис. 16

Остальные графики производных строятся аналогично.

Удельная стоимость линии связи от УСД к ЭВМ (магистраль связи) накладывает ограничение на область поиска оптимальных координат. Эта удельная стоимость приводится в П4 таблица 3. Она зависит от типа выбранного интерфейса. В нашем случае она равна 5. Поэтому область поиска оптимальных координат уменьшается до диапазона изменения производной [-5; +5]. Таким образом область поиска координаты X для УСД1 равна (60; 72) (заштрихованная область на рис. 16).

₁₁ и ₁₂ - оконечные точки области поиска оптимальных координат.

Уточним, что в нашем примере Q₀₁ характеризует связь УСД1 с датчиками 1-ой группы, Q₀₂ - связь УСД2 с датчиками 1-ой, 2-ой и 3-ей групп, Q₀₃ - связь УСД3 с датчиками 4-ой группы, Q₀₄ -связ УСД4 с датчиками 4-ой группы, Q₀₅ - связь ЭВМ с двумя МР и одним ПУ.

Приведем диаграмму координатных осей, о которых говорилось выше, для Q₀₁^x, Q₀₂^x, Q₀₃^x ,Q₀₄^x ,Q₀₅^x (рис.17).

Под этими координатными осями располагаются координатные оси для Q₁₁, Q₁₂, Q₁₃, Q₁₄, Q₁₅, которые характеризуют те же связи УСД с датчиками и ЭВМ с ВУ, что и Q уровня 0, но с учетом ограничений накладываемых на область поиска стоимостью магистрали. На этих координатных осях размечаются только те узловые точки, которые принадлежат области поиска. Значения производной копируются с уровня 0, за исключением того, что за областью поиска значение производной принимается равным удельной стоимости магистрали. В нашем случае этими значениями будут +5 и -5.

Под этими линиями располагается еще одна линия для Q₂₁, которая характеризует стоимость связи между объектами 1-го и 2-го уровня. В нашем случае это будет связь УСД1, УСД2,УСД3,УСД4 и ЭВМ. На эту линию переносятся все узловые точки с линий Q₁₁^х, Q₁₂^х, Q₁₃^х, Q₁₄^x ,Q₁₅^x , а значения производной от Q₂₁^х в получившихся отрезках образуется как сумма значений производных в соответствующих отрезках на линиях Q₁₁, Q₁₂, Q₁₃ , Q₁₄, Q₁₅

Рис.17

В первую очередь определяется координата единственного объекта на самом нижнем уровне иерархии. В качестве значения координаты выбирается любое значение из отрезка, где производная ближе всего к 0 (чаще это концевые точки такого отрезка). В нашем случае это будет координата разъема, примем ее равной: x₂₁=27. Остальные координаты (всех объектов уровня 1) определяются по правилу:

Где ₁₁ и ₁₂ - оконечные точки области поиска оптимальных координат для соответствующих Q_1i. Определим оптимальные координаты для УСД1, УСД2, УСД3,УСД4 и ЭВМ:

УСД1

УСД2

УСД3

УСД4

ЭВМ

Рассмотрим координату У

Определим диапазон изменения производной. Примем, что у>у_i, тогда

Т.е. диапазон изменения производной [-40; +40].

Пусть теперь у<у_i, тогда

Т.е. диапазон изменения производной [-36; +36].

Пусть теперь x<x_i, тогда

Т.е. диапазон изменения производной [-32; +32].

Пусть теперь x<x_i, тогда

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

Рис.18

Рис.19

В первую очередь определяется координата единственного объекта на самом нижнем уровне иерархии. В качестве значения координаты выбирается любое значение из отрезка, где производная ближе всего к 0 (чаще это концевые точки такого отрезка). В нашем случае это будет координата разъема, примем ее равной: y₂₁=40.

Остальные координаты (всех объектов уровня 1) определяются по правилу:

Где ₁₁ и ₁₂ - оконечные точки области поиска оптимальных координат для соответствующих Q_1i.

Определим оптимальные координаты УСД1, УСД2, УСД3,УСД4 и ЭВМ:

УСД1

УСД2

УСД3

УСД4

ЭВМ

После того как найдены оптимальные координаты и по оси X и по оси Y для всех объектов АСНИ, данные о них сводятся в таблицу 11:

Таблица 11

	X	Y
УСД1	60	75
УСД2	71	35
УСД3	27	40
УСД4	27	35
ЭВМ	27	50

ЭТАП 4. РАСЧЕТ СТОИМОСТИ СЕТИ СВЯЗИ

Приводится ортоганальный план размещения АСНИ в монтажном пространстве:

Рис.20

В соответствии с выбранной метрикой рассчитывается длина магистрали l_маг., проверяется: является ли она допустимой для выбранного интерфейса. Рассчитывается стоимость магистрали Q_маг..

В нашем примере выбран интерфейс с кольцевой магистралью и соединение между устройствами следующее: УСД1-УСД2-УСД3- УСД4-ЭВМ- УСД1 (см. рис.25). Расчет длины магистрали выполняется следующим образом:

Для нашего интерфейса длина магистрали не должна превышать 150 м, проверим это условие:

l_маг.=149,36 м > 150 м

Т.е. данное условие выполняется.

Рассчитаем стоимость магистрали:

Q_маг.=5*149,36=746,8руб.

Qсовм=2087,935 руб

Рассчитывается длина линии связи между датчиками и УСД, а также между ВУ и ЭВМ, и ее стоимость Q_л.св

Совместим расчет длины линии связи и ее стоимости в одном выражении:

Вычисляется конечная стоимость сети связи в АСНИ:

Q_сети=Q_маг.+Q_л.св. [руб]

АНАЛИЗ РЕЗУЛЬТАТОВ И ОЦЕНКА ПОКАЗАТЕЛЕЙ ЭФФЕКТИВНОСТИ АСНИ

Указывается перечень аппаратных и программных средств оптимальной АСНИ и стоимость последней:

Q_s=Q_усд+Q_эвм+Q_бп+Q_ос+Q_сети

Здесь Q_усд - стоимость всех УСД в АСНИ;

Q_эвм, Q_бп - стоимость ЭВМ и блоков оперативной памяти соответственно;

Q_ос - стоимость операционной системы;

Q_сети - стоимость сети связи.

Перечень аппаратно-программного обеспечения:

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

1-й вариант ЭВМ

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

1-й вариант алгоритма ЗСД _по=5,4545 мкс

1-ой вариант ОС _п(С)=,706336

3 блоков памяти по 16 кб

Расчет стоимости АСНИ:

Q_усд=4*500=2000 у.е.

Q_эвм=810³=8000 у.е.

Q_озу=10³h=10³3 =3000 у.е.

Q_ос=1000 у.е. (с настройкой)

Q_сети=2087,935+746,8=2834,735 у.е.

Таким образом, стоимость АСНИ:

Q_s=2000+8000+3000+1000+2834,735 =16834,735 у.е.

Параметры временной диаграммы работы АСНИ на этапе сбора данных сводятся в Таблицу 12

Таблица 12

Параметры временной диаграммы
№ задачи	Параметры циклограммы	Частота вызова, Гц	Избыточность по загрузке
I	ri	i	foi	fi	i=(fi-foi)
11	32	3	44,9428	44,9428	0
12	32	4	44,9428	44,9428	0
13	32	7	44,9428	44,9428	0
14	32	8	44,9428	44,9428	0
15	32	11	44,9428	44,9428	0
16	32	12	44,9428	44,9428	0
17	32	15	44,9428	44,9428	0
18	32	16	44,9428	44,9428	0
19	32	19	44,9428	44,9428	0
21	4	1	13,26425	719,0860	0,159897
22	4	2	13,26425	719,0860	0,159897
31	128	20	611,45	22,44144	-0,13343
32	128	23	611,45	22,44144	-0,13343
33	128	24	611,45	22,44144	-0,13343
34	128	27	611,45	22,44144	-0,13343
35	128	28	611,45	22,44144	-0,13343
41	512	31	3,415579	5,61786	0,000499
42	512	32	3,415579	5,61786	0,000499
43	512	52	3,415579	5,61786	0,000499
44	512	55	3,415579	5,61786	0,000499
45	512	56	3,415579	5,61786	0,000499
46	512	59	3,415579	5,61786	0,000499
47	512	60	3,415579	5,61786	0,000499
48	512	63	3,415579	5,61786	0,000499
49	512	64	3,415579	5,61786	0,000499
410	512	68	3,415579	5,61786	0,000499
411	512	71	3,415579	5,61786	0,000499

Параметр характеризует избыточную загрузку процессора выполнением i-ой задачи в следствии запуска последней с избыточной частотой f_i.

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

К_ц=₀ / _рт

Здесь - нижний теоретический предел загрузки процессора в РТ.

Информационная избыточность, вносимая при работе по временной диаграмме, оценивается коэффициентом:

Здесь - суммарная частота запуска прикладных задач на этапе сбора данных.

Производительность АСНИ в РТ вычисляется по формуле:

Резерв АСНИ по нагрузке в рабочей точке на этапе сбора данных:

,

где С_s - производительность системы в РТ.

Информационная производительность АСНИ в РТ на этапе сбора данных (скорость записи информации в ОЗУ) не превышает величины:

B_s=8n₀C / K_сж [бит/сек]

Здесь n₀ - коэффициент определяемый по формуле n₀=[n_ацп / 8], где = []-значение выражения в скобках, дополненное до большего целого;

К_сж - средний коэффициент сжатия данных в системе (с учетом датирования).

Объем оперативной памяти, требуемый для регистрации измерительной информации на этапе сбора данных, равен:

V_озу=B_sT_сб / 8 [байт]

Информационная избыточность, вносимая при записи данных в ОЗУ, оценивается коэффициентом:

_д= 8n₀ / n_ацп - 1

Информационная избыточность АСНИ в целом равна:

_s=(_усд+1)( _вд+1)( _д+1)-1

Стоимостной дисбаланс АСНИ характеризует асимметрию загрузки системы и рассчитывается по следующей формуле:

Здесь Q_i - стоимость i-ой компоненты АСНИ; _i - загрузка (или коэффициент использования) i-ой компоненты. Компонентами АСНИ являются УСД, ЭВМ, блоки памяти, сеть связи, ОС. Их стоимости соответственно: Q_усд, Q_эвм, Q_бп, Q_сети, Q_ос. Загрузка компонент:

_усд = С_рт / С_усд

_эвм = С_рт / С_s

_бп = h без приведения к большему целому [h]

_сети = К_з.ц.

_ос = 1

₀=С₀(_по+_инт+_усд)=1731,291 (2710^-6+10010^-6+10010^-6)=0,3930

_рт = 0,3967

К_ц = _0/рт = 0,3930 / 0,3967=0,99

С =44,942889+719,086082+22,45+5,6211=2016,64

С₀ = 1731,291

_вд = (С - С₀) / С₀ =(2016,64-1731,291) / 1731,291 = 0,1648

_рт / С = 0,3967/2016,64= ,000196 сек = 196 мксек

Что при длительности такта:

₀ = 1 / f₀ = 1 / 2016,64 =495 мксек

составляет более половины процессорного времени, затрачиваемого на выполнение полезной работы.