Background Image
Previous Page  5 / 11 Next Page
Information
Show Menu
Previous Page 5 / 11 Next Page
Page Background

Для решения операторного уравнения используем трехслойный

итерационный метод сопряженных градиентов, записанный в кано-

ническом виде [8]

ϕ

k

+1

=

α

k

+1

ϕ

k

+ (1

α

k

+1

)

ϕ

k

1

α

k

+1

β

k

+1

r

k

, k

= 0

,

1

, . . . ,

(12)

где

r

k

— невязка,

r

k

= A

ϕ

k

ϕ

. Для вычисления итерационных пара-

метров

β

k

+1

, α

k

+1

используются следующие рекуррентные формулы

β

k

+1

=

(

r

k

, r

k

)

(A

r

k

, r

k

)

,

k

= 0

,

1

, . . .

;

(13)

α

k

+1

= 1

β

k

+1

β

k

(

r

k

, r

k

)

(

r

k

1

, r

k

1

)

1

α

k

1

, k

= 1

,

2

, . . . , α

1

= 1

.

(14)

В этом итерационном методе вычисление вектора

y

n

k

осуществля-

ется численной реализацией прямой задачи

ˆ

y

k

2

y

k

+ ˘

y

k

τ

2

+ A(

σ

ˆ

y

k

+ (1

2

σ

)

y

k

+

σ

˘

y

k

) = 0

, t

ω

τ

,

(15)

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

y

0

k

=

ϕ

k

, y

1

k

=

ϕ

k

+

τν, x

ω

h

.

(16)

Следовательно, невязка определяется по формуле

r

k

=

y

n

k

ϕ

,

x

ω

h

.

Аналогично выполняется вычисление вектора

z

k

= A

r

k

:

ˆ

z

k

2

z

k

+ ˘

z

k

τ

2

+

A

(

σ

ˆ

z

k

+ (1

2

σ

)

z

k

+

σ

˘

z

k

) = 0

,

(

x, t

)

ω

,

(17)

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

z

0

k

=

r

k

, z

1

k

=

r

k

, x

ω

h

.

(18)

Таким образом, в случае неклассической задачи для уравнения коле-

баний струны итерационный метод сопряженных градиентов реализу-

ется в следующем порядке.

1. Полагаем

k

= 0

и задаем начальное приближение искомой функ-

ции

ϕ

k

, x

ω

h

.

2. Запускаем счетчик итераций

k

=

k

+ 1

, последовательно решая

прямые задачи (15), (16) и (17), (18), определяем невязку

r

k

=

y

m

k

ϕ

,

x

ω

h

, и вектор

z

k

= A

r

k

,

x

ω

h

.

3. По рекуррентным формулам (13), (14) находим значения итера-

ционных параметров

β

k

+1

, α

k

+1

.

4. По формуле (12) рассчитываем очередное приближение искомого

начального условия

ϕ

(

x

)

k

+1

,

x

ω

h

.

5. Процесс продолжается до тех пор, пока не выполнится критерий

выхода из итерационного цикла

k

r

k

k

< ε

, иначе продолжаем процесс,

возвращаясь к пункту 2.

ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2015. № 3

81