Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
2343.pdf
Скачиваний:
43
Добавлен:
07.01.2021
Размер:
5.89 Mб
Скачать

 

 

 

7. Задачи для самостоятельного решения

 

 

 

 

 

 

 

1.

Найти

методом

сеток

приближенное

решение

уравнения

 

U

2U

, удовлетворяющего условиям;

 

 

 

 

 

 

 

 

 

 

 

 

t

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, 0) sin x при 0 x 1;

 

 

 

 

 

 

 

 

 

 

 

 

 

U (0, t) U (1, t) 0 , если

0 t 0,025 .

 

 

 

 

 

 

2. Методом сеток найти решения задачи

 

 

 

 

 

 

 

 

 

 

 

 

С

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U U ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, 0) 0,2 sin (1 x)x;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ut (x, 0) 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (1, t) 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, t)

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Методом сеток найти с точностью приближенное решение урав-

 

 

 

 

 

 

 

 

 

 

 

, где U (x, y) определена внутри прямо-

 

нения Лапласа U xx

U yy 0

 

 

 

 

 

б

y d , если на границе Г

 

угольника GABCD :

G : (x, y), a

x b, c

 

области

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г x a ;

x b ; c y d ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аy c ; y d ; a x b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

известны граничные условия U (x, y) f (x, y) ,

(x, y) Г .

 

 

 

 

 

 

Положить M 6; N 4 ;

0,001 . Значения функций и парамет-

 

ров взять из табл. 33.

 

 

 

 

Д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 33

 

 

а

 

b

 

c

 

 

d

 

U

 

AC

 

U

 

BD

 

U

 

AB

 

U

 

CD

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п/п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

4

 

10

 

–1

 

1

 

3(1–у)

 

x–2y

 

 

x

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

2

 

–1

 

2

 

0

 

 

2

 

5siny

 

xy

 

 

0

 

x+y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

9

 

6

 

–1

 

1

 

4 y 2

 

3 y

 

 

1

 

2x

 

 

4

 

–2

 

1

 

–1

 

0

 

5y

 

xy

 

 

0

 

x 2

 

 

5

 

–1

 

2

 

1

 

 

3

 

x+y

 

y 2

 

 

1

 

1–x

 

196

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 33.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

b

 

c

 

d

 

U

 

AC

 

U

 

BD

 

U

 

AB

U

 

CD

 

 

 

 

 

 

 

 

 

 

 

 

п/п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

0

1

–1

 

1

 

 

2y

 

y2 1

 

 

0

 

xy

 

 

7

–1

0

 

0

 

4

 

xy

 

1+x

 

 

1

 

x+y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 1

 

8

–1

2

 

0

 

2

 

siny

 

 

0

 

 

 

1

 

 

9

0

6

–1

 

1

 

1+y

 

1 y2

 

 

x

 

 

0

 

10

1

4

 

0

 

1

 

2y–1

 

xy

 

1–x

 

 

1

 

 

и

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

11

2

3

 

1

2

 

y–1

 

 

 

 

 

 

0

 

xy

 

12

4

7

–1

0

 

 

 

y

 

1–x

 

 

0

 

 

x

 

13

0

3

 

0

1

 

2x+y

 

 

0

 

 

sinx

x2 1

 

 

 

 

б

 

y2 x

 

 

 

 

 

 

 

 

 

14

0

1

–1

1

 

2+y

 

 

cosx

 

1–x

 

15

–1

2

 

0

1

 

 

 

0

 

x+cos y

 

 

1

x2 y2

 

16

–1

0

 

0

2

 

 

 

1

 

xy 2

 

 

sinx

 

yх

 

 

 

 

 

 

 

А

 

 

 

 

 

 

 

 

 

17

–1

2

 

0

2

 

siny

 

1+y

 

 

x

 

 

1

 

18

0

3

–4

0

 

1–y

 

1 y2

 

 

x

 

 

0

 

19

1

4

 

0

4

 

 

 

y

 

1+x

 

 

0

 

 

x

 

20

–1

2

 

0

4

 

 

 

y

 

 

0

 

 

1+x

ysin x

 

 

 

 

 

 

 

 

 

 

 

 

 

Д

 

 

 

 

 

21

0

2

–1

3

 

 

2y

 

1+x

x2

 

x+y

 

22

0,5

1

 

0

2

 

y(1–y)

 

sin y

 

1

 

 

x

 

23

–2

0

 

1

5

 

 

2y

 

2 cos y

 

0

 

x+y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

24

0

1

 

0

2

 

 

 

0

 

2xy 2

 

x3

 

 

y

 

25

–1

2

 

0

2

 

3sin y

5y 1 y2

 

0

 

 

1

 

26

–2

1

–1

0

 

xy

 

sin

 

 

x

 

2

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

27

–1

1

 

0

3

 

sin

 

x

 

x–2

 

3

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

28

2

4

–1

2

cos

x 1

 

x+2

 

x

 

x+4

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

29

0

3

 

2

4

 

1 y

 

x y

2

 

 

x

 

x y

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

–2

0

–2

0

 

x y2

 

y2 x

 

0

x2 y2

 

197

4. Применяя метод сеток с шагом h 0,1, найти решение уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Utt

 

U xx , удовлетворяющее условиям

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, 0) f (x) ;

Ut (x, 0)

Ф(x) ;

U (0, t) (t) ; U (1, t) Ψ (t)

для

 

0 t 0,5; 0 x 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вариантов 1– 10:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) ax2

1,1 sin x ; Ф(x) 0,

(t) (t) 0.

 

 

 

 

 

 

 

 

1.

a 1,1;

 

2.

a 1,6;

3.

a 1,2;

4.

a 1,7;

 

5. a 1,3;

 

 

 

 

 

6.

a 1,8;

7.

 

a 1,4;

8.

a 1,9;

9.

a 1,5;

 

10. a 2,0.

 

 

 

 

 

Для вар антов 11– 20 f (x) задана табл. 34.

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 34

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х

 

 

0

 

0,1

 

0,2

 

 

 

0,3

 

0,4

 

0,5

 

 

 

0,6

 

 

0,7

 

 

0,8

0,9

 

1,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x)

0

 

0,0145

0,0511

0,0921а

0,1114

0,1825

 

0,1905

 

0,148

 

 

0,1025

0,0502

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11.

 

a 1,02;

 

 

12.

a 0,975;

13.

a 1,0;

14.

 

a 1,025;

 

 

 

 

 

15.

 

a 1,5;

 

 

 

16.

a 1,03;

17.

a 1,045;

18.

 

a 1,06;

 

 

 

 

 

19.

 

a 0,9;

 

 

 

20.

a 0,995.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

А

 

 

 

 

 

 

5. Методом прогонки решить уравнение теплопроводности Ut U xx ,

 

удовлетворяющее

 

 

условиям

 

 

 

U (x, 0) f (x) ;

 

U (0, t) (t) ;

 

U (a, t) (t).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для вариантов 1–15:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) ax

 

 

b sin x ;

(t) (t) 0;

T 0,02 .

 

 

 

 

 

 

 

 

1. а 1,1;

 

b 1,1;

2. а 1,3;

 

b 1,4;

3.

а 1,1;

b 1,2;

 

 

 

 

4. а 1,3;

 

b 1,5;

5. а 1,1;

 

b 1,3;

6.

а 1,5;

b 1,1;

 

 

 

 

7. а 1,1;

 

b 1,4;

8. а 1,5;

 

b 1,2;

9.

а 1,1;

b 1,5;

 

 

 

 

10. а 1,5;

 

b 1,3;

11. аД1,3; b 1,1; 12. а 1,5; b 1,4;

 

 

 

 

13. а 1,3;

 

b 1,2;

14. а 1,5;

 

b 1,5;

 

15. а 1,3; b 1,3.

 

 

 

 

Для вариантов 16–30:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) e bx sin ax;

(1) 0; (t) e b sin at ;

 

T 0,02 .

 

 

 

 

 

 

 

 

16. а 12 ;

b 0,1;

 

17. а 4 ;

Иb 0,4;

 

 

 

 

18. а 12 ;

b 0,2;

 

19. а 4 ;

b 0,5;

 

 

 

 

 

 

 

 

 

 

20. а 12 ;

b 0,3;

 

21. а 3 ;

b 0,1;

 

 

 

 

 

 

 

 

 

 

22. а 12 ;

b 0,4;

 

23. а 3 ;

b 0,2;

 

 

 

 

 

 

198

24.

а 12 ;

b 0,5;

25.

а 3 ;

b 0,3;

26.

а 4 ;

b 0,1;

27.

а 3 ;

b 0,4;

28.

а 4 ;

b 0,2;

29.

а 3 ;

b 0,5;

30.

а 4 ;

b 0,3.

 

 

 

Вопросы

задания для самопроверки [1, 2, 3, 4, 5, 6, 7]

производных

 

 

 

С1. Как выглядят сеточные аналоги первой и второй частных

 

?

 

2.

Объясн те,

как получаются сеточные аналоги первой и вто-

 

б

рой частных

.

3.

В чем заключается метод сеток численного решения диффе-

ренциальных уравнен й с частными производными?

4.

Как решается задача Дирихле в прямоугольнике для уравне-

 

 

А

ния Лапласа методом сеток?

5.

Как решается уравнение математической физики гиперболи-

ческого типа методом сеток?

 

 

6. Как решается уравнение математической физики параболиче-

ского типа методом сеток?

Д

 

 

 

7. Как выглядит метод прогонки для уравнения теплопроводно-

сти?

 

 

 

 

И

199

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]