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

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

СОДЕРЖАНИЕ

Задание

Аннотация

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

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

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

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

Литература

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

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

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

Таблица 1

Таблица 2

Таблица 3

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

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

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

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

Подпись_____________

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

Подпись_____________

АННОТАЦИЯ

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

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

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

Информационный расчет (для интерполяции 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^opt=5

f^opt=50,51661

B^opt=252,5831

С₁=454,6495

Рис.4

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

=80 [1/сек].

m=2

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

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

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

²(n,t) ²

Таблица 5

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^opt = 5

f^opt = 13,37803

B^opt = 66,89014

С₁= 66,89014

Рис. 6

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

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

m=11;

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

R()=.

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

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

²(n,t) ²

Таблица 7

n^opt = 4

f^opt = 3,454699

B^opt = 17,2735

С₁= 38,00169

Рис. 7

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

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

Таблица8

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

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

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

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

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

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

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

Снятие показания с датчика и его первичную обработку будем считать прикладной задачей. Для опроса 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

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

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

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

Таблица 10