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

Дипломы / Бакалаврская_работа_(Вешкурцев)_2

.pdf
Скачиваний:
101
Добавлен:
30.03.2015
Размер:
2 Мб
Скачать
☆

Данная программа была написана на языке высокого уровня С# 4.0 с применением .NET Framework 4.0. Вид окна программы приведен на рисунке 10. Листинг программы приведен в приложении A.

Рисунок 10 – Вид окна программы фильтрации и усреднения

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

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

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

На рисунке 11 приведен сигнал, получаемый на выходе случайно-неоднородной среды. Данный сигнал повторяет форму квазидетерминированного сигнала, подаваемого с генератора нижних частот на источник излучения. Сигнал на выходе случайно-неоднородной среды содержит помеху промышленной частоты 50 Гц, а также высокочастотные составляющие. Данная помеха вносит определенную погрешность при измерении оценок вероятностных характеристик сигнала.

41

Рисунок 11 – График сигнала, полученного на выходе случайно-неоднородной среды

На рисунке 12 представлен график оценки плотности вероятности данного сигнала, полученный при помощи характериометра.

Рисунок 12 – График оценки плотности вероятности с помехой

На рисунке 13 показан график оценки плотности вероятности данного сигнала, полученный по результатам пяти экспериментов после проведения фильтрации и усреднения.

Рисунок 13 – График оценки плотности вероятности без помехи и с усреднением

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

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

42

Таблица 2 – Относительная разность начальных и центральных моментных функций сигнала без обработки и с обработкой

Момент

Сигнал без фильтрации и

Отфильтрованный и усредненный

Относительная

усреднения

сигнал

разность

 

 

 

 

 

m1

0,0272

0,0394

30,99%

 

 

 

 

m2

1,5020

1,5120

0,66%

 

 

 

 

m3

0,0969

0,1350

28,24%

 

 

 

 

m4

7,6307

7,6754

0,58%

 

 

 

 

M1

0,0000

0,0000

-

 

 

 

 

M2

1,5012

1,5105

0,61%

 

 

 

 

M3

-0,0258

-0,0438

41,21%

 

 

 

 

M4

7,6268

7,6682

0,54%

 

 

 

 

По таблице можно количественно оценить, насколько отфильтрованные и усредненные результаты отличаются от результатов без фильтрации. Увеличение точности измерения достигает

1/0,4 = 2,5 раза.

Таким образом, можно сделать вывод, что внедрение средств автоматизации позволяет повысить точность измерений до 2-2,5 раз по сравнению с данными без фильтрации и усреднения.

6 Анализ результатов исследования

6.1 Анализ результатов экспериментальных исследований

Нами в одинаковых лабораторных условиях (Т °С = 23±2 °С, давление 760 мм рт. ст.) с постоянной погрешностью 8 % исследованы случайно-неоднородные среды следующих видов:

воздушная смесь среды проживания человека;

красное сухое вино «Каберне»;

красное сухое вино «Chateau Potensac Medoc 2004 года»;

красное сухое вино «Chateau Grand Monteil»;

красное сухое вино «Chateau Cadillac 2005 года»;

фальсифицированное красное сухое вино «Каберне».

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

Эта гипотеза проверялась многократно при разных начальных условиях эксперимента, в число которых входили следующие переменные: частота амплитудной модуляции электромагнитного поля (ЭМП); энергия ЭМП; амплитуда сигнала на выходе приемника излучений при исследовании воздушной среды. С помощью усилителя амплитуда сигнала устанавливалась равной 2,8 В. Значения переменных сведены в таблице 3, а результаты экспериментальных исследований приведены в Приложении C.

43

Таблица 3 – Значения переменных, характеризующих эксперимент

Наименование переменной величины

Значение величины

 

 

Частота модуляции

60 Гц; 120 Гц

 

 

Энергия ЭМП

1,28 мВт; 2,86 мВт; 5,53 мВт

 

 

 

Амплитуда сигнала на выходе приемника

2,8

В

излучений при исследовании воздушной

 

 

среды

 

 

 

 

 

Геометрические размеры или объем среды

0,4

дм3

 

 

 

Вкаждом эксперименте зафиксирована потеря энергии ЭМП в среде, т. к. амплитуда сигнала на выходе приемника излучений получилась меньше установленной ранее 2,8 В. Объяснение этому явлению дает молекулярно-кинетическая теория, согласно которой ЭМП взаимодействует с мельчайшими частицами вещества. При этом энергия ЭМП расходуется на поглощение, на рассеяние волны ЭМП, на отражение волн случайно-неоднородной средой. Приемник излучений нашей лабораторной установки фиксировал суммарную потерю энергии ЭМП, не разделяя еѐ на три составляющих: отражение, рассеяние, поглощение. Других физических эффектов взаимодействия ЭМП со случайно-неоднородной средой нами не наблюдалось. Полагаем, индуцированного излучения радиоволн самой средой не существовало, поскольку эксперименты проводились с использованием слабых электромагнитных полей. Такие поля не разрушают среду и обеспечивают неразрушающий контроль вещества.

Вприложении C вероятностные характеристики сигнала на выходе приѐмника излучений представлены с помощью функций и числовых значений. Любая из них получена в результате обработки 4096 значений реализации сигнала длительностью 0,165 с. Функции повторяют форму

кривых, рассчитанных по известным формулам 33 при нахождении вероятностных характеристик квазидетерминированного сигнала. Такое может быть, когда изучаемый четырехполюсник является линейным. При этом он не искажает форму вероятностных характеристик, а делает другими только значения числовых параметров этих характеристик. Следовательно, каждая изученная нами случайно-неоднородная среда является линейной в диапазоне частот 0…120 Гц. Этот диапазон частот получается при работе радиопреобразователя нашей лабораторной установки, в которой источником ЭМП является инерционная лампа накаливания типа МС 6,3×0,3, подключенная к генератору квазидетерминированного сигнала с частотой 30 Гц, 60 Гц. По-видимому, инерционность лампы накаливания и приемника излучений (фотодиод ФД256) не позволила нам увидеть быстрые флуктуации параметров сигнала, вызванные взаимодействием молекул вещества с волнами электромагнитного поля. Такое несовершенство радиопреобразователя в будущем следует устранить.

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

Рассмотрим это явление на примере исследования красного сухого вина «Каберне» (смотри приложение C). Увеличение частоты модуляции ЭМП с 60 Гц до 120 Гц приводит к возрастанию потерь энергии поля, при этом плотность вероятности мгновенных значений сигнала на выходе приемника излучений сосредоточена в пределах от -1,5 В до +1,5 В при частоте 60 Гц и в пределах от -1,4 В до +1,4 В при частоте 120 Гц. Такая же ситуация наблюдается при анализе функции

44

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

Таблица 4 - Начальные и центральные моментные функции красного вина «Каберне» при энергии ЭМП, равной 2,86 мВт

Частота

m1

m2

m3

m4

M1

M2

M3

M4

модуляции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60 Гц

0,051

1,092

0,106

5,26

0

1,08

-0,06

5,26

 

 

 

 

 

 

 

 

 

120 Гц

0,037

0,951

0,077

4,5

0

0,95

-0,03

4,5

 

 

 

 

 

 

 

 

 

Для проведения анализа данных таблицы 4 нами предполагается следующая гипотеза относительно физического смысла моментных функций высокого порядка. Наряду с центрированными и нецентрированными сигналами вводим термин энергетический и «неэнергетический» сигналы. При этом считаем m3 математическим ожиданием «неэнергетического» нецентрированного сигнала, m4 - дисперсия «неэнергетического» нецентрированного сигнала. Тогда М3 - математическое ожидание центрированного «неэнергетического» сигнала, а М4 – дисперсия центрированного «неэнергетического» сигнала. На наш взгляд, моментные функции высокого порядка сохраняют физический смысл моментных функций 1-го и 2-го порядков на уровне нелинейностей и вариаций величин второго уровня значимости. В соответствии с нашей гипотезой нелинейность второго уровня значимости случайно-неоднородной среды изменяется с плюса для «неэнергетического» нецентрированного сигнала на минус для «неэнергетического» центрированного сигнала, тогда как дисперсии названных сигналов совпадают для обеих частот модуляции. Эта нелинейность среды зависит от энергии ЭМП. В таблице 5 показаны значения моментных функций при исследовании красного вина «Каберне» полем с энергией 5,53 мВт.

Таблица 5 - Начальные и центральные моментные функции красного вина «Каберне» при энергии ЭМП, равной 5,53 мВт

Частота

m1

m2

m3

m4

M1

M2

M3

M4

модуляции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60 Гц

0,031

1,092

0,051

4,88

0

1,27

-0,04

4,88

 

 

 

 

 

 

 

 

 

120 Гц

0,031

1,2

0,078

5,9

0

1,2

-0,03

5,5

 

 

 

 

 

 

 

 

 

По-видимому, увеличение энергии ЭМП приводит к массовой поляризации молекул вещества, при этом случайно-неоднородная среда становится изотропной и линейной при модуляции радиоволн сигналом с частотой 60 Гц. Начальная моментная функция m3 уменьшилась в два раза по сравнению с m3 =0,106 (смотри таблицу 4). Для частоты модуляции 120 Гц нелинейность среды одинакова, а повышение энергии ЭМП привело к росту дисперсии m4, M4.

45

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

Таблица 6 - Начальные и центральные моментные функции разного вещества при энергии ЭМП, равной 2,86 мВт

Моментальная

Частота модуля 60 Гц

 

Частота модуляции 120 Гц

 

функция

 

 

 

 

 

 

 

 

вино

вино

вино

вино

вино

вино

вино

вино

 

«Каберне»

«Cha au

«Chateau

«Chateau

«Каберне»

«Cha au

«Chateau

«Chateau

 

 

Potensac

Grand

Cadillac

 

Potensac

Grand

Cadillac

 

 

Medoc

Monteil»

2005

 

Medoc

Monteil»

2005

 

 

2004

 

года»

 

2004

 

года»

 

 

года»

 

 

 

года»

 

 

 

 

 

 

 

 

 

 

 

m2

1, 09

2,3

1,39

1,67

0,95

2,7

1,13

1,51

m3

0,106

0,101

0,11

0,13

0,077

-0,65

-0,04

0,135

 

 

 

 

 

 

 

 

 

m4

5, 26

12,85

6,98

8,59

4,5

15,71

5,47

7,65

M2

1, 09

2,3

1,39

1,67

0,95

2,7

1,13

1,51

 

 

 

 

 

 

 

 

 

M3

-0,06

0,017

-0,058

-0,09

-0,03

0,519

0,07

-0,04

 

 

 

 

 

 

 

 

 

M4

5,26

12,85

6,98

8,59

4,5

15,7

5,47

7,67

Значения оценок вероятностных характеристик, представленных в таблице 6, разные у каждого вида красного сухого вина. Показательно выглядят характеристики сухого вина «Chateau Potensac Medoc 2004 года». Для этого вещества характеристики отличаются широким диапазоном изменений и зависимостью от частоты модуляции амплитуды радиоволны. Так, например, отношение начальных моментов m4/m2=5,6, а смешанных моментов m3/М3=5,9. Данный признак обязательно следует учитывать при разработке интегрированного показателя качества вещества. Для других веществ таблицы 6 эти отношения меньше, однако они не опускаются ниже 1,4. Следовательно, вариации вероятностных характеристик сигнала, полученного в результате взаимодействия ЭМП с веществом, таковы, что могут служить материалом для контроля интегрального показателя качества продукта. Забегая вперед, скажем, что предложенный нами признак для определения качества вещества использует все моментные функции таблицы 6.

Обратим внимание на следующий интересный факт. Ценовая шкала исследованных нами красных сухих вин начиналась от 0,2 тыс. руб. и заканчивалась 1,85 тыс. руб. Начальная и центральная моментные функции 2-го порядка точно распределили сухие вина по ценовой шкале. Красное сухое вино «Chateau Potensac Medoc 2004 года» по всем показателям занимает 1 место,

второе место - вино «Chateau Cadillac 2005 года», третье место - вино «Chateau Grand Monteil»,

четвертое место – вино «Каберне», причем это вино было упаковано в картонную тару. Моментные функции 2-го порядка уменьшали свои значения в соответствии с качеством красного вина. Следовательно, можно прогнозировать контроль качества вещества по значениям этих моментных функций.

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

46

Возможность выявления фальсифицированного вещества проиллюстрируем с помощью данных таблицы 7, подготовленной с помощью материалов приложения C.

Таблица 7 - Начальные и центральные моментные функции фальсифицированного красного сухого вина «Каберне» при энергии ЭМП, равной 2, 86 мВт

Моменталь

Частота модуля 60 Гц

 

 

Частота модуляции 120 Гц

 

ная

 

 

 

 

 

 

воздушная

вино

фальсифицирова

воздушная

вино

фальсифицирова

функция

среда

«Каберне»

нное

вино

среда

«Каберне»

нное

вино

 

 

 

 

 

 

 

 

 

 

 

«Каберне»

 

 

 

«Каберне»

 

 

 

 

 

 

 

 

 

 

m2

3,21

1, 09

2,63

 

3,22

0,95

2,86

 

m3

-0,41

0,106

-0,39

 

-0,985

0,077

-0,48

 

m4

19,8

5, 26

15,17

 

19,8

4,5

16,9

 

 

 

 

 

 

 

 

 

 

M2

3,21

1, 09

2,63

 

3,19

0,95

2,86

 

M3

0,344

-0,06

0,25

 

0,76

-0,03

0,41

 

 

 

 

 

 

 

 

 

 

M4

19,8

5,26

15,15

 

19,7

4,5

16,9

 

 

 

 

 

 

 

 

 

 

Для исследования фальсификат вина «Каберне» приготовлен нами путем добавления в сухое красное вино «Каберне» 40 % воды из водопровода, выдержанной в сосуде в течение суток.

Моментные функции фальсификата вина «Каберне» существенно отличаются от таких же функций при исследовании нефальсифицированного вина «Каберне». Значения функций отличаются в 2 и более раз. Этого вполне достаточно для разработки показателя, характеризующего фальсификат вещества. Помимо моментных функций, другие вероятностные характеристики сигнала на выходе приемника излучений изменяются при взаимодействии ЭМП с фальсифицированным веществом.

Кроме этого, в таблице 7 приведены вероятностные характеристики сигнала приемника излучений при взаимодействии ЭМП с воздушной средой. Отличие характеристик воздушной среды и вина «Каберне» значительное. Это ещѐ раз подтверждает, что радио-статистический метод исследования случайно-неоднородной среды показывает однозначный результат.

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

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

1)вещество представляет собой среду с детерминированными параметрами и линейной передаточной функцией, которая не изменяет законы распределения и параметры сигнала, а изменяет только количественные параметры законов распределения;

2)на интервале времени 0,165 с, равном длительности реализации сигнала на выходе источника излучений, вещество представляет собой случайно-неоднородную среду с постоянными параметрами;

47

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

(флуктуации) возникают на интервале периода колебания электромагнитной волны, равном 25×10-12 с. Если фотодиод ФД 256 имеет полосу пропускания 12 кГц, то его постоянная времени

больше 80 мкс. Следовательно, имеется усреднение флуктуаций параметров и мгновенных значений сигнала в 80×10-6/25×10-12 = 3,2×106 раз (в миллион и более раз);

4)среда является случайно-неоднородной и взаимодействует с электромагнитной

волной, имеющей амплитудную модуляцию с частотой 30 Гц и 60 Гц. Поскольку волна взаимодействует с молекулами вещества, которые свое положение меняют за время 10-10 с, то на

огибающей электромагнитной

волны с периодом Т = 1/30+3,33×10-2 с будет

присутствовать

N = 3,33×10-2/10-10 = 3,33×108

раз соприкосновение вектора напряженности

молекулы с

вектором напряженности внешнего поля . Для изотропного вещества вектор и вектор

направлены в одну сторону. В результате этого огибающей электромагнитной волны будут возникать впадины глубиной . Эти впадины можно увидеть, если рассматривать период колебания сигнала приемника излучений через 0,1 пс. Это очень высокая частота просмотра, равная 10 ГГц. Имеющаяся в нашем распоряжении аппаратура не позволяет сделать вышеописанное, она имеет полосу пропускания 1 МГц.

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

6.2 Разработка алгоритмов вычисления интегрального показателя качества вещества

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

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

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

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

48

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

в прямоугольной трехмерной системе координат по осям откладываются значения функций M2, M3, M4;

находится площадь треугольника, построенного на данных вершинах;

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

√ (

)(

)(

 

),

 

 

 

 

где

 

(

 

)

 

– стороны треугольника. Эти стороны равны:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

√

 

 

√

 

 

√

.

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

Пример построения треугольника по значениям центральных моментных функций приведен на рисунке 14. Блок-схема алгоритма представлена в приложении B.

z

М2

0

М3

М4

у

х

Рисунок 14 – Построение треугольника по центральным моментным функциям

49

6.3 Расчет интегрального показателя качества

По разработанному алгоритму, блок-схема которого представлена в приложении B, был произведен расчет интегральных показателей качества для четырех видов вина. Результаты расчетов для частоты модуляции 60 Гц и энергии ЭМП 2,86 мВт приведены в таблице 8.

Таблица 8 – Расчет интегральных показателей качества четырех видов вина для частоты модуляции 60 Гц и энергии ЭМП 2,86 мВт

 

M2

M3

M4

a

b

c

S

стоимость, руб

Chateau Potensac Medoc

2,30

0,02

12,85

2,30

12,85

13,05

14,78

1850

 

 

 

 

 

 

 

 

 

Сha au Cad lla 2005

1,67

-0,09

8,59

1,67

8,59

8,75

7,18

675

 

 

 

 

 

 

 

 

 

Chateau Grand Monteil

1,39

-0,06

6,98

1,39

6,98

7,12

4,86

505

 

 

 

 

 

 

 

 

 

Каберне

1,09

-0,06

5,26

1,09

5,26

5,37

2,87

200

 

 

 

 

 

 

 

 

 

Столбцы M2, M3, M4 содержат соответственно значения центральных моментных функций второго, третьего и четвертого порядков, столбцы a, b, c – длины сторон треугольника, S – площадь треугольника, т. е. интегральный показатель качества.

Примечательным является тот факт, что интегральные показатели качества вин выстроились в соответствии с их стоимостью.

Расчет интегрального показателя качества воздушной среды при тех же начальных условиях дает значение 31,99, что показывает минимум двукратное отличие интегрального показателя воздушной среды от вещества. Данный факт подтверждает различие между вероятностными характеристиками воздушной среды и вещества.

50

Соседние файлы в папке Дипломы