Научная статья на тему 'Приближенное решение эллиптических уравнений на нейронных сетях Хопфилда'

Приближенное решение эллиптических уравнений на нейронных сетях Хопфилда Текст научной статьи по специальности «Математика»

CC BY
504
126
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕЙРОННАЯ СЕТЬ ХОПФИЛДА / КРАЕВЫЕ ЗАДАЧИ ДЛЯ ЭЛЛИПТИЧЕСКИХ УРАВНЕНИЙ / УСТОЙЧИВЫЕ МЕТОДЫ РЕШЕНИЯ / HOPFIELD NEURAL NETWORK / BOUNDARY PROBLEMS FOR ELLIPTIC EQUATIONS / STABLE SOLUTION METHODS

Аннотация научной статьи по математике, автор научной работы — Бойков Илья Владимирович, Баулина Ольга Александровна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

APPROXIMATE SOLUTION OF ELLIPTIC EQUATIONS ON HOPFIELD NEURAL NETWORKS

Background. Solution of mathematical physics’ problems on artificial neural networks is an actively developing concept combining methods of calculus mathematics and computer science. Application of neural networks is especially effective in solution of reverse and incorrect problems and equations with inaccurately set parameters. At the present time the main method of solution of mathematical physics’ problems on artificial neural networks is minimization of functional error. The study is aimed at development of a stable and quick method of solving mathematical physics’ problems on artificial neural networks, based on the theory of differential equation solution stability. Materials and methods. The article describes an approximate method of elliptic equations’ solution on Hopfield neural networks. The method consists in approximation of the source boundary problem of difference scheme and formation of a system of regular differential equations, the solution of which is reduced to solution of the difference scheme. Results. The authors suggest a method of boundary problem solution for linear and non-linear elliptic equations, based on the methods of the stability theory. Effectivenes of the method is demonstrated by the model examples. Conclusions. The results of the study may be used for solution of a wide class of boundary problems for linear and non-linear elliptic equations, determined in sectionally smooth areas.

Текст научной работы на тему «Приближенное решение эллиптических уравнений на нейронных сетях Хопфилда»

УДК 004.02

И. В. Бойков, О. А. Баулина

ПРИБЛИЖЕННОЕ РЕШЕНИЕ ЭЛЛИПТИЧЕСКИХ УРАВНЕНИЙ НА НЕЙРОННЫХ СЕТЯХ ХОПФИЛДА

Аннотация.

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

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

Результаты. Предложен метод решения краевых задач для линейных и нелинейных эллиптических уравнений, основанный на методах теории устойчивости. Эффективность метода проиллюстрирована модельными примерами.

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

Ключевые слова: нейронная сеть Хопфилда, краевые задачи для эллиптических уравнений, устойчивые методы решения.

I. V. Boykov, O. A. Baulina

APPROXIMATE SOLUTION OF ELLIPTIC EQUATIONS ON HOPFIELD NEURAL NETWORKS

Abstract.

Background. Solution of mathematical physics’ problems on artificial neural networks is an actively developing concept combining methods of calculus mathematics and computer science. Application of neural networks is especially effective in solution of reverse and incorrect problems and equations with inaccurately set parameters. At the present time the main method of solution of mathematical physics’ problems on artificial neural networks is minimization of functional error. The study is aimed at development of a stable and quick method of solving mathematical physics’ problems on artificial neural networks, based on the theory of differential equation solution stability.

Materials and methods. The article describes an approximate method of elliptic equations’ solution on Hopfield neural networks. The method consists in approximation of the source boundary problem of difference scheme and formation of a system of regular differential equations, the solution of which is reduced to solution of the difference scheme.

Results. The authors suggest a method of boundary problem solution for linear and non-linear elliptic equations, based on the methods of the stability theory. Effec-tivenes of the method is demonstrated by the model examples.

Conclusions. The results of the study may be used for solution of a wide class of boundary problems for linear and non-linear elliptic equations, determined in sectionally smooth areas.

Key words: Hopfield neural network, boundary problems for elliptic equations, stable solution methods.

Применение нейронных сетей Хопфилда [1, 2] для решения задач математической физики основано на возможности представления нейрона в виде электронной схемы, описываемой нелинейным обыкновенным дифференциальным уравнением. Согласно этому представлению i -й нейрон, соединенный с N нейронами сети (включая самого себя), описывается системой уравнений

du^ u^ v ^

Ci~dt = - R+(u j)+^, i = 1,2,---,N,

i j=i

где Wj - синаптические веса нейронов сети; Ij - ток, представляющий внешнее смещение; ui - индуцированное локальное поле на входе функции активации f (ui); f (ui) - нелинейные функции активации; Ri и Q -сопротивление утечки и емкость утечки соответственно.

Опишем архитектуру нейронной сети Хопфилда, используемой в данной работе. Предлагаемая сеть, состоящая из n нейронов, показана на рис. 1.

В нейронную сеть входят нелинейные устройства, реализующие нелинейные функции f (xi, Х2,..., xn), i = 1, 2,.., m.

Выходные сигналы f■ (xi,X2,...,xn), i = 1,2,..,n, суммируются в блоке сумматора с синаптическими коэффициентами w-, i = 1,2, ..,n, j = 1,2,..,m, и после прохождения RC цепочки подаются на вход устройства активации, реализующего функцию x = 9(u).

В данной архитектуре используется функция 9(u) = au. Таким образом, представленная на рис. 1 нейронная сеть Хопфилда реализует систему нелинейных обыкновенных дифференциальных уравнений

Ci dxi (t) + xi (t) = ^

j, и Lwv4-Hj a dt aR,■

+ ~br = ^WUfU(xb • • •, xn) + Ii, i = 1,2,., n. 'г' j=1

Отметим, что в качестве сопротивлений Я могут браться достаточно большие значения, а также Iу могут полагаться равными нулю. В результате нейронные сети Хопфилда могут моделировать системы уравнений вида

Су ах( (о т г ( ч ■ 1 2

----(x1, -, хп X ■ =1,2, -,п

}=1

Рис. 1. Архитектура нейронной сети

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

Ниже используются следующие обозначения:

Я(а,г) = {г є В: ||г - а|| < г}, ^(а,г) = {г є В : ||г - а|| = г},

Л(К) = Нш(( + ИК\\ -1)-1.

НІ0

Здесь В - банахово пространство; а є В; К - линейный оператор, действующий из В в В; Л(К) - логарифмическая норма [6] оператора К; через I обозначен тождественный оператор.

im- |л(('(ф(т))) т = -а, а > 0.)

1. Непрерывные методы решения операторных уравнений

Рассмотрим нелинейное операторное уравнение

A( х) = 0, (1)

действующее из банахова пространства B в B.

Рассмотрим в банаховом пространстве B задачу Коши

^ = A( x(t)); (2)

dt

х(0) = хо. (3)

Пусть оператор A имеет непрерывную производную Гато; A(0) = 0. Теорема 1 [7, 8]. Пусть на любой дифференцируемой кривой ф(0, расположенной в шаре B(0, r) достаточно малого радиуса г, интеграл

t

|л(('(ф(т))) т

0

не положителен (отрицателен и

1t

lim

t ——^ t 0

Тогда тривиальное решение уравнения (2) устойчиво (асимптотически устойчиво).

Замечание. Теорема справедлива и при г = «>.

Из теоремы 1 следует, что если для любой непрерывно дифференцируемой функции g (t), определенной в банаховом пространстве B, выполняется неравенство

lim - ]Л(A'(g (T))dт<-а, а >0, (4)

t—” 0

*

то задача Коши (2)-(3) сходится к решению х уравнения (1).

Таким образом, справедливо следующее утверждение.

-1—г * -1—г

Теорема 2 [9]. Пусть уравнение (1) имеет решение х . Пусть на любой дифференцируемой кривой g (t), расположенной в банаховом пространстве B, справедливо неравенство (4). Тогда решение задачи Коши (2)-(3)

*

сходится к решению х уравнения (1) при любом начальном приближении.

Замечание. Из неравенства (4) следует, что логарифмическая норма Л( A'( х)) может обращаться в нуль или принимать положительные значения в конечном или счетном числе точек пространства B.

*

Теорема 3 [9]. Пусть уравнение (1) имеет решение х . Пусть на любой

*

дифференцируемой кривой g(t), расположенной в шаре B(х , г)

выполняются следующие условия:

I

1) при любом t (t > 0) выполняется неравенство Ja(A'g(T))dT< 0;

0

2) справедливо равенство lim 1 Jл(A'(g(x))dт = -а, а >0.

t^^10

*

Тогда решение задачи Коши (2)-(3) сходится к решению x уравнения (1).

Решение линейных эллиптических уравнений на нейронных сетях Хопфилда

Пусть D = [0,1; 0,1]. В области D рассмотрим задачу Дирихле

a( xi, Х2 )Au (xi, Х2) + b( Xi, X2)du( x1, x2) +

dxi

+c( x1, x2)du (xb x2) + d (x1, x2 )u( x1, x2) = f (x1, x2) (5)

dx2

с граничными условиями

u( xb x2)|r = U0(x1, x2). (6)

Здесь Г = dD.

Замечание 1. Двумерные уравнения рассматриваются для простоты обозначений.

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

Введем сетку узлов (у^ , vi), k, l = 0,1,..., N, где = к / N,

к = 0,1,..,N.

Во внутренних узлах у, V[), к, l = 0,1,., N -1, оператор Au

аппроксимируется пятиточечной разностной схемой

u(vk-ь vi) + u (vk+Ь vl) - 2u (vk, vl) + u (vk, vl+1) + u(vk, vl-1) - 2u (vk, vl) h2 h2 ’

du du

а операторы ------ и ----- аппроксимируются двухточечными разностными

3x1 3x2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

схемами

Фк+Ъ vl)- u (vk, vl) u (vk, v/+l)- u(vk, vl)

---------------------- и------------------------,

где h = 1 / N.

h

h

В результате задача Дирихле (5), (6) аппроксимируется разностной схемой

(

a (vk, vl)

и (vk-1, Vl) + и (vk+1, Vl) - 2u (vk, Vl) +

h2

+

u (Vk , Vl+1) + U(Vk , Vl-1) - 2u(Vk , Vl ) 1 + w v ) u(vk+1, Vl ) - U (Vk , Vl )

2

+ b(vk, vl)-------------1;----------2— + c(vk, vl) x

х и (Ч,^+0 и(Ч,У1) + а( ,У/)ы^^) - / V,VI) = 0, к,/ = 1,2,...,N -1. (7) п

Замечание. Из условия (6) следует, что

и ^0, V/ )= моО>о, V/), и (V N, V/ )= ыоО>м, V/),

и ^к ,г;0)= ы0^к ,v0), Фк, VN )= и0 0к, VN).

Представим систему уравнений (7) в виде

Фк, V)(ы (г;к+1, V) + ы (г;к-1, V) + Фк, V+1) + и ^к, V -1) - 4ы ^к, V)) + +пФк, V)(ы ^к+1, V) - и ^к, V)) + пФк, V)(ы (г;к, V+1) - и ^к, V)) +

+И2МV, V/)ыV, V/) - Н2/^к, V/) = 0, к,/ = 1,2,.,N -1. (8)

Системе уравнений (8) поставим в соответствие систему диффе-

ренциальных уравнений

--= ак/ (Фк, V )(ык+1,/(^) + ик-1,/(^) + ик,/+1(Ґ) + ик,/-1() - 4ык,/(ґ)) +

Ш

+пЪ^к, V)(ик+1,/^) - ик/(()) + ПФк, V )(u(vk, V+1) - и(vk, V )) +

+П2МV, V/)и^к, V/) - Н2/^к, V/)), к,/ = 1,2,.,N -1. (9)

Здесь ак/ - константы, выбор которых будет описан ниже,

ик,/ ^ ) = и^к, V, ').

Пусть система уравнений (8) однозначно разрешима. Обозначим через

*

и (Ук, VI), к, / = 0,1,2,., Ы, решение этой системы и введем обозначения

wк/(0 = ик/(0 -и*/, и*/ = и0к,VX k,/ Ы.

Тогда система уравнений (9) принимает вид

--= ак/(а^к, V)(^+1,/() + wk-1,/(() + wk,/+1(() + wk,/-1(() - 4^к,/(()) +

М

+ЪЪ (vk, V )(wk+1,/ ^) - wk/ (()) + пфк, V )(w(vk, V+1) - Фк, V )) +

+п2МV,V/)wk,/ (Г)), к,/ = 1,2,.,N -1. (10)

Замечание. Очевидно, что

*0,1(0 = *N,1 = *к,0(0 = *кN(0 = 0, к ^Л..^N.

Теорема 4. Пусть существует функция а(^,¿2) такая, что а(*1, ¿2 )а(Ь, ¿2) ^ 0, а(^, ^2)* (¿1, ¿2) ^ 0, а(^,¿2)Фъ¿2) ^ 0, а(^,^)ё(¿1,¿2) < 0 (11)

при всех ¿1, ¿2 £ В. Тогда система уравнений (10) устойчива в целом.

Доказательство. Представим систему уравнений (10) в операторной форме:

&^ = АМ,), (12)

где м!^) - вектор функция,

* (Г) = { *11 (Г),..., -1 (Г),..., ^-l,N -1 ^);

матрица А строится очевидным образом.

Нетрудно видеть, что при выполнении перечисленных условий (11)

2

логарифмическая норма матрицы А равна Л(А) = тах кёу,у)<0.

1<к ,1 < N-1

Пусть функция ё(Х1, Х2) < 0 в области В. Следовательно

2

ё (Х1, Х2) <- <0 в этой области, тогда Л(А) < -вк .

Из этого условия следует, что решение системы уравнений (12) при , абсолютно устойчиво при любых начальных условиях м>к1 (0),

к, 1 = 1,2,...,N -1, и при единственном ограничении:

*0,1 (0) = WN,/ (0) = *к ,0(0) = *к, N (0) = 0, к, 1 = 0,1,..., N.

Таким образом, при указанном ограничении система уравнений (12) устойчива в целом. Отсюда следует, что система уравнений (9) устойчива в целом при любых начальных возмущениях, не нарушающих граничных значений.

Поэтому решение системы уравнений (9) сходится при í

*

к решению и (ук,VI), к, 1 = 1,2,.,N -1, системы уравнений (7).

Теорема доказана.

Из доказательства теоремы 4 следует, что при í решение системы дифференциальных уравнений (9) стремится к решению системы алгебраических уравнений (7).

Решение нелинейных эллиптических уравнений на нейронных сетях Хопфилда

2

Рассмотрим в области В = [0,1] систему нелинейных уравнений

a( X, »2) Si (Аи (xl, х2>) + b( X, »2) §2

Эх

+

/

+c( хъ Х2) §3

du( xj, Х2) дх2

\

+ d (xj, Х2) §4 (u( xj, x2)) = f( xj, Х2) (13)

V 2 /

с граничными условиями

и ( *1, Х2)|г = «о( *1, *2).

(14)

Здесь Г = ЭВ; gi(*), i = 1,2,3,4 - непрерывно-дифференцируемые функции на Я.

По аналогии с линейным случаем введем сетку узлов V, V/), к,/ = 0,1,., N и поставим граничной задаче (13), (14) в соответствие систему разностных уравнений:

a(vk, vl) §1

u (vk -l, vi) + u (vk+1, vl) - 2u (vk, vl) h2

+

+

u (vk, vl+1) + u (vk, vl-1) - 2u (vk, vl) h2

+

+4(Vk, „ )g2 I u(Vk+1'Vl)~ u(Vk •Vl) 1 + c(Vk, v, )§3 i u(Vk •Vl+1) - u(Vk •Vl) ' I +

+d (vk, v, ) g4(u (vk, v, )) = f(vk, vl), k, l = 1,2,..., N -1.

(15)

Системе уравнений (15) поставим в соответствие систему нелинейных дифференциальных уравнений

^ ик-1,1(() + ик+1,1(() + ик,/+1(() + ик,/ -1(() - 4ик/(() ^

dt

+

+b(vk, v, ) § 2

uk+1,l(t) - ukl(t)

+ c(vk, vl) §3

ukl+1(t) - ukl(t)

+

Здесь

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+^V,VI),§4 (ик,/ (0) - /^к,VI)), к,/ = 1,2,.,N -1. (16)

ик1 (0 = иV,V/, 0; к,/ = 0,1,., N; ио1 (0 = ио(0, V/), / = 0,1,., N,

UkN (0 = и0 V ,1), к = 0,1,., N,

UN/(0 = и0(1,V/), I = 0,1,.,N, ик0(0 = и^к ,0), к = 0,1,., N.

*

Предположим, что система уравнений (15) имеет решение и (ук,V/), к, / = 0,1,., N. Введем новую функцию w(vk, VI, ^), положив

и(Ук,VI,0 = ,VI, О + и (у*,VI).

Тогда систему уравнений (16) можно представить в виде

dwk/(t) _

dt

= ak/ (a (Vk , V/) g

wk+1,/(t) + wk-1,/(t) + wk,/+1(t) + wk ,/-1(t) - 4 wk/(t) h2

+

* * * *

+

uk+1,/ + uk-1,/ + uk,/+1 + uk ,/-1 - 4uk/

- a(vk, v/) g1

( 4+11 + u*-1,/ + ukj+1 + ukj-1(t) - 4u*/ ^

+

+b(vk, v/) g2

( wk+1,/(t) - wk,/(t) . uk+1,/ - uk,/ ^

+

-b(vk, v/) g2

/ * * \ uk+1,/- uk,/

+

+c(vk, v/) g3

k

wk,/+1(t) - wk,/(t) . uk,/+1- uk,/

+

■c(vk, v/) g3

f * * \ uk,/+1- uk,/

+

+й(vk,V)^4(^*,/ (ї) + и*,/) - й(vk,v/)£4(и*1)), к,I = 1,2,...,N -1 (17)

В первом приближении система уравнений (17) имеет вид

dwk/ (t) , , ч ,

—— = аы (а (vk, v/) g1 dt

( u*+1,/ + u*-1,/ + u*,/+1 + u*,/-1 - 4u*/ ^

x

x

wk+1,/(t) + wk-1,/(t) + wk,/+1(t) + wk,/-1(t) - 4 wk/(t) ^

+

+b(vk, v/) g2

+c(vk, v/) g3

( * * \

uk+1,/- uk,/

h

V У

/ * * \

uk,/+1- uk,/

wk+1,/(t) - wk ,/(t)

+

wk,/+1(t) - wk ,/(t)

+

+й (V*, V/) £4 (и*,[ )щ ,/ (ї)), к ,1 = 1,2,..., N -1.

Исследуем устойчивость решения системы уравнений (18).

п2

(18)

Теорема 5. Пусть В = [0,1] . Пусть а(л^, Х2) - функция, определенная на В, а(Х1, Х2) Ф 0 при (Х1, Х2) є В, и такая, что при любых (Х1, Х2) є В, и любом натуральном N, N > 2, выполняются условия:

1) а( xb x2)a( xi, Х2) g[

( uk+U + uk-U + ukJ+1 + ukJ-1 - 4ukl ^

>0;

2) а(xi,X2)b(xi,Х2)g2

( k k \ uk+1,l - uk ,l

> 0;

3) а(x1, x2)c(x1, x2)g3

( k k \ uk ,l+1 - uk ,l

> 0;

4) а(X!, Х2 )й(X!, Х2)g4 (и*,/) < “Р <0

при 1 < к, I < N -1.

Тогда решение системы уравнений (18) устойчиво в целом. Доказательство. Представим систему уравнений (18) в операторной форме:

^ = Bw(,),

dt

(19)

где

^(г),..., ^1,N-1 (г), М>2^),..., W2,N-1 (г),..., ^N-1,1^),•••, ^N-1,N-1 (')};

построение матрицы В очевидно.

Нетрудно видеть, что при выполнении условий теоремы логарифмическая норма матрицы В отрицательна. Тогда из неравенства Винтнера (см. [6; 8, с. 51] следует устойчивость в целом системы уравнений (18). Теорема доказана.

Устойчивость в целом решения системы уравнений (18) дает основание полагать, что решение системы уравнений (16) будет устойчиво относительно

*

вектора Ык1, к, I = 1,2,., N -1. Для строгого доказательства устойчивости

*

решения системы уравнений (16) относительно вектора и^, к, I = 1,2,., N -1, нужно использовать критерии, предложенные в работах [7, 8]. В применении к системе уравнений (16) эти критерии оказываются трудно обозримыми. Поэтому при решении нелинейных эллиптических уравнений на нейронных сетях Хопфилда достаточно ограничиться критериями теоремы 5. Выполнение этих критериев дает достаточные основания для проведения вычислений.

Эффективность предложенного метода проиллюстрирована на следующих примерах.

2

Пример 1. Рассмотрим уравнение Ли = 0 в области Б = [0,1] с граничными условиями и0(х, у):

прямая АВ = {у = 0, 0 < х < 1} иАв = х; прямая АБ = {х = 0, 0 < у < 1} иАБ = 0; прямая СБ = {у = 1, 0 < х < 1} исБ = х; прямая ВС = {х = 1, 0 < у < 1} ивс = 1.

Точное решение и(х, у) = х.

Задача Дирихле была приведена к системе уравнений (9), которая решалась методом Эйлера.

Результаты решения задачи Дирихле приведены на рис. 2-5. Здесь т -число итераций метода Эйлера; N - число узлов разностной схемы по одной переменной; И - шаг метода Эйлера.

т, итер.

Рис. 2. Результаты решения задачи с шагом в методе Эйлера И = 0,1

т, итер

Рис. 3. Результаты решения задачи с шагом в методе Эйлера И = 0,01

Рис. 4. Результаты решения задачи с шагом в методе Эйлера И = 0,001

т,итер

Рис. 5. Результаты решения задачи с шагом в методе Эйлера И = 0,0001

2

Пример 2. Рассмотрим уравнение Л^ - ^ = 0 в области Б = [0,1] при граничных условиях ^0(х, у):

прямая АВ = {у = 0, 0 < х < 1} wAB = х + 1; прямая АБ = {х = 0, 0 < у < 1} wAD = еЬ у + у; прямая СБ = {у = 1, 0 < х < 1} wCD = (х + 1)(еЬ(1) + 8Ь(1)); прямая ВС = {х = 1, 0 <у < 1} wвс = 2(еЬ у + у);

Точное решение - w(x, у) = (х + 1)(еЬ у + у).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Краевая задача была приведена к системе уравнений (9), которая решалась методом Эйлера.

Результаты решения уравнения приведены на рис. 6-8. Здесь т - число итераций метода Эйлера; N - число узлов разностной схемы по одной переменной; И - шаг метода Эйлера.

пп., итер

Рис. 6. Результаты решения задачи с шагом в методе Эйлера И = 0,1

пп., итер

Рис. 7. Результаты решения задачи с шагом в методе Эйлера И = 0,01

Заключение

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

Результаты работы могут быть использованы при решении широкого класса краевых задач для линейных и нелинейных эллиптических уравнений, определенных в кусочно-гладких областях.

од

0,01

| 0,001 X ]

0,0001

о

Е

0,00001

0,000001

пп, итер

Рис. 8. Результаты решения задачи с шагом в методе Эйлера h = 0,001 Список литературы

1. Hopfield, J. J. Neural Networks and Physical Systems with Emergent Collective Computational Abilities / J. J. Hopfield // Proc. Natl. Acad. Sci. USA. - 1982, April. -Vol. 79. - P. 2554-2558.

2. Hopfield, J. J. Neurons with Graded Response have Collective Computational Properties like those of Two-State Neurons / . J. Hopfield // Proc. Natl. Acad. Sci. USA. - 1984, May. - Vol. 81. - P. 3088-3092.

3. Горбаченко, В. И. Нейрокомпьютеры в решении краевых задач теории поля / В. И. Горбаченко. - М. : Радиотехника, 2003. - 336 с.

4. Бойков, И. В. Приближенное решение задач математической физики на нейронных сетях Хопфилда / И. В. Бойков, В. А. Руднев, А. И. Бойкова // Нейрокомпьютеры, разработка, применение. - 2013. - № 10. - С. 13-22.

5. Бойков, И. В. Приближенное решение гиперсингулярных интегральных уравнений на нейронных сетях Хопфилда / И. В. Бойков, О. А. Баулина // Журнал Средневолжского математического общества. - 2013. - Т. 15, № 1. - С. 41-51.

6. Далецкий, Ю. Л. Устойчивость решений дифференциальных уравнений в банаховом пространстве / Ю. Л. Далецкий, М. Г. Крейн. - М. : Наука, 1970. -534 с.

7. Бойков, И. В. Об устойчивости решений дифференциальных и разностных уравнений / И. В. Бойков // ДАН СССР. - 1990. - Т. 314, № 6. - С. 1298-1300.

8. Бойков, И. В. Устойчивость решений дифференциальных уравнений / И. В. Бойков. - Пенза : Изд-во ПГУ, 2008. - 244 с.

9. Бойков, И. В. Об одном непрерывном методе решения нелинейных операторных уравнений / И. В. Бойков // Дифференциальные уравнения. - 2012. -Т. 48, № 9. - С. 1308-1314.

10. Гюнтер, Н. М. Теория потенциала и ее применение к основным задачам математической физики / Н. М. Гюнтер. - М. : ГИТТЛ, 1953. - 415 с.

References

1. Hopfield J. J. Proc. Natl. Acad. Sci. USA. 1982, April, vol. 79, pp. 2554-2558.

2. Hopfield J. J. Proc. Natl. Acad. Sci. USA. 1984, May, vol. 81, pp. 3088-3092.

-N=10

■-N =30

3. Gorbachenko V. I. Neyrokomp’yutery v reshenii kraevykh zadach teorii polya [Neurocomputers in solution of boundary problems of the field theory]. Moscow: Radio-tekhnika, 2003, 336 p.

4. Boykov I. V., Rudnev V. A., Boykova A. I. Neyrokomp’yutery, razrabotka, primenenie [Neurocomputers, development, application]. 2013, no. 10, pp. 13-22.

5. Boykov I. V., Baulina O. A. Zhurnal Srednevolzhskogo matematicheskogo obshchestva [Journal of Middle Volga mathematical society]. 2013, vol. 15, no. 1, pp. 41-51.

6. Daletskiy Yu. L., Kreyn M. G. Ustoychivost’ resheniy differentsial’nykh uravneniy v banakhovomprostranstve [Stability of differential equations solution in Banach space]. Moscow: Nauka, 1970, 534 p.

7. Boykov I. V. DANSSSR [Reports of the USSR Academy of Sciences]. 1990, vol. 314, no. 6, pp. 1298-1300.

8. Boykov I. V. Ustoychivost’ resheniy differentsial’nykh uravneniy [Stability of differential equations solutions]. Penza: Izd-vo PGU, 2008, 244 p.

9. Boykov I. V. Differentsial’nye uravneniya [Differential equations]. 2012, vol. 48, no. 9, pp. 1308-1314.

10. Gyunter N. M. Teoriya potentsiala i ee primenenie k osnovnym zadacham matematich-eskoy fiziki [Theory of potentials and application thereof to main problems of mathematical physics]. Moscow: GITTL, 1953, 415 p.

Бойков Илья Владимирович

доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)

E-mail: math@pnzgu.ru

Баулина Ольга Александровна аспирант, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)

E-mail: math@pnzgu.ru

Boykov Il'ya Vladimirovich Doctor of physical and mathematical sciences, head of sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street,

Penza, Russia)

Baulina Ol'ga Aleksandrovna Postgraduate student, Penza State University (40 Krasnaya street, Penza, Russia)

УДК 004.02 Бойков, И. В.

Приближенное решение эллиптических уравнений на нейронных сетях Хопфилда / И. В. Бойков, О. А. Баулина // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. -№ 1 (29). - С. 39-53.

i Надоели баннеры? Вы всегда можете отключить рекламу.