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

Учебное пособие 1121

.pdf
Скачиваний:
4
Добавлен:
30.04.2022
Размер:
791.65 Кб
Скачать

Для этого функцию f (x) продолжают с отрезка 0; l на отрезок l ; l и продолженную функцию разлагают в ряд Фурье на отрезке l ; l стандартным способом.

Ряд Фурье для функции f (x) называют ещё разложением f (x) в ряд Фу-

рье на отрезке 0; l .

Рис. 8. Продолжение функции на отрезок l; l

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

Рис. 9. Продолжение функции чётным и нечётным способами

Аналитически продолжение задаётся следующим образом. Если черезf (x) обозначить продолженную функцию, то чётное продолжение будет задаваться соотношением

f (x) , 0 x l ,

f (x)

f ( x), l x 0;

под нечётным продолжением понимается функция f (x) , определяемая соотношением

71

f (x) , 0 x l ,

f (x)

f ( x), l x 0.

Заметим, что если f (x) есть чётное продолжение исходной функции, то коэффициенты Фурье для f (x) могут быть найдены по формулам (2.36), (2.37). Если же f (x) есть нечётное продолжение функции f (x) , то для на-

хождения коэффициентов Фурье функции f (x) можно воспользоваться

формулой (2.42). Во всех этих формулах фигурирует только исходная функция f (x) .

П р и м е р Разложить в ряд Фурье по косинусам на отрезке [0; ] функцию

f (x)

 

 

x

.

4

2

 

 

 

В рассматриваемом случае l . Продолжим данную функцию на отрезок

[ ;0] чётным образом и обозначим продолженную функцию через

f (x) .

Функция

f (x) чётная, поэтому её ряд Фурье состоит только из косинусов.

Найдём коэффициенты an , n 0,1, 2,... по формулам (2.36), (2.37).

 

 

 

 

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

x

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

2

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

 

 

f (x) dx

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

4

2

 

 

 

4

4

 

4

 

4

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

an

 

 

f (x) cos nx dx

 

 

 

 

 

 

 

 

 

 

 

 

cos nx dx

 

 

 

cos nx dx

 

 

 

x cos nx dx

 

 

 

 

4

2

 

2

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

x sin nx

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin nx

 

 

 

 

x d (sin nx )

 

 

 

 

 

 

 

 

 

 

sin nx dx

 

cos nx2

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

2n

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, n 2k ,

 

 

 

 

 

(cos n cos0)

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(( 1)

 

1)

 

 

 

(1 ( 1)

 

)

2

 

 

 

 

 

 

 

 

n

2

 

n

2

 

n

2

 

 

 

,

n 2k 1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

k N .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, a0 0, a2k

0, a2k 1

 

 

 

 

 

 

 

 

2

 

 

 

 

, k N

. Поэтому

 

 

 

 

 

(2k 1)2

 

 

 

 

 

72

 

2

 

cos(2k 1)x

 

f (x)

 

 

(2k 1)

2

.

 

 

k 1

 

 

Заметим, что функция f (x)

является

непрерывной и кусочно-

дифференцируемой на отрезке [ ; ] , удовлетворяющей условию f ( )f ( ). Поэтому

 

 

 

 

2

 

 

 

 

cos(2k 1)x

 

 

 

 

 

f (x)

 

 

 

 

 

 

(2k 1)

2

 

 

,

x .

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

Поскольку f (x) f (x) , 0

x , то получаем требуемое разложение

 

 

x

 

 

2

 

 

 

 

cos(2k 1)x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2k 1)

2

 

, 0 x .

4

2

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

Положив в последней формуле x 0, получим равенство

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

k 1

 

(2k 1)

 

 

 

 

из которого следует:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

1

 

...

2

 

 

 

 

 

 

 

1

 

 

 

.

 

(2k 1)

2

2

 

2

 

 

k 1

 

 

 

 

 

3

 

 

 

5

8

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

1.Дайте определения функционального ряда и области сходимости такого ряда.

2.Дайте определение степенного ряда. Сформулируйте теорему Абеля и

следствие из неё. Будет ли сходиться степенной ряд вида an xn в точке

n 0

x 5 , если он сходится в точке x 7 ?

3.Дайте определения радиуса сходимости и интервала сходимости степенного ряда.

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

5. Дайте определение разложимости функция f (x) на интервале ( R; R) , R 0 в степенной ряд.

73

6.Сформулируйте свойства сумм степенных рядов.

7.Дайте определение ряда Маклорена функции f (x) .

8.Сформулируйте необходимое и достаточное условие разложимости функции f (x) в ряд Маклорена.

9.Выпишите ряды Маклорена для функций ex , sin x , cos x , (1 x) , 0,N, ln(1 x), arctg x и для каждой из этих функций укажите промежуток, на

котором она является суммой ряда.

10.Дайте определение ряда Фурье функции f (x) , интегрируемой на отрез-

ке[ l ;l] .

11.Дайте определение функции f (x) , кусочно-дифференцируемой на отрезке [a ;b]. Сформулируйте теорему о сходимости рядов Фурье.

12.Какой вид имеет ряд Фурье чётной интегрируемой функции? Нечётной функции? Как преобразуются формулы для коэффициентов Фурье чётной (нечётной) функции?

§3. Дифференциальные уравнения в частных производных

3.1.Основные понятия

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

двумя независимыми переменными x, y называется уравнение вида

 

u ,

u

 

 

2

u

 

 

2

u

 

 

2

u

 

 

F x, y, u(x , y),

,

 

,

 

,

 

, ... 0,

(3.1)

y

 

 

2

 

 

2

x y

 

x

 

x

y

 

 

которое связывает неизвестную функцию u(x , y) , её аргументы и частные производные (по x и по y ). Здесь F заданная функция своих аргументов.

Предполагается, что равенство (3.1) выполняется тождественно в некоторой области E2 .

Аналогично определяется уравнение и для большего числа независимых переменных.

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

тельно непрерывны в области и которая при подстановке её в уравнение обращает его в тождество.

74

П р и м е р Решить дифференциальное уравнение в частных производных второго по-

рядка 2u 0.

x y

Для решения уравнения перепишем его в виде

u 0.y x

Обозначив ux v , получим

v 0.

y

Из последнего уравнения следует, что v (x), или ux (x), где(x) произвольная функция аргумента x. Поэтому

u(x, y) (x)dx ( y) (x) ( y) ,

где ( y) - произвольная функция переменной y . Таким образом, u(x, y) (x) ( y).

Отметим, что здесь (x) и ( y) - произвольные функции своих перемен-

ных, а полученное решение зависит от двух произвольных функций. Легко видеть, что если (x) и ( y) есть произвольные дважды непрерывно диф-

ференцируемые функции своих аргументов, то функция u(x, y) действительно является решением рассматриваемого уравнения.

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

Приведём некоторые примеры дифференциальных уравнений в частных производных, которые наиболее часто встречаются при описании физиче-

ских явлений и которые называются уравнениями математической физики.

Уравнения математической физики подразделяют на уравнения гиперболического, параболического и эллиптического типов (см. [3, с. 19]).

1. Простейшее уравнение гиперболического типа имеет вид

2u a2 2u , u u(x,t).

t2 x2

75

Здесь x пространственная переменная, t время. Это дифференциальное уравнение второго порядка называется одномерным волновым уравнением. Оно описывает колебательные процессы при поперечных колебаниях упругой струны, продольных колебаниях упругого стержня и т.д.

Рассматривают также уравнения вида

2u

a

2

 

2u

 

2u

, u u(x, y,t).

t

2

 

 

x

2

y

2

 

 

 

 

 

 

 

 

 

 

Это так называемые двумерные волновые уравнения, часто используемые в акустике.

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

2. Простейшее уравнение параболического типа имеет вид

u a2 2u , u u(x,t).

t x2

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

3. Простейшее уравнение эллиптического типа имеет вид

2u 2u 0, u u(x, y).

x2 y2

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

3.3. Вывод дифференциального уравнения колебаний струны

Рис. 10. Профиль колеблющейся упругой струны длины l в момент времени t

76

Рассмотрим натянутую упругую струну, закреплённую на концах (в точках x 0 и x l). Пусть в положении равновесия струна направлена по оси

Ox. Если её отклонить от положения равновесия, то она начнёт совершать

колебания.

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

жению (оси Ox ).

Обозначим через u x,t смещение точек струны в момент времени t от

положения равновесия,

0 x l . При

каждом фиксированном

значении

t график функции u x,t будет давать,

очевидно, форму струны в этот мо-

мент времени (рис. 10).

 

 

 

Обозначим через

угол, который

составляет касательная

к графику

функции u x,t с осью Ox . При этом x,t , т.е. это также будет функ-

ция координаты и времени.

 

 

 

 

Отметим, что

u x,t

 

 

tg

;

 

 

 

 

x

 

 

u

v x,t скорость точек струны;

2u

w x,t ускорение точек стру-

t

 

t2

 

 

ны.

 

 

 

Будем изучать малые колебания струны, т.е. будем считать, что угол мал, так что его квадратом и более высокими степенями можно пренебречь,

2 1.

Всилу сказанного будем считать, что

sin

3

 

5

 

... ; cos 1

2

 

4

... 1;

3!

5!

 

2!

4!

 

 

 

 

 

 

 

 

 

 

 

 

tg

 

sin

 

 

 

 

sin

sin , т.е. tg ;

cos

 

 

1 sin 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

2

 

u 2

 

 

 

 

 

tg

 

0 tg

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

Запишем баланс сил для любого участка струны вдоль осей координат.

Для этого сначала найдём изменение длины участка M1M2 струны при её колебаниях (рис. 11)

77

 

 

x2

 

 

u

2

x2

M1M2

 

 

1

 

dx dx x2 x1 .

 

 

 

 

 

 

 

 

x1

 

 

x

 

x1

 

 

 

 

 

 

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

любой её участок сохраняет свою длину.

Пусть T T x,t

величина силы натяжения (x,t) , действующей в

 

 

 

 

каждой точке струны. Тогда уравнение равновесия для участка

M1M 2

по оси

Ox имеет вид

T2 cos 2 T1 cos 1 0 T1 T2 T .

Таким образом, величина силы натяжения при малых колебаниях не изменяется в процессе колебаний.

Рис. 11. К выводу уравнений равновесия при колебаниях струны

Пусть const

стка струны M1M 2

или

линейная плотность струны. Тогда баланс сил для учапо оси Oy будет иметь вид

T

sin

2

T sin

1

x

2u

,

2

 

 

 

1

 

 

t

2

 

 

 

 

 

 

 

 

 

x2 x1

 

 

T tg

 

T tg x

2u

.

 

(3.2)

2

 

 

 

 

 

 

1

 

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим величину tg 2 tg 1 :

78

tg 2 tg 1 u x2 ,t u x1,t 2u x .

x x x2

Поэтому соотношение (3.2) можно записать в виде

T 2u x x 2u .

x2 t2

Разделив полученное равенство на x , получим соотношение

T2u 2u .

x2 t2

Вводя обозначение a

T

, окончательно получаем уравнение колебаний

 

 

 

 

 

 

 

упругой струны в виде

 

 

 

 

 

 

 

 

2u

a2 2u

, a2

T .

(3.3)

 

 

t2

x2

 

 

 

Уравнение (3.3) является уравнением свободных колебаний упругой струны. Если струна колеблется в какой-то среде или на неё действуют некоторые силы, например, вес, электромагнитные силы и т.п., то уравнение примет вид

2u a2 2u g x,t ,

t2 x2

где g x,t плотность внешних сил. Такое неоднородное уравнение назы-

вают уравнением вынужденных колебаний струны. В дальнейшем мы будем изучать только свободные колебания.

Уравнение колебаний струны описывает все процессы, происходящие при колебаниях, связывает колебания струн различной длины в различных условиях. В дальнейшем будем рассматривать только колебания ограниченной струны длиной l , закреплённой на концах.

Чтобы из множества всех решений уравнения с частными производными выбрать определённое решение, необходимо задать дополнительные условия, которым должно удовлетворять решение уравнения. Укажем такие условия для уравнения колебаний струны.

79

Будем предполагать, что в начальный момент времени t 0 заданы начальные смещения и скорости всех точек струны, т.е.

u x, t

 

t 0 u x, 0 f x ,

u x,t

 

t 0

u x,0

F x , 0 x l ,

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

t

 

 

t

 

 

 

 

 

 

u x,0

 

 

 

 

 

 

 

 

 

u x, 0 f x ,

 

 

F x , 0 x l ,

(3.4)

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

где f x ,

F x

известные функции,

заданные на отрезке 0;l . Сформу-

лированные условия называются начальными условиями для уравнения колебаний струны.

Условие закрепления струны на концах приводит к следующим соотно-

шениям:

 

u 0, t 0 , u l, t 0 , t 0 ,

(3.5)

которые означают, что перемещения в концевых точках струны отсутствуют (равны нулю). Последние два условия – это так называемые краевые (граничные) условия. Заметим, что функции f (x) и F(x) , определённые на от-

резке [0;l], из-за условий закрепления струны на концах должны удовлетво-

рять условиям f (0) f (l) 0, F(0) F(l) 0 .

В общем случае концевые точки струны могут двигаться, и тогда гранич-

ные условия принимают вид

 

u 0, t 1 t , u l, t 2 t , t 0 ,

(3.6)

где 1 и 2 заданные функции времени.

Итак, решение задачи о свободных колебаниях упругой струны с закреплёнными концами сводится к решению уравнения (3.3) при начальных условиях (3.4) и граничных условиях (3.5).

3.4.Метод Фурье решения задачи о колебаниях ограниченной струны с закреплёнными концами

Найдём сначала решения уравнения (3.3), не равные тождественно нулю, удовлетворяющие однородным (т.е. нулевым) граничным условиям (3.5)

u(0,t) u(l,t) 0

80