CC BY
1Естественные науки Нагуслаева И.Б. Старший научный сотрудник
Область научных интересов - распространение радиоволн, сейсмоэлектромагнетизм Ангархаева Л.Х. Кандидат физико-математических наук, доцент
Область научных интересов - математическое моделирование в геофизике, прямые и обратные задачи, интерпретация радиоимпедансных зондирований Башкуев Ю.Б. Доктор технических наук, профессор.
Основной круг научных интересов - электродинамика слоистых структур, радиоинтроскопия горного массива, излучение и распространение радиоволн. УДК 621.371+550.837.7
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СПЕКТРАЛЬНО-КОРРЕЛЯЦИОННЫХ, СТАТИСТИЧЕСКИХ И ЭЛЕКТРОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ФИЗИКО-
ТЕХНИЧЕСКИХ СИСТЕМ
Аннотация. Рассмотрен лабораторный практикум для моделирования физико-технических систем, имеющий широкие функциональные возможности. Программный комплекс использован при изучении курсов «Теория передачи сигналов», «Основы геоэлектромагнетизма», выполнении курсовых, дипломных проектов и УИРС, а также дистанционного обучения.
Ключевые слова: Моделирование, обработка сигналов, интерпретация, прямая и обратная задачи.
MATHEMATICAL MODELING OF THE SPECTRAL CORRELATING, STATISTICAL AND ELECTRODYNAMIC CHARACTERISTICS OF THE APPIED PHYSICS SYSTEMS Naguslaeva I.B., Angarhaeva L.H., Bashkuev Yu.B.
Abstract. The practical laboratory work for applied physics systems modeling with the broad functional specter is observed in the article. The program complex is used in studying the courses "the theory of transferring the signals" and "the basis of geomagnetism", in writing the course research papers, degree research papers and UIRS, and also in the remote educational process.
Key words: Modeling, processing of signals, interpretation, direct and inverse problems.
Введение
Повышение уровня фундаментальной математической и физико-технической подготовки актуально в связи с существенным повышением требований к подготовке молодых специалистов по современным телекоммуникационным и геоинформационным технологиям со стороны работодателей. Студенты 1-2 курсов достаточно уверенно владеют основами компьютерных и телекоммуникационных технологий благодаря массовому внедрению сотовой связи и Интернета. Однако одновременно с этим следует отметить снижение общей математической и физической «культуры мышления» основной массы студентов, истоки которой закладываются и формируются еще в начальной школе. В высшей школе существует серьезное противоречие между математической и физико-технической подготовкой студентов. Причина этого явления заключается в отсутствии тесных связей между учебными курсами, а также в том, что, как правило, преподаватели кафедр математического цикла глубоко не вовлечены в научную работу современного физико-технического уровня из-за большой загруженности учебной работой.
Цель статьи - показать связь и внутреннее единство математического и физико-технического подходов при компьютерном моделировании электромагнитных процессов в линейных и нелинейных системах на примере двухсеместрового лабораторного практикума. В статье развивается подход к определению «оптимальной» структуры лабораторных работ, описаны соответствующие математические (цифровые) модели и результаты их практиче-
ского применения по дисциплинам «Теория передачи сигналов» и «Основы геоэлектромагнетизма» для студентов физико-технических и телекоммуникационных специальностей Бурятского государственного университета и Иркутского государственного университета путей сообщения. Следует отметить, что профессиональная научная деятельность авторов доклада связана с Сибирским отделением РАН. Преподавание в ВУЗе позволяет нам отчетливо увидеть достоинства и недостатки в подготовке молодых специалистов, дать им правильную практическую оценку и тем самым повысить эффективность обучения.
В нашем подходе учебный курс рассматривается как множество элементов знаний (наиболее значимых математических и физическо-технических понятий, свойств, утверждений, законов и правил). Под типовым лабораторным заданием нами понимается задание, которое направлено на выполнение и проверку некоторого подмножества элементов знаний за академическое занятие. Структурой лабораторного занятия будем называть набор типовых заданий, включенных в конкретную лабораторную работу, а оптимальной структурой - набор заданий, который позволяет студенту понять сущность физических процессов и явлений и их математическое описание адекватной моделью, а преподавателю - сделать объективный вывод о степени усвоения рассматриваемой темы, удовлетворяющей определенному критерию. Лабораторный практикум включает необходимые и традиционно сложившиеся темы (810 тем) по теории сигналов и электродинамике. Он базируется на пакете программ «Редактирование и обработка сигналов» [1,2] и программном комплексе «Импеданс» [3,4]. Перечень работ по разделам курса высшей математики:
1. « Ряды и интеграл Фурье. Корреляционный анализ»
- Исследование спектров сигналов.
- Гармонический синтез и анализ периодических сигналов.
- Корреляционный анализ детерминированных сигналов.
- Исследование сигналов с линейной частотной модуляцией (ЛЧМ).
2. «Теория вероятностей и математическая статистика»
- Исследование законов распределения случайных сигналов.
- Прохождение случайных сигналов через линейные и нелинейные цепи.
- Исследование корреляционных функций и энергетических спектров случайных процессов.
3. «Интерполяция функций»
- Дискретизация непрерывных сигналов во времени (теорема Котельникова).
4. «Прямые и обратные задачи»
- Интерпретация радиоимпедансных зондирований методом подбора.
- Интерпретация радиоимпедансных зондирований методом регуляризации А.Н.Тихонова.
1. Спектрально-корреляционные и статистические характеристики сигналов Пособия [1,2] содержат описания и порядок выполнения лабораторных работ. Для удобства приведены методические рекомендации по использованию пакета программ ЕОБ'^ Разработанный лабораторный практикум с использованием пакета программ «Редактирование и обработка сигналов» позволяет применять любую операционную среду, начиная с Windows-95. Программа моделирует лабораторный стенд с полем для графиков исследуемых сигналов и необходимым набором контрольно-измерительной аппаратуры. Общий вид главного меню программы EDSW представлен на рис. 1.
Цифрован Обработка Сигналов
Файл Редактор Обработка Запись Параметры Ог.но Подсказка
шш
'■W ellcome to есЫ.еке
Free Size: [СЛИБ7 МБ
Ä
Рис.1. Главное меню программы EDSW «Генератор сигналов» формирует тестовые сигналы с хорошо известными математическими моделями, например:
1) гармонический сигнал s(t) = Aq coS(sin)(wt + pjwt + рн),
2) периодическая последовательность прямоугольных импульсов с амплитудой U, длительностью t и периодом повторения T s(t) = s(t ± nT);
3) гауссовский шум - случайный процесс с плотностью вероятности распределения амплитуд
p(u) =
1
- u
/(2s2)
где s - дисперсия, s - СКО (рис.2).
Рис.2. Гауссовский шум и плотность вероятности распределения амплитуд В качестве примера рассмотрим тему «Ряды и интеграл Фурье», которая раскрывается в лабораторной работе «Исследование спектров сигналов. Гармонический синтез сигналов». Цель работы - развить у студентов спектрально-временное представление в системе «амплитуда-частота» (спектроанализатор) и. «амплитуда-время» (осциллограф). Здесь главное -показать, что периодические сигналы полностью описываются рядом Фурье по гармоникам основной частоты f = 1/T, где T - период повторения. А непериодические (однократные) сигналы представляются интегралом Фурье. Так, для импульсный периодической последовательности униполярных импульсов с амплитудой U = 30000 у.е., t = 1 мсек, T = 5 мсек, получаем дискретный (линейчатый) спектр с гармониками на частотах 200, 400, 600, 800, 1200 Гц, ..., а частота 1000 Гц «зануляется», т. к. при скважности T/1 = 5 пятая, 10-ая, 15-ая и другие кратные 5 гармоники равны 0 (рис.3).
40000 30000 20000 10000
4.053 4.054
4.053 4.059
4.062 4.064 4.Q65 4.066 4.067
4.070 4.071
1
i 1 г,..;
щ t j ■ * ШШ 1. 1 .j. : i : A
268.555 537.109 8G5.664 1074.21 Э 1342.773 1611.328 187Э.883 2148.438 2416.ЭЭ2 2Б85.547 2954.102 3222.656 3491.211 3759.766 4028.320 4296.875 4565.430 4833
Рис.3. Импульсная периодическая последовательность и амплитудно-частотный спектр
При увеличении периода повторения T (стремлении T к бесконечности, например, взяв T = 1 сек) получаем сплошной спектр, описываемый функцией «координатный синус» Sinex = (Sinx)/x (рис.4). Здесь первый нуль находится на частоте 1000 Гц = 1/1, второй нуль - 2000 Гц = 2/1 и т.д.
Рис. 4. Импульсная последовательность (T ® ¥ ) и амплитудно-частотный спектр Теперь рассмотрим гармонический синтез сигналов. Пусть необходимо синтезировать пилообразное напряжение:
S(t) = 2A0 (sin Wt -1/2 sin 2Wt +1/3 sin3 Wt -1/4 sin 4Wt +1/5 sin 5 Wt -1/6 sin 6Wt +...) Зададим 2Ao = 10000 = 10000 у.е., W = 200, 400, 600 Гц и т.д. В результате суммирования гармонических колебаний с кратными частотами получаем периодический сигнал, аппроксимирующий пилообразное напряжение. Его спектр показывает пять использованных при синтезе гармоник, амплитуды которых убывают обратно пропорционально номеру гармоники (рис.5).
Описания лабораторных работ в методических пособиях [1,2] включают краткую теорию, цель работы, лабораторное задание, программу работ, требования к оформлению отчета и контрольные вопросы.
1 0Ш 266.555 537.109 №664. 1074.219 1342.773 1611.328 1579:883 2149.438.2416.992 2685.547 2954.102 3222.656 3431.211 3759.768 4026,320 4296.875 4565.430 4933.994.51.02.539 5371.094
Рис.5. Пилообразное напряжение и его амплитудно-частотный спектр 2. Прямые и обратные задачи радиоимпедансного зондирования
В лаборатории геоэлектромагнетизма Отдела физических проблем при Президиуме Бурятского научного центра СО РАН развита теория, разработаны численные методы и создано программное обеспечение решения прямых и обратных задач радиоимпедансного зондирования для вертикально-неоднородных (дискретно-слоистых и градиентных) сред с учетом эффектов частотной дисперсии основных электрических параметров земной коры.
Создан пакет прикладных программ «Импеданс» [3]. Он предназначен для решения прямой и обратной задач распространения электромагнитной волны в слоистой среде:
1) расчет поверхностного импеданса по известному геоэлектрическому разрезу;
2) восстановление параметров геоэлектрического разреза слоистой среды (sj,sj,hj)
по измерениям ее импеданса.
Методика расчетов прямой задачи основана на понятии приведенного поверхностного импеданса плоской волны на поверхности горизонтально-слоистой структуры, с помощью которого можно учесть влияние подстилающей среды на распространение радиоволн:
5= ЕТ/(ИТ2 о ),
где Ет, Ит - горизонтальные взаимно перпендикулярные компоненты электрического и магнитного полей на границе раздела «воздух-земля»; 2о = ^ то /80 = 377 Ом - характеристический импеданс свободного пространства.
Для решения обратной задачи электромагнитного зондирования для слоистой полупроводящей среды, удовлетворяющей импедансным граничным условиям, разработан алгоритм метода регуляризации на основе гипотетического разреза. Сглаживающий функционал
ь 2
А.Н. Тихонова имеет вид: Ма [§■ ] = I[§■ ] + аО[^], где I] = -3[ - функционал не-
I=1 Э
вязки;
N ,2 N п 2 N-1 п 2
n=1
g ] =Z Pn en+Z qn sn-s" +Z
n=1
■n -s"
n=1
hn hn
- стабилизирующий функционал;
g = (s,s,h) - параметры многослойной полупроводящей среды; a - параметр регуляризации; (e ",s",h") - гипотетическое распределение параметров разреза; pn,qn,rn - веса соответствующих параметров; 5\ - экспериментальные значения импеданса на l -ой частоте.
Минимизация рассматриваемого функционала производилась с помощью метода Ньютона-Канторовича (метода линеаризации). Осуществляя линеаризацию оператора прямой задачи
A в окрестности имеющегося приближения g(n), получили формулу:
Dg = (BTБ+al)~l ■ (БТдэ + alg° ), где B - матрица Якоби для A с элементами Bj = (dAj / dgj) ^. Для решения системы линейных уравнений на каждом шаге итерационного процесса использовался метод псевдообращения. С использованием сингулярного разложения B(gn )= ULVT получена окончательная формула для решения:
Dg = Л2 + alr)-1 U Tдэ + (Vn1V T +alr ag".
Критерием выбора «наилучшего» геоэлектрического разреза служит минимум сред-неквадратического отклонения экспериментальных значений импеданса от рассчитанных для модели n -слойной среды.
Программы пакета созданы с помощью объектно-ориентированного языка программирования Visual Basic for Applications и функционируют в среде Excel. Дружественный графический интерфейс с системой меню облегчает работу пользователям с различным уровнем компьютерной подготовки. Поэтапно осуществляется ввод, просмотр, обработка измеренных данных, вывод результатов в наглядной форме в виде таблиц и диаграмм, сохранение результатов в файлы. Интерактивный режим позволяет контролировать процесс интерпретации на каждом этапе работы. При этом существенно снижается погрешность интерпретации и повышается надежность получаемых результатов. Все программы прошли тестирование, апробированы на модельных примерах. Проведена оценка погрешности определения параметров геоэлектрического разреза. Кроме указанных выше программ, пользователю предоставлены также следующие возможности:
- статистические расчеты импеданса в заданном диапазоне частот или получение статистических оценок модуля и фазы поверхностного импеданса однородной в геологическом отношении области (комплекса горных пород);
- учет частотной дисперсии электрических параметров среды;
n
- расчеты для градиентных полупроводящих сред путем аппроксимации большим числом плоскопараллельных однородных слоев с электрическими параметрами, соответствующими профилю неоднородного слоя (кусочно-градиентная модель наилучшим образом отражает реальное строение многолетнемерзлого грунта);
- получение пространственного распределения параметров геоэлектрического разреза при радиоимпедансном зондировании при площадной съемке или на профиле.
В лабораторном практикуме для студентов предложены различные модели природных слоистых сред типа «лед-соленая (морская) вода», «лед-пресная вода», «лес на горном хребте» и другие. Особый интерес представляет «вечная мерзлота», занимающая около 25% площади суши Земного шара. Одной из особенностей мерзлых пород является наличие значительной частотной дисперсии электрических свойств малопроводящих мерзлых рыхлых отложений. Слоистые геоэлектрические структуры, встречающиеся на мерзлоте, дают сильный эффект осцилляции модуля || и фазы ср§ импеданса в индуктивной и емкостной областях
(рис. 6). Учет частотной дисперсии электрических свойств мерзлых рыхлых отложений приводит к почти монотонной частотной зависимости импеданса. На рис.6 принято, что в диапа-2 6
зоне частот 10 ^ 10 Гц сопротивление Р2 геоэлектрического разреза (ГЭР) изменяется по формуле ^ Р2 = 4,8 - 0,166 ^ / .
40 ч 20 0 --20 --40 --60 --80
фа
1д р2=4,8-0,166 1д {
1
0,1
0,01 а
0,001 .
0.0001
\
с дисперсией р2
р1=500 Ом-м, с-| = 10, м; 1\.=63000 Ом-м, «5-10, 160 м; рз=500 ОМ-М, 1&5И0, 1Пз=оо
т кГц
0,01 0,3 1 10 100 -1000
Рис. 6. Осцилляции модуля и фазы импеданса для модели мерзлоты.
Предложенный регуляризирующий алгоритм продемонстрировал свою эффективность в процессе обработки и анализа экспериментальных данных и позволяет восстанавливать структуру достаточно сложных геоэлектрических разрезов. Были проинтерпретированы результаты зондирований в диапазоне 17-529 кГц в 23 типичных пунктах, характеризующих мерзлые осадочные и кристаллические породы, а также акваторию реки Колыма и тундровые озера (р-н п. Черский). Установлено, что электрическое сопротивление рыхлых мерзлых пород во много раз меньше известных по ВЭЗ значений, а зоны с аномально низким сопротивлением (10-60 Ом.м) имеют большое территориальное распространение. Полученные результаты важны при прогнозировании параметров ГЭР криолитозоны.
Создана программа графической интерпретации радиоимпедансных зондирований [4], предназначенная для определения параметров (1 ^ 4)-слойного геоэлектрического разреза (диэлектрических проницаемостей, удельных сопротивлений и толщин слоев) по частотной зависимости поверхностного импеданса. В ней реализован графический вариант метода подбора. Искомые параметры разреза выбираются из критерия наилучшего совмещения теоретических кривых и экспериментальных значений импеданса. Достоинства: интерактивный режим, наглядность, быстрота обработки. Программа позволяет выбрать модель среды, выявить эквивалентные разрезы, получить начальные параметры разреза, которые затем можно уточнить методом регуляризации. На рис. 7 представлен вид рабочего окна программы для 3-слойного разреза.
В целом лабораторные работы проводились с интересом, хорошей адаптацией и усвоением изучаемого материала. Разработанный цифровой программный комплекс характеризуется многофункциональностью и универсальностью.
Рис. 7. Вид рабочего окна программы графической интерпретации данных радиоимпедансных зондирований для 3-слойного разреза.
Заключение
Разработан комплекс лабораторных работ для моделирования спектрально-корреляционных, статистических и электродинамических характеристик физико-технических систем, имеющий широкие функциональные возможности. Цифровой программный комплекс может быть использован при УИРС, выполнении курсовых, дипломных работ и проектов.
Комплекс может быть полезен при изучении курсов «Основы радиоэлектроники», «Статистическая радиофизика», «Радиотехнические цепи и сигналы». В связи с развитием и использованием дистанционного образования открываются новые направления в применении данного электронного программно-методического лабораторного практикума как в целях обучения, так и переподготовки кадров в различных учебных заведениях Восточной Сибири.
Библиография
1. Нагуслаева И.Б. Статистические характеристики случайных электромагнитных процессов: Учеб.-метод. пособие. Часть I./ Ответ.ред. Ю.Б.Башкуев. Улан-Удэ: Изд-во Бурятского госуниверситета, 2005. 42с.
2. Нагуслаева И.Б. Спектральные и корреляционные характеристики электромагнитных процессов: Учеб.-метод. пособие. Часть II./ Ответ.ред. Ю.Б.Башкуев. Улан-Удэ: Изд-во Бурятского госуниверситета, 2005. 38с.
3. Ангархаева Л.Х. Свидетельство об официальной регистрации программы для ЭВМ № 2002610893. Пакет программ «Импеданс» для решения задач радиоимпедансного зондирования. М.: РОСПАТЕНТ, 06.06.2002.
4. Ангархаева Л.Х. Свидетельство об официальной регистрации программы для ЭВМ № 2006610089. Графическая интерпретация радиоимпедансных зондирований. М.: РОСПАТЕНТ, 10.01.2006.
1. Naguslaeva I. B. Statistical characteristics of random electromagnetic processes: methodical textbook, part 1./Editor-in-chief Y.B. Bashkuev. Ulan-Ude: Buryat State University publishing, 2005, p.42
2. Naguslaeva I. B. Spectral and correlating characteristics of electromagnetic processes: methodical textbook, part 2./Editor-in-chief Y.B. Bashkuev. Ulan-Ude: Buryat State University publishing, 2005, p.38
3. Angarhaeva L. H. The official registration of the computer software program certificate 2002610893. The Impedance software package for solving the radio-impedance probe problems. M.: ROSPATENT, 06.06.2002.
4. Angarhaeva L. H. The official registration of the computer software program certificate 2006610089. Graphic interpretation of the radio-impedance probe. М: ROSPATENT, 10.01.2006.
