CC BY
21
Выход из программы может быть осуществлен досрочно в том случае, если полученное на ранних итерациях решение не улучшается, в том числе не снижается степень риска выбранного портфеля ценных бумаг.
Таким образом, в данной работе предложен подход к решению класса задач, определенных нелинейной целевой функцией и нелинейными ограничениями на принципах генетических алгоритмов оптимизации.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Дружинин М.А. Сборник материалов тринадцатой Междунар. науч.-практ. конф. «Управление качеством» // ПР0бЕл-2000, МАТИ, 2014. - С. 109-1l0.
2. Гладков Л.А., Курейчик В.В., Курейчик В.М. Генетические алгоритмы. - М.: Физматлит, 2006. - 318 с.
3. Балыбердин В.А., Белевцев А.М., Дружинин М.А. Генетические алгоритмы поиска в задачах оптимизации систем сетецентрического управления специального назначения // Известия ЮФУ. Технические науки. - 2011. - № 5 (118). - С. 153-159.
4. Белевцев А.М., Дружинин М.А. Разработка и исследование адаптивного поискового алгоритма для решения многоэкстремальных задач оптимизации информационных процессов в информационных системах с распределенной обработкой данных // Известия ЮФУ. Технические науки. - 2013. - № 5 (130). - С. 162-165.
Статью рекомендовал к опубликованию д.т.н., профессор В.А. Петраков.
Белевцев Андрей Михайлович - Научно-исследовательский институт технологий нового поколения; e-mail: [email protected]; 121522, г. Москва, ул. Оршавская, 3; тел.: 89037691788; директор; д.т.н.; профессор.
Дружинин Михаил Александрович - 3 ЦНИИ Минобороны России (3 Центральный научно-исследовательский институт Минобороны России); e-mail: [email protected]; 107564, Москва, Погонный пр-д, 10; тел.: 89261334779; заместитель начальника отдела.
Belevtsev Andrey Michailovich - Research institute of the new generation of Technologies; e-mail: [email protected]; 3, Orshavckay street, Moscow, 121522, Russia; phone: +79037691788; director; dr. of eng. sc; professor.
Druzhinin Mihail Aleksandrovich - 3 the Central Scientific Research Institute of the Ministry of Defence of the Russian Federation; e-mail: [email protected]; 10, Pogonny’j, Mockow, 107564, Russia; phone: +79261334779; deputy division chief.
УДК 681.3.06:681.323(519.6)
А.Н. Голиков
КУСОЧНО-ИНТЕРПОЛЯЦИОННОЕ НЬЮТОНОВСКОЕ ПРИБЛИЖЕНИЕ ФУНКЦИЙ ДВУХ ПЕРЕМЕННЫХ, ЧАСТНЫХ ПРОИЗВОДНЫХ И ДВОЙНЫХ ИНТЕГРАЛОВ ПО КРУГОВОМУ СЕКТОРУ
Излагается кусочно-интерполяционная схема приближения действительных функций двух действительных переменных на круговом секторе интерполяционными многочленами Ньютона с прямоугольной системой узлов в полярной системе координат с началом в центре сектора. Приближение строится на равномерном разбиении исходной области двумя семействами параллельных прямых таких, что прямые из разных семейств взаимно перпендикулярны. В каждой такой подобласти в узлах, расположенных с постоянным шагом вдоль координатных осей, строится интерполяционный многочлен Ньютона. Многочлен преобразуется к каноническому виду, после чего используется для численного дифференцирования и интегрирования приближаемой функции. Степень многочлена и мелкость разбиения исходной области выбираются алгоритмически минимальными так, чтобы дос-
l3l
тичь априори заданной границы абсолютной погрешности приближения. Схема является базой программного комплекса для уточнённого моделирования электрон-фононного транспорта в квантовых проволоках. Схема отличается научной новизной по построению на основе интерполяционных многочленов Ньютона с приведением к каноническому виду и минимизацией временной сложности и памяти.
Кусочно-интерполяционное приближение; интерполяция по Ньютону; приближение частных производных; приближенное вычисление двойных интегралов.
A.N. Golikov
PIECWISE INTERPOLATION FOR BIVARIATE FUNCTIONS, PARTIAL DERIVATIVES, AND DOUBLE INTEGRALS OVER A CIRCLE SECTOR
The piecewise interpolation for real bivariate functions, partial derivatives, and double integrals is presented for a circle domain. The computer approximation use Newton’s interpolation on a rectangle. The polar coordinate system is defined with the center at the center of the circle. The approximation is constructed using a uniform sub-division of the domain by families ofparallel lines. In each sub-domain the polynomial is constructed at uniformly distributed nodes; thereafter it is transformed to a canonical form. At this form, the interpolation polynomial is used for numerical differentiation and integration. Degree ofpolynomial and size of the sub-domain are varied by an algorithm; at that, an a priori given boundary of approximation error is reached. The presented numerical scheme is a base of program for a computer simulation of electron-phonon scattering with a high accuracy. The scheme differs from analogues by a time complexity minimization, by a construction on a base of Newton’s interpolation with a transformation to a canonical form.
Piecewise interpolation; Newton’s interpolation; approximation of partial derivatives; approximation of double integrals.
Введение и постановка задачи. Задачи управления [1] обусловливают необходимость создания быстродействующей и энергоэффективной электроники. Одним из перспективных направлений технологического решения таких задач на элементном уровне является разработка электронных устройств на основе квантовых проволок [2].
В процессе моделирования наноэлементов на основе квантовых проволок требуется приближать, а также численно дифференцировать и интегрировать функции вида
f: G^R, R2 зG = {а(г,^)\ р(С0,А)<ртах, ^min <^<^max}, (1)
где f е СNx+n+1(g) , СN+n+1 (g) - класс функций, непрерывных вместе со своими частными производными до (N + Ny +l) -го порядка включительно на круговом секторе G; (pmin и i^max - углы между прямолинейными границами сектора и осью абсцисс; р(С0, а) - евклидово расстояние между центром кругового сектора С0 и текущей точкой A. Для простоты в дальнейших рассуждениях считается заданной полярная система координат с центром в центре сектора G из (1).
В статье решается задача построения компьютерного приближения функции f в замкнутой области G, f и G из (1), с достижением априори заданной границы абсолютной погрешности приближения, а также с учётом ограничений памяти и временной сложности. С этой целью приближающий ньютонов многочлен строится на равномерном разбиении области (1) с преобразованием к каноническому виду с постоянными числовыми коэффициентами [3, 4]. Такое преобразование позволяет использовать полученный многочлен для численного дифференцирования и интегрирования функции (1), используя реализованные в виде стандартной подпрограммы простые правила вычисления частных производных и первообразных от степенной функции двух переменных.
Ниже описывается алгоритм разбиения области (1) с минимизацией степени приближающего многочлена и числа подобластей.
Минимизация числа подобластей и степени многочлена при кусочноинтерполяционном приближении на равномерном разбиении исходной области. Работа алгоритма осуществляется в ограничениях памяти (числа подобластей), и временной сложности приближения (степени интерполяционного многочлена). Названные ограничения формализованы в виде условий
[N < N , N < N ,
I r max’ p max’
Ik < k k < k (2)
^ r — map p — max’
где N и N - степени интерполяционного многочлена относительно переменных r и ф соответственно; 2kr и 2kp суть число подобластей вдоль соответствующих координатных осей, значения Nmax и kmax задаются априори.
Область (1) разбивается на прямоугольные подобласти G~ прямыми
r = І • hr , (pj = Фшіп + j • К , (3)
где hr =^ , hp=Pmax -Pmin , Pmax , Фшіп и p„,ax из (1), І = 0, 2k -1,
r)kr т r>fkp
] = 0, 2кг -1.
Построение приближающего многочлена Р~ ы ы (г,р) начинается при минимальных значениях N г = = 1, кг = к^= 0. После преобразования интерпо-
ляционного многочлена к каноническому виду
N N
Р N Ы=ХХ аїі/ф3 , (4)
1=0 >=0
в проверочных точках, расположенных вдоль координатных осей с постоянным априори заданным шагом, меньшим шага интерполяции, проверяется условие
P~N (r,P)- f (r,P)
<є, (5)
где граница е задаётся априори.
Если условие (5) выполнено во всех проверочных точках, то приближающий многочлен считается построенным. Если условие (5) не выполнено хотя бы в одной точке, то значения Nг и Nу увеличиваются на единицу и построение повторяется до удовлетворения условия (5) либо до исчерпания пределов (2) по Nг и
—у . В последнем случае устанавливаются единичные значения Nг и , а кг и к увеличиваются на единицу, и так далее.
Таким образом, если существует набор некоторых Nг, N у , кг, ку , при которых условие выполняется во всех проверочных точках, то значения степени и числа подобластей минимальны по построению, что обеспечивает минимизацию временной сложности приближения и памяти, требуемой для хранения коэффициентов а ~ из (4).
После построения многочлена (4) коэффициенты а~ записываются в память компьютера для дальнейшего использования. В таком случае приближение функции с использованием схемы Горнера займёт время ф) = ^хЫу (^ + tc), где
^ и tc - время одного умножения и сложения соответственно; при этом индекс к подобласти 0~ вычисляется при помощи соотношений
к = j • 2кр + i +1. i =
r +1. j = Р Qmin
К _ К Р
+1. (б)
где [а] - целая часть числа а, Пг и Пр из (3).
Ниже описывается преобразование интерполяционного многочлена Ньютона в текущей подобласти О ~ к виду (3).
Кусочно-интерполяционное приближение функций двух переменных на основе многочленов Ньютона. Преобразование интерполяционного многочлена к виду (4) инвариантно относительно степени, числа подобластей и номера подобласти, поэтому значения Nг, , кг, кр и к считаются определёнными по опи-
санному в предыдущем параграфе алгоритму и фиксированными.
В подобласти О~ задаётся прямоугольная система равноотстоящих узлов
О~ = {( Г',/ ’Рт )
к _ К
и g — - шаги интерполяции; пг и
' г р
добласти (3). В узлах (7) строится интерполяционный многочлен Ньютона
r, = ri + Рj,m = Pj + mg, , = 0 Nr, m = 0 Np Г . (7)
где П —~^ и ^ — ~^ - шаги интерполяции; Пг и Пр - линейные размеры по-
Nr Np АП+\ f (r Р ) m , к
-r- Aг~рк/ Vi,o,pj,o)т—г/ \т—і
WM-iІ '/ТГП(r-r,)П(p-p„.). (8)
т!Птк! рк АА ‘,{/АА
т=0 к=0 т!п к! 6 £=0 ?=0
где Ат+кк/(.,0,9-,0) - двойные конечные разности; Д0^/(рЛ0) = I(.,0, 9Л0)• После преобразований многочлен (8) запишется в виде
N. • 9 т к
Рк,мг ,мг ^ и(9 £ йт,к I! (t - г ЯП( и - -), (9)
9 т=0 к=0 г=0 -=0
, Дтгт9к1 (гг,0, 9,0) Г - Г,.0 9 - 9,0
где Ьт,к =----------171-------; t =—Г~ , 9 =------------— •
т! к! п 6
В процессе приведения многочлена (9) к виду (4) по корням многочленов
т т к к ^
Кт ^ )=П( t - г )=£ ак/ , Ок (и )=П( и - -)=£ ~т-к-и
г=0 г=0 у=0 у=0
вычисляются их коэффициенты при помощи матричной схемы из [5]. Для равномерных сеток узлов интерполяции схема, используя только операции сложения и умножения, применяется к целочисленным данным и поэтому не копит погрешности округления.
Многочлены Ят (t) и Qk (и) зависят только от степени интерполяционного многочлена и инвариантны относительно приближаемой функции (1), поэтому коэффициенты dk,г и dm_k,^ после однократного вычисления и занесения в память компьютера считаются априори известными.
Вычислив коэффициенты dk г и dm-k,j, получим
к т-к к т-к к т-к
П(t - /)П(и - ^ХХ dk,idm-k,jtU =ХХ °т,к,и?и , (10)
где ^т,к,1^ = dk,idm-k, j .
Учитывая (10), раскрыв скобки и приведя подобные члены в (9), окончательно запишем
N
РКМ((ги((Р)) = XXак,ц , (11)
N
где а~,у = XXЬткВт,к,г^ , * = 1 Nr , 7 = 1 ^ .
т=у к=г
Изложенная схема нацелена на компьютерное приближение функций (1) из стандартной библиотеки, например при моделировании электронного транспорта в
квантовых проволоках [2, 6]. В таком случае коэффициенты а~ ^ заносятся в память компьютера, и приближение функции (1) сводится к вычислению многочлена
(11) по схеме Горнера с предварительной дешифрацией (6) индекса к подобласти.
Схема (11) обладает параллелизмом и инвариантностью относительно размера области (1) [7].
Ниже в рамках излагаемого метода приближаются частные производные функции (1) порядка не выше степени многочлена (11) по соответствующим переменным.
Кусочно-интерполяционное приближение частных производных от функции двух действительных переменных. В предположении, что многочлен (11) построен, частные производные от функции (1) приближаются частными производными от многочлена (11), т.е. имеет место приближённое равенство
а"У(х,у) * к-I„-к'Хр’Х'—_____________— а~ (12)
алу к я х х (г - £)! (/ - к)! а~ ^ и , (12)
где 0 < £ < Ыг, 0 < к < Nу .
В (12) многочлен (11) продифференцирован как сложная функция от х и у ,
, \ , ч т!
а также учтено, что т (т -1)... (т - п +1) = . Последнее позволяет, предва-
(т - п)!
рительно записав в память компьютера члены последовательности факториалов: 0!, 1!, ... N, вычислять произведения вида т(т-1)...(т-п + 1) за время одного деления, что снижает временную сложность приближённого дифференцирования.
При запомненных предварительно коэффициентах (11) приближение производных сводится к дешифрации (6) и вычислению значения многочлена (12) по схеме Горнера за время t (1)= N - к)(^ -£)(2^ + 2^ + ^ ), где ^ ^, tc - время одного умножения, деления и сложения соответственно [7].
г=0
г=0 у=0
Кусочно-интерполяционное приближение двойных интегралов по круговому сектору. Пусть многочлен (11) построен, тогда в рамках излагаемого метода двойной интеграл от функции f по круговому сектору G, f и G из (1) приближается соответствующим интегралом от многочлена (11). В силу аддитивности интеграла по области интегрирования имеет место приближённое равенство
2kr+кф
jj f (x y)dxdy * 2 If P~N,Nv (t(r) u(<p)) drdv . (13)
G ~=0 G~
k
С учётом (11), а также с учётом замен переменных t(r) и u(<p) из (9) соотношение (13) запишется в виде
2kr +ky-1 Nr Nr Nr Nr
If f (x , y) dx dy * 2 hg f dt f a~ t u^ du ,
G k~=0 0 0 *'=0 J=0
откуда окончательно получим
2kr+^-1 Nr N 1
x,y) dxdy *hg 2 2 2{;+lw/+,)a~,j tMui++. (14)
G k=0 i=0 i=0 (i + V(i + V
Схема (14) ориентирована на компьютерное вычисление двойных интегралов по круговому сектору (1). Схема по построению инвариантна относительно размеров области и обладает естественным параллелизмом, как показано в [7].
Численный эксперимент. Для численного эксперимента использовался компьютер на базе процессора Intel Core 2 Quad Q6600 с 8 Гб RAM.
Погрешность приближения функций на изложенной основе не превосходит по порядку величину 10-18, частных производных первого порядка - 10-15, прямых производных второго порядка -10-11, смешанной производной второго порядка -10-9. Абсолютная погрешность приближённого вычисления двойных интегралов не более 10-18.
Таким образом, в статье изложены компьютерные схемы приближённого вычисления функций двух переменных, их частных производных и двойных интегралов по круговому сектору. Схемы инвариантны относительно размера области определения функции, обладают естественным параллелизмом [7]. Подход нацелен на использование для приближения функций, производных и интегралов из стандартной библиотеки. За счёт алгоритмического подбора наименьшей степени приближающего многочлена для данного числа подобластей удаётся снизить временную сложность вычислений при достижении априори заданной границы абсолютной погрешности.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Гайдук А.Р., Плаксиенко Е.А. Управление нелинейными объектами с компенсацией неопределённого возмущения // Мехатроника, автоматизация, управление. - 2013. - № 1.
- С. 2-8.
2. Борздов А.В., Поздняков Д.В., Борздов В.М., Орликовский А.А., Вьюрков В.В. Моделирование влияния поперечного электрического поля на дрейфовую скорость электронов в GaAs квантовой проволоке // Микроэлектроника. - 2010. - Т. 39, № 6. - C. 436-442.
3. Голиков А.Н. Кусочно-полиномиальная схема аппроксимации функций двух переменных, частных производных и двойных интегралов с повышенной точностью // Известия ЮФУ. Технические науки. - 2011. - № 5 (118). - С. 179-186.
4. Ромм Я. Е., Фирсова С.А. Минимизация временной сложности вычисления функций с приложением к цифровой обработке сигналов. - Таганрог: Изд-во Таганрог. гос. пед. ин-та, 2008. - 125 с.
5. Ромм Я.Е. Локализация и устойчивое вычисление нулей многочлена на основе сортировки. II // Кибернетика и системный анализ. - 2007. - № 2. - С. 161-174.
6. Голиков А.Н. Самосогласованный расчёт электрон-фононного рассеяния в GaAs нанопроволоках на основе кусочно-полиномиальных схем. - Таганрог: ТГПП, 2011. - 123 с.
7. Голиков А.Н. Кусочно-полиномиальные схемы вычисления функций двух переменных, частных производных и двойных интегралов на основе интерполяционного полинома Ньютона. - Таганрог: ТГПИ., 2010. - 150 с.
Статью рекомендовал к опубликованию д.т.н., профессор Л.П. Фельдман.
Голиков Александр Николаевич - НОУ ВПО «Таганрогский институт управления и экономики»; e-mail: [email protected]; 347900, г. Таганрог, ул. Петровская, 45; кафедра прикладной математики и информационных технологий; к.т.н.; старший преподаватель.
Golikov Alexander Nikolayevich - Taganrog Management and Economics Institute; e-mail: [email protected]; 45, Petrovskaya street, Taganrog, 347900, Russia; the department of applied mathematics and information technologies; cand. of eng. sc.; lecturer.
УДК 519.876.5
Е.Н. Бородулина
ГЕНЕТИЧЕСКИЙ АЛГОРИТМ РЕШЕНИЯ ЗАДАЧ СИТУАЦИННОГО
УПРАВЛЕНИЯ
Рассмотрена задача разработки модели поддержки жизненного цикла системы в условиях многовариантного развития событий и инструментария для ее реализации. Предложена адаптивная модель, включающая три компоненты: структурную модель управляемой системы, сценарно-вероятностную модель внешней среды и модель управляющей системы. Отмечается, что система, проектированная на основе данной модели, является достаточно гибкой и позволяет решать задачи прогнозирования перехода системы из одного состояния в другое путем предварительного задания параметров внешней среды, управления изучаемой системой с помощью выбранного управляющего воздействия, управления внешней средой с целью превращения воздействия в воздействие, улучающее параметры функционирования системы, изучения параметров среды при разных ситуациях и синтеза эффективной управляемой системы. В качестве ядра инструментария предложен генетический алгоритм с динамическим выбором генетических операторов в группе, позволяющий решить задачу выбора оптимального управляющего воздействия при условии предварительного формирования целевой функции на основе анализа внешнего воздействия, параметров внутренней среды и рисков.
Генетический алгоритм; ситуационное управление; социотехническая система; сценарно-вероятностная модель; управляемая система; управляющая система.
E.N. Borodulina
THE GENETIC ALGORITM OF SITUATION MANAGEMENT PROBLEM
SOLUTION
The problem of construction of life cycle support system model under multi-version event development and instruments for implementation is considered. The adaptive model including three components such as structural model of controlled system, scenary probabilistic model of surroundings and controlling system model and allowing solving the problems of sophisticated social
