Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Ankilov

.pdf
Скачиваний:
36
Добавлен:
07.02.2015
Размер:
2.14 Mб
Скачать

10. Задача y // 2 y / y 0 , y(0) 1, y / (0) 0 называется:

а) краевой задачей; б) задачей Вронского; в) задачей Коши;

г) фундаментальной задачей.

11.Общее решение уравнения y // 2 y / y 0 имеет вид:

а) y C1ex C2 ex ; б) y C1 x C2 ex ; в) y C1ex C2 xex ; г) y C1ex C2 e x ; д) y C1ex .

12.Какая из следующих функций является решением уравнения y / xy x :

а) y x ;

б) y x2 5 ; в) y x 5 ;

г) y x2 x3 ;

д)

y x2 .

 

 

13. Какие две

следующие функции являются решением системы уравнений

x/ 3x 4 y , y / 2x 5y :

а)

x et ,

y e t ;

б)

x e t ,

y et ;

в)

x 2e t , y e t ;

г)

x 2e t , y et ;

д)

x 2e t ,

y e t .

14. Какое из следующих уравнений является характеристическим для линейной

однородной системы x/ y ,

y / x :

а)

r 2

1 0 ;

 

б)

r 2

1 0 ;

 

в)

r 3

1;

 

г)

r 2

r 0 ;

 

д)

r 2

r 0 .

 

15. Какие две следующие функции описывают общее решение системы x/ y ,

y / x :

а)

x C1 cos t C2 sin t , y C1 sin t C2 cos t ;

 

б)

x C1 cos t , y C1 sin t ;

 

в)

x cos t sin t , y sin t cost ;

 

г)

x sin t ,

y sin t ;

 

д)

x cos t ,

y 0 .

 

181

6.8.2. Задачи Образцы решения задач

Задача 1. Найти общий интеграл уравнения y'sin x y cos x 0 .

Решение. Дано уравнение с разделяющимися переменными. Заменяем у' по формуле y' dydx и разделяем переменные:

dydx sin x y cos x , dy sin x y cos x dx , dyy cossin xx dx

(таким образом, левая часть уравнения зависит только от у, а правая зависит только от х). Интегрируем:

dy

 

cos x

dx ;

dy

ln

 

y

 

,

cos x

dx

d sin x

ln

 

sin x

 

ln

 

C

 

ln

 

C sin x

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

sin x

y

 

 

 

 

 

 

 

 

 

sin x

 

sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем общий интеграл ln

 

y

 

ln

 

C sin x

 

.

Выражая у, находим общее решение y C sin x ,

 

 

 

 

где С – произвольная константа. В процессе решения мы делили на у, поэтому могли потерять решение у = 0 (это легко проверить, подставляя у = 0 в исходное уравнение). Однако решение у = 0 входит в общее решение при С = 0.

Задача 2. Найти решение задачи Коши y' 2xy 2xe x2

, у(0) = 2.

 

 

 

 

 

Решение. Дано линейное уравнение 1-го порядка. Решение уравнения ищем в виде

y U V ,

где

U U (x) ,

V V (x)

некоторые

функции.

Находим

производную

y' (UV )' U 'V UV '

и

подставляем

у,

 

у'

в

исходное

уравнение:

(U 'V UV ') 2xUV 2xe x2 ,

или

U 'V U (V ' 2xV ) 2xe x2 .

Приравнивая

к

нулю

выражение в скобках, получаем два уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V ' 2xV 0 ,

 

 

 

 

 

 

 

(6.30)

 

 

 

 

 

U 'V 2xe x2 .

 

 

 

 

 

 

(6.31)

Уравнение (6.30) – это уравнение с разделяющимися переменными. Подставляя в (6.30)

V ' dV и разделяя переменные, получаем

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

2xV 0 ,

dV 2xV , dV 2xdx .

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

dx

 

 

V

 

 

 

 

 

 

 

Интегрируем: dV 2xdx ,

ln

 

V

 

x2 C .

Положим С = 0 и

найдем

V e x2 .

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя полученное значение V в уравнение (6.31), находим U: U ' e x2

 

2xe x2

, U ' 2x ,

U 2xdx x2

C . В

итоге

общее

решение

задачи

есть

y UV (x2 C)e x2 ,

где

С

произвольное число. Потребуем, чтобы

функция

y

удовлетворяла

начальному

условию

y (0) = 2.

Подставим в уравнение y (x2

C) e x2

значения х = 0,

у = 2:

2 = (0+С)е0, или

2 = С. Значит, решение задачи Коши получается из общего решения при С = 2:

 

 

y (x2 2)e x2 .

182

Задача 3. Решить уравнение (sin xy xy cos xy)dx x2 cos xydy 0 .

Решение. Проверим, что данное уравнение является уравнением в полных

дифференциалах.

Обозначим M (x, y) sin xy xy cos xy ,

N (x, y) x2 cos xy

и найдем

частные производные:

 

 

 

 

 

 

 

M

x cos xy x cos xy x2 y sin xy 2x cos xy x2 y sin xy ,

 

 

 

 

 

y

 

 

 

 

 

 

 

 

N 2x cos xy x2 y sin xy ,

 

 

 

 

 

 

 

x

 

 

так что

M

 

N

при любых значениях х и у. Значит, условие (6.4) выполняется. Таким

y

x

 

 

 

 

 

 

образом, данное уравнение есть уравнение в полных дифференциалах и согласно (6.5),

 

 

 

 

 

F(x, y) M (x, y) sin xy xy cos xy ,

(6.32)

 

 

 

 

 

x

 

 

 

 

 

 

 

F(x, y) N (x, y) x2 cos xy .

 

(6.33)

 

 

 

 

 

y

 

 

Из соотношения (6.32) находим интегрированием по переменной х (у принимается за константу):

F(x, y) (sin xy xy cos xy)dx x sin xy ( y) ,

где φ(у) – пока неопределенная функция. Подставляя F(x,y) в уравнение (6.33), получаем:

 

 

x sin xy ( y) x2 cos xy ,

 

y

 

 

 

 

откуда x2 cos xy '( y) x2 cos xy ,

'( y) 0 , так что

( y) C сonst . Таким образом,

F(x, y) x sin xy C . Общий интеграл исходного уравнения есть: x sin xy C 0 .

Задача 4.1. Найти общее решение однородного уравнения y 4 y 5y 0 .

Решение. Характеристическое уравнение K 3 4K 2 5K 0 имеет три простых корня: K1 0, K2 2 i, K3 2 i . Значит,

y y0.0. C1 C2e2 x cos x C3e2 x sin x .

Задача 4.2. Найти общее решение линейного однородного уравнения y(4) 4 y 4 y 0 .

Решение. Характеристическое уравнение K 4 4K 2

4 0 имеет корни

K1,2 i 2

кратности 2. Следовательно, общее решение уравнения

 

 

y y0.0. C1 cos 2x C2 x cos 2x C3 sin 2x C4 x sin 2x .

Задача 4.3. Найти общее решение дифференциального уравнения: y 5y 6 y x2 x 2 .

183

 

 

 

 

Решение.

Характеристическое уравнение K 2 5K 6 0

имеет корни кратности 1:

K

1

2, K

2

3.

 

Значит,

y

0.0.

C e2 x

C

2

e3x . Так

как i 0

 

не

совпадает

 

с

корнями,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частное решение следует искать в виде

 

y ч.н.

Ax2

Bx C . Подставляя

y ч.н.

в исходное

уравнение, получаем тождество:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2A 5 2Ax B 6 Ax2 Bx C x2 x 2 .

 

 

 

 

 

 

 

 

 

Приравниваем

 

коэффициенты

при

одинаковых

степенях

 

x в левой

и

правой

частях

тождества: x : 10A 6B 1;

 

x2

: 6A 1; свободные члены:

2A 5B 6C 2 .

 

 

 

 

 

 

 

 

 

Решая систему, находим:

A

 

1

,

B 1

C

10

. Значит, y

ч.н.

1 x2

 

1 x 10

; общее

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

9

 

27

 

 

 

 

 

 

 

 

 

6

 

 

9

 

27

 

 

решение уравнения имеет вид y y

0.0.

y ч.н. C e2 x C

e3x

1 x2

1 x 10 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

6

 

 

9

27

 

 

 

 

 

 

 

 

 

 

Задача 4.4. Найти общее решение уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 5y 6 y 18x2

6x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

 

Характеристическое

 

уравнение

имеет

 

вид:

 

 

K 3 5K 2 6K 0 ,

или

K K 2 5K 6 0 ;

его

 

 

корни

 

 

 

K1

0, K2

2,

K3 3

 

 

(все

 

корни

 

простые).

Тогда

y

0.0.

C e0 x C

2

e2 x C

e3x C C

e2 x

C

e3x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

3

 

 

 

 

1

 

2

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как среди корней характеристического уравнения есть корень K 0 , совпадающий

с числом i 0 i0 0 , то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ч.н. x Ax2 Bx C Ax3 Bx2 Cx .

 

 

 

 

 

 

 

 

 

 

 

Подставляя y ч.н. в исходное уравнение, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6A 5 6Ax 2B 6 3Ax2 2Bx C 18x2 6x .

 

 

 

 

 

 

 

 

 

Приравниваем коэффициенты при одинаковых степенях x:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18A 18, 30A 12B 6, 6A 10B 6C 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая систему, находим:

A 1,

 

B 3,

C 4 ,

тогда

 

y ч.н. x3 3x2

4x ,

и

общее

решение имеет вид y C

C

e2 x

C

e3x

x3 3x2

4x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача 4.5. Найти общее решение уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(4)

5y(3) 6 y(2) x2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

Имеем:

 

K 4 5K 3

 

6K 2

0 ,

или

K 2 K 2

5K 6 0 .

 

Корни:

 

K 1

0

(кратности 2), K

2

2

и

K

3

3 (простые). Получаем

y

0.0.

C C

2

x C

e2 x

C

e3x

. Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

3

 

 

4

 

 

 

 

 

i 0

 

является

 

корнем

кратности

2

характеристического

 

уравнения,

то

y ч.н. x2 Ax2 Bx C . Далее решение строится как в предыдущей задаче.

 

 

 

 

 

 

 

 

 

 

 

Находим:

 

 

y ч.н

 

= Ax4 Bx3 Cx2 ,

 

y ч.н= 4Ax3 3Bx2

2Cx,

y ч.н=12Ax2 6Bx 2C,

y ч.н= 24Ax 6B,

 

y(4) ч.н= 24A. Подставляя

 

функцию

yч.н.

 

и

ее

 

производные

 

в

исходное

уравнение, получаем тождество:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24A 5 24Ax 6B 6 12Ax2 6Bx 2C x2 .

 

 

 

 

 

 

 

 

 

184

Приравнивая коэффициенты при одинаковых степенях x , получаем систему линейных

уравнений относительно неизвестных А, В, С:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

72A 1,

 

 

 

 

 

 

72A 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

120A 36B 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10A 3B 0 .

 

 

 

 

 

 

 

1

24A 30B 12C 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4A 5B 2C 0

 

 

 

 

 

 

 

Решая систему, находим: A

1

,

B

 

5

 

,

 

C

19

, y

ч.н.=

1

 

 

x4

 

 

5

x3

19

x2 .

 

108

216

 

 

108

216

 

 

72

 

 

 

 

 

 

 

 

72

 

 

 

 

 

Общее решение уравнения есть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y C

C

x C

e2 x

C

e3x

 

1

x4

5

x3

19

x2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

3

 

 

4

 

 

 

72

 

108

 

216

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача 4.6. Найти общее решение уравнения y 3y 2 y xex .

 

 

 

 

 

 

 

Решение. Характеристическое уравнение K 2 3K 2 0

имеет простые корни K1 1 и

K

2

2 . Значит,

y

0.0.

C ex C

e2 x

. Переходим к отысканию

y ч.н..

Правая часть уравнения

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеет вид f (x) xex . Здесь 1,

0,

s 1. Число i 1 является корнем кратности

1 характеристического уравнения, значит k =1, y ч.н. Ax B xex .

 

 

 

 

 

 

Подставив

 

 

y ч.н.

 

и

 

ее

производные

 

в

исходное

уравнение, получим

ex Ax2 4Ax Bx 2A 2B 3ex Ax2 2Ax Bx B 2ex Ax2 Bx xex .

 

 

 

 

Сокращаем на ex

 

и приводим подобные члены:

2Ax 2A B x .

 

 

 

 

Приравнивая коэффициенты при одинаковых степенях

x

 

в левой и правой частях

тождества, находим: A 0,5,

B 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем y ч.н.

x

 

 

x ex .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение уравнения имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y y0.0. y ч.н. C e

x

 

 

2 x

x2

 

 

x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

e

 

 

 

 

x

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

Задача 4.7. Найти общее решение уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 2 y 5y 17 cos 2x 16e x .

 

 

 

 

 

 

 

 

Решение.

K 2

2K 5 0,

 

K1,2 1

4 1 2i,

1,

 

2,

y

0.0.

e x C cos 2x C

2

sin 2x . Так как правая часть уравнения – сумма слагаемых вида

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.16), то согласно замечанию 2 пункта 6.3.2 частное решение ищем в виде суммы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ч.н. Acos 2x B sin 2x Сex .

 

 

 

 

 

 

 

 

Подставляя y ч.н. в исходное уравнение, получаем тождество

 

 

 

 

 

 

 

4Acos2x 4Bsin2x Cex

2 2Asin2x 2Bcos2x Cex

 

 

 

 

 

5 Acos2x Bsin2x Cex 17cos2x 16ex

0.

 

 

 

 

Приравнивая к нулю суммарные коэффициенты при cos 2x , sin 2x и ex, получаем

 

 

 

 

cos 2x : 4A 4B 5A 17;

 

 

sin 2x : 4B 4A 5B 0 ;

 

 

ex: C 2C 5C 16.

Решая систему, находим: A 1,

 

B 4 , С = 2, тогда

y ч.н. cos 2x 4sin 2x 2e x . Для общего

решения получаем выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y C e x

cos 2x C

e x sin 2x cos 2x 4sin 2x 2ex .

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

185

Задача 5. Найти общее решение линейной неоднородной системы д. у. 2-го порядка методом исключения неизвестных:

 

 

 

 

 

 

 

 

 

 

 

y' x y t ,

x' 3x 4 y 2t .

 

 

 

 

(6.34)

 

 

 

Решение. Здесь t – независимая переменная, х = x(t), y = y(t) – неизвестные функции.

Из

первого

уравнения

системы (6.34)

 

выражаем

x y' y t и

 

находим

производную

x' ( y' y t)' y'' y' 1 .

Подставляем

 

х

и

 

х' во

второе уравнение (6.34): y'' y' 1

3( y' y t) 4 y 2t .

Приводя подобные слагаемые, получаем уравнение относительно

функции

у:

y'' 2 y' y 5t 1. Общее

 

решение этого уравнения

y (C C

t)e t 5t 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

( y

o.o.

(C

C

t)e t ,

y

ч.н.

5t 9 , y y

o.o.

y

ч.н.

). Находим

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x y' y t (C

e t (C

C

t)e t 5) (C e t C

te t 5t 9) t 2C

te t (C

2

2C )e t

6t 14.

 

 

2

 

 

1

2

 

 

1

 

 

 

2

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение системы есть:

y (C1 C2t)e t 5t 9 , x (C2 2C1 2C2t)e t 6t 14 .

Задача 6.1. Найти общее решение системы д. у. 2-го порядка методом Эйлера: x' x y , y' 2x 4 y .

Решение. Согласно методу Эйлера составляем характеристическое уравнение:

 

0

 

1

 

(1 r)(4 r) 2 r 2 5r 6 .

 

1 r

 

 

 

2

4 r

 

 

 

 

 

 

 

Корни уравнения = 0

есть r1 = 2, r2 = 3.

Система (6.26) имеет вид: (1 r) 0 ,

2 (4 r) 0 , причем

второе

уравнение

лишнее.

Имеем:

(1 r) .

Полагая

1 1, r r1 2 , получаем 1 1. Полагая 2

1, r r2

3 , получаем

2 2 .

Подставляя найденные значения в (6.29), находим общее решение системы:

 

x C e2t

C

e3t , y C e2t 2C

e3t .

 

 

 

 

1

2

 

1

2

 

 

 

Задача 6.2. Решить линейную однородную систему д. у. методом Эйлера: x' 2x y, y' x 4 y .

Решение. Составляем характеристическое уравнение:

 

0

2 r

1

 

 

(2 r)(4 r) 1 r 2 6r 9 .

 

 

 

1

4 r

 

 

 

 

 

 

 

 

Его корень r r1 r2

3

(кратный корень). Решение

 

системы ищем в виде

x ( 1 1t)e3t ,

y ( 2 2t)e3t . Подставляя х, у в первое уравнение системы получаем:

 

 

e3t 3(

t)e3t

2(

t)e3t (

2

 

2

t)e3t ,

 

1

 

1

1

 

 

1

1

 

 

 

1 3( 1 1t) 2( 1 1t) ( 2 2t) ,

1 3 1 3 1t 2 1 2 (2 1 2 )t .

Приравнивая

коэффициенты

при

 

одинаковых

степенях t слева и справа,

имеем:

1 3 1 2 1 2 ,

3 1 2 1 2 ,

откуда 2 1 1 ,

2

 

1 . Обозначая 1

С1 ,

1 С2 ,

получаем общее решение системы

 

 

 

 

 

 

 

 

 

 

x (C C

t)e3t , y (C C

2

C

t)e3t .

 

 

 

 

1

2

1

 

2

 

 

 

186

 

Задача 6.3. Решить систему

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x' x 5y, y' 2x y .

 

 

 

 

 

Решение. Характеристическое уравнение

 

 

 

 

 

 

 

 

 

0

 

1 r

5

 

(1 r)( 1 r) 10 r 2 9 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1 r

 

 

 

 

 

 

 

 

Имеем комплексные корни: r 2

9,

 

r 9,

r

3i (т. е. 0, 3 ). Система

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

 

 

 

(6.26)

есть:

(1 r) 5 0, 2 ( 1 r) 0 .

Из

первого

уравнения

находим

(1 r) / 5 .

 

Полагаем

1

5, r 3i ,

тогда

1 1 3i ;

значит,

x

er1t 5e3it ,

y

 

er1t

(1 3i)e3it

. Выделяем действительную и мнимую части х1, у1 с

1

1

 

1

1

 

x1 5(cos 3t i sin 3t) , поэтому Re x1 5cos 3t, Im x1 5sin 3t;

помощью формул Эйлера:

y

(1 3i)e3it (1 3i)(cos3t i sin 3t) cos 3t 3i cos 3t i sin 3t 3i2 sin 3t

cos 3t 3sin 3t

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i( 3cos 3t sin 3t) , так как i2 1 , поэтому Re y cos 3t 3sin 3t,

Im y

3cos 3t sin 3t .

 

Общим решением системы будет

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

x C1 Rex1 C2 Im x1

5C1 cos 3t 5C2 sin 3t ,

 

 

 

 

 

y C1 Re y1 C2

Im y1 C1 (cos3t 3sin 3t) C2 ( 3cos 3t sin 3t) (C1

3C2 ) cos 3t

 

(3C1

C2 ) sin 3t.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчетное задание

Задача 1. Найти общий интеграл дифференциального уравнения.

1. yy' x 1.

3.1 1 y' e y 0 .

5.y2 x2 y' 0 .

7. 2 1 ex yy' ex . 9. 2x2 yy' y2 1 .

Задача 2. Найти решение задачи Коши.

1.

xy' 2 y 2x3ex , y 1 0 ;

2.

y' y sin x xe cos x , y 0 0 ;

3.

y' y ex x , y 1 0 ;

4.xy' y 4x3 2x , y 1 0 ;

5.xy' lnyx 2x2 ln x , y e e2 ;

6.xy' y 4x3 , y 1 1;

7. y' 2 y 3ex 2 , y 0 1;

8.xy' y 3x3 2 , y 1 0 ;

9.y' y e2 x 1 , y 0 1;

10.y' xy 3x3 3x , y 1 0 .

2. x xy yy' 1 x 0 . 4. tg y dx tg x dy 0 . 6. ex 1 y2 2 y 1 ex y' . 8. ex 1 y 1 ex y' 0 .

10. e y 1 x2 y' 2x 1 e y 0 .

187

Задача 3. Найти общий интеграл уравнения.

1.x 2 9xy2 dx 4 y2 6x3 ydy 0 .

2.x 2x2 y2 y x2 2 y2 y' 0 .

3.x 2 y4 xdx 1 4x2 y2 ydy 0 .

4.x3 x y4 y2 x4 y 2 y' 0 ..

5.x 3y 2 x dx 3x2 y 6 y 2 1 dy 0 .

6.cos x y2 sin x dx 2 y cos x sin y dy 0 .

7.sin y 1 y cos x dx 1 x cos y sin x dy 0 .

8.3x2 y y3 cos x x3 3xy2 sin y y' 0 .

9.sin y dx 1 x cos y dy 0 .

10.x cos y sin x dx y sin y cos x dy 0 .

Задача 4. Решить уравнение.

1.

y''' y'' 2 y 2 4x 6e x .

2.

 

y'' y' 20e2 x 2x .

3.

y'' y x 2 2e x 2ex .

4.

 

y''' y'' 2 y' 4x sin x 3cos x .

5.

y'' 2 y' 2 y 5x 4 ex e x .

6.

y''' 2 y'' y' 2x 2e2 x .

7.

y'' 3y' 18x cos x 3sin x .

8.

 

y'' 2 y' 8y 9ex 4cos 2x 12sin 2x .

9.

y''' y' 6e2 x 2x .

10.

y'' 4 y 8e 2 x 3sin x 3cos x .

Задача 5. Решить систему методом исключения неизвестных.

1.x' y t .y' x t

3.x' 2x 4 y cos t .

y' x 2 y sin t

 

 

 

 

2t

 

5.

x' y 2x e

2t .

6e

 

y' 2 y 3x

 

7.

x' x

2 y

 

.

 

 

5sin t

 

 

y' x

 

 

9.

x' 3

2 y

.

 

 

 

 

 

 

 

y' 2x 2t

 

 

 

2.x' 3x y et .y' x 3y et

4.

x' y x t

 

.

 

4x 2t

 

y' 3y

 

6.

x' 4x y

 

t .

 

y' 2x y 2e

 

 

8.x' x y et .y' 3x y

10. x' t 2 y

y' t x .

Задача 6. Решить однородную систему дифференциальных уравнений методом Эйлера.

1.

x' 3x

5y

 

 

.

 

 

8y

 

y' 2x

 

3.

x' 3x y

.

 

 

 

 

 

 

 

 

y' x y

 

 

 

 

5.

x' 3x

8y

 

.

 

 

 

 

 

 

 

y' 3y x

 

 

7.

x' x 5y

 

.

 

 

 

 

 

 

y' 3y x

 

 

 

9.

x' 2x

y

 

.

 

 

 

2 y

 

 

 

y' 3x

 

 

 

2.x' x 4 y .y' x y

4.

x' 3x

y

.

 

 

 

y' 4x y

 

6.x' 2x y .y' 4 y x

8.

x' x y

 

.

 

2x

 

y' 4 y

 

10.

x' 4x 2 y

 

 

y' y x .

188

ГЛАВА 7. ЧИСЛЕННЫЕ МЕТОДЫ И ИХ РЕАЛИЗАЦИЯ В СИСТЕМЕ

MATHCAD

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

Для ознакомления студентов с основными численными методами, предусмотренными программой курса высшей математики для инженерных специальностей, их особенностями и возможностями разработан цикл лабораторных работ, который включает в себя следующие работы:

решение систем линейных уравнений методом Гаусса;

решение нелинейных уравнений;

вычисление определенных интегралов;

численное интегрирование обыкновенных дифференциальных уравнений первого порядка;

аппроксимация функций на основании экспериментальных данных по методу наименьших квадратов.

Весь теоретический материал, необходимый для выполнения работ, можно почерпнуть из конспекта лекций и из учебников [6]–[9].

В данном учебном пособии для каждой лабораторной работы (подразделы 7.1 – 7.5) приведены теоретические сведения, постановки задач и основные численные методы их решений, примеры выполнения работ, варианты заданий и вопросы по прочитанному материалу.

Необходимый для выполнения лабораторных работ минимум знаний о математической системе Mathcad содержится в подразделе 7.6.

Для выполнения лабораторных работ разработан специально подготовленный в прикладной системе Mathcad пакет программ Lab, включающий пять файлов:

Lab1.mcd – решение систем линейных уравнений методом Гаусса; Lab2.mcd – решение нелинейных уравнений;

Lab3.mcd – вычисление определенных интегралов;

Lab4.mcd – численное интегрирование обыкновенных дифференциальных уравнений первого порядка;

Lab5.mcd – аппроксимация функций на основании экспериментальных данных по методу наименьших квадратов.

Для удобства реализации численных методов в системе Mathcad, в данном учебном пособии приведены тексты программ.

Прежде чем приступать к выполнению своего задания, необходимо в каждом файле рассмотреть пример, для которого исследование уже проведено. Далее нужно руководствоваться подсказками, указаниями и заданиями, выделенными в тексте программ жирным шрифтом. В каждом файле помимо сведений, необходимых для выполнения работы, представлены также некоторые дополнительные возможности системы Mathcad по соответствующему разделу математики.

189

7.1. Решение систем линейных уравнений методом Гаусса

7.1.1. Постановка задачи

Пусть задана система n линейных алгебраических уравнений (СЛАУ) с n неизвестными

x1, x2,...,xn :

 

 

 

 

 

 

 

 

a(1) x

a(1) x

2

... a(1) x

n

b(1)

,

 

 

 

 

 

 

 

 

11 1

12

1n

1

 

 

 

 

 

 

 

 

 

a(1) x

a(1) x

2

... a(1) x

n

b(1)

,

 

 

 

 

 

 

 

 

21 1

22

2n

2

 

 

 

 

 

 

 

 

 

...................................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a(1) x

a(1) x

2

... a(1) x

n

b(1)

,

 

 

 

 

 

 

 

 

n1 1

n2

nn

n

 

или в матричной форме: A X B , где

 

 

 

 

 

a(1)

a(1)

 

....

a(1)

 

 

 

 

 

 

 

 

 

11

12

 

 

 

1n

 

 

 

 

 

 

 

 

a(1)

a(1)

 

....

a(1)

 

основная матрица системы,

 

A

21

22

 

 

 

2n

 

 

 

.... .... .... ....

 

 

 

 

 

 

 

 

 

(1)

(1)

 

....

(1)

 

 

 

 

 

 

 

 

an1

an2

 

ann

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

элементов,

X

x2

 

столбец неизвестных.

 

 

 

 

....

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

(7.1)

b(1)

 

 

 

1

 

 

 

(1)

 

столбец свободных

B b2

 

 

....

 

 

b(1)

 

 

 

n

 

 

Решением системы (7.1) называется совокупность значений неизвестных x1, x2 ,..., xn ,

удовлетворяющая одновременно каждому уравнению из системы (7.1). Система решена полностью, если все решения найдены.

Теорема Кронекера-Капелли: Для того, чтобы система (7.1) была совместна (имела хотя бы одно решение), необходимо и достаточно, чтобы ранг основной матрицы A и ранг расширенной матрицы системы (основная матрица системы с добавлением справа столбца свободных элементов)

 

a(1)

 

 

11

 

a(1)

 

A

21

 

 

....

 

 

(1)

an1

были равны: rang A rang A r .

a12(1) a22(1)

....

an(12)

.... a1(1n)

.... a2(1n)

.... ....

.... ann(1)

|

(1)

 

 

b1

 

 

|

b(1)

 

(7.2)

 

2

 

| ....

 

 

|

b(1)

 

 

 

n

 

 

При этом, если ранг равен числу неизвестных r n , то система (7.1) имеет единственное решение, т. е. определена. Если r n , то система (7.1) имеет бесконечное множество решений, зависящих от (nr) произвольных параметров, т. е. неопределена.

Существует много методов решения таких систем [6]–[8]. В данной лабораторной работе будем решать ее методом Гаусса с частичным выбором ведущего элемента. Суть

этого метода состоит в последовательном исключении неизвестной x1 из 2, 3,…, n-го уравнений, x2 – из 3, 4,…, n-го уравнений и т. д. Для этого на каждом шаге преобразования сначала выбираем так называемое «ведущее уравнение». На i-м шаге (т. е. при исключении неизвестного xi , i=1, 2,…, n–1) в качестве ведущего уравнения нужно взять из i-го, (i+1)-

го,…, n-го уравнений то уравнение, в котором коэффициент перед xi имеет наибольшую абсолютную величину. Ведущее уравнение ставим на место i-го уравнения, и во всех ниже

190

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