WWW.NET.KNIGI-X.RU
БЕСПЛАТНАЯ  ИНТЕРНЕТ  БИБЛИОТЕКА - Интернет ресурсы
 

«ПРИКЛАДНАЯ МАТЕМАТИКА ПОМЕХОУСТОЙЧИВОСТЬ СИСТЕМЫ СВЯЗИ С БИФАЗОВОЙ МАНИПУЛЯЦИЕЙ А.Т.Гаврилин Нижегородский госуниверситет Для определения ...»

Труды Научной конференции по радиофизике, ННГУ, 2004

ПРИКЛАДНАЯ МАТЕМАТИКА

ПОМЕХОУСТОЙЧИВОСТЬ СИСТЕМЫ СВЯЗИ

С БИФАЗОВОЙ МАНИПУЛЯЦИЕЙ

А.Т.Гаврилин

Нижегородский госуниверситет

Для определения вероятности ошибочного приема в системах передачи дискретных сообщений необходимо знать вероятностное распределение статистики (функционала от принятого сигнала), используемой в решающей процедуре. В рассматриваемой системе связи [1] при поэлементном приеме такой статистикой является аргумент (1, 2) нормированного выборочного биспектра &&& ) ) ) n n n & & & b = (mod b ) eiY = 1 2 3 = 1 + 1 1 + 2 1 + 3, (1) s3 s1 s2 s3 s1 s2 &&& & & & где (1, 2, 3), (n1, n2, n3), (s1, s2, s3) с надстрочными символами – это комплексные амплитуды (при частотах 1, 2 опорного дуплета и частоте 3= 1+2 информативной компоненты манипулируемого бота) соответственно: принятого сигнала (t), гауссовского шума n(t) (с дисперсиями n2), составляющих бота sk = Akexp[ik] (k=1,2,3); 1+2–3=. Здесь – информационный символ, постоянный на тактовом интервале. Далее для определенности он полагается равным нулю.

Попытка прямого вычисления распределения статистики (1) приводит к вероятностной “задаче трех тел”: необходимости находить распределение кубичного функционала от компонент гауссовского случайного вектора. Выходом из положения может служить аппроксимация статистики (1) внешне более сложной статистикой, основывающаяся на приближенном стохастическом равенстве 1 + z e z z, (2) где – положительная константа, выбираемая, однако, не по формуле Тейлора, а из условия совпадения дисперсий левой и правой частей равенства (2).



Этот прием, с одной стороны, напоминает, а с другой, противоположен известному методу статистической линеаризации [2] и может быть назван статистическим потенцированием. В надпороговом режиме (|z|2 0,5) величина близка к 1/4. Обозначив отношения комплексных амплитуд nk sk числами zk = xk + iyk, в силу (2) имеем z k z k x k x k + y k i y k (1 2 x k ) ) b e k =1 = e k =1 e k =1.

Труды Научной конференции по радиофизике, ННГУ, 2004 Таким образом, оценка бифазы бота, осуществляемая биспектральном приемником, м

–  –  –

где xk, yk – условно (при фиксированных фазах k) гауссовские случайные величины с дисперсиями (k Ak)2. Наличие в статистике (3) только парных произведений случайных величин позволяет выразить ее характеристическую функцию в достаточно простом виде

–  –  –

Сравнение вероятностей ошибок классификации, вычисленных, с одной стороны, по данной формуле с использованием (4) и, с другой, при помощи прямого компьютерного моделирования, показало хорошее совпадение для широкого диапазона отношений сигнал/шум (k Ak, k=1,2,3). Вместе с тем вряд ли стоит ожидать столь же точной аппроксимации ошибок посредством статистического потенцирования в случае, когда передаваемая информация закладывается в модуль биспектра передаваемого сигнала. Это сомнение продиктовано тем обстоятельством, что при энергетическом информационном параметре “хвосты” распределений мгновенных значений процессов более продуктивны в части воздействия на решающую схему.

[1] Гаврилин А.T. //В кн.: Тр. 6-й научн. конф. по радиофизике. 7 мая 2002 г. /Ред.

А.В.Якимов. -Н.Новгород: ТАЛАМ, 2002, с.280.

[2] Свешников А.А. Прикладные методы теории случайных функций. -М.: Судпромгиз, 1963, 512 с.





–  –  –

Как видно из (4), для отыскания x необходимо знать не саму функцию F(t), а ее лаплас-образ. Среднее в формуле (5) можно найти через так называемую производящую функцию моментов аддитивного функционала марковского процесса x0(t) [1], которая для нашего случая имеет вид следующего условного среднего

–  –  –

[1] Теория связи /Пер. с англ. под ред. Б.Р.Левина. -М.: Связь, 1972, с.89.

[2] Fauve S., Heslot F. //Phys. Lett. A 1983. V.97, №1-2. P.5.

[3] McNamara B., Wiesenfeld K., Roy R. //Phys. Rev. Lett. 1988. V.60, №25. P.2626.

Труды Научной конференции по радиофизике, ННГУ, 2004

ЭФФЕКТ УСКОРЕНИЯ ДИФФУЗИИ В СЛУЧАЙНО ПЕРЕКЛЮЧАЮЩЕМСЯ

ПИЛООБРАЗНОМ ПОТЕНЦИАЛЕ

А.А.Дубков Нижегородский госуниверситет Ранее в работах [1,2] изучалась диффузия броуновских частиц в симметричном периодическом потенциале, модулируемом гауссовым белым шумом, и было обнаружено ускорение диффузии по сравнению со случаем свободной диффузии при любой форме потенциального профиля. В настоящем докладе рассматривается более сложная задача, связанная с определением условий возникновения данного эффекта при модуляции потенциала марковским дихотомическим шумом.

Рассмотрим одномерное передемпфированное броуновское движение в случайно переключающемся симметричном периодическом потенциале U(x), описываемое ланжевеновским уравнением dx dU ( x) (t ) + (t ).

= (1) dt dx Здесь: x(t) – координата броуновской частицы, (t) –гауссов белый шум с интенсивностью 2D и нулевым средним, (t) – марковский дихотомический шум со значениями ±1 и средней частотой переключений. Как видно из рис.1, потенциал переключается между двумя противоположными конфигурациями в случайные моменты времени, причем в “перевернутой” конфигурации максимумы и минимумы потенциала меняются местами. Будем далее полагать выполненным условие суперсимметрии потенциала с пространственным периодом L [3] E U (x) = U (x L / 2), (2) Рис. 1 при котором отсутствует направленный поток частиц (рэтчет-эффект).

Пусть в начальный момент времени реализуется верхняя конфигурация потенциала на рис.1, а броуновская частица находится в потенциальной яме в начале координат. Тогда процесс диффузии можно грубо описать в виде последовательных скачкообразных переходов частицы из точек минимума потенциала xm=mL в соседние минимумы xm±1. В этой скачкообразной модели диффузии время перехода из xm в xm±1 представляет собой время первого достижения частицей, стартующей из точки x=xm, поглощающих границ x=xm±1, а длина скачка является случайной величиной, принимающей с равной вероятностью значения +L и –L.

Описанная математическая модель случайного процесса x(t) позволяет просто и в то же время строго рассчитать асимптотическое значение эффективного коэффициента диффузии Dэфф. Несложные вычисления в соответствии с определением Dэфф дают Труды Научной конференции по радиофизике, ННГУ, 2004

–  –  –

с граничными условиями (0)=0, (L)=0.

Найти решение уравнения (5) для произвольного потенциала U(x) не представляется возможным, поэтому далее подробно рассматривался пилообразный периодический потенциальный профиль, удовлетворяющий условию суперсимметрии (2), 0 x L/2 kx, U ( x) =, (6) k ( L x), L / 2 x L где k=2E/L. В результате подстановки (6) в (4), (5) и довольно громоздких расчетов было получено точное выражение для эффективного коэффициента диффузии, справедливое для любых значений высоты потенциальных барьеров E, интенсивности шума D и средней частоты переключений.

Рис. 2 На рис.2 изображена область ускоренной (по сравнению со случаем U(x)=0) диффузии на плоскости двух безразмерных параметров задачи: относительной высоты потенциальных барьеров =E/D и относительного темпа переключений потенциала =L2/(2D). Как видно из рис.2, ускорение наблюдается лишь для достаточно высоких барьеров, когда броуновские частицы при очередном переключении достаточно быстро скатываются в ближайшие минимумы потенциала, и для достаточно частых переключений.

Работа выполнена при поддержке РФФИ (грант № 02-02-17517), гранта Ведущие Научные Школы НШ-1729.2003.2 и научной программы “Университеты России” (проект УР.01.01.020).

–  –  –

[1] Малахов А.Н. //Письма в ЖТФ. 1998. Т.24, вып.21. С.9.

[2] Дубков А.А. //Письма в ЖТФ. 2003. Т.29, вып.3. С.18.

[3] Reimann P. //Phys. Rev. Lett. 2001. V.86, №22. P.4992.

РАСЧЕТ ЧИСЛА СИЛЬНОСВЯЗНЫХ ФЕЙНМАНОВСКИХ ДИАГРАММ

МЕТОДОМ ПРОИЗВОДЯЩИХ ФУНКЦИЙ

А.А.Дубков, В.Ю.Семенов Нижегородский госуниверситет Фейнмановские диаграммы находят широкое применение во многих областях физики: в квантовой теории поля, в теории распространения волн в турбулентных средах, в физике высоких энергий, в эйнштейновской теории гравитации. Число сильносвязных фейнмановских диаграмм может служить оценкой скорости сходимости рядов теории возмущений.

В настоящей работе рассматриваются только одномерные диаграммы, получаемые всевозможными попарными замыканиями 2n точек на прямой. Их можно разбить на два вида: сильносвязные и слабосвязные. На рис.1 представлена слабосвязная диаграмма. На ней между точками 2 и 3 можно “разломить” прямую, не повредив связей. На рис.2, в свою очередь, приведена сильносвязная диаграмма.

Здесь нет места, где можно было бы “разломить” прямую, не повредив связей.

Рис.1 Рис.2

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

Если обозначить общее число фейнмановских диаграмм для 2n точек на прямой через V2n, а число сильносвязных через S2n, то для первого положения границы получим согласно правилам комбинаторики S2V2n-2 слабосвязных диаграмм, для следующего положения – S4V2n-4 диаграмм и т.д. Для последнего мыслимого положения границы разрыва (перед предпоследней точкой на прямой) имеем S2n-2V2 слабосвязных диаграмм. Общее число слабосвязных диаграмм находим суммированием всех вариантов.

Учитывая то, что сумма сильносвязных и слабосвязных диаграмм составляет общее число фейнмановских диаграмм V2n=(2n-1)!!, приходим окончательно к следующему рекуррентному соотношению для числа сильносвязных диаграмм S2n n 1 S 2 n = (2n 1)!! S 2 k [2(n k ) 1] !!, n 2. (1) k =1 Труды Научной конференции по радиофизике, ННГУ, 2004 Для решения соотношения (1) введем производящую функцию вида

–  –  –

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

–  –  –

[1] Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. -М.: Лань, 2003, 832с.

[2] Справочник по специальным функциям. -М.: Наука, 1979, 832с.

–  –  –

где: j – номер гидроопоры; i – номер частоты возбуждения гидроопоры, при i=0,1,2…n; индекс н перед индексами ji в обозначении функции H(s) указывает, что она нормирована; a2 – нормировочный коэффициент, равный a2 по абсолютной величине и является безразмерной величиной; a2 = m – масса нагрузки; a1i = bzi – гидравлический коэффициент трения гидроопоры; a0ji = cji – динамическая жесткость данного номера гидроопоры на данной частоте возбуждения; b0 = Fz – сила воздействия на гидроопору; b1i = Bzi = z0iвibzi – многомерный параметр гидроопоры, входящий в определение гидравлической составляющей гидроопоры; pi Sпор= z0iвibzi cos вit – сила, действующая на жидкость при pi – перепаде давления между рабочей и компенсационными камерами гидроопоры и Sпор – поршневая площадь опорной платы.

Нормированные системные функции H(s)нji вида (1) полностью характеризуют собственное движение элементов гидроопоры и определяют ее динамическую жесткость и, являясь рациональными функциями, полностью определяются корнями полиномов числителя и знаменателя – нулями Szi и полюсами Spi. Расположение Труды Научной конференции по радиофизике, ННГУ, 2004 нулей и полюсов на комплексной плоскости S обуславливается как типом системной функции H(s)нji, так и элементами, входящими в состав гидроопоры. Положение нулей и полюсов содержит значительную информацию – как качественного, так и количественного характера – относительно свойств гидравлической виброопоры.

Тип используемых элементов гидравлической виброопоры, ее конструкция или ее динамическая жесткость ограничивают область комплексной S-плоскости, в которой могут располагаться нули и полюса функций H(s)нji. Расположение полюсов и нулей функций H(s)нji на частотах возбуждения гидроопоры характеризуется “диаграммами полюсов-нулей”.

Ввиду того, что расположение полюсов и нулей функций H(s)нji гидропоры на каждой частоте возбуждения вi различно, то для настройки гидроопоры на заданную частоту вибрации вi необходимы следующие требования (условия):

- для перестройки гидроопоры на заданные резонансные частоты требуется учитывать определенную экспериментом зависимость динамической жесткости от частоты возбуждения;

- необходимо, чтобы расположение полюсов и нулей на заданной частоте вибрации вi, до настройки гидроопоры, перешло при ее настройке на расположение полюсов и нулей, соответствующее собственной частоте гидроопоры в0 при ее минимальной динамической жесткости C0;

- динамическая жесткость настроенной гидроопоры на заданной частоте вi должна соответствовать динамической жесткости C0 ненастроенной гидроопоры с собственной частотой в0.

Работа выполняется при финансовой поддержке гранта РФФИ 03-02-16924 и ФПЦ "Интеграция" гранта Президента РФ для поддержки ведущих научных школ.

[1] Ложкин Ф.В. Исследование гидроопор для виброзащиты транспортных средств.

Автореферат диссертации на соискание ученой степени кандидата технических наук

. -Н. Новгород, 2002.

[2] Макаров И.М., Менский Б.М. Линейные автоматические системы (элементы теории, методы расчета и справочный материал). -М.: Машиностроение, 1982, 504с.

[3] Справочное пособие по теории систем автоматического регулирования и управления. /Под общей ред. Е.А. Сашковского. -Минск: Вышэйш. школа, 1973, 584с.

[4] Сибирт У.М. Цепи, сигналы, системы. Ч.1. -М.: Мир, 1988, 336с.

[5] Введение в теорию и расчет активных фильтров. -М.: Радио и связь, 1984, 384с.

[6] Гордеев Б.А., Охулков С.Н. Прикладная механика и технологии машиностроения: сборник научных трудов /Ред. В.И. Ерофеева, С.И. Смирнова и Г.К. Сорокина. -Н. Новгород: Интелсервис, 2003, вып.2(6), с.112.

Труды Научной конференции по радиофизике, ННГУ, 2004

ВЛИЯНИЕ ВНЕШНЕГО МАГНИТНОГО ПОЛЯ НА ПЛОТНОСТЬ ПОТОКА

ЭНЕРГИИ В МАГНИТОРЕОЛОГИЧЕСКОМ ТРАНСФОРМАТОРЕ

Б.А.Гордеев, П.Н.Морозов Нижегородский филиал института машиноведения РАН Сложность и нелинейность поведения электро- и магнитореологических жидкостей (ЭРЖ и МРЖ), а также относительно высокие напряженности полей, необходимых для управления вязкостью этих составов сыграли свою роль в том, что на протяжении почти полувека это явление не получило распространения на практике.

Только в последнее время, используя новые технологии, начинают находить применение ЭРЖ и МРЖ в системах гашения вибраций и демпфирования ударов [1,2].

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

Одной из возможных модернизаций гидроопоры является применение электрореологических жидкостей, заполняющих рабочую и компенсационную камеры, впервые предложенная в 1989 г. в Нф ИМАШ РАН [1].

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

–  –  –

Из (8) следует, что плотность потока энергии зависит от Н0 как функция пятой степени. Следует отметить, что даже в простейшем случае совпадения частот внешних вибрационных и электромагнитных полей плотность переноса энергии существенно нелинейна.

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

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

Работа выполняется при финансовой поддержке Гранта РФФИ 03-02-16924 и ФЦП “Интеграция” гранта Президента РФ для поддержки ведущих научных школ.

[1] А.С. № 17799843 (СССР) Виброизолирующее устройство. //Гордеев Б.А., Образцов Д.И., Юдин В.А., Поташев О.А. кл. F 16 F 6/00. 07.12.92. Бюл. № 45.

[2] Гордеев Б.А., Ерофеев В.И., Синев А.В. //Проблемы машиностроения и надежности машин. 2003. №6. С.22.

Труды Научной конференции по радиофизике, ННГУ, 2004

ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ ДВУМЕРНОЙ УПРУГОЙ СИСТЕМЫ

С ДВИЖУЩЕЙСЯ НАГРУЗКОЙ

О.А.Волкова, Е.Е.Лисенкова Нижегородский госуниверситет Рассматривается упругая система, представляющая собой пластину модели Кирхгофа-Лява на упругом основании. Вдоль нее безотрывно по некоторому закону x=l(y,t) движется объект (нагрузка), совершающий изгибно-крутильные колебания.

Предполагается, что на объект действует источник сил, изменяющийся по гармоническому закону P(y,t)=P0·eit·cos kyy.

Исходя из вариационного принципа Гамильтона-Остроградского, на основе разработанной в [1,2] процедуры постановок краевых задач о согласованном динамическом поведении упругих систем с движущимися нагрузками получим, что изгибные колебания пластины определяются решением уравнения U tt + D(U xxxx + 2U xxyy + U yyyy ) + hU = 0, (1) удовлетворяющим при x=l(y,t) условиям непрерывности

–  –  –

0ltt + E0 Jl yyyy = [ F ] + ql, (6) U (y,t) - изгибные, (y,t) – крутильные колебания нагрузки, l(y,t) – закон движения нагрузки, [N]=N(l(y,t)+0,y,t)-N(l(y,t)-0,y,t), F – сила давления волн, ql – плотность внешней обобщенной силы, – поверхностная плотность, 0 – погонная плотность, h – толщина пластины, J0, J – полярный момент инерции и момент инерции попеТруды Научной конференции по радиофизике, ННГУ, 2004 речного сечения относительно осей, перпендикулярных оси стержня,, h’ – коэффициенты упругости нагрузки и основания пластины соответственно.

Ограничиваясь случаем равномерного движения объекта x=Vt и рассматривая лишь вынужденные установившиеся колебания в данной системе, численно и аналитически решались задачи кинематики и динамики излучения волн.

Найдены частоты и волновые вектора излучения волн источником в диапазоне частот 0100Гц, движущимся по стальной пластине толщины 1,5мм, ширины 1м, коэффициент жесткости упругого основания которой 5107 мн3. Определены критические скорости движения нагрузки, при переходе через которые качественно меняется картина волнообразования. На плоскости параметров (,V) им соответствуют кривые (рис.), разбивающие ее на области, в которых характер возбуждения волн различен.

Так если частота внешней силы не превосходит наименьшую частоту возбуждаемых в пластине волн (), а скорость движения объекта меньше первой критической (VV1*), то движущаяся нагрузка волн не излучает. В сопутствующей системе координат прогиб пластины экспоненциально спадает по обе стороны от нагрузки. С увеличением скорости движения (V1*VV2*) наряду с собственным полем, представляющим собой спадающие по экспоненте осцилляции, излучается по одной волне справа и слева от объекта. При переходе через вторую критическую скорость, (VV2*) источник собственного поля не создает, но излучаются слева и справа по две волны. Исследован был также случай при.

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

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

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 03-01-00664).

[1] Болдин В.П., Весницкий А.И. //Машиноведение. 1989. №1. С.70.

[2] Болдин В.П., Маланов С.Б., Уткин Г.А. //ПММ. 1992. Т.56, вып.1. С.34.

–  –  –

МАТЕМАТИЧЕСКОЕ ПРОГНОЗИРОВАНИЕ СТРУКТУРЫ ВАЛОВОГО

РЕГИОНАЛЬНОГО ПРОДУКТА

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

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

В настоящей работе решается следующая задача. Известна динамика и структура валового регионального продукта (ВРП) Нижегородской области (без налогов) в период с 1994 г. по 2000 г. (в млн. руб.) (до 1998 г. – в млрд. руб.) (табл.1) [1].

Зная структуру ВРП (без налогов) в период с 1994 г. по 1998 г., оценим, как она изменится в последующие два года, и сравним прогнозные оценки с имеющимися реальными данными.

–  –  –

Рассматривая каждую из представленных в табл. 1 отраслей как инвестиционную альтернативу и проводя сравнение альтернатив согласно методике [2], использующей критерии стохастического доминирования и элементы теории голосований, получаем следующие предпочтения инвестора: 1 f 2 f 3 f 4 f 5. Это означает, что в ближайшие год-два капитал, задействованный во всей экономике Нижегородской области, перераспределится таким образом, что наибольший приток инвестиционных средств будет наблюдаться в промышленности, а наибольший отток средств – в сфере нерыночных услуг. При этом уже в 1999 г. и 2000 г. в наиболее перспекТруды Научной конференции по радиофизике, ННГУ, 2004 тивных с точки зрения инвестора отраслях будет наблюдаться отдача от инвестиций, т.к. инвестор уже в 1998 г. реагирует на изменения в рыночной конъюнктуре.

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

и 2000 г. Для этого на основе реальных данных 1999 г. и 2000 г. рассчитывается удельный вес валовой добавленной стоимости (ВДС) каждой отрасли в общей величине ВРП (без налогов), затем производим подобные расчеты для периода с 1994 г. по 1998 г. После этого, используя метод наименьших квадратов, можем установить зависимость удельного веса ВДС каждой отрасли от порядкового номера года в виде уравнения линейной регрессии y= +x. Далее сравнивается реальный удельный вес ВДС каждой отрасли с прогнозным, т.е. рассчитанным по уравнению регрессии (табл.2).

–  –  –

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

Прогноз на период в два года тоже достаточно точен. Искажения оценок проявляются лишь по отраслям 3 и 4, причем искажения незначительны. Поэтому в целях повышения точности прогноза перелива капитала между отраслями необходимо каждый период (в нашем случае каждый год) проводить анализ, представленный в работе [2] и в настоящей статье.

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

[1] Статистические материалы “Нижегородская область 1991-2001 гг.” (электронный вариант) -Н.Новгород: Нижегородский областной комитет государственной статистики, 2002, с.105.

[2] Данилова Т.Н., Кошелев Е.В. //Ученые записки. Т.6. -Н.Новгород: Изд-во ВВАГС, 2003, с.74.

Труды Научной конференции по радиофизике, ННГУ, 2004

БИФУРКАЦИЯ РОЖДЕНИЯ ИНВАРИАНТНОГО ТОРА В СЛУЧАЕ

ГОМОКЛИНИЧЕСКОГО КАСАНИЯ МНОГООБРАЗИЙ СЕДЛО-ФОКУСА

И.И.Овсянников Нижегородский госуниверситет В данной работе рассматривается трехмерный диффеоморфизм, имеющий неподвижную точку типа седло-фокус, устойчивое и неустойчивое многообразия которой касаются квадратичным образом. При этом, якобиан отображения J=2 (где и – абсолютные величины устойчивых и неустойчивого мультипликаторов) равен 1. Целью работы является изучение основных бифуркаций однообходных периодических траекторий в системах, близких к исходной особое внимание следует уделить бифуркации рождения инвариантного тора из периодической траектории.

Наличие негрубой гомоклинической траектории влечет за собой нетривиальную динамику, в том числе, возникновение странных аттракторов, квазиаттракторов и др. [1,2]. Переход через границу J=1 означает смену глобальной устойчивости (сжатия) на неустойчивость и, возможно, рождение инвариантного тора. В физической системе рождение тора – это срыв с периодического режима на биения.

Исследование проводится в трехпараметрическом семействе диффеоморфизмов. Один из параметров отвечает за расщепление многообразий седла, второй - за отклонение величины J от 1, а третий – это аргумент устойчивых мультипликаторов.

С помощью перемасштабирования координат и параметров отображение первого возвращения можно привести к следующему виду:

X = Y, Y = Z, (1) Z = M BY CX Z 2.

Здесь M, B и C – перемасштабированные исходные параметры.

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

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

–  –  –

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

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

Работа поддержана грантами “Ведущие научные школы” №НШ-838.2003.2;

РФФИ №02-01-00273.

[1] Gonchenko S.V., Shilnikov L.P., Turaev D.V. //Computers Math. Applic. 1997. V.34, №2-4. P.195.

[2] Гаврилов Н.К., Шильников Л.П. //Матем. сб. 1972. Т.88, №4. C.475.

–  –  –

Так как exp{t}1, то LT(,z) – оператор сжатия. По теореме о сжатых отображениях для случайных операторов [2] уравнение LT(,z)=z имеет единственное решение, и x(t)= Lt(,z*) – решение исходной СКЗ.

[1] Репин О.Н., Саичев А.И. //Дифференциальные и интегральные уравнения. 1986.

С.82.

[2] Кириллов П.В. Случайные уравнения. -Кишинев: Штиинца, 1982.

ПРОМЕЖУТОЧНАЯ АНОМАЛЬНАЯ ДИФФУЗИЯ.

МНОГОМЕРНЫЙ СЛУЧАЙ

А.И.Саичев, С.Г.Уткин Нижегородский госуниверситет Аномальной диффузии, отличие которой от линейной заключается в нелинейном росте среднего квадрата диффузионного процесса со временем, в последние годы проявляется особое внимание [1]. Явления аномальной диффузии затрагивают широкий спектр физических проблем, таких как перенос заряда в аморфных полупроводниках [2-4], геометрия фракталов [5], квантовая оптика [6], турбулентная диффузия [7], хаотическая динамика гамильтоновых систем [8] и т.д.

Чаще всего зависимость среднего квадрата удовлетворяет степенному закону X2(t) ~ t, возникающему из-за нарушения закона больших чисел (ЗБЧ) и (или) центральной предельной теоремы (ЦПТ). Для описания асимптотического (на больших временах и в больших пространственных масштабах) поведения подобных стохастических процессов применяют аппарат дробных производных. Зачастую решение уравнений Труды Научной конференции по радиофизике, ННГУ, 2004 в частных дробных производных затруднено, или возможно только в частных случаях.

Кроме “чисто” аномально-диффузионных процессов, для многочисленных приложений представляют интерес, до сих пор практически неизученные, “суб-аномальные” случайные процессы, до некоторого момента подчиняющиеся аномальнодиффузионным закономерностям, а на очень больших временах удовлетворяющие ЗБЧ и ЦПТ. В данной работе исследуются подобные процессы случайных блужданий, проявляющие промежуточный аномально-диффузионный характер, и подчиняющиеся линейному закону диффузии при t.

Был рассмотрен процесс многомерных случайных блужданий, подчиняющийся стохастическому уравнению dX(t ) = h k (t t k ), dt k

–  –  –

Работа выполнена при частичной финансовой поддержке грантов Ведущие Научные Школы НШ-838.2003.2, РФФИ 03-02-16680, а также гранта Минобразования РФ Е02-3.5-232.

[1] Metzler R., Klafter J. //Phys. Rep. 2000. V.339, №1. P.1.

[2] Учайкин В.В. //ЖТФ. 1998. T.68, вып.1. С.138.

[3] Учайкин В.В. //ТМФ. 1998. Т.115, №1. С.154.

[4] Учайкин В.В., Саенко В.В. //ЖТФ. 2001. T.71, вып.2. С.8.

[5] Rinn B., Dieteich W., Maass P. //Phil. Mag. 1998. V.77, №.5. P.1283.

[6] Schaufler M.F., Schleich W.P., Yakovlev V.P. //Phys. Rev. Lett. 1997. V.39. P.383.

[7] Zaslavsky G.M., Benkanda S. Chaos, Kinetics and Nonlinear Dynamics in Fluids and Plasmas. -Berlin: Springer, 1998.

[8] Zaslavsky G.M., Edelman M., Niyazov B.A. //Chaos. 1997. V.7, №1. P.159.

[9] Саичев А.И., Уткин С.Г. //Актуальные проблемы статистической радиофизики (малаховский сборник). 2002. Т.1, №1. С.5.

[10] Saichev A.I., Woyczynski W.A. Distributions in the Physical and Engineering Sciences. Volume I. –Boston: Birkhuser, 1997, 336p.

–  –  –

где число N определяет размерность базиса волновых функций фi, Сj – константы разложения собственных функций гамильтониана H по этому базису.

Поверхность потенциальной энергии – это значительный фактор, определяющий точность вариационных расчетов колебательно-вращательных спектров. При расчетах, в основе потенциальной поверхности обычно используется высококачественная ab initio поверхность, рассчитанная методами квантовой химии. Для повыТруды Научной конференции по радиофизике, ННГУ, 2004 шения точности расчетов проводят оптимизацию поверхности потенциальной энергии. Обычно она осуществляется умножением потенциальной функции на некий морфический полином, параметры которого подбираются из условия минимизации стандартного отклонения. Интраполяционные свойства оптимизированного потенциала, как правило, получаются довольно высокими как по энергии, так и по квантовым числам. Экстраполяционные – несколько хуже.

Авторами были проведены совместные оптимизации поверхностей потенциальной энергии для молекул – H216O, H217O и H218O. При оптимизации в основу была положена высококачественная, недавно рассчитанная нами ab initio поверхность потенциальной энергии молекулы воды [1]. К ней были добавлены поправки к приближению Борна-Оппенгеймера: адиабатическая [2], релятивистская [3,4], поправка, учитывающая лэмбовский сдвиг [5] и неадиабатическая поправка [2].

Оптимизации производилась до энергий 25000 см-1.

Первая оптимизация проводилась для J=0, 2 и 5. Она включала 1788 уровней энергии. При оптимизации было использовано 27 параметров. Получено стандартное отклонение рассчитанных уровней от экспериментальных 0,071 см-1. Во вторую оптимизацию были добавлены уровни энергии с J=10. Было получено стандартное отклонение 0,077 см-1. При третьей оптимизации были взяты лишь те экспериментальные уровни энергии, которые дают не большую ошибку по сравнению с расчетом. Было использовано 1004 уровня, что позволило достигнуть точности расчетов 0,028 см-1.

В таблице приведены результаты расчета уровней энергии с полученными поверхностями потенциальной энергии. N – это количество уровней, использованное при вычислении стандартного отклонения. При J=20 и N=334 все уровни с v22 и Ka15 были исключены.

J N стандартное отклонение опт. 1 опт. 2 опт. 3 20 352 0.312 0.206 0.37 20 334 0.204 0.159 0.23 39 10 0.087 0.047 – 40 6 0.055 0.063 – Работа выполнена при поддержке РФФИ (проект № 03-02-16125).

[1] Polyansy O., Csaszar A., Shirin S., Zobov N., Barletta P., Tennyson J., Schwenke D., Knowles P. //Science. 2003. V.299. P.539.

[2] Schwenke D. //J. Phys. Chem. A. 2001. V.105. P.2352.

[3] Csaszar A., Kain J., Polyansky O., Zobov N., Tennyson J. //Chem. Phys. Let. 1998.

V.293. P.317.

[4] Quiney H., Barletta P., Tarczay G., Csaszar A., Polyansky O., Tennyson J. //Chem.

Phys. Let. 2001. V.344. P.413.

[5] Pyykko P., Dyall K.G., Csaszar A., Tarczay G., Polyansky O., Tennyson J. //Phys.

Rev. A. 2001. V.63. P.024502.



Похожие работы:

«mini-doctor.com Инструкция Коджинейт Фс лиофилизат 250 мо ВНИМАНИЕ! Вся информация взята из открытых источников и предоставляется исключительно в ознакомительных целях. Коджинейт Фс лиофилизат 250 мо Действующее вещество: Фактор свертывания крови VIII Лекарственная форма: Лио...»

«Новиков Сергей Витальевич Энергетический беспорядок и транспорт носителей заряда в неупорядоченных органических материалах 02.00.04 – физическая химия АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математи...»

«ОБЩАЯ ХАРАКТЕРИСТИКА Практика является необходимой составляющей производственного процесса студентов по направлению подготовки 02.03.03 "Математическое обеспечение и администрирование информационных систем" и проводится в соответствии с учебным планом. Форма проведения Курс Название практик...»

«1981 г. Ноябрь Том 135, вып. 3 УСПЕ ХИ ФИЗНЧ ЕС КИХ НАУК -53(092) ВКЛАД МАКСВЕЛЛА В РАЗВИТИЕ МОЛЕКУЛЯРНОЙ ФИЗИКИ И СТАТИСТИЧЕСКИХ МЕТОДОВ М. А. Елъяшевич, Т. С. Протъко СОДЕРЖАНИЕ 1. Введение 381 2. Истоки исследований Максвелла по молекулярно-кинетической теории газов 383 3. Общий обзор...»

«1 А. Скляров Приложения к трактату „Основы физики духа“ Приложение N 1: Единая Физика Духа и Материи как новая научная парадигма (Стр. 2 28) Приложение N 2: Возможна ли у нас партия нового типа?. (Стр. 29 74) Приложение N 3: Идеологический манифест Па...»

«Учреждение образования "БЕЛОРУССКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ" БИБЛИОГРАФИЯ Рекомендовано учебно-методическим объединением высших учебных заведений Республики Беларусь по химико-технологическому образованию в качестве пособия для студентов высших учебных заведений, обучающихся по специальности 1-47 01 01 "Издательское дел...»

«В. И. Светцов ОПТИЧЕСКАЯ И КВАНТОВАЯ ЭЛЕКТРОНИКА Учебное пособие Иваново МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ Ивановский государственный химико-технологический университет В. И. Светцов ОПТИЧЕСКАЯ И КВАНТОВАЯ ЭЛЕКТРОНИКА...»








 
2017 www.net.knigi-x.ru - «Бесплатная электронная библиотека - электронные матриалы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.