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

5561

.pdf
Скачиваний:
2
Добавлен:
13.11.2022
Размер:
1.87 Mб
Скачать

За базисную переменную можно взять α1; тогда α2 и α3 – свободные переменные. Полагаем α2 = С2, α3 = С3 2, С3 – произвольные постоянные); тогда α1 = С2 + С3. Кроме того, так как С = 0, то β1 = β2 = β3 = 0. Подставляем найденные коэффициенты в (11.12), получим частные решения, соответствующие двукратному корню k = 3.

Таким образом, общее решение исходной системы таково:

x C et

C

e3t

C

e3t ,

 

1

2

 

 

3

 

y

C et

C

e3t

 

 

,

 

1

2

 

 

 

 

z

C et

 

 

C

e3t ,

 

1

 

 

3

 

 

где C1 , C2 и C3 – произвольные постоянные.

Пример 7. Найти общее решение системы дифференциальных уравнений

dx

 

3x

y

z,

 

 

dt

 

 

 

 

dy

 

x

y

z,

 

 

dt

 

 

 

 

 

dz

 

2x

 

2z.

 

 

 

dt

 

 

 

 

 

 

Решение. Характеристическое уравнение для этой системы таково:

 

 

3

k

1

1

 

 

 

 

 

 

 

 

1

1 k

1

 

 

0.

 

 

 

 

 

 

2

0

2

k

 

 

 

 

 

Корнями характеристического уравнения

k 3

 

6k 2 12k

8

0 являются

действительные числа

k1

k2

k3

2 ,

т.е.

корень

k

=

2 является

трёхкратным.

 

 

 

 

 

 

 

 

 

 

 

Соответствующие этому трёхкратному корню решения надо искать в виде

 

 

 

x = (α1 + β1t + γ1t2)e2t,

 

 

 

 

 

 

 

y = (α2 + β2t + γ2t2)e2t,

 

 

 

(11.13)

 

 

 

z = (α3 + β3t + γ3t2)e2t.

 

 

 

 

Для определения

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

функции

(11.13)

и

их

производные

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

β1 = α1 + α2 – α3,

1 = β1 + β2 – β3,

γ1 + γ2 – γ3 = 0,

β2 = α1 – α2 + α3,

2 = β1 – β2 + β3,

γ1 – γ2 + γ3 = 0, (11.14)

β3 = 2α1,

γ3 = β1,

γ1 = 0.

 

61

 

Первая группа уравнений (11.14) выражает коэффициенты β1, β2, β3 через α1, α2, α3. Подставляя их во вторую группу уравнений, получим выражения γ1, γ2 , γ3 через α1, α2, α3:

γ1 = 0, γ2 = α1 + α2 – α3, γ3 = α1 + α2 – α3.

Третья группа уравнений (11.14) автоматически обращается в тождество.

Ранг системы алгебраических уравнений (11.14) равен 6, а свободных переменных – три. В качестве свободных неизвестных можно взять α1, α2, α3. Пусть α1 = С1, α2 = С2, α3 = С3 1, С2, С3 – произвольные постоянные). Тогда базисные переменные выражаются через них следующим образом:

β1 = С1 + С2 – С3,

γ1 = 0,

β2 = С1 – С2 + С3,

γ2 = С1 – С2 + С3,

β3 = 2С1,

γ2 = С1 – С2 + С3.

Подставляя решения системы (11.14) в (11.13) и преобразуя, получим

следующее общее решение исходной системы примера 7:

x C (1 t) C

t C

3

( t) e2t ,

1

 

 

 

2

 

 

 

 

y C (1 t 2 ) C

2

(1 t t 2 ) C

3

(t t 2 ) e2t ,

1

 

 

 

 

 

 

 

 

z C (2t t 2 ) C

t

2 C

3

(1 t 2 ) e2t .

1

 

 

2

 

 

 

 

 

 

По виду общего решения

легко

 

выписать три линейно независимых

решений заданной системы дифференциальных уравнений и её фундаментальную систему решений (предоставляем это сделать читателю).

Контрольные задания

В соответствующем номере варианта требуется найти:

1.Общее решение (общий интеграл) дифференциального уравнения с разделёнными переменными, записанного в дифференциальной форме.

2.Общее решение дифференциального уравнения с разделяющимися переменными, записанного в дифференциальной форме.

3.Общий интеграл однородного дифференциального уравнения первого

порядка.

4.Общий интеграл уравнения в полных дифференциалах.

5.Общее решение линейного дифференциального уравнения первого порядка методом Бернулли (методом подстановки) и методом Лагранжа (методом вариации произвольной постоянной).

62

6. Частное решение

линейного уравнения первого порядка, из

предыдущего номера, удовлетворяющее данному начальному условию у(х0) = у0.

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

8.Частное решение уравнения из предыдущего номера, удовлетворяющее заданным начальным условиям у(х0) = у00, y (х0) = у01.

9.Общее решение линейного однородного дифференциального уравнения второго порядка с правой частью f(x) при каждой из трёх её указанных видов.

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

11.Общее решение однородной системы с постоянными коэффициентами, состоящей из двух дифференциальных уравнений первого порядка с неизвестными функциями x(t), y(t).

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 1

 

 

 

 

 

 

 

 

 

 

 

 

1.

 

ydy

2x dx .

 

 

 

2.

y(1

 

x2 )dy

 

x(1 y2 )dx

0.

3.

y

 

y

 

tg

 

y

.

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

4.

(2xy

 

 

y

1)dx

(x2

x

2y2 )dy 0.

5.

y

y

x2 .

 

 

 

 

 

6.

у(0) = 1.

 

 

 

 

 

 

 

7.

y

6 y

10 y

0.

 

 

8.

у(0) = 5, y (0)

2.

 

 

 

9.

y

4 y

f (x);

 

 

 

1) f(x) = x2 + 1; 2) f(x) = (x + 2)e-2x; 3) f(x) = e3xcos 2x.

10. Скорость изменения количества населения есть первая производная по времени от количества населения. Предполагается, что скорость прироста населения прямо пропорциональна количеству населения. Найти зависимость между количеством населения N и временем t , если известно, что в некоторый

63

момент, принимаемый за начальный, количество населения составляло N0 , а

через год население увеличилось на n %. Найти предполагаемое на этой основе население города S на 1 сентября 2010 года, если известно, что 1 сентября 2008 года оно составляло 5,342 миллиона человек. Годовой прирост принять равным

1,5%.

 

dx

 

x 3y;

 

 

 

11.

dt

 

 

dy

 

 

 

 

2x 2 y.

 

 

 

 

dt

 

 

 

 

Вариант 2

1.y dy xdx .

2.ey(1 + x2)dy – 2x(1+ey)dx = 0.

3. y

x2

y2

xy

.

 

x2

 

 

 

 

 

4.

(2x 1

 

y

)dx (

1

2 y)dy 0.

 

 

 

 

 

 

 

 

x2

 

x

 

5.

y

 

(ctgx) y

sin x.

 

6.

y(

 

)

.

 

 

 

 

2

 

 

 

 

7.

y

 

y

 

6 y

0.

 

 

8.

y(0) = 2,

 

y (0)

5 .

 

 

9.

y 10 y 25y

f (x) ;

1) f(x) = 2x + 5; 2) f(x) = xe5x; 3) f(x) = ex(x cos 3x + sin 3x).

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

тела (согласно

механическому

смыслу

производной

 

dT

) прямо

 

 

 

 

 

 

 

 

 

 

dx

пропорциональна

 

разности

температур

T t

с

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

пропорциональности

k , где

T

температура

тела,

t температура

 

 

 

 

64

 

 

 

 

 

окружающего среды, x искомое время охлаждения. Из экспериментального опыта известно, что тело в комнате с температурой 240 С в течении 10 минут охлаждается с 900 С до 400 С. Через сколько времени от момента начала охлаждения тело, нагретое до 900 С, охладится до 270 С?

 

dx

 

4x 6 y;

 

 

 

11.

dt

 

 

dy

 

 

 

 

2x 5 y.

 

 

 

 

dt

 

 

 

 

Вариант 3

1.y dy – x dx = 0.

2.sin x dy = y ln y dx.

3. 3xy x2 dy 2xy y2 dx 0 .

4. (exy + 5y)dx + (ex + 5x + 1)dy = 0.

5. y

2

y

1

1

.

 

 

 

x

 

x

6.

y(1)

 

1

.

 

 

 

2

 

 

 

 

 

 

 

 

 

7.

y

2 y

y

0 .

 

8.

y(0) = 6,

y (0)

5 .

 

9.

y

6 y 13y

f (x) ;

 

 

1) f(x) = x2;

2) f(x) = xe3x;

3) f(x) = e3x sin 2x.

10. Пусть

y t

выпуск продукции некоторого предприятия. Скорость

увеличения

выпуска продукции

прямо пропорциональна с коэффициентом

пропорциональности k доходу ( p y y доход от продажи выпуска продукции).

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

y t

есть

dy

. Составить уравнение,

связывающее скорость изменения выпуска

 

 

 

dt

 

 

 

 

 

продукции и доход от продажи выпуска y t

по цене p y . Предполагается, что

с увеличением выпуска будет проходить насыщение рынка и цена товара

p y

будет падать. Известно,

что p y

7

3y .

Составить дифференциальное

уравнение для функции y t

. Найти функцию y t

при условии, что y 0

10.

 

 

 

 

 

65

 

 

 

l , где
3) f(x) = e6x(cos 2x + x sin 2x).

 

 

 

 

dx

3x

4 y;

 

 

 

 

 

 

 

 

 

 

 

 

11.

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

2x

6 y.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 4

1.

 

dy

 

 

 

cos x dx .

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

 

1 x2 dy

1 y2 dx 0 .

3.

 

y

 

 

y

ln

 

y

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x

 

 

 

 

4. (x + 2xy2)dx + (2x2y + y2)dy = 0.

5.

 

y

2

y

 

x2ex .

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

6.

 

y(1)

 

 

 

 

e .

 

 

 

 

 

 

 

7.

 

y

 

 

2 y

 

 

2 y

0 .

 

8.

y(0) = 8,

 

y (0)

2.

 

9.

 

y

 

36 y

 

 

f (x) ;

1) f(x) = 2x2 + 3x + 5; 2) f(x) = xe6x;

10. Пусть y t объём продукции некоторой отрасли, реализованной к моменту времени t . В экономической модели естественного роста полагают, что скорость выпуска продукции (акселерация) пропорциональна величине инвестиций I t с

коэффициентом пропорциональности 1l норма акселерации. Модель предполагает, что вся производимая отраслью продукция реализуется по

некоторой фиксированной цене p .

При таких условиях очевидно, что доход к

моменту времени t

составит Y t

py t . Величина

I t

в модели составляет

фиксированную

часть дохода

I t

m Y t

mpy t ,

где

коэффициент

пропорциональности m является нормой инвестиций

0

m

1 . Найти

выражение для объёма реализованной продукции

y

y t

, если известно, что

кривая спроса

p y

задаётся уравнением p y

5

y ,

норма

акселерации

1 l 2 , норма инвестиций m 0,6 ,

y 0

0,4 .

 

 

 

 

 

 

 

 

66

 

 

 

 

 

 

 

dx

 

2x 4 y;

 

 

 

11.

dt

 

 

dy

 

 

 

 

3x 3y.

 

 

 

 

dt

 

 

 

 

Вариант 5

1. y2 dy = x2 dx.

2. y1 x2 dy x1 y2 dx 0 .

 

y

 

 

2xy

 

 

 

 

 

3.

 

 

 

 

 

.

 

 

 

 

 

 

x2

 

y2

 

 

 

4.

(sin y + y + x)dx + (x cos y + x + y)dy = 0.

 

 

 

 

1 x

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

5.

y

 

y

e x .

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

6. y(1) = e.

 

 

 

 

 

 

 

 

7.

y

 

 

y

12 y

0 .

 

 

8.

y(0) = 3,

y (0)

5.

 

 

9.

y 12 y 36 y

f (x) ;

 

 

 

1) f(x) = 3x + 2;

2) f(x) = xe6x;

3) f(x) = ex cos 3x.

 

10.

Пусть

 

выпуск

продукции y t

является функцией времени.

Из

механического смысла производной скорость выпуска продукции представляет собой производную dydt . Пусть с y издержки производства, а цена единицы продукции есть p y . Рассматривается модель, в которой скорость выпуска

продукции прямо пропорциональна прибыли предприятия с коэффициентом пропорциональности k ( p y y c y представляет собой прибыль производства). Составить дифференциальное уравнение для нахождения

выпуска продукции

y t ,

если

цена единицы

продукции

задана

функцией

p y 4,1

y . Найти

выпуск продукции y t

при условии,

что

y 0 0,5.

Издержки

производства выражает функция с y 0,1y 3,

а

цену единицы

продукции – функция p y

4,1

y .

 

 

 

 

67

 

 

 

 

dx

 

 

2x

3y;

 

 

 

 

 

 

 

 

11.

 

dt

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

5x

4 y.

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 6

1.

 

 

dy

 

 

 

dx

 

0.

tgy

 

 

ctgx

 

 

 

 

2. dy = (y2 – 4)dx.

3.

x dy

 

 

x2

 

y2

y dx 0 .

 

4.

(3x2 y2 + 1)dx + (y – 2xy)dy = 0.

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

5.

y

y

e x .

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

6. y(1) = e.

 

 

 

 

 

7.

y

3y

 

10 y

0.

 

8.

y(0) = 10, y (0)

6 .

 

9.

y

4 y 13y

f (x) ;

 

 

1) f(x) = x3 + 5;

2) f(x) = xe2x;

3) f(x) = xe2x cos 3x.

10. Пусть

x t

количество радиоактивного вещества в момент времени t . Из

эксперимента известно, что скорость радиоактивного распада в каждый момент времени (единицу времени) прямо пропорциональна наличной его массе x t с

коэффициентом пропорциональности k k 0 . Предполагается, что происходит убывание массы вещества. Из механического смысла производной скорость радиоактивного распада представляет собой производную функции x t . За 100 лет распалась половина первоначального количества радиоактивного вещества (период полураспада составляет 1600 лет). Определить временной промежуток, по истечении которого останется 6% первоначального количества радиоактивного вещества.

 

dx

 

x 2 y;

 

 

 

11.

dt

 

 

 

 

 

 

dy

 

4x 3y.

 

 

 

 

dt

 

 

 

 

 

 

68

 

 

 

 

 

 

 

 

 

Вариант 7

1.

 

1

dy

 

2x dx .

 

 

 

 

 

 

y

 

 

 

 

2.

tg x dy y dx = 0.

 

 

 

 

 

 

 

 

 

3.

 

y

 

x2

y2 dx x dy .

4. (x + 4y2)dx + (8xy + ey)dy = 0.

5.

 

y

y sin x

sin x .

6.

 

y

 

 

2.

 

 

 

2

 

 

7.

 

y

2 y

y

0 .

8.

y(0) = 8, y (0)

12 .

9.

 

y

5y 6 y

f (x) ;

1) f(x) = 3x2 + 2x + 10; 2) f(x) = (x + 1)e3x; 3) f(x) = xe2x sin 3x.

10. В рамках некоторого экономического процесса предполагают, что доход

Y t , полученный к

моменту

времени

t определённой отраслью,

является

суммой инвестиций

I t и

величины

потребления C t : Y t = I

t + C t .

Предполагают, что произведение скорости увеличения дохода на коэффициент капиталоёмкости b прироста дохода пропорционально величине инвестиций. Напомним, что под капиталоёмкостью понимают показатель, характеризующий отношение основного капитала к чистому доходу (прибыли). Найти функцию дохода y y t , если известно, что величина потребления задаётся функцией

C 3t ; коэффициент капиталоёмкости прироста дохода b 1 3 , а Y (0) 4 .

 

 

 

dx

 

 

4x

5 y;

 

 

 

 

 

 

 

11.

 

dt

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

3x

6 y.

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вариант 8

1.

dy

 

 

 

dx

 

0 .

 

y

 

 

x

 

 

 

 

 

2. (1 + ex)y dy = ex dx.

69

 

y

x2

y2

 

3. y

 

 

 

0 .

 

x

 

 

 

 

 

 

4. (3x2y2 + 4xy + 2)dx + (2x3y + 2x2 – 3)dy = 0.

5.

y

y tg x

cos x .

6.

y(0) = 4.

 

 

 

7.

y

y

6 y

0.

8.

y(0) = 8,

y (0)

6 .

9.

y

10 y 25y

f (x) ;

1) f(x) = 2x3 + 5; 2) f(x) = xe-5x; 3) f(x) = xe2x cos x.

10. Из эксперимента известно, что скорость радиоактивного распада некоторого вещества в каждый момент времени (единицу времени) прямо пропорциональна наличной его массе m t с коэффициентом пропорциональности k k 0 . Полагается, что происходит убывание массы вещества. Из механического смысла производной известно, что скорость радиоактивного распада представляет собой производную массы. За 300 лет распалась половина первоначального количества (период полураспада составляет 1300 лет). Определить временной промежуток, по истечении которого останется 8% первоначального количества вещества.

 

dx

 

7x 14 y;

 

 

 

11.

dt

 

 

dy

 

 

 

 

3x 8 y.

 

 

 

 

dt

 

 

 

 

Вариант 9

1.y dy = – x dx.

2.(1 + x2)dy + xy dx = 0.

3. y

y2

2

y

.

x2

x

 

 

 

4.(2xy – ln y)dx + (x2 + 1 – xy )dy = 0.

5.y y e x .

70

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