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

Для численной реализации системы линейных алгебраических

уравнений (6), (7) воспользуемся итерационным методом сопряжен-

ных градиентов, позволяющим уточнить смещение струны в началь-

ный момент времени. На каждой итерации решается прямая задача

для уравнения колебаний струны с помощью стандартной симметрич-

ной трехслойной разностной схемы, обладающей вторым порядком

аппроксимации по

h

и

τ

[7]. Итак, пусть вместо обратной задачи

рассматривается прямая задача для этого же уравнения, когда второе

условие в (5) заменено начальным условием

y

(

x,

0) =

ϕ

(

x

)

,

x

ω

h

.

Дискретный аналог прямой задачи имеет вид

ˆ

y

2

y

+ ˘

y

τ

2

+

A

(

σ

ˆ

y

+ (1

2

σ

)

y

+

σ

˘

y

) = 0

, t

ω

τ

;

(8)

y

(0) =

ϕ, y

(

τ

) =

y

(0) +

τν, x

ω

h

.

(9)

Разностная схема (8), (9) устойчива при выполнении условия [7]

σ

1

4

1

τ

2

k

A

k

(10)

и справедлива следующая априорная оценка:

ˆ

Y

(

t

) =

k

Y

(

t

+

τ

)

k

≤ k

Y

(

t

)

k

. . .

≤ k

Y

(

τ

)

k

, t

ω

τ

,

(11)

где

ˆ

Y

— норма,

ˆ

Y

2

=

1

4

k

ˆ

y

+

y

k

2

A

+

τ

2

k

y

t

k

2

R

(1

/

4)

A

,

R

=

1

4

E

+

σA.

Тем самым норма решения задачи Коши со временем убывает. Отме-

тим, что для явной разностной схемы (

σ

= 0

) условие устойчивости

(10) принимает вид

τ

h/

k

2

.

Итерационный метод.

Для нахождения решения сеточной задачи

(6), (7) используем наиболее быстро сходящийся итерационный ме-

тод вариационного типа — метод сопряженных градиентов [8], осно-

ванный на последовательном уточнении искомого начального условия

ν

(

x

)

с дальнейшим решением на каждой итерации дискретного ана-

лога прямой задачи (8), (9). Придадим этой задаче соответствующую

операторную формулировку. Последовательно исключая промежуточ-

ные значения сеточной функции

y

(

t

)

, t

ω

τ

при заданных значениях

ϕ

и

ν

на конечный момент времени, получаем

y

(

t

n

) = A

ϕ

+B

ν,

где

A

,

B

— операторные полиномы

n

-й и

(

n

1)

-й степеней от положитель-

но определенного и самосопряженного оператора

A

. С учетом этого

приближенному решению обратной задачи естественно сопоставить

решение следующего операторного уравнения:

A

ϕ

=

ϕ

(

x

)

B

ν

= ˉ

ϕ.

В силу самосопряженности оператора

A

самосопряженным является

и операторный полином

A

.

80

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