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

Точность представления спектральной характеристики определяется разрешением по частоте

(2.3)

Обратное дискретное преобразование Фурье определяется соотношением (2.4):

(2.4)

Из сравнения формул (2.2) и (2.4) следует, что они отличаются знаком показателя экспоненты, множителем перед знаком суммы, а также переменной суммирования. Это позволяет строить единые программы для прямого и обратного преобразований Фурье.

Применяя формулу Эйлера, выражение (2.2) можно привести к виду (2.5):

(2.5)

где

(2.6)

Оценивание спектральной плотности мощности (СПМ) с помощью дискретного преобразования Фурье осуществляется по формуле (2.7):

(2.7)

Где X(k) - дискретное преобразование Фурье (спектральная характеристика) временного ряда

, соответствующего процессу x(t);

T - период дискретизации процесса x(t);

N - длина временного ряда.

Черта в правой части формулы (2.7) означает операцию осреднения. Применение формулы (2.7) без операции осреднения приводит к получению "грубой" оценки СПМ. Формула (2.7) позволяет вычислить оценку СПМ посредством статистического осреднения модуля спектральной характеристики совокупности данных, поделенного на длину записи данных. Статистическое осреднение необходимо здесь потому, что ординаты спектральной характеристики являются случайными величинами изменяющимися для каждой используемой реализации случайного временного ряда

.

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

и длина записи .

Можно показать, что оценки СПМ приближенно имеют распределение

с n степенями свободы, где . Более того, для достаточно больших n, например, , распределение аппроксимируется гауссовским (нормальным) распределением. В этом случае нормированное стандартное отклонение (стандартное отклонение, связанное с оцениваемой величиной, т.е. процентная ошибка, или, в статистической терминологии, "коэффициент разброса") определяется соотношением (2.8):

(2.8)

Величину

называют стандартной ошибкой.

Если

, то .

Последний результат означает, что вычисление оценки СПМ с использованием полной длины временного ряда имеет стандартную ошибку, равную 100%.

Если отрезок T_p поделить на m участков, то в этом случае

.

Подставляя полученный результат в (2.8), найдем

.

Таким образом, для повышения точности оценивания СПМ необходимо исходный временной ряд длины N разбить на m участков длины N_у, вычислить для каждого i-го участка

по формуле (1), а затем найти осредненную оценку по формуле

.

Следует иметь в виду, что разрешение по частоте в рассмотренном случае определяется из соотношения

. Число степеней свободы для найденной оценки СПМ можно найти следующим образом

.

Следовательно, для повышения степеней свободы и, соответственно, статистической устойчивости оценок СПМ необходимо увеличивать число участков для осреднения.

Повышение числа степеней свободы можно достичь другим способом – осреднением по частотам.

Сглаженная оценка

,

(2.9)

полученная осреднением l соседних оценок спектральной характеристики, имеет распределение

с числом степеней свободы, равным примерно 2l. Это следует из теории о сложении величин, имеющих распределение .

Следует отметить, что разрешение по частоте в данном случае определится из соотношения

.

Поскольку операция осреднения линейная, оценку СПМ можно найти, комбинируя осреднение по участкам с осреднением по частотам. При этом сначала выполняется осреднение по участкам, а затем – по частотам. При осреднении по m участкам с последующим осреднением l соседних спектральных оценок в итоге получаются оценки, число степеней, свободы которых равно

. Разрешение в этом случае равно .

3 РАЗРАБОТКА ТЕХНИЧЕСКОГО ОБЕСПЕЧЕНИЯ

В состав технического обеспечения включены информационно-измерительный канал и персональная ЭВМ.

3.1 Структура АСНИ

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

· подсистема измерений – информационно-измерительный канал (ИИК);

· подсистема передачи данных;

· подсистема обработки данных;

· подсистема визуализации и документирования результатов.

Обобщенная структура АСНИ представлена на рисунке 3.1

Рисунок 3.1 – Обобщенная структура АСНИ

На рисунке 3.1 приняты следующие обозначения:

ИИК- информационно-измерительный канал;

УСО – устройство связи с объектом;

БД - база данных.

Состав и структура ИИК приведены на рисунке 3.2.

Рисунок 3.2 – Структура информационно-измерительного канала

На рис.3.2 приняты следующие обозначения:

Д - датчик;

ДУ – дифференциальный электронный усилитель;

ФНЧ - фильтр нижних частот;

НУ – нормирующий усилитель;

АЦП - аналого-цифровой преобразователь.

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

ФНЧ выполняет функции противомаскировочного фильтра и служит для подавления высокочастотных составляющих сигнала в целях исключения ошибок аналого-цифрового преобразования.

НУ предназначен для согласования входного напряжения АЦП с динамическим диапазоном изменения преобразуемого аналогового сигнала.

3.2 Расчёт дифференциального усилителя

В простейшем случае в качестве измерительного усилителя может быть использован операционный усилитель (ОУ) в дифференциальном включении (рисунок 3.2). При выполнении условия R₁/R₂=R₃/R₄ усиление дифференциального сигнала намного больше усиления синфазного сигнала и коэффициент ослабления синфазного сигнала (КОСС) будет максимальным.

Рисунок 3.3 - Схема простейшего измерительного усилителя

В соответствии с техническим заданием:

U_c.max = 3 мВ

U_сф = 1,2 В – синфазная помеха

d = 3 % – допустимый процент подавления синфазной помехи

R_c = 300 Ом – внутреннее сопротивление датчика

,

(3.1)

Определяется коэффициент усиления:

.

(3.2)

К_синф определяется из следующих соображений: если резисторы имеют допуск