
Вариационные методы Старинова / Приближенные методы-Буханько АА
.pdf
Связь |
между |
решением |
y(x) интегрального уравнения (31) и решением |
||||||
u(x) |
x2 2 ln x |
краевой задачи (32) определена выражениями (II.3), (II.5), |
|||||||
см. приложение II. Следовательно, точным |
|||||||||
решением уравнения (31) будет |
|
|
|
|
|||||
|
y |
x |
|
d 2u x |
2 |
2 |
|
|
|
|
|
dx2 |
x2 |
||||||
|
|
|
|
|
|||||
и, обратно, приближенное решение задачи |
|||||||||
(32) будет иметь вид: |
|
|
|
|
|
||||
|
xi |
|
|
|
|
|
|
|
|
|
ui |
xi |
|
s yi ds 1 i 0,5 . |
1
На рис. 3.1 представлено графическое сравнение точных и приближенных решений
уравнений (31), (32); в табл. 3.2 – их значения.
Рис. 3.1
Таблица 3.2 – |
Сравнение приближенных |
||||
|
|
|
и точных решений |
||
xi |
y(x) |
|
yi |
u(x) |
ui |
1,0 |
4,0000 |
|
4,0000 |
1,0000 |
1,0000 |
1,2 |
3,3889 |
|
3,3825 |
1,0754 |
1,0676 |
1,4 |
3,0204 |
|
3,0118 |
1,2871 |
1,2409 |
1,6 |
2,7813 |
|
2,7719 |
1,6200 |
1,4989 |
1,8 |
2,6173 |
|
2,6077 |
2,0644 |
1,8345 |
2,0 |
2,5000 |
|
2,4906 |
2,6137 |
2,2453 |
3.3 Метод вырожденного ядра
Ядро K x, s называется вырожденным, если оно может быть представлено в виде конечной суммы парных произведений:
|
n |
|
K x, s |
i x i s , |
(3.16) |
i |
1 |
|
где функции i x , i s можно считать линейно независимыми. Для таких ядер интегральное уравнение Фредгольма второго рода (3.2):
|
|
b |
y x |
f x |
K x, s y(s)ds |
|
|
a |
51

решается весьма просто. Подставляя выражение (3.16) в уравнение (3.2), будем иметь
|
|
|
n |
|
|
|
|
y x |
f |
x |
ci |
i (x) , |
(3.17) |
||
|
|
|
i 1 |
|
|
|
|
где |
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
ci |
i |
s y |
s ds i |
1, n |
(3.18) |
||
|
a |
|
|
|
|
|
|
некоторые постоянные коэффициенты. Если в выражения (3.18) подставить
функцию (3.17), то для определения коэффициентов ci |
получим алгебраиче- |
||||||||||||||
скую систему линейных уравнений |
|
|
|
|
|
|
|
|
|
|
|
||||
|
b |
|
|
b |
|
|
n |
|
|
|
|
|
|
|
|
ci |
i |
s f s ds |
|
i |
s |
|
ci |
1 s ds |
i 1, n |
|
|||||
|
a |
|
|
a |
|
|
i |
1 |
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ci |
|
cj |
ji |
|
fi , |
|
|
|
|
|
|
(3.19) |
|
|
|
|
j |
1 |
|
|
|
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
fi |
i |
s f s |
ds, |
ij |
|
i |
s j s |
ds . |
(3.20) |
|||||
|
|
a |
|
|
|
|
|
a |
|
|
|
|
|
|
|
Запишем систему (3.19) в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ij |
ji |
cj |
|
fi |
i |
1, n . |
|
|
|
|
||
|
|
j |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
Обозначим определитель этой системы
1
det |
ij |
ji |
|
11 |
|
21 |
|
n1 |
12 |
1 |
22 |
|
n2 |
|
|
|||
1n |
|
2n |
1 |
nn |
|
|
и через |
|
ij |
- алгебраические дополнения соответствующих элементов |
|||||
|
|
|
|
|
|
|
|
|
ij |
ij |
этого определителя. Если |
|
0 , то по правилу Крамера нахо- |
||||
|
|
|
|
|
|
|
||
дим |
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
ji ( |
) f j |
|
|
|
|
|
|
ci |
j 1 |
|
i 1, n . |
||
|
|
|
|
|
||||
|
|
|
|
|
Следовательно, в силу (3.17) интегральное уравнение (3.2) имеет единственное решение:
52

y x |
f x |
n |
n |
|
|
|
|
|
ji |
f j i |
x . |
i 1 |
j 1 |
Отсюда, подставляя вместо |
fi |
соответствующее выражение (3.20) и заменяя |
||||||||
сумму интегралов интегралом суммы, получим |
|
|
|
|||||||
|
|
|
|
|
|
b |
x, s, |
|
|
|
|
|
y |
x |
f |
x |
|
f s ds |
, |
(3.21) |
|
|
|
a |
|
|||||||
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
n |
n |
|
|
|
|
|
|
|
|
где |
x, s, |
i |
x |
j |
s ji |
|
. Из формулы |
(3.21) |
вытекает, что |
|
|
i 1 |
j 1 |
|
|
|
|
|
|
|
|
функция
x, s,
R(x, s, )
n |
n |
|
ji |
|
i x |
j s |
|
|
|
||
|
|
||
i 1 |
j 1 |
|
|
есть резольвента интегрального уравнения (3.2).
Собственные значения ядра K x, s определяются из уравнения |
|
||
0 . |
(3.22) |
||
|
|
|
|
Если k k 1, m; m n есть корень уравнения (3.22) (очевидно, |
k 0 ), |
то соответствующие собственные функции k x ядра K x, s , т.е. нетривиальные решения однородного уравнения
y x
b
k K x, s y(s)ds ,
a
будут иметь вид
k x
|
n |
|
|
k |
c(k ) |
i |
x , |
i |
|
||
i |
1 |
|
|
где c(k ) |
- ненулевые решения линейной однородной системы |
||||||
i |
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
ij |
ji cjk 0 i 1, n . |
|
||
|
|
|
j 1 |
|
|
|
|
Если |
k |
есть собственное значение ядра K x, s |
, то неоднородное |
||||
|
|
|
|
|
|
|
уравнение (3.2) или не имеет решений, или имеет бесконечное множество решений.
Пример 3.3. Найти приближенное решение уравнения
|
1 |
|
y x |
sh(xs) y s ds 1 x2 . |
(33) |
|
0 |
|
53

Заменяем ядро K(x, s) sh(xs) суммой первых трех членов ряда Тейлора:
|
|
|
|
|
(xs)3 |
|
(xs)5 |
|
|
||
|
sh(xs) |
xs |
|
|
|
|
|
. |
|
|
|
|
3! |
|
|
5! |
|
|
|||||
|
|
|
|
|
|
|
|
|
|||
Тогда решение уравнения (33) будем искать в виде |
|
|
|||||||||
|
y(x) 1 x2 |
c x c x3 |
c x5 . |
|
|
||||||
|
|
|
|
1 |
2 |
|
3 |
|
|
|
|
Обозначим f (x) 1 x2 , |
1 |
x, |
2 |
|
x3 , |
3 |
x5 , |
1 |
(s) |
||
|
|
|
|
|
|
|
s5
3 и найдем по формулам (3.20) коэффициенты системы
5!
1 |
|
|
1 |
1 |
|
f |
|
(s) f (s)ds |
(s s3 )ds |
, |
|
1 |
|
||||
1 |
|
|
4 |
|
|
0 |
|
|
0 |
|
|
|
|
|
|
s3 s, 2 (s) 3! ,
(3.19):
|
|
|
|
|
|
|
f |
|
|
|
1 |
|
|
|
(s) f (s)ds |
|
1 |
1 |
|
|
(s3 |
s5 )ds |
|
|
|
1 |
|
, |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 3! |
|
|
|
|
|
|
|
|
|
|
|
72 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
f |
|
1 |
|
|
(s) f (s)ds |
1 |
1 |
(s5 |
|
s7 )ds |
1 |
|
|
, |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
3 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 5! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2880 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
1 |
|
(s) |
|
(s)ds |
1 |
|
s2 ds |
|
1 |
, |
|
|
|
|
|
|
|
|
1 |
1 |
|
s |
4 ds |
|
1 |
|
, |
|
|
|
|
|
|
|
|
1 |
1 |
s6 ds |
|
1 |
, |
||||||||||||||||||||||||
11 |
1 |
1 |
|
|
|
|
|
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
13 |
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
0 3! |
|
|
|
30 |
|
|
|
|
|
|
|
5! |
|
|
840 |
|
||||||||||||||||||||||||||||
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|||||||||||||||||||||||||||
|
|
|
1 |
s4 ds |
1 |
|
, |
|
|
|
|
|
|
|
|
1 |
|
1 |
s6 ds |
|
1 |
|
, |
|
|
|
|
|
|
1 |
|
1 |
|
s8 ds |
|
1 |
, |
|
|
|
|||||||||||||||||||||||||
|
21 |
|
|
|
|
|
|
22 |
|
|
|
|
|
|
|
|
|
|
23 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
3! |
|
|
|
|
|
|
|
|
|
|
42 |
|
|
|
0 5! |
|
|
|
1080 |
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
1 |
|
6 |
|
|
1 |
|
|
|
|
|
|
|
1 |
1 |
8 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
1 |
1 |
|
|
|
10 |
|
1 |
|
|
|
|
|
||||||||||||||||||
|
31 |
|
s |
|
ds |
|
|
|
|
, |
|
|
|
32 |
|
|
|
|
|
|
|
|
|
s |
ds |
|
|
|
|
|
, |
|
33 |
|
|
|
|
|
|
|
|
|
s |
|
ds |
|
|
|
|
|
. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
7 |
|
|
|
|
|
|
|
|
|
|
3! |
|
|
54 |
|
|
|
|
|
|
|
5! |
|
|
|
|
1320 |
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
Таким образом, получаем систему |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
c |
|
1 |
c |
|
|
1 |
|
c |
|
1 |
c |
1 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
1 |
|
3 |
1 |
5 |
2 |
|
7 |
3 |
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
c |
|
|
1 |
|
c |
|
1 |
c |
|
1 |
|
|
c |
|
1 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
30 |
|
1 |
|
|
|
|
42 |
2 |
|
54 |
|
3 |
72 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
c |
|
|
1 |
|
c |
|
|
|
|
1 |
|
c |
|
|
|
1 |
|
c |
|
|
|
1 |
. |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
3 |
|
840 |
1 |
|
|
1080 |
2 |
|
|
|
|
1320 |
3 |
|
2880 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
Решив ее методом итерации, получаем c1 |
0,3833, |
c2 |
|
|
|
|
|
0, 0273, c3 |
|
0, 0008. |
Таким образом, приближенное решение уравнения (33) можно записать в виде
y(x) 1 x2 0,3833x 0,0273x3 0,0008x5.
54

4 ЗАДАНИЕ НА КУРСОВУЮ РАБОТУ
I Краевые задачи для обыкновенных дифференциальных уравнений
I.1. Найти решение краевой задачи
y p(x) y q(x) y f (x), x a, b
0 y(a) |
1 y (a) |
A1 , |
0 y(b) |
1 y (b) |
B1 |
с шагом h 0,1 :
а) методом конечных разностей, используя конечно-разностные и цен- трально-разностные отношения; сравнить результаты вычислений с точным решением; сделать вывод;
б) методом прогонки, используя разностную схему, давшую в п. а) лучший результат;
в) сравнить результаты вычислений методом конечных разностей, методом прогонки с точным решением; сделать вывод.
I.2. Найти решение краевой задачи:
y p(x) y q(x) y f (x), x a, b
y(a) A2 , y(b) B2
с шагом h 0,1 методом Галеркина, методом коллокаций и методом Ритца. Привести сравнение результатов вычислений с точным решением. Сделать вывод.
Примечание: Для составления краевых задач использовать функции y(x) , p(x) , q(x) и коэффициенты 0 , 1 , 0 , 1 , приведенные в табл. 4.1.
II Уравнения в частных производных
II.1. Методом сеток найти приближенное решение первой смешанной задачи для уравнения теплопроводности на отрезке:
u |
a2u |
xx |
, x |
0, l , t 0; |
||||||||
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2x2 |
|
|
|
l |
|
|||
|
|
|
|
|
|
, |
0 |
x |
|
|
, |
|
|
|
|
|
|
|
|
|
|||||
u |
x, 0 |
|
|
|
l |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
||
|
|
|
l |
x, |
|
x |
l; |
|
||||
|
|
|
2 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
u 0,t |
u l,t |
0. |
|
|
|
Сравнить результаты решений с аналитическим решением.
55

№ |
|
1 |
2 |
3 |
|
4 |
|
|
5 |
|
|
6 |
|
7 |
8 |
9 |
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
4 |
5 |
3 |
|
5 |
|
|
5 |
|
|
3 |
|
2 |
3 |
2 |
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
3 |
5 |
2 |
|
8 |
|
|
6 |
|
|
10 |
|
5 |
3 |
4 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
№ |
|
11 |
12 |
13 |
|
14 |
|
|
15 |
|
|
16 |
|
17 |
18 |
19 |
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
2 |
3 |
4 |
|
3 |
|
|
4 |
|
|
5 |
|
4 |
4 |
6 |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l |
|
7 |
4 |
8 |
|
5 |
|
|
2 |
|
|
4 |
|
6 |
4 |
3 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
II.2. Методом сеток найти приближенное решение первой смешанной за- |
||||||||||||||||
дачи для волнового уравнения на отрезке [0, l] : |
|
|
|
|
|||||||||||||
|
|
|
|
u |
a2u |
xx |
, x |
0,l , t 0; |
|
|
|
|
|||||
|
|
|
|
tt |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
u x, 0 |
x x l , ut |
x, 0 |
0; |
|
|
|
|||||||
|
|
|
|
|
|
u 0,t |
|
u |
l,t |
0. |
|
|
|
|
|
Сравнить результаты вычислений с аналитическим решением.
№ |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
a |
3 |
3/2 |
2 |
2/3 |
4 |
2 |
1/2 |
1/3 |
1/3 |
2/3 |
|
|
|
|
|
|
|
|
|
|
|
l |
3 |
3 |
1 |
2/3 |
2 |
2 |
1 |
1/2 |
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
№ |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
|
|
|
|
|
|
|
|
|
|
|
a |
4 |
1/3 |
3 |
3/2 |
3 |
2 |
1/2 |
1/2 |
3/2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
l |
3 |
2 |
3/2 |
3/2 |
2 |
1/2 |
3/2 |
1/2 |
1 |
3 |
|
|
|
|
|
|
|
|
|
|
|
II.3. Найти приближенное решение задачи Дирихле в квадрате с верши-
нами A(0, 0), B(0,1), C(1,1), |
D(1, 0) : |
|
|
|
|
|
|
|
|||
|
|
|
|
u |
0, |
|
|
||||
|
|
|
|
|
|
|
|
|
|
||
u |
|
AB |
1 |
y , |
u |
|
|
BC |
1 |
x , |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||
u |
|
CD |
2 |
y , |
u |
|
AD |
2 |
x |
||
|
|
|
|
||||||||
с шагом h 1 3 , применяя: |
|
|
|
|
|
|
|
|
|
а) метод сеток; б) метод усреднения Либмана с точностью 0,1.
56

№
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
1 y |
|
1 x |
2 y |
|
|
2 x |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30 y |
|
30 1 |
x2 |
0 |
|
|
|
|
|
0 |
|
|
|
|||
30 y |
|
30 cos |
|
|
x |
|
30 cos |
|
|
y |
|
30x |
|
|
||
|
2 |
|
|
|
2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
50 y 1 |
y2 |
0 |
|
|
|
|
0 |
|
|
|
|
|
50sin |
x |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
20 y |
|
20 |
|
|
|
|
20 y2 |
|
|
50x 1 |
x |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
0 |
|
50x 1 x |
50 y 1 y2 |
50x 1 x |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30sin |
y |
20x |
20 y |
|
|
30x 1 |
x |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30 1 |
y |
|
|
|
|
|
20 y |
|
|
30 1 |
|
x |
||||
20 |
x |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
50sin |
y |
|
|
|
|
|
30 y2 |
|
|
50sin |
x |
|||||
30 |
x |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
40 y2 |
40 |
|
|
|
|
40 |
|
|
|
|
|
40 sin |
|
x |
||
|
|
|
|
|
|
|
|
|
2 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
50 y |
|
50 1 |
|
x |
0 |
|
|
|
|
|
0 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
40 1 |
y |
30x |
30 y |
|
|
40 1 |
|
x2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
10 y2 1 |
y2 |
0 |
|
|
|
|
0 |
|
|
|
|
|
10sin |
x |
||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
50sin |
y |
20x |
20 y2 |
|
|
50sin |
x |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
30 y2 |
30 1 |
|
x |
0 |
|
|
|
|
|
0 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
20 y |
|
20 cos |
|
|
x |
|
20 cos |
|
y |
|
20x |
|||||
|
2 |
|
|
2 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
40sin |
y |
|
|
|
|
|
50 y2 |
|
|
40sin |
x |
|||||
50 |
x |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
10 y |
|
10 |
|
|
|
|
10 y2 |
|
|
30x 1 |
x |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
15 y |
|
15 cos |
|
|
x |
|
15 cos |
|
|
|
y |
|
15x |
|
|
|
|
2 |
|
|
2 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
10 y 1 |
y2 |
0 |
|
|
|
|
0 |
|
|
|
|
|
20sin |
x |
||
|
|
|
|
|
||||||||||||
0 |
|
30x 1 x |
30 y 1 y2 |
30x 1 x |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
57

III Интегральные уравнения
III.1. Для краевой задачи
y p(x) y q(x) y f (x), x a, b
y(a) A3 , y (a) B3
получить соответствующее уравнение Вольтерра второго рода:
|
x |
y x |
K x, s y s ds f x |
a
и найти его приближенное решение. Результат сравнить с точным решением.
III.2. Составить интегральное уравнение Фредгольма второго рода:
|
b |
y x |
K x, s y s ds f x |
a
и найти его приближенное решение. Результат сравнить с точным решением.
Примечания: 1) Для составления уравнений использовать функции y(x) , K(x, s) , приведенные в табл. 4.1.
2) При составлении уравнения Фредгольма необходимо уточнить сходимость будущего решения, т.е. определить диапазон значений величины
1 .
M b a
Требования к оформлению курсовой работы
Оформление курсовой работы (структура, оформление текста, нумерация страниц и разделов, графический материал, таблицы, формулы и уравнения, список использованной литературы, приложения) производится согласно [7].
При выполнении работы предлагается использовать программные про-
дукты MathCAD, Maple и др.
58

№
y(x)
[a,b]
p(x)
q(x)
0,1, 0,1
K(x, s)
|
|
|
1 |
|
|
|
|
|
2 |
|
3 |
|
|
|
|
|
4 |
|
|
|
|||||||
|
|
x3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
x |
|
|
ln Nx |
0,1N |
|
|
0,1N |
|||||||||||||
|
|
|
|
|
0,1N |
e N 0,1N |
|
|
|
|
|
|
|
|
|||||||||||||
|
N |
1 |
|
|
N |
|
x |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
[N |
1, N |
2] |
[N |
|
2, N |
3] |
[1, 2] |
|
[N 1, N 2] |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
(N 1)x 1 |
(N 1)x 1 |
(N 2)x 1 |
(N 2)x 1 |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
0,1N |
|
|
|
|
|
0,1Nx |
|
|
0,1N 1 |
|
|
0,1N |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
1 |
|
|
0 |
|
|
3 |
-2 |
3 |
-2 |
1 |
0 |
3 |
2 |
1 |
0 |
1 |
0 |
3 |
-2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
1 |
|
|
|
(1 |
|
s)(e0,3xs |
1) |
|
1 x |
s |
|
|
|
x |
|
s |
|
|||||||
1 |
x2 |
s2 |
2 |
|
xs |
1 |
|
xs |
|||||||||||||||||||
|
|
|
|
|
|
|
|
№ |
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
y(x) |
|
|
0,1Nx4 |
|
|
|
x2 |
0,1N |
(N |
1)e 0,1Nx |
0,1N ln Nx |
|||||||||||||||||||||||||||||||||||
|
|
|
N |
1 |
|
|
|
|
|
|
|
|
|
N |
1 |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
[a,b] |
[N |
|
2, N |
1] |
[N |
|
|
|
3, N |
2] |
[N |
|
|
4, N |
3] |
|
|
|
[1, 2] |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
p(x) |
(2N 1)x 1 |
(1 2N )x 1 |
(1 N)x 1 |
|
(2 N )x 1 |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
q(x) |
|
|
|
0,3N |
|
|
|
|
|
0,3Nx |
|
|
|
0,3Nx2 |
|
|
|
|
|
0,3N |
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,1, 0,1 |
1 |
|
3 |
|
0 |
|
|
1 |
2 |
|
|
-1 |
|
1 |
0 |
1 |
0 |
|
|
1 |
|
2 |
3 |
|
|
-1 |
|
1 |
|
0 |
||||||||||||||||
K(x, s) |
|
|
1 |
|
s |
|
|
|
|
|
|
|
x |
|
s |
|
|
|
|
x s |
|
|
|
|
|
1 |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
xs |
|
|
|
|
|
|
|
x |
|
s |
|
|
|
|
|
|
|
s |
|
|
1 |
|
|
x2 |
|
s2 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
№ |
|
|
|
|
9 |
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
11 |
|
|
|
|
|
|
|
12 |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y(x) |
|
|
|
N |
2 |
|
|
|
|
cos(Nx) 0,1N |
sin(Nx) |
0,1N |
(N |
|
1)x2 |
0,1N |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
x 0, 2N |
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
[a,b] |
[N |
|
5, N |
4] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[N |
|
|
6, N |
5] |
|||||||||||||
|
|
|
|
[ 2 , ] |
|
|
|
[ 2 , ] |
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
p(x) |
|
(N 3)x 1 |
(N 3)x 1 |
(3 N )x 1 |
(3N 1)x 1 |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
q(x) |
|
0,3N 1 |
|
|
|
|
|
N |
|
1 |
x |
|
|
|
N |
|
|
1 |
x2 |
|
|
|
|
|
N |
1 |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
10 |
|
|
|
10 |
|
|
|
|
|
|
10x |
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
0,1, 0,1 |
2 |
|
-1 |
|
1 |
|
|
0 |
1 |
|
|
0 |
|
3 |
2 |
0 |
1 |
|
2 |
|
-1 |
2 |
|
|
-1 |
|
0 |
|
1 |
|||||||||||||||||
K(x, s) |
(1 |
|
s)(e0,2 xs |
1) |
|
|
|
|
2 |
|
s |
|
|
|
|
2 |
|
|
|
x |
s |
|
|
|
|
|
|
x |
s |
|
|
|
|
|||||||||||||
|
|
|
|
1 |
|
xs |
|
|
|
1 |
|
|
xs |
|
|
|
|
|
|
|
x |
s |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
59

№ |
|
|
13 |
|
|
|
|
|
|
|
|
|
14 |
|
|
|
|
|
|
15 |
|
|
|
|
|
|
16 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
x3 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
x |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
0, 2N |
|
|
|
|
|
|
||||||||
y(x) |
|
0, 2N |
|
|
|
|
|
e N |
|
|
|
ln Nx |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
N |
1 |
|
|
|
0,1N |
|
|
0,1N |
|
|
|
N |
|
|
x |
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
[a,b] |
|
|
|
|
[N, |
|
|
|
|
|
[N |
1, |
|
|
|
|
[1, 2] |
|
|
|
|
|
[N |
|
|
2, |
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
N |
1] |
|
|
|
|
|
|
|
N ] |
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
1] |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
p(x) |
|
(1 3N )x 1 |
|
(0,1N |
|
1)x 1 |
(1 |
|
0,1N)x 1 |
(0, 2N |
1)x 1 |
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
q(x) |
|
|
|
|
N |
2 |
x |
|
|
|
|
|
0,1Nx2 |
|
|
|
|
0,1N |
|
|
|
|
0,3N |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
x |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
0,1, 0,1 |
0 |
|
1 |
|
-2 |
|
|
3 |
1 |
|
|
0 |
|
|
1 |
|
|
2 |
0 |
|
1 |
|
-2 |
|
|
1 |
-2 |
|
3 |
|
|
0 |
|
|
1 |
|
||||||||||||||
K(x, s) |
(1 |
|
s)(e 0,1xs |
1) |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
s |
|
|
1 |
x |
|
s |
|||||||||||||||||||||||
|
|
2 |
|
x2 |
|
|
|
s2 |
|
|
|
|
2 |
|
xs |
|
|
|
|
2 |
|
|
xs |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
№ |
|
|
17 |
|
|
|
|
|
|
|
|
|
18 |
|
|
|
|
|
|
19 |
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
y(x) |
|
|
(N 1)x4 |
|
|
1 |
|
cos(0, 2Nx) |
sin(0,1Nx) |
N |
(0,1N |
1)x2 |
|
|
|
3N |
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
N |
3 |
|
|
|
|
N |
10 |
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
[a,b] |
|
|
|
|
[N |
1, |
|
|
|
|
|
|
|
[ |
|
, ] |
|
|
|
|
[ |
|
, ] |
|
|
|
|
|
[N |
|
|
2, |
|
|
|
|
|
|||||||||||||
|
|
|
|
N 2] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N 1] |
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
p(x) |
|
(0,1N |
1)x 1 |
(1 |
|
0, 2N)x 1 |
(0, 2N |
|
1)x 1 |
(N |
|
5)x 1 |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
q(x) |
|
|
1 |
|
N |
x |
|
|
|
|
|
2 N |
x |
|
|
|
2 |
|
|
N |
x2 |
|
|
0,3Nx |
1 |
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
0,1, 0,1 |
0 |
|
1 |
|
2 |
|
|
3 |
-3 |
|
|
2 |
|
|
0 |
|
1 |
0 |
|
1 |
|
2 |
|
|
-3 |
1 |
|
0 |
|
|
-3 |
|
2 |
|
||||||||||||||||
K(x, s) |
|
|
|
1 |
x |
|
|
(1 |
|
|
s)(e0,3xs |
1) |
|
|
|
1 |
|
x |
|
|
1 |
x |
|
s |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
1 |
|
xs |
|
|
|
|
|
2 |
|
xs |
|
|
|
1 |
|
|
xs |
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
60