Численная эффективность одной модификации r-алгоритма
Рассматривается модификация r-алгоритма – алгоритма минимизации с использованием операции растяжения пространства в направлении разности двух последовательных субградиентов. В отличие от r-алгоритма, значения коэффициентов растяжения в предложенной модификации программно рассчитываются в процессе ра...
Збережено в:
Дата: | 2017 |
---|---|
Автор: | |
Формат: | Стаття |
Мова: | Russian |
Опубліковано: |
Інститут кібернетики ім. В.М. Глушкова НАН України
2017
|
Назва видання: | Теорія оптимальних рішень |
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/131435 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | Численная эффективность одной модификации r-алгоритма / Н.Г. Журбенко // Теорія оптимальних рішень: Зб. наук. пр. — 2017. — № 2017. — С. 33-38. — Бібліогр.: 4 назв. — рос. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-131435 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1314352018-03-24T03:03:05Z Численная эффективность одной модификации r-алгоритма Журбенко, Н.Г. Рассматривается модификация r-алгоритма – алгоритма минимизации с использованием операции растяжения пространства в направлении разности двух последовательных субградиентов. В отличие от r-алгоритма, значения коэффициентов растяжения в предложенной модификации программно рассчитываются в процессе работы алгоритма. Алгоритм может использоваться с постоянным шагом. Приводятся результаты исследования численной эффективности алгоритма. Розглядається модифікація r-алгоритму – алгоритму мінімізації з використанням операції розтягування простору в напрямку різниці двох послідовних субградієнтів. На відміну від r-алгоритму, значення коефіцієнтів розтягування у запропонованій модифікації програмно розраховуються в процесі роботи алгоритму. Алгоритм може використовуватися з постійним кроком. Наводяться результати дослідження чисельної ефективності алгоритму. The modification of r-algorithm, the minimization algorithm using space dilation operation along the direction of the difference of two successive subgradients, is considered. In contrast to r-algorithm, in the proposed modification the values of dilation coefficients are calculated during the execution of the algorithm. The algorithm can be used with a constant step size. The results of the study of the numerical efficiency of the algorithm are presented. 2017 Article Численная эффективность одной модификации r-алгоритма / Н.Г. Журбенко // Теорія оптимальних рішень: Зб. наук. пр. — 2017. — № 2017. — С. 33-38. — Бібліогр.: 4 назв. — рос. 2616-5619 http://dspace.nbuv.gov.ua/handle/123456789/131435 519.8 ru Теорія оптимальних рішень Інститут кібернетики ім. В.М. Глушкова НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
Russian |
description |
Рассматривается модификация r-алгоритма – алгоритма минимизации с использованием операции растяжения пространства в направлении разности двух последовательных субградиентов. В отличие от r-алгоритма, значения коэффициентов растяжения в предложенной модификации программно рассчитываются в процессе работы алгоритма. Алгоритм может использоваться с постоянным шагом. Приводятся результаты исследования численной эффективности алгоритма. |
format |
Article |
author |
Журбенко, Н.Г. |
spellingShingle |
Журбенко, Н.Г. Численная эффективность одной модификации r-алгоритма Теорія оптимальних рішень |
author_facet |
Журбенко, Н.Г. |
author_sort |
Журбенко, Н.Г. |
title |
Численная эффективность одной модификации r-алгоритма |
title_short |
Численная эффективность одной модификации r-алгоритма |
title_full |
Численная эффективность одной модификации r-алгоритма |
title_fullStr |
Численная эффективность одной модификации r-алгоритма |
title_full_unstemmed |
Численная эффективность одной модификации r-алгоритма |
title_sort |
численная эффективность одной модификации r-алгоритма |
publisher |
Інститут кібернетики ім. В.М. Глушкова НАН України |
publishDate |
2017 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/131435 |
citation_txt |
Численная эффективность одной модификации r-алгоритма / Н.Г. Журбенко // Теорія оптимальних рішень: Зб. наук. пр. — 2017. — № 2017. — С. 33-38. — Бібліогр.: 4 назв. — рос. |
series |
Теорія оптимальних рішень |
work_keys_str_mv |
AT žurbenkong čislennaâéffektivnostʹodnojmodifikaciiralgoritma |
first_indexed |
2025-07-09T15:26:14Z |
last_indexed |
2025-07-09T15:26:14Z |
_version_ |
1837183562029203456 |
fulltext |
Теорія оптимальних рішень. 2017 33
ТЕОРІЯ
ОПТИМАЛЬНИХ
РІШЕНЬ
Рассматривается модификация
r-алгоритма – алгоритма мини-
мизации с использованием опера-
ции растяжения пространства в
направлении разности двух по-
следовательных субградиентов. В
отличие от r-алгоритма, значе-
ния коэффициентов растяжения
в предложенной модификации
программно рассчитываются в
процессе работы алгоритма. Ал-
горитм может использоваться с
постоянным шагом. Приводятся
результаты исследования числен-
ной эффективности алгоритма.
Н.Г. Журбенко, 2017
УДК 519.8
Н.Г. ЖУРБЕНКО
ЧИСЛЕННАЯ ЭФФЕКТИВНОСТЬ
ОДНОЙ МОДИФИКАЦИИ
r -АЛГОРИТМА
Введение. Более 40 лет назад разработан
субградиентный алгоритм минимизации с
растяжением пространства в направлении
разности двух последовательных градиентов
– r-алгоритм [1]. Практика использования
r-алгоритма показывает, что до настоящего
времени – это один из наиболее эффектив-
ных алгоритмов негладкой оптимизации. Од-
нако теоретическое исследование эффектив-
ности алгоритма далеко не закончено. Ос-
новная проблема теоретического обоснова-
ния r-алгоритма состоит в согласованном
выборе значений коэффициента растяжения
пространства и шагового множителя. В рабо-
те [2] предложено семейство модификаций
r-алгоритма – ( )r -алгоритмы. Величины ко-
эффициентов растяжения пространства на
итерациях ( )r -алгоритмов не постоянны,
они вычисляются в процессе его работы.
Алгоритмы не требуют использования проце-
дуры одномерного спуска по направлению. В
работе [3] приведены результаты исследова-
ния численной эффективности конкретных
вариантов ( )r -алгоритмов. В данной работе
предлагается новая модификация r-алгоритма
из этого семейства и приводятся результаты
численного исследования ее
эффективности.
Численная схема r-алгоритма. Рассмат-
ривается задача безусловной минимизации
субдифференцируемой функции ( )f x в .nR
*Работа выполнена при частичной поддержке
Volkswagen Foundation (грант № 90 306).
Н.Г. ЖУРБЕНКО
34 Теорія оптимальних рішень. 2017
( )f x обозначим множество субградиентов функции ( )f x в точке .x В r -алго-
ритме используется оператор растяжения пространства [4] ( ) ( 1) ,TR I
где ,nR – направление и коэффициент растяжения пространства,
| | 1, 0. Вычислительная схема r -алгоритма применительно к задаче отыс-
кания безусловного минимума функции ( )f x состоит в следующем.
0-ой шаг алгоритма (инициализация). Выбираем начальное приближение
0x и невырожденное линейное преобразование 0.B Вычисляем: 0 0( ) ( )g x f x ;
* *
0 0 0( )g B g x (субградиент в преобразованном пространстве 1
0 0 0 ,Y B X A X X
– исходное пространство). Пусть на шаге k алгоритма ( 0,1,2, )k получены
определенные значения векторов ,kx
*
kg (субградиент в преобразованном про-
странстве) и матрицы .kB ( 1
k kA B – матрица преобразования пространства).
(k + 1)-й шаг алгоритма ( 0,1,2, ).k
Вычисляем:
1) 1kh шаговый множитель, 1 0;kh
2) * *
1 1 / | |;k k k k k kx x h B g g
3) 1( ) ( ),k kg x f x (субградиент в точке 1kx );
4)
* *
1 1( )k k kg B g x (субградиент в преобразованном пространстве k kY A X );
5)
* * * *
1 1 1( )/ | |k k k k kg g g g (направление растяжения пространства kY );
6) 1 1 11, 1/k k k ( k коэффициент растяжения пространства kY );
7)
11 ( ),
kk k kB B R
(обратный оператор преобразования пространства
1
1 1 1k k kY A X B X
);
8)
1
*
1 1 1( )
kk k kg R g
* *
1 1( ( ),k k kg B g x (субградиент в преобразованном
пространстве
1
1 1 1k k kY A X B X
).
Переходим к ( 2)k -му шагу алгоритма или прекращаем работу при
выполнении критерия останова.
В r -алгоритме значения коэффициентов растяжения пространства (пара-
метр r -алгоритма) выбираются одинаковыми на всех итерациях: 1.k На
практике рекомендуется это значение выбирать порядка 2.0. Величина шагового
множителя определяется процедурой минимизации («спуска») по направлению
* */ | | .k k kB g g Обычно применяемые процедуры «спуска» являются достаточно
грубой реализацией алгоритма локализации минимума по направлению. Наибо-
лее часто используется процедура адаптивной регулировки шаговых множителей
[1, 4]. Эта процедура определяется параметрами 10 1,q 2 1q и целым чис-
лом 2.L Эти величины являются параметрами (константами) алгоритма.
ЧИСЛЕННАЯ ЭФФЕКТИВНОСТЬ ОДНОЙ МОДИФИКАЦИИ r-АЛГОРИТМА
Теорія оптимальних рішень. 2017 35
Вычислительная схема r ( )-алгоритмов соответствует приведенной схе-
ме r -алгоритма. Отличие состоит лишь в следующем. Вместо оператора растя-
жения ( )R используем оператор
( ) ( 1) ,TR I (1)
где ,nR – нормирующий множитель, 1,R 0. В отличие от операто-
ра ( ),R вектор не нормирован, т. е. выполнение условия | | 1 не требует-
ся. Различные варианты алгоритма будут определяться выбором нормирующего
множителя .
Отметим простейшие свойства оператора ( ) :R * ;R R 1( )R
2/ (1 ) ;T I (0)R I . Пусть | | 0 и / | | . Тогда
( ) ( ) ,R R I где
21 | | . (2)
Таким образом, если 0, то оператор ( )R является оператором растя-
жения по направлению / | | . Значение коэффициента растяжения определяет-
ся формулой (2). Значение нормирующего множителя будет определяться на
основании субградиентов (в преобразованном пространстве)
* *
1, :k kg g
* *
1 1( , ).k k kg g Естественным требованием на функцию 1 2( , )g g будет
выполнение условия (условия «однородности»)
2
1 2 1 2( , ) ( , ) / ,g g g g
где
1, 0.R Это условие обеспечивает независимость работы алгоритма
от множителя на целевую функцию.
Легко видеть, что алгоритм 0( ) :r
2
0 1 2 2 1( , ) 1/ | |g g g g фактически
является r -алгоритмом с коэффициентом растяжения равным 2. Отметим, что
именно это значение рекомендуется на практике использования r -алгоритма.
Выбирая различные нормирующие множители , мы будем получать различные
алгоритмы рассматриваемого класса.
В данной работе рассматривается 1( )r -алгоритм с нормирующим множите-
лем 1 1 2 2 1( , ) 1/ (| || |).g g g g Основная особенность этого нормирующего мно-
жителя состоит в следующем. Согласно (2) значение коэффициента растяжения
определяется формулой
2 1 1 21 | | / | | | | / | | 2cos( ),g g g g (3)
где – угол между субградиентами 2 1, .g g Отсюда следует, что значение коэф-
фициента растяжения не ограничено сверху, оно может быть сколь угодно боль-
Н.Г. ЖУРБЕНКО
36 Теорія оптимальних рішень. 2017
шим. Коэффициент может принимать большие значения в случае когда норма
одного из субградиентов существенно больше нормы другого субградиента.
Пусть, например, 2 1| | | | .g g Тогда направление растяжения пространства будет
близким к направлению субградиента 2.g Поэтому операция растяжения про-
странства будет существенно уменьшать значение 2 1| | / | | .g g То есть операция
преобразования пространства будет направлена не только на уменьшение угла
между субградиентами, но и на выравнивание их норм. Следует отметить, что
это свойство может быть полезным в ситуации, когда по одним направлениям в
точке минимума производная по направлению равна нулю, а по другим направ-
лениям терпит разрыв (образно выражаясь, по одним направлениям функция
гладкая, а по другим существенно не дифференцируемая).
Численная эффективности 1( )r -алгоритма. Приведем результаты
численных исследований эффективности 1( )r -алгоритма в сравнении с
r -алгоритмом. Результаты численных исследований эффективности других
вариантов алгоритм ( )r приведены в [1, 2].
В качестве тестовых задач рассматривались задачи минимизации двух сле-
дующих функций: 1 2
1
1
( ) ,
n
i
n i
i
f x x
1
2
1
( ) | |,
n
i
n i
i
f x x
где параметр n выби-
рался в зависимости от размерности задачи n по формуле
6/( 1)10 .n
n
Таким
образом, степень вытянутости линий уровня (степень «овражности») функций
определяется значением параметра
1 610 ,n
n
она одинакова для всех функций
независимо от числа переменных. Начальная точка 1.0, 1,2,... .ix i n Критерий
останова:
610 ,kf
где kf – значение функции на итерации останова.
Результаты решения тестовых задач минимизации функций 1 2( ), ( )f x f x
приведены в табл. 1 и 2 соответственно, где приняты следующие обозначения:
1( )r – алгоритм с регулировкой шаговых множителей r-алгоритма;
*
1( )r –
алгоритм с постоянным шагом; n – размерность пространства переменных; k –
номер итерации, на которой алгоритм прекратил работу; gk – количество вы-
числений субградиента; max – максимальное значение коэффициента растяже-
ния; avrg – среднее значение коэффициента растяжения; r – результаты работы
r -алгоритма с параметрами 1 2 1.q q В алгоритмах отмеченных символом
“
**” используется постоянная величина шага (в преобразованном простран-
стве). В алгоритмах, не отмеченных символом “
**”, используется адаптивная ре-
гулировка шаговых множителей r-алгоритма.
ЧИСЛЕННАЯ ЭФФЕКТИВНОСТЬ ОДНОЙ МОДИФИКАЦИИ r-АЛГОРИТМА
Теорія оптимальних рішень. 2017 37
ТАБЛИЦА 1
Параметры
Алгоритм
n k gk
max avrg
1( )r 10 73 106 9.95 4.42
*
1( )r 10 78 79 9.39 3.94
r 10 135 168 2.00 2.00
1( )r 100 603 659 14.43 4.92
*
1( )r 100 629 630 12.36 4.83
r 100 1197 1382 2.00 2.00
1( )r 300 1775 1951 15.03 5.08
*
1( )r 300 1841 1842 15.03 4.98
r 300 3416 3898 2.00 2.00
1( )r 500 2972 3306 22.99 5.15
*
1( )r 500 3041 3042 27.43 5.00
r 500 6275 6946 2.00 2.00
ТАБЛИЦА 2
Параметры
Алгоритм
n k gk
max avrg
1( )r 10 158 305 7.86 4.15
*
1( )r 10 181 182 7.52 3.64
r 10 304 350 2.00 2.00
1( )r 100 1536 1554 5.38 4.49
*
1( )r 100 1540 1541 5.76 4.49
r 100 3256 3264 2.00 2.00
1( )r 300 4669 4678 5.29 4.57
*
1( )r 300 4672 4673 5.22 4.57
r 300 10015 10023 2.00 2.00
1( )r 500 7847 7873 5.19 4.58
*
1( )r 500 7850 7851 5.19 4.58
r 500 16850 16864 2.00 2.00
Н.Г. ЖУРБЕНКО
38 Теорія оптимальних рішень. 2017
Выводы. 1( )r -алгоритм – это модификация r -алгоритма. Вычислитель-
ная схема 1( )r -алгоритма с постоянным шагом существенно проще схемы
r -алгоритма. Величины коэффициентов растяжения пространства на итерациях
1( )r -алгоритма не постоянны, они вычисляются в процессе его работы. Алго-
ритм может использоваться с постоянным шаговым множителем. Численные
эксперименты показали достаточно высокую эффективность 1( )r -алгоритма.
М.Г. Журбенко
ЧИСЕЛЬНА ЕФЕКТИВНІСТЬ ОДНІЄЇ МОДИФІКАЦІЇ r-АЛГОРИТМУ
Розглядається модифікація r-алгоритму – алгоритму мінімізації з використанням операції роз-
тягування простору в напрямку різниці двох послідовних субградієнтів. На відміну від
r-алгоритму, значення коефіцієнтів розтягування у запропонованій модифікації програмно ро-
зраховуються в процесі роботи алгоритму. Алгоритм може використовуватися з постійним
кроком. Наводяться результати дослідження чисельної ефективності алгоритму.
N.G. Zhurbenko
NUMERICAL EFFICIENCY OF ONE MODIFICATION OF R-ALGORITHM
The modification of r-algorithm, the minimization algorithm using space dilation operation along the direc-
tion of the difference of two successive subgradients, is considered. In contrast to
r-algorithm, in the proposed modification the values of dilation coefficients are calculated during the execu-
tion of the algorithm. The algorithm can be used with a constant step size. The results of the study of the
numerical efficiency of the algorithm are presented.
1. Шор Н.З., Журбенко Н.Г. Метод минимизации, использующий операцию растяжения
пространства в направлении разности двух последовательных градиентов. Кибернети-
ка. 1971. № 3. С. 51 – 59.
2. Журбенко Н.Г. Об одной модификации r-алгоритма. Материалы III Международной
конференции. Математическое моделирование, оптимизация и информационные тех-
нологии. Кишинеу: Эврика, 2012. С. 355 – 361.
3. Журбенко Н.Г., Чумаков Б.М. Программное управление коэффициентами растяжения
r-алгоритма. Теорія оптимальних рішень. К.: Ін-т кібернетики імені В.М. Глушкова НАН
України, 2012. С. 113 – 118.
4. Шор Н.З. Методы минимизации недифференцируемых функций и их применение. Киев:
Наук. думка, 1979. 200 с.
Получено 10.03.2017
|