Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Исследование средней длительности выбросов огибающей нормального шума. Методические указания к лабораторной работе №4 по курсу «Автоматизация обработк

.pdf
Скачиваний:
0
Добавлен:
09.08.2026
Размер:
409 Кб
Скачать

3037

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования

«ЛИПЕЦКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

Кафедра физики и биомедицинской техники

ИССЛЕДОВАНИЕ СРЕДНЕЙ ДЛИТЕЛЬНОСТИ ВЫБРОСОВ ОГИБАЮЩЕЙ НОРМАЛЬНОГО ШУМА

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

к лабораторной работе № 4

по курсу «Автоматизация обработки экспериментальных данных»

Липецк Липецкий государственный технический университет

2015

УДК 537.86 (07) И 889

Составители: В.Ф. Осинин, Д.А. Подлесных, С.И. Шарапов, А.А. Демидова, С.Ф. Четвериков

Рецензент – доктор технических наук, профессор В.Н. Малыш

И 889 Исследование средней длительности выбросов огибающей

нормального шума [Текст]: методические указания к лабораторной работе №4 по курсу «Автоматизация обработки экспериментальных данных» / сост.: В.Ф. Осинин [и др.]. – Липецк: Изд-во Липецкого государственного технического университета, 2015. – 12 с.

Методические указания предназначены для студентов 3,4-го курсов направления «201000.62 Биотехнические системы и технологии» очной формы обучения.

Табл.2 . Ил.2. Библиогр.: 5 назв.

© ФГБОУ ВПО «Липецкий государственный технический университет», 2015

ЛАБОРАТОРНАЯ РАБОТА № 4

ИССЛЕДОВАНИЕ СРЕДНЕЙ ДЛИТЕЛЬНОСТИ ВЫБРОСОВ ОГИБАЮЩЕЙ НОРМАЛЬНОГО ШУМА

Цель работы

1.Получить функцию распределения вероятностей и распределение сред-

него числа выбросов для нормального шума.

2.Найти распределение средней длительности выбросов огибающей нор-

мального шума

3.Проверить работу статистического анализатора в режиме одновремен-

ного измерения распределения вероятностей и среднего числа превы-

шения пороговых уровней.

4.Определить эффективную полосу пропускания узкополосного прием-

ника.

Блок-схема аппаратурного комплекса

Аппаратурный комплекс для исследования средней длительности выбро-

сов огибающей нормального шума показан на рис.1.

Рис. 1. Функциональная схема аппаратурного комплекса

3

Комплекс состоит из четырех основных блоков:

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

тенный усилитель с полосой пропускания частот от 3 до 30 кГц и коэффи-

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

2. Приемно-регистрирующий блок. Он состоит из: узкополосного прием-

ника с диапазоном рабочих частот 2 150 кГц и набором эффективных полос пропускания от 10 до 3000 Гц; низкочастотного фильтра для выделения оги-

бающей узкополосного процесса; многоканального регистратора с динамиче-

ским диапазоном измеряемых напряжений порядка 80 дБ.

3.Блок калибровки и настройки содержит генератор нормального шума с эффективной полосой генерации 50 кГц и коммутирующее устройство.

4.Блок обработки информации и управления экспериментом (МЭВМ).

Для пересчета опорных уровней анализатора помех в эффективные зна-

чения напряжения шумового генератора используется генератор широкопо-

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

теля при отключенной антенне. Это также необходимо для последующей кали-

бровки уровней напряженности поля в абсолютных единицах, исходя из извес-

тных параметров системы.

Нормальный шум, создаваемый генератором, через антенный усилитель,

приемник и статистический анализатор поступает в решающий блок. Таким образом, мы получим статистические данные шума в полосе Δfэ на частоте приема f. Время анализа выбирается порядка 1 минуты, поскольку процесс

4

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

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

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

мальный фон, который в координатах Релея описывается прямой линией с угло-

вым коэффициентом наклона равным двум.

Нахождение длительности пересечений некоторого порогового

уровня огибающей напряжения случайного процесса

Будем считать, что стационарный случайный процесс задан своим зако-

ном распределения и корреляционной функцией. Зафиксируем уровень a и бу-

дем искать вероятностные характеристики выбросов над уровнем а. Найдем,

прежде всего, вероятность Р(х>а, t), т.е. вероятность того, что за малый интер-

вал времени t будет иметь место одно пересечение уровня а снизу вверх.

Очевидно, для этого должны быть одновременно выполнены два условия x(t)<a и x(t + t)>a. Если t мало, то

P x > a, Δt Δt

 

yW2 a, y dy = ΔtN a ,

(1)

 

>y>0

 

 

5

где W2 a, y - плотность вероятности.

Величину

yW2 a, y dy = N a можно рассматривать как среднюю

>y>0

 

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

Рис.2

Теперь рассмотрим задачу о среднем времени Та, соответствующем выполнению условия х>а. Очевидно, время Та является частью всего рассматриваемого интервала Т. Удобно ввести единичную импульсную случайную функцию la(t) по условию

Тогда

a

 

 

 

 

 

 

 

 

 

 

 

 

dx .

(2)

P l t

 

= 1 = p x

 

t

 

> a

=

 

xW

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>x>a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ta

=

 

xW x dx .

 

 

(3)

Очевидно, Ta = θi ;

 

 

 

 

 

 

 

 

 

 

i=1

 

T

 

 

 

>x>aт

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

Теперь найдем среднее число выбросов за время Т (рис. 3). Так как за время t может быть только один выброс с вероятностью р(х>а, t), то за время

Т будет nа,τ

 

выбросов, где na,τ =

T

P x > a, Δt , т.е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Δt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yW2 a, y dy = TN a

 

 

 

(4)

 

 

 

 

 

 

 

na,τ = T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>y>0т

 

 

 

 

 

 

 

 

 

 

 

Отсюда видно, что N(а) представляет собой среднюю частоту выбросов,

т.е. среднее число выбросов в секунду.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N(a)= nа, T / Т

 

 

 

 

 

 

 

 

 

 

(5)

Средняя продолжительность одного выброса равна

 

 

 

 

 

 

 

 

 

 

 

 

 

T

xW x dx

 

 

 

 

 

xW x dx

 

 

 

 

 

Na,τ

=

 

>x>aт

 

=

>x>aт

 

 

 

 

 

(6)

θ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

T

 

yW2 a, y dy

N

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>y>0т

 

 

 

 

 

 

 

 

 

 

 

 

а= Р(х а)/ N(a).

 

 

 

 

 

 

 

 

 

 

(7)

Рис. 3. К вычислению длительности выбросов случайного процесса x(t), где i – продолжительность одного i-го выброса

7

Распределение средней длительности импульсов огибающей

напряжения нормального шума

Итак, для стационарного эргодического процесса, когда усреднение по ансамблю заменяется усреднением по времени, при известных распределениях вероятности Р(a) и среднего числа выбросов N(a) средняя длительность выбро-

сов огибающей напряженности поля может быть найдена из выражения

(a)= P(a)/ N (a) , (8)

где (a) средняя длительность выбросов на уровне a, P(a)- вероятность превышения уровня a, N(a) среднее число выбросов огибающей с положительной производной в единицу времени на уровне a.

Для огибающей напряжения V(t) нормально-флюктуационного процесса распределения P(V) и N(V) известны, поэтому формула для длительности импульсов принимает вид

θ =

P(V >V0 )

 

=

 

 

σ

,

(9)

 

 

 

 

 

 

N(V >V0 )

 

δω

 

V0

 

где V - напряжение огибающей

нормального шума,

- эффективное

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

Исследование распределения средней длительности выбросов

огибающей нормального шума удобно проводить в системе координат

(y=lg (V), x=lg(V/ )), где данная зависимость представляется прямой линией с угловым коэффициентом, равным -1 , т.е.

lg (V)=lg θ =

P(U >U0 )

 

=

 

 

σ

-lg(V/ ).

(10)

N(U >U0 )

δω

U0

Уклонение показаний счетчиков анализатора для нормально-

флюктуирующего напряжения V от представленной линейной зависимости говорит о нарушении работы установки в этом режиме.

8

Отметим, что из этой формулы вытекает важный практический результат: при известных распределениях P(V) и N(V) огибающей нормального шума по формуле легко рассчитывается функция (V), поведение которой в координатах lg (V), lg(V/ ) при V= после простейших вычислений сразу же дает значение среднеквадратичной (эффективной) полосы спектра флюктуаций приемника, вычисление которой другими методами связано с некоторыми трудностями.

ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

1.

Задать продолжительность измерения с помощью пульта управления на

 

статистическом анализаторе. Время анализа выбирается порядка 1

 

минуты, так как процесс стационарный и эргодический.

2.

Дождаться окончания анализа и, используя пульт управления, получить

 

данные по wi - времени и ni - количеству превышений огибающей

 

нормального шума для каждого порогового уровня Vi статистического

 

анализатора.

 

 

 

3.

Полученные данные записать в таблицу 1.

4.

Записать общее время эксперимента

T

w

 

i .

 

 

 

 

i 1

5.

По данным таблицы 1 рассчитать амплитудное распределение

 

вероятностей гауссова шума ( P(V> V0)

и lg (-ln P(V> V0))) по формуле

 

n

 

 

 

 

w (V )

/T, значения которых записать в таблицу 2.

 

Pn(V>V0)= i i

 

i 1

 

 

 

6.Рассчитать распределение среднего числа выбросов огибающей нормального шума N(V> V0) и lg (N((V> V0)) по формуле ni/T=Ni(V>V0),

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

7.Рассчитать распределение средней длительности выбросов огибающей нормального шума θ(V>V0) и lg (θ(V>V0) по формуле θ(Vi)=P(Vi)/N(Vi),

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

9

8.По данным таблицы 2 построить график зависимости распределения средней длительности выбросов огибающей нормального шума в координатах lg θ(V>V0) от lg V0.

9.По методу наименьших квадратов найти коэффициенты k и b для формулы y=kx+b, где y=lg θ(V), x=lgV.

10. Найти

эффективную

 

полосу пропускания, используя формулу

 

P(U >U0 )

 

 

 

 

 

 

 

 

 

 

σ

b=lg θ =

 

 

=

 

 

 

 

.

N(U >U0 )

δω

U0

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

Таблица 1

Измеренные значения

i

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

Рассчитанные данные

 

 

 

 

 

 

 

 

 

 

i

P(V>V0)

N(V>V0)

 

θ(V>V0)

lg θ(V>V0)

V0

lg V0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

………

………

 

………

………

………

………

 

 

 

 

 

 

 

 

12

 

 

 

 

 

4

0,602

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

 

 

 

 

 

2

0,301

 

 

 

 

 

 

 

 

14

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

15

 

 

 

 

 

0

 

 

 

 

 

 

 

 

10

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]