Решение систем полиномиальных уравнений на ЭВМ

ления менее релевантных ИО. Второй тип штрафа менее предпочтителен, так как приводит к тому, что актор может покинуть площадку.

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

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

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

доставляющего максимальному количеству сотрудников предприятия доступ в единое информационное пространство.

Литература

1. Колчин А.Ф., Овсянников М.В., Стрекалов А.Ф., Сумароков С.В. Управление жизненным циклом продукции. М.: Ахарсис, 2002. 304 с.

2. Leitäo P. Holonic rationale and bio-inspiration on design of complex emergent and evolvable systems // Trans. on large-scale data and knowledge centered systems I, LNCS 5740, 2009. SpringerVerlag Berlin Heidelberg, pp. 243-266.

3. Бурков В.Н., Коргин Н.А., Новиков Д.А. Введение в теорию управления организационными системами. М.: Либро-ком, 2009. 264 с.

4. Каляев И.А. Стратегии группового управления в распределенных системах // Управление в распределенных сете-центрических и мультиагентных системах: матер. науч.-технич. сем. СПб: Концерн ЦНИИ «Электроприбор», 2010. C. 8-9.

5. Прикладной анализ случайных процессов; [под ред. С.А. Прохорова]. Самара: Изд-во СНЦ РАН, 2007. 582 с.

6. Иващенко А.В. Модель многоакторной интегрированной информационной среды предприятия // Вестн. СамГУПС. Самара: СамГУПС, 2012. № 1 (15). C. 103-109.

7. Иващенко А.В. Управление согласованным взаимодействием пользователей интегрированной информационной среды предприятия. Самара: Изд-во СНЦ РАН, 2011. 100 с.

References

1. Kolchin А.Б., Ovsyannikov М^., Strekalov А.Б., Sumaro-kov S.V., Upravlenie zhiznennym tsiklom produktsii (Product Lifecycle Management), Moscow, Akharsis, 2002, 304 p.

2. Leitäo P., Holonic rationale and bio-inspiration on design of complex emergent and evolvable systems. Trans. on large-scale data and knowledge centered systems I, LNCS 5740, 2009, SpringerVerlag Berlin Heidelberg, pp. 243-266.

3. Burkov V.N., Korgin NA, Novikov DA., Vvedenie v teo-riyu upravleniya organizatsyonnymi sistemami (Introduction in Organizational System Management Theory), Moscow, Librokom, 2009, 264 p.

4. Kalyaev 1.А., Materialy nauchno-tekhnicheskogo seminara «Upravlenie v raspredelennykh setetsentricheskikh i multiagentnykh sustemakh» (Materials of Scientific and Technical Seminar «Management in Distributed Network-Centric and Multiagent Systems»), St. Petersburg, Kontsern «Elektropribor», 2010, pp. 8-9.

5. Prikladnoy analiz sluchainykh protsessov (Application Analysis of Random Processes), [ed. by SA. Prokhorov], Samara, Izd. Samara Research Center RAS, 2007, 582 p.

6. Ivashchenko А^., VestnikSamGUPS (Journ. of the Samara State Univ. of Railway Transport), 2012, no. 1 (15), pp. 103-109.

7. Ivashchenko А^., Upravlenie soglasovannym vzaimo-deystviyem polzovateley integrirovannoy informatsionnoy sredy predpriyatiya (Management of Consistent Interaction of Enterprise Integrated Information Environment Users), Samara, Izd. Samara Research Center RAS, 2011, 100 p.

УДК 681.3

РЕШЕНИЕ СИСТЕМ ПОЛИНОМИАЛЬНЫХ УРАВНЕНИЙ НА ЭВМ

И.А. Лёзин, кт.н..

Предложен ускоренный алгоритм поиска базисов Грёбнера, используемых при решении систем полиномиальных уравнений. Данный алгоритм рассматривает проблему переполнения разрядной сетки и некорректных вычислений при проведении расчетов на ЭВМ. Классический алгоритм Бухбергера и обновленный алгоритм Фужера требуют

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

Ключевые слова: системы полиномиальных уравнений, базис Грёбнера, алгоритм Фужера, алгоритм Бухберге-

ра.

SOLUTION OF SIMULTANEOUS POLYNOMIAL EQUATIONS ON A COMPUTER Lyozin I.A., Ph.D.

Abstract. The article offers the accelerated algorithm of Groebner bases search used in the course of solution of simultaneous polynomial equations. The proposed algorithm considers the issue of word size overflow and incorrect computations in the course of calculating on the computer. The classical Buchberger's algorithm and the updated Faugere's algorithm require a lot of steps; meanwhile they do not contain any additional actions on correction of floating point numbers computations that causes some errors and incorrect computation of simultaneous equations roots and irrational factors. The new algorithm by several times reduces the number of operations of polynom quantity reduction, does not store all possible pairs of polynoms and does not check them for the fact that one can be expressed in terms of the other. The algorithm proposed in the article operates only with that number of polynoms which are now a basis. The basis is expanded if none of the polynoms represented in it can be represented in the form of combination of the other polynoms of this basis. Otherwise, the extra polynom shall be removed from the basis that helps avoid overgrowth of the basis and excess operations with this polynom. In order to avoid the problem of incorrect computations because of the word size of the computer the algorithm envisages post-check. If any of the polynoms is reduced, then the initial polynom should be represented by the combination of the reduced polynom and the rest of the basis polynoms. If the values of coefficients of the equal monoms do not coincide after that, then the coefficients of the reduced polynom require correction by the error size.

Keywords: systems of polynomial equations, Grobner basis, Faugere algorithm, Buchberger algorithm.

Системы полиномиальных уравнений в общем случае решаются построением базиса Грёбнера из базиса, составленного многочленами рассматриваемой системы, и решением системы уравнений, левые части которых являются многочленами, составляющими базис Грёбнера [1]. В числе самых широко распространенных методов для построения базиса Грёбнера - алгоритм Бухбергера и алгоритм Фужера [2-4]. Однако оба они при относительной простоте выполняемых преобразований имеют очень высокую операционную сложность, что затрудняет, а зачастую делает невозможным вычисление базиса.

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

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

Итак, пусть K - поле вещественных чисел, К[х, у] - кольцо многочленов от двух переменных с коэффициентами из К. Рассмотрим некоторый идеал I над кольцом многочленов, в который входят полиномы из системы. Множество H={g, К} является базисом идеала I. Зададим лексикографическое упорядочение мономов в многочленах:

Хрхуру > XV о

Рх > Чх М Рх = Чх Ру > Чу )•

Используя операции сложения и умножения над многочленами, входящими в базис Н, необходимо редуцировать его так, чтобы построить базис Грёбнера идеала I.

Введем понятие ^-полинома [2]:

ЬСМ

5

Здесь ЬТ(/) представляет старший терм многочлена /(.х, у), а ЬСМ/ /) равен наименьшему моному, который делится на старшие мономы многочленов/(х, у) и/(х, у). Операцией редуцирования произвольного полинома /(х, у) из К[х, у] относительно Н является нахождение такого минимального многочлена г(х, у), который выражается в виде

f

где а(х, у) еК[х, у] и /(х, у) еН.

Базис Н можно назвать базисом Грёбнера, если для любых двух полиномов /(х, у) и /(х, у) из Н выполняется условие

я

При заданном лексикографическом упорядочении (1) базис Грёбнера будет составлен парой многочленов, удовлетворяющих условию (4), вида

Мх1 Мху

/ (х,у )=х*пх) х-+5>ту) у-,

п=1 т=0

Т

/у (х у ЬХ^' ут.

т=0

Множество решений системы, составленной из полиномов (5), будет соответствовать множеству решений системы.

Алгоритм нахождения базиса Грёбнера представлен в виде блок-схемы (см. рис.).

На вход алгоритма редукции подается базис Н идеала I, составленный из множества полиномов исходной системы уравнений. Каждая итерация алгоритма редуцирования базиса начинается с того, что в базисе Н ищется полином /(х, у) с минимальным старшим мономом, который может быть редуцирован относительно множества полиномов, состоящего из всех полиномов базиса, кроме выбранного полинома. Дальнейшее выполнение алгоритма разбивается на две ветви.

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

ства полиномов, если в этом множестве есть многочлены, старшие мономы которых являются делителями старшего монома редуцируемого полинома. Редукция полинома - итеративный процесс. Если перед редукцией остаток г(х, у) принять равным исходному редуцируемому полиному /(х, у), то на каждом шаге уменьшения остатка осуществляется поиск максимального полинома /(х, у) в множестве Н, который может редуцировать остаток г(х, у), и выполняется редукция

ЬСМ

г

где ЬС(г) - коэффициент перед старшим мономом многочлена г(х, у).

Редуцирование по формуле (6) следует продолжать до тех пор, пока остаток г(х, у) не станет нередуцируемым. Когда редуцирование завершено и остаток г(х, у) не является нулевым полиномом, он добавляется в базис Н.

Перед добавлением остатка в базис следует провести корректировку коэффициентов для г(х, у). В кольце К[х, у] необходимо найти такие многочлены а&1(х, у), аь,г(х, у), а&г(х, у) и аИ,г(х, у), чтобы можно было представить полиномы g(x, у) и И(х, у) в базисе Н+{г}:

8

Ь

/

Если решение системы (7) не находится, значит, при вычислениях в коэффициентах остатка накопилась методическая погрешность из-за способа представления вещественных чисел в памяти ЭВМ. Так как уже имеющиеся в базисе полиномы имеют правильные значения коэффициентов, необходимо откорректировать коэффициенты остатка г(х, у) в соответствии со значениями погрешности для каждого из коэффициентов в системе (7). Эта задача решается подобно задаче редуцирования: сначала определяются значения коэффициентов для полиномов а&(х, у) и ак(х, у), затем вычисляются разности:

г8 (х у) = 8 (х у)- X а8,1 (х у)л (х у) =

= а8 ,г ( х, у ) г ( х, у ) ,

гь (х у) = И (х у)- X аК1 (х у)л (х у) =

= акг

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

/ H={g, И} /

/х, у)= шт (/¡^Н),/(х,у)^ *Н-{Л}г(х, у)<>0

Да

/(х,у)<>0

Нет

Н=Н-{/}

Нет

г(х, у)<>0

\ Г

5(/, /)^>(х, у)<>0, ш1п (ЬСЫ/ /)) 11

Да

Н=Н+{г}

/ Н={/х, /у} /

Блок-схема алгоритма редукции системы уравнений

базис Н и алгоритм редуцирования базиса идет на очередную итерацию.

По второй ветви алгоритм редуцирования базиса развивается, если в базисе Н не найден полином, который может быть редуцирован. В таком случае набор базисных функций должен быть расширен для продолжения редукции. Из всех возможных пар полиномов базиса /(х, у) и /(х, у), для которых не выполняется условие (4), нужно найти пару с минимальным ЬСЫ(/и /).

Если такая пара найдена, редуцированный остаток ^-полинома этой пары включается в базис Н, предварительно пройдя корректировку коэффициентов, как это было описано ранее. Если же такой пары многочленов в базисе Н нет, значит, этот базис полностью редуцирован и является базисом Грёбнера.

Таким образом, в результате работы алгоритма получим базис, составленный из пары многочленов (5). Теперь для решения исходной системы нужно вычислить все корни у полиномиального уравнения одной переменной/у(х, у)=0, подставить по очереди эти значения в многочлен /х(х, у), превращая его в полином одной переменной, и рассчитать соответствующие корни ; уравнения

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

Литература

1. Прасолов В.В. Многочлены. М.: МЦНМО, 2001.

2. Buchberger B. Grobner Bases: an Algorithmic Method in Polynomial Ideal Theory / Recent trends in multidimensional system theory, U. Reidel Publishing Company, 1985.

3. Buchberger's algorithm. URL: http://www.geocities.com/ famancin/buchberger.html (дата обращения: 01.06.12).

4. Faugere J.C. A new efficient algorithm for computing Grobner bases (F4) // Journal of Pure and Applied Algebra, 1999, no. 139, pp. 61-88.

References

1. Prasolov V.V., Mnogochleny (Polynomials), Moskovsky Tsentr Nepreryvnogo Matematicheskogo Obrazovaniya, 2001.

2. Buchberger B. Grobner Bases: an Algorithmic Method in Polynomial Ideal Theory. Recent trends in multidimensional system theory, U. Reidel Publishing Company, 1985.

3. Buchberger's algorithm, available at: www.geocities.com/ famancin/buchberger.html (accessed 01.06.12).

4. Faugere J.C., A new efficient algorithm for computing Grobner bases (F4), Journal of Pure and Applied Algebra, 1999, no. 139, pp. 61-88.

УДК 681.3

ИССЛЕДОВАНИЕ АППРОКСИМАТИВНЫХ ВОЗМОЖНОСТЕЙ РАДИАЛЬНО-БАЗИСНОЙ СЕТИ С ОРТОГОНАЛЬНЫМИ ПОЛИНОМАМИ

И.В. Лёзина, кт.н.

Описывается постановка задачи аппроксимации, обосновывается возможность использования в качестве аппрок-симатора плотности распределения вероятности радиально-базисной нейронной сети, приводятся аппроксимирующее выражение для данной сети и выражение для целевой функции, с помощью которой происходит подбор параметров базисных функций и значений весов. Рассматривается возможность использования при аппроксимации плотности распределения вероятности радиально-базисной сетью не только традиционных функций Гаусса, но и сигмои-дальных и степенных функций и ортогональных полиномов Лежандра, Чебышева I и II рода, Лагерра и Эрмита. Приводятся соответствующие формулы. Сравниваются погрешности аппроксимации путем вычисления среднего квадратического отклонения. В качестве примеров приводится аппроксимация плотности вероятности Симпсона и Рэлея радиально-базисной сетью c сигмоидальными, степенными функциями, а также полиномами Лежандра, Чебышева I и II рода, Лагерра и Эрмита. Дается рекомендация по использованию радиально-базисной сети с полиномами Лежандра, Чебышева I и II рода в качестве базисных функций при увеличении числа нейронов в скрытом слое, так как такая сеть позволяет достичь более низких значений среднего квадратического отклонения, чем сеть с традиционными функциями Гаусса.

Ключевые слова: аппроксимация, радиально-базисная сеть, ортогональные полиномы, Лежандр, Чебышев, Ла-герр, Эрмит.

STUDY OF APPROXIMATIVE POSSIBILITIES OF RADIAL-BASIC NETWORK WITH ORTHOGONAL POLYNOMS Lyozina I.V., Ph.D. (Samara State Aerospace University, chuchyck@yandex.ru) Abstract. The article describes set up of the problem of approximation, provides rationalization for possibility to use the radial-basic neural network as an approximator of probability density function, gives the approximating statement for this