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

21

Метод Эйлера. Значения искомой функции находят по формуле y(xk+1 ) = yk+1 = yk + h f (xk , yk ) ;

Метод разложения в ряд Тейлора. Решение представляют в виде не-

скольких первых членов ряда:

 

 

 

 

 

y(x) = y(x0 ) + y(x0 )(x x0 ) +

 

y′′(x0 )

(x x0 )2 KK=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2!

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

x0 ) +

 

 

x0 )

K

 

 

 

 

 

= y(x0 ) + f (x0 , y0 )(x

2!

f (x0 , y0 )(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задания контрольной работы

 

Проинтегрируйте дифференциальные уравнения:

 

 

1.1

(1+ y)(ex dx e2 y dy) (1+ y2 )dy = 0;

2)

 

 

 

1.2.

y′ = 4x + 2 y 1

 

 

 

 

 

 

 

 

1.3.

xy′− y = x2 y2

;

 

 

 

 

 

 

 

 

x( y′+ e

y

 

 

 

 

 

 

 

 

1.4.

 

) = y ;

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

1.5.

( y2 2x2 )dy + 2xydx = 0, y(1) =1;

 

 

 

 

1.6.

dy = ( ytgx 1)dx,

y(0) = 3;

 

 

 

 

 

 

 

1.7.

(1+ y2 )dx = (arctgy x)dy ;

 

 

 

 

 

 

 

 

 

 

 

 

y3

 

 

 

 

 

 

 

 

 

1.8.

y

= yctgx + sin x ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.9.(x +1)(y′+ y2 ) = −y

1.10.yx3 sin y = xy′−2y ;

1.11.

y

′′

= xe

x

= 0

 

 

, y(0) =1, y (0)

1.12.

y′′+ (tgx) y′ = sin 2x ;

 

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

Задача 1. Поглощение света при прохождении через воду Поглощение светового потока тонким слоем воды пропорционально толщине слоя и потоку, падающему на его поверхность. Зная, что при прохождении

через слой толщиной 2м поглощается первоначального светового потока, определить, какой процент его дойдет до глубины 12м ?

 

 

 

22

 

Решение.

 

 

 

Составим дифференциальное уравнение. Обозначим через

световой

поток,

падающий на поверхность на глубине . При прохождении через

слой воды толщиной

поглощенный световой поток

равен диф-

ференциалу

,

где – коэффициент пропорциональности

(

).

 

 

 

Общее решение дифференциального уравнения получаем путем разделения переменных . В результате общее решение имеет вид:

 

 

 

 

.

 

 

 

 

 

 

 

По условию задачи при

имеем

 

 

поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

и

 

 

,

 

До глубины

 

 

 

 

 

м дойдет световой поток

 

 

 

 

 

 

 

 

 

 

 

 

что составляет 8,78 первоначального светового потока.

Задача 2

Найти кривую, у которой сумма длин касательной (точнее длины её отрезка от точки касания до точки пересечения с осью абсцисс) и подкасательной в любой её точке равна произведению координат точки касания.

Решение.

Пусть

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

вольной точке

кривой

Согласно условию задачи

 

.

 

Из прямоугольного треугольника

Тогда дифференциальное уравнение примет вид

 

 

 

,

 

 

23

Умножая обе части полученного уравнения на дробь

 

 

 

получим

 

 

.

Преобразуем его. Возводим обе части в квадрат.

Разделим обе части на (при условим, что ). откуда

Разделяя переменные и интегрируя, получим общее решение

Задания контрольной работы

По заданным условиям составить дифференциальное уравнение и решить его.

1.13. Скорость остывания пропорциональна разности температур тела (T) и окружающей среды Тс = 20o С. За 10 минут тело остыло со 100oC до

60oC . За какое время тело остынет до 25oC ?

1.14. Найдите ток I (t)в электрической цепи, заданной уравнением

LI′+ RI = E sinωt, I (0) = 0, E = const, ω = const

1.15.В начале семестра один день тянется как два. Через три месяца два дня тянутся как один. К началу сессии времени ни на что не хватает. Найдите

время начала сессии t0 , если скорость течения времени обратно пропорциональна квадратному корню из времени, оставшегося до сессии.

1.16.Найти кривую, проходящую через точку , зная, что угловой коэффициент касательной в любой точке кривой в три раза больше углового коэффициента прямой, соединяющей эту же точку с началом координат.

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

Теоретическое введение

Общее решение линейного однородного дифференциального уравнения порядка n с постоянными коэффициентами

y(n) + an1 y(n1) + an2 y(n2) +KK+a1 y′+ a0 y = 0

24

находят как линейную комбинацию фундаментальных (линейно независимых или базисных) решений y1 (x), y2 (x),KKyn (x) :

y00 =C1 y1 (x) + C2 y2 (x) +KKCn yn (x) .

Базисные решения определяются корнями характеристического уравнения, которое получается из исходного уравнения заменой производных

y(k ) на λk :

λ(n) + an1λn1 + an2λn2 +KKa1λ +a0 = 0.

Каждому действительному корню λ кратности r в общем решении

соответствует решение:

(C0 +C1x +C2 x2 +KKCr1xr1 ) eλx ;

Каждой паре комплексно сопряженных корней α ±i β кратности r

соответствует решение:

e

α x

[(C0

+C1x +KKCr1x

r1

ˆ

ˆ

ˆ

r1

)sin βx].

 

 

)cos βx + (C0

+ C1xKKCr1x

 

Согласно этому правилу решения однородных уравнений 2-го порядка

записываются следующим образом:

Если корни характеристического уравнения λ2 + a1λ +a0 = 0 действительные и различные λ1 λ2 , то y00 = C1eλ1 x +C2 eλ2 x ;

Если корни характеристического уравнения λ2 + a1λ +a0 = 0

действительные и одинаковые (кратность к = 2) λ1 = λ2 = λ0 , то

y00 = ( C0 +C1x) eλ0 x ;

Если корни характеристического уравнения λ2 + a1λ +a0 = 0

комплексные λ1,2 =α ±i β , то y00 = eα x (C1 cos βx +C2 sin βx) .

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

y(n) + an1 y(n1) + an2 y(n2) +KK+a1 y′+ a0 y = f (x)

находят как сумму y(x) = y00 (x) + y(x) , где y(x) - любое частное решение исходного неоднородного уравнения, y00 - общее решение соответствующего однородного уравнения. При этом частное решение y(x) определяется правой

25

частью уравнения f (x) и может быть найдено методом вариации произвольной постоянной или методом подбора по правой части

специального вида

f (x) = (an xn +an1xn1 +an2 xn2 +KK+a1x + a0 ) e(α+iω) x .

Если число α ± iω является корнем характеристического уравнения

(прочитайте в пояснениях к заданию 2, что такое мнимая единица, комплексное число, формула Эйлера), соответствующего данному неоднородному уравнению, кратности r , то частное решение подбирается в

виде:

 

 

 

 

 

 

) e(α+iω) x .

 

 

y = xr (b xn +b

xn1 KKb x +b

 

 

 

n

n1

 

 

1

0

 

 

 

 

 

При этом рассматриваются частные случаи:

 

 

 

 

 

 

Правая часть – многочлен

f (x) = a

n

xn +a

n1

xn1

+KK+a x +a

0

, то и

 

 

 

 

 

 

1

 

решение подбирается в виде многочлена той же степени, но с произвольными коэффициентами

y = (bn xn +bn1xn1 KKb1x +b0 ) xr ,

где r - кратность корня λ =α ±iω = 0 характеристического уравнения;

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

f (x) = (an xn +an1xn1 +KK+a1x +a0 ) eα x ,

то и частное решение подбирают в виде экспоненты, умноженной на многочлен той же степени:

y = eα x (bn xn +bn1xn1 KKb1x +b0 ) xr ,

где r - кратность корня λ =α ±iω =α характеристического уравнения;

Правая часть-многочлен, умноженный на тригонометрическую функцию

1)f (x) = (an xn +an1xn1 +KK+a1x +a0 ) соsω x

2)f (x) = (an xn +an1xn1 +KK+a1x +a0 ) sinω x

или сумма этих выражений. Тогда частное решение подбирается в виде

z = eiω x (bn xn + bn1xn1 KKb1x + b0 ) xr ,

26

где r - кратность корня λ =α ±iω = ± iω характеристического уравнения.

При этом в первом случае y = Re z , а во втором случае y = Im z .

Замечание. Также в случае, когда правая часть уравнения имеет вид

f (x) = (ak xk + ak 1xk 1..... + a1x + a0 ) cosωx +(lm xm +lm1xm1....l1x +l0 )sinωx

и, ± iω не являются корнями характеристического уравнения, частное решение в этом случае может быть подобрано в виде

y = (cn xn + cn1xn1 +KK+ c1x + c0 ) соs ω x + (bn xn + bn1xn1.....b1x + b0 ) sin ωx

,где max(k, m) .

Втом случае ± iω являются корнями характеристического уравнения

кратности r , частное решение в этом случае может быть подобрано в виде

y = xr ((cn xn + cn 1xn 1 +KK+ c1x + c0 ) соsω x +(bn xn +bn 1xn 1.....b1x +b0 )sinωx)

Задания контрольной работы

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

коэффициентов′′

1.17. y 4 y = sin x

1.18. y′′+ 3y′+ 2 y = exp 2x

1.19. y′′+ 2 y′+ y = cos x

1.20. y′′+ 4 y′+5y = x

1.21. y′′−4y = e2 x

1.22. y′′+ 2y′+ y = x2 1

1.23. y′′+ 4y′+5y = e3 x

27

1.24.y′′+ 9 y = sin 3x

1.25.y(4) + y′′ = x,

1.26.y(4) + y′′ = x2

1.27.y′′−6y′+ 25y = e4x (17x2 + 38x + 40)

1.28.y′′−9 y = sin 3x

Контрольные вопросы

1.Что является решением дифференциального уравнения? Дайте определение общего и частного решений дифференциального уравнения

2.Определите тип дифференциального уравнения:

 

a) x dy y = y2ex ; b) dy

=

y

+

y3

, c) dy 3x2 + y = 0; d) y dy + x3

= 0

 

x

x3

 

 

dx

dx

 

 

dx

 

dx

 

3.

Дифференциальное уравнение семейства кривых

y = (C +C x)ex

имеет вид:

 

 

a) y′′− y′−2y = 0 ;

 

b) y′′−2y′+ y = 0

1

2

 

 

 

 

 

 

 

4.

Частное решение дифференциального уравнения y′′−3y′−4y = xex имеет

 

вид:

a)

y = xex (ax +b) ;

 

 

 

b)

y = (ax + b)ex

 

 

 

5. Общее решение дифференциального уравнения имеет вид …

ЗАДАНИЕ 2. Функции комплексной переменной

1.Комплексные числа. Алгебраическая, тригонометрическая, показательная форма записи. Действия с комплексными числами

2.Функции комплексной переменной. Аналитические функции

3.Вычет аналитической функции в изолированной особой точке. Вычисление контурных и несобственных интегралов при помощи вычетов.

4.Преобразования Лапласа.

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

28

6.Решение линейных дифференциальных уравнений методом свертки (формулы Дюамеля, Грина).

Литература [1,2, 3, 4, 5]

Тема 1. Комплексные числа. Алгебраическая, тригонометрическая, показательная форма записи. Действия с комплексными числами

Теоретическое введение

Комплексные числа

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

х2 +1 = 0(поскольку квадрат любого действительного числа неотрицателен). Поэтому возникла потребность расширить множество действительных чисел таким образом, чтобы новое множество содержало корни всех алгебраических уравнений. Введение комплексных чисел позволяет достигнуть этой цели. Прежде всего введем новый символ i , который называют мнимой едини-

цей таким образом, что i2 = −1. Тогда корни уравнения x2 +1 = 0 запишутся как x = ±1 = ±i .

Алгебраическая форма записи комплексного числа

Комплексное число в алгебраической форме записывается как z = x + iy .

Здесь x – действительное число, называемое реальной или действительной частью комплексного числа. Обозначают: x = Re z . Действительное число y называют мнимой частью комплексного числа. Обозначают: y = Im z .

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

тельных чисел (х; у). Если y = 0 , то комплексное число совпадает с действительным и изображается точкой на действительной оси ОХ. При х = 0 получаются чисто мнимые числа (0; у)= 0 + iy , которые изображаются точкой на мнимой оси OY. Комплексное число z = (x; у)= x + iy можно отождествить с

точкой плоскости OXY или радиусом – вектором rr = (x; у)= xi + yj . Плоскость OXY будем называть комплексной плоскостью (рис. 1)

29

Рис. 1

Два комплексных числа равны, если равны их действительные и мнимые части:

x1 + iy1 = x2 + iy2 x1 = x2 ; y1 = y2 .

Два комплексных числа называют сопряженными, если действительные

части этих

чисел

 

равны, а мнимые

отличаются знаком. Обозначают:

z = x iy =

(x,y).

 

 

Над комплексными числами в алгебраической форме определены сле-

дующие операции:

 

 

 

а) сложение

 

 

 

 

z1 + z2 = x1 + iy1 + x2 + iy2 = (x1 + x2 )+ i(y1 + y2 ).

Сумма

z + z

 

изображается вектором – суммой векторов (x , у ) и

(х2 , y2 ).

1

 

2

 

1 1

б) вычитание

 

 

 

 

z1 z2 = x1 + iy1 (x2 + iy2 )= (x1 x2 )+ i(y1 y2 ).

Разность z1

z2 изображается вектором – разностью векторов (x1 , у1 ) и

(х2 , y2 ).

 

 

 

 

 

 

 

 

 

z1–z2

z1+ z2

 

 

 

 

z1

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

в) умножение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z1

z2 =

(x1 + iy1 ) (x2 + iy2 )= x1 x2 + ix1 y2 + iy1 x2 y1 y2 =

 

 

 

 

=

(x1 x2 y1 y2 )+ i(x1 y2 + y1 x2 ).

 

 

 

 

 

 

 

 

 

 

 

Особо отметим произведение комплексно-сопряженных чисел, ко-

торое является числом действительным:

 

 

 

 

 

 

 

 

 

 

 

 

z z = (x + iy) (x iy)= x2 (iy)2 = x2 + y2 .

 

 

г) деление

 

 

 

(x1 + iy2 )(x2 iy2 )

(x1 x2 + y1 y2 )+ i(y1 x2 x1 y2 )

 

 

z1

 

x1 + iy1

 

 

 

 

 

 

=

 

 

 

 

=

 

 

 

 

 

 

 

 

 

=

 

x 2

+ y 2

 

.

 

z

2

x

+ iy

2

 

(x

2

+ iy

2

)(x

2

iy

2

)

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

2

 

 

Пример 1.

Найти значение функции

f (z)= z4 +

2 +i

(3 + 2i)

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

z =1 2i .

f (1 2i)= (1 2i)4 +12+2ii (3 + 2i).

Для вычисления первого слагаемого воспользуемся формулой бинома Ньютона:

(1 2i)4 =1 + 4(2i)+ 6(2i)2 + 4(2i)3 + (2i)4 =

=1 8i 24 + 32i +16 = −7 + 24i

Второе слагаемое есть частное от деления двух комплексных чисел. Ис-

пользуя правило деления комплексных чисел в алгебраической форме, получим:

 

2 + i

 

(2 + i)(1 + 2i)

5i

 

=

 

=

 

= i .

 

1 2i

(1 2i)(1 + 2i)

5

Итак,

f (1 2i)= (7 + 24i)+ i (3 + 2i)= −4 + 23i .

Пример 2.

Решить уравнение (2 +i)z2 (5 i)z + (2 2i)= 0 .

По формуле для корней квадратного уравнения имеем:

z = 5 i ± (5 i)2 4(2 + i)(2 2i) = 5 i ± − 2i .

1,2

 

2(2 + i)

 

4 + 2i

 

 

 

 

Извлекая корень квадратный из числа 2i , получим:

2i = 1 2i 1 = (1 i)2 = ±(1 i).

Следовательно,

Соседние файлы в папке Методические указания для контрольной работы