Исследование средней длительности выбросов огибающей нормального шума. Методические указания к лабораторной работе №4 по курсу «Автоматизация обработк
.pdf
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 ) |
|
= |
|
2π |
|
σ |
, |
(9) |
|
|
|
|
|
|||||
|
N(V >V0 ) |
|
δω |
|
V0 |
|
|||
где V - напряжение огибающей |
нормального шума, |
- эффективное |
|||||||
напряжение нормального шума, - среднеквадратичная полоса спектра флюктуаций на выходе приемника, V0 - пороговый уровень напряжения.
Исследование распределения средней длительности выбросов
огибающей нормального шума удобно проводить в системе координат
(y=lg (V), x=lg(V/ )), где данная зависимость представляется прямой линией с угловым коэффициентом, равным -1 , т.е.
lg (V)=lg θ = |
P(U >U0 ) |
|
= |
|
2π |
|
σ |
-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 ) |
|
|
|
|
|
|
|
|
|
|
2π |
|
σ |
|||
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
