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

3447

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

 

 

 

 

 

(x)

 

exp

dx

 

eln x x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножив обе части уравнения на x , получим уравнение в

полных дифференциалах:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3x2

 

 

6x2 y 3xy2 )dx (2x3

3x2 y)dy 0.

 

 

Общий интеграл этого уравнения найдем по формуле (2.74),

положив x0

y0

0 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

(3x

2

6x

2

y 3xy

2

)dx C ,

x

3

2x

3

y

 

x

2

y

2

 

C.

 

 

 

 

 

 

2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общий

интеграл

 

получен.

Поскольку

 

1

 

 

1

 

0 ,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решения, не вытекающего из общего нет.

Решение

x

0 ,

входящее в состав общего интеграла (при C

0) ,

не является

посторонним, что показывает его подстановка в исходное уравнение.

Пример. Найти общий интеграл уравнения y(1 xy)dx xdy 0.

Решение. Это уравнение является уравнением Бернулли, но мы построим решение, приведя его к виду уравнения в полных дифференциалах.

P(x, y)

 

 

 

( y

xy2 ) 1

2xy;

 

Q(x, y)

1;

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

 

y

 

 

 

 

 

 

1

 

P(x, y)

 

Q(x, y)

1

(2 2xy)

(x).

 

Q(x, y)

 

 

y

 

x

 

x

 

 

 

 

 

 

 

 

Таким образом, уравнение не имеет интегрирующего множителя, зависящего только от x . Проверим выполнение условия (2.89):

1 Q(x, y)

P(x, y) x

P(x, y)

2(1

xy)

2

( y).

 

 

 

 

 

 

y

 

y(1

xy)

 

y

 

 

 

70

Уравнение имеет интегрирующий множитель, зависящий от y . Находим его по формуле (2.90):

 

 

 

( y)

exp(

( y)dy) exp

 

 

2

dy

 

 

e

 

2 ln y

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

y 2

 

 

Умножив, данное уравнение на 1 y2 , превращаем его в

уравнение в полных дифференциалах:

 

 

 

 

 

 

 

 

 

 

 

 

 

1

x dx

 

x

 

dy 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Его общий интеграл найдем по формуле (2.74), полагая

 

x0 0,

y0

1 . Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

1

 

 

 

y

 

 

 

 

 

x

 

x 2

 

 

 

 

 

 

 

 

 

x dx

0dy

 

C,

 

 

 

 

C.

 

 

 

 

 

 

 

 

y

 

 

 

 

y

2

 

 

 

 

 

 

 

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кроме общего интеграла получаем дополнительное

решение,

не

 

вытекающее

из

 

общего,

 

из

уравнения

1

y 2

0,

y

 

0.

Поскольку

 

 

 

1

 

0, то

 

в общем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интеграле не содержится посторонних решений.

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

Теорема. Если

(x, y) - интегрирующий множитель

уравнения (2.78), т.е.

(x, y)P(x, y)dx (x, y)Q(x, y)dy

dU(x, y) ,

то интегрирующим множителем является функция

(U) ,

где (z) - произвольная непрерывная функция, не равная нулю.

Доказательство.

Пусть

F(z)

- первообразная для

функции

(z) , т.е. F (z)

(z) . Тогда имеем:

(U ) Pdx Qdy

(U )dU F (U )dU dF(U ) .

Это

и доказывает,

что

(U)

- интегрирующий

множитель уравнения (2.78).

71

 

 

 

 

 

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

 

 

 

 

 

 

Найти интегрирующий множитель

и проинтегрировать

уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

(x2

 

y2

2x)dx

 

2ydy

 

0. Ответ:

ex ;

ex (x2 y2 )

C.

2.

 

dy

 

(3x2 cos y

sin y) cos y.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

 

1

 

 

; x(x2

tg y)

C;

y

 

 

k .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

2

 

 

3.

(x2

sin2 y)dx xsin 2ydy

0.

 

Ответ:

 

1

;

x2

sin2 y

Cx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

4.

1

 

x

dx

2xy

 

x

 

 

x 2

 

dy

0.

 

 

 

 

 

 

 

 

 

 

 

y

 

y

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

1

 

 

 

 

 

 

 

y 2

x

 

 

 

 

 

 

 

 

 

 

; ln

x

 

ln

y

 

 

C.

 

 

 

 

 

 

 

 

x

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.8*. Пфаффовы формы, полный дифференциал

итермодинамика

1.Дифференциальные соотношения, используемые в термодинамике, такие как уравнение первого закона термодинамики для простых систем

Q

dU

A,

(2.91)

где, Q - количества тепла подводимого (или отводимого) к

рассматриваемой термодинамической системе;

U - внутренняя

энергия системы; A -

работа,

производимая

системой (или

производимая над системой); соотношение для работы термодинамической системы,

связанной с увеличением

объема V системы

против сил

внешнего давления p

 

 

A

pdV;

(2.92)

уравнение второго закона термодинамики

 

Q

TdS,

(2.93)

72

где, S и T - энтропия и температура системы;

объединенное уравнение первого и второго законов термодинамики (для простой системы)

TdS dU pdV;

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

Z A1(x1, , xn )dx1 An (x1, , xn )dxn ,

где x1, , xn - независимые переменные.

Выражение вида

 

n

 

đ Z

Ak (x1, , xn )dxk

(2.94)

k

1

 

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

(Пфафф Иоганн Фридрих 1765-1825 г.г. – немецкий математик и астроном, один из учителей К. Гаусса. Наибольшие достижения Пфаффа относятся к области исследования теории важного класса дифференциальных уравнений). В формуле (2.94) перечеркнутая буква đ означает, что в правой части равенства стоит, вообще говоря, дифференциальное выражение, которое может быть и полным дифференциалом. Таким образом, введенный символ “đ” объединяет две возможности: đ d , когда правая часть этого равенства представляет собой полный дифференциал; đ - в том случае, когда правая часть не является полным дифференциалом.

Сравнивая, например, уравнение первого закона термодинамики для простых систем (2.91) с выражением (2.94) видим, что первый закон термодинамики представляет собой пфаффову форму двух переменных; для термодинамических систем, совершающих, помимо работы расширения, еще один вид работы - пфаффову форму трех переменных

 

Q dU pdV

W ,

где

- любая обобщенная сила, за исключением давления, W -

любая

обобщенная координата,

за исключением объема.

73

Выражения (2.92) и (2.93) очевидно являются пфаффовыми формами одного переменного.

2. При рассмотрении дифференциальных выражений вида

(2.94) важнейшим вопросом для теории дифференциальных уравнений является вопрос о том, будет ли это выражение полным дифференциалом. Этот вопрос в частном случае двух переменных мы изучали в §2.6 и §2.7.

Различие в поведении пфаффовых форм, представляющих собой полный дифференциал

 

n

 

 

 

 

 

 

n

Z (x1

, , xn )

 

 

 

dZ

A (x , , x

n

)dx

k

 

 

dx

k

,

 

 

 

 

 

 

k

1

 

 

 

 

 

 

xk

 

k

1

 

 

 

 

 

 

k 1

 

 

 

 

Ak (x1, , xn )

 

 

Z (x1, , xn )

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

xk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

от пфаффовых

форм,

не

являющихся

таковыми, весьма

существенна и для термодинамики. Дело в том, что вычисленный по любой линии, соединяющий точки M1 и M 2 , интеграл от полных дифференциалов

dZ Z (M 2 ) Z (M1)

(2.95)

M1M 2

равен разности значений функции Z в этих точках (не зависит

от вида линии). Если же

đ Z dZ , то

величина интеграла

зависит не только от координат точек M1

и M 2 , но и от вида

кривой, соединяющей эти точки:

 

đ Z

Z (M1M 2 ).

(2.95)

M1M 2

В термодинамике при отсутствии внешних воздействий на простую однокомпонентную систему ее состояние определяется заданием двух независимых интенсивных параметров (Напомним: термодинамические параметры, не зависящие от массы или числа частиц в системе, называются интенсивными). Любой другой параметр состояния такой

74

системы является функцией двух заданных параметров. Иными словами – любые три параметра состояния (например, давление p , удельный объем v V m или температура T )

однозначно связаны между собой. Уравнение, связывающее эти параметры, называется уравнением состояния вещества.

Например, уравнения состояния идеального газа –

уравнение Клапейрона

- Менделеева, имеет вид: pV vRT ,

где v m

- число молей газа; m - масса газа;

- молярная

масса;

R 8314,4

Дж (киломоль.K ) -

универсальная

(молярная) газовая постоянная.

Всистеме координат, по осям которой отложены величины этих параметров, связь между параметрами состояния представляется в виде поверхности состояния. Изменение состояние термодинамической системы приводит к перемещению по поверхности состояния изображающей состояние точки; кривые на поверхности состояния – траектории этой точки.

Втермодинамике разделяют все термодинамические функции на два вида – функции состояния и функции процесса.

Первые из них характеризуются тем, что при малом изменении состоянии системы изменение этой функции является полным дифференциалом, т.е. при конечном перемещении точки, изображающей состояние системы, имеет место формула (2.95). Для функций процесса – имеет место формула (2.96).

В§§2.6 и 2.7 мы рассмотрели два класса пфаффовых форм, зависящих от двух переменных: пфаффовы формы, которые не являются полным дифференциалом, но имеют интегрирующий множитель, умножение на который превращает дифференциальное выражение в полный дифференциал некоторой функции.

Например, второй закон термодинамики в формулировке Р. Клаузиуса (Рудольф Юлиус Эмануэль Клаузиус 1822 - 1888 г.г. – выдающийся немецкий физик) вводит понятие энтропии

соотношением dS

1

Q . Энтропия, в отличии от количества

T

 

 

75

тепла, является функцией состояния и потому

dS

- есть

полный дифференциал. Последнее означает, что

1

 

является

T

 

 

интегрирующим множителем для пфаффовой формы đ Q Q .

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

Голономными пфаффовыми формами называются такие линейные дифференциальные формы (2.94), которые имеют интегрирующий множитель (очевидно, что формы 1-го класса

– частный случай 2-го класса).

Пфаффовы формы третьего класса (не имеющие интегрирующего множителя) называются неголономными.

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

P(x, y, z)dx Q(x, y, z)dy R(x, y, x)dz

голономна в том случае, если выполняется соотношение:

P

Q

R

Q

R

Q

R

P

Q

0.

 

 

 

 

 

 

 

 

 

z x, y

y

x y, z

z

y x, z

x

 

x, z

x, y

y, z

Приведенный в этом параграфе материал указывает перспективу развития теории, изложенной в §§2.6 и 2.7 и связь формальной математической теории с физикой тепловых процессов. В третьем семестре курса “Высшая математика” мы покажем связь этой теории с теорией потенциала, имеющей важнейшее значение для гидро- и аэродинамики, термодинамики и теории электричества.

76

2.9. Применение дифференциальных уравнений первого порядка к задачам физики и экологии

З а д а ч а о р а д и о а к т и в н о м р а с п а д е. Экспериментальным путем установлено, что скорость распада радиоактивного вещества, т. е. скорость изменения его массы в зависимости от времени, прямо пропорциональна его количеству. Установим закон изменения массы m радиоактивного вещества в зависимости от времени t , считая,

что начальная масса вещества при

t 0

была

m0 .

Пусть в момент времени t

масса вещества есть m , в

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

t масса составляет m

m . За время t

распадается масса

m .

Отношение

 

 

m

средняя скорость

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

распада за время

t , a

lim

m

 

dm

 

мгновенная скорость

 

 

 

 

 

 

 

t

 

dt

 

 

 

 

t

0

 

 

 

 

распада в момент времени

t .

 

 

 

 

 

 

 

 

 

Согласно условию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dm

km,

 

 

 

 

 

 

(2.97)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где k

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

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

надо найти зависимость массы m

от времени t . Решая

уравнение, находим

dm

 

kdt, ln m

kt

ln C,

откуда

m

 

 

 

 

 

 

 

 

 

 

 

 

m

Ce k t .

 

 

 

(2.98)

Формула (2.98) дает зависимость массы вещества как

функции времени. В данной задаче

постоянная

C имеет

определенное

значение,

а именно:

при t

0

получаем

m Ce0

C.

Подставляя это значение

C в формулу (2.98),

0

 

 

 

 

 

 

 

 

 

получаем искомую зависимость массы радиоактивного вещества от времени:

77

 

 

m m e kt.

(2.99)

 

 

0

 

Равенство (2.98) представляет собой общее решение

дифференциального уравнения, а равенство (2.99)

частное

решение, отвечающее начальному условию данной задачи.

Коэффициент

k

определяется экспериментально.

Например, для радия k

0,000447. Промежуток времени T , за

который распадается половина первоначальной массы радиоактивного вещества, называют периодом полураспада этого вещества. Подставляя в формулу (2.99) вместо m

значение

m0 / 2 ,

вместо

k

значение 0,000447,

получаем

уравнение для определения периода полураспада T

радия:

 

 

m0

 

m e 0,000447 T ,

0,000447T

ln 2,

 

 

 

 

 

 

 

2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда T

 

 

 

ln 2

 

лет.

 

 

 

 

 

; T

1550

 

 

 

 

0,000447

 

 

З а д а ч а о з а к о н е «е с т е с т в е н н о г о р о с т а 1 ».

Закон «естественного роста»

это закон, согласно которому

скорость «роста» вещества прямо пропорциональна его количеству. Найдем формулу для определения изменения

количества вещества

y в зависимости от времени t , считая,

что в начальный момент t

0 количество вещества было равно

y0 .

 

 

Здесь независимой

переменной является время t , а

искомой величиной

количество вещества в любой момент

времени. Скорость «роста» вещества есть скорость изменения величины y в зависимости от переменной t . Используя, как и

в предыдущей задаче, физический смысл производной, можно записать закон «естественного роста» следующим образом:

dy

ky,

(2.100)

 

dt

 

 

_________________________________________________________________________

1) Примеры заимствованы из книги проф. МГУ Ризниченко Г.Ю. Лекции по математическим моделям в биологии. Ч.1. – Ижевск: НИЦ “Регулярная и хаотичная динамика”, 2002. 232 с.

78

где k 0 коэффициент пропорциональности. Уравнение (2.100), отличающееся только знаком правой части от уравнения (2.97), описывает многие процессы «размножения».

Решение уравнения (2.100), удовлетворяющее заданному

начальному условию t 0 , y

y0 ,

имеет вид

y

y ekt.

(2.101)

 

0

 

Формула (2.101) и выражает закон «естественного роста». Согласно этому закону, например, происходит «размножение» нейтронов в ядерных реакциях, размножение бактерий, рост кристаллов и т. п.

Отметим, что при линейной зависимости скорости роста (2.100) от численности популяции в данный момент времени, сама численность популяции испытывает экспоненциальный рост (рис. 7).

dy

 

y

dt

 

 

k tg

 

y0

 

0

 

y

0

t

 

 

a)

 

 

б)

Рис. 7 Зависимость скорости роста от численности

популяции (а) и численности от времени (б)

Рассматриваемая модель о динамике численности популяции была предложена Мальтусом в широко известной работе “О законе роста народонаселения”. Томас Роберт Мальтус (1766 – 1834 г.г., английский экономист и демограф,

79

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