Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
высшая математика теоретический материал часть 3.pdf
Скачиваний:
100
Добавлен:
26.07.2016
Размер:
1.43 Mб
Скачать

1.3. Вычисление двойных интегралов в полярных координатах

Для вычисления двойного интеграла I = ∫∫ f (x, y)dσ мы пользовались до сих пор

D

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

рис.15

x = r cos φ, у = r sin φ.

Разобьём область интегрирования D на частичные области σi двумя системами координатных линий: r = const, φ = const. Этими линиями будут соответственно концентрические окружности с центром в полюсе и лучи, исходящие из полюса (рис. 15). При этом частичными областями σi будут криволинейные четырехугольники, ограниченные дугами концентрических окружностей и их радиусами. Площадь σi области σi будет

 

σ i =

1

 

(ri +

2

ϕi

1

ri2

 

 

r

ri

ϕi

 

2

 

ri )

2

ϕi = ri +

 

i

 

 

 

 

 

 

 

 

 

 

2

 

 

 

Или

σ

i

= r'

r ϕ

i

где

 

 

 

 

 

 

 

 

 

 

 

 

i

i

 

 

 

 

 

 

 

 

 

r

' = r +

 

 

ri

есть средний радиус между r

, и r

+

r

.

 

 

 

i

i

 

 

2

 

 

 

 

 

 

 

i

i

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть дана функция f(x, у), непрерывная в области D. Составим для неё интегральную сумму, разбивая область D на частичные области σi, и выбирая в качестве произвольных точек Рi (xi , yi) областей σi точки, лежащие на средних окружностях радиуса ri т.е. полагая

x

= r'

cosϕ

i

y

= r' sinϕ

i

.

i

i

 

i

i

 

Тогда

 

 

 

 

 

 

 

n

 

 

 

n

 

 

In = f (xi , yi )

σ i = f (ri' cosϕi ,ri' sinϕi )ri' ri ϕi

 

i=1

 

 

 

i=1

 

 

Так как в правой части стоит интегральная сумма для функции f(r cos φ, r sinφ) r по переменным r и φ, то, переходя к пределу, получим

I = ∫∫ f (x, y)dxdy = ∫∫ f (r cosϕ , r sinϕ )rdrdϕ .

D D

Это равенство является формулой преобразования двойного интеграла от декартовых координат к полярным. Выражение dσ = r dr dϕ называется элементом площади в полярных координатах.

Правило преобразования двойного интеграла к полярным координатам.

Для того чтобы преобразовать двойной интеграл в декартовых координатах в двойной интеграл в полярных координатах, нужно x и у в подынтегральной функции заменить соответственно через r cos φ и r sin φ, а произведение dx dy заменить произведением r dr dφ.

49

Вычисление двойного интеграла в полярной системе координат, также как и в декартовой, сводится к последовательному интегрированию по переменным r и φ. Укажем правило расстановки пределов.

Пусть полюс не содержится внутри области интегрирования D, заключенной между лучами φ=φ1 и φ = φ2, и координатные линии φ = const встречают её границу не более чем в двух точках (рис. 16а). Область может также иметь вид, изображенный на рис. 16б.

Полярные уравнения кривых ADC и АВС пусть будут r = r1 (φ) и r = r2 (φ).

Интегрируя сначала по r в пределах его изменения при постоянном φ, т.е. от r1 (φ) до r2 (φ), а затем по φ от φ1 до φ2, получим

ϕ

r

(ϕ )

I = 2 dϕ

2

f (r cosϕ , r sinϕ )rdr .

ϕ1

r1

(ϕ )

Интегрирование в обратном порядке, т.е. сначала по φ, а потом по r, обычно не встречается.

Рис.16.

В частном случае, когда областью интегрирования служит часть кругового кольца r1 r r2 ,ϕ1 ϕ ϕ2 , пределы интегрирования постоянны по обеим переменным:

ϕ

r

 

I = 2 dϕ 2

f (r cosϕ , r sin ϕ )rdr.

ϕ1

r1

 

Пример 7. Расставим пределы интегрирования в полярных координатах, если область D – круг x2 + y2 ax (рис. 17). Переходя к полярным координатам, получим уравнение окружности в виде r = a cos φ. Здесь r1 (φ) = 0 и r2 (φ) = a cos φ. Пределы изменения по φ от π2

до π2 . Поэтому

π

 

a cosϕ

 

I = 2

 

 

dϕ

f (r cosϕ , r sinϕ )rdr .

 

π

 

0

 

2

 

 

 

Пример 8. Вычислим объём V общей части шара радиуса а и кругового цилиндра радиуса

а

2

 

 

 

при условии, что центр шара лежит на

Рис. 17

Рис. 18.

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

50

Имеем

1 V = ∫∫ a2 x2 y2 dxdy ,

 

4

D

где D – полукруг, являющийся половиной основания цилиндра.

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

соответствии с правилом преобразования имеем

1 V = ∫∫ a2 r2 rdrdϕ

 

4

D

Так как полярное уравнение полуокружности, ограничивающей область D, есть r = a cos φ (см. пример 7), то, интегрируя сначала по r, а затем по φ, найдём

 

 

 

 

 

π

 

 

 

 

a cosϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V =

2 dϕ

 

 

a2 r2 rdr .

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находя внутренний интеграл, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

3

 

 

a cosϕ

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

(a2

r2 )

 

 

 

 

a3 2

 

 

a3 π

 

2

 

 

 

 

 

2

3

 

 

 

 

V

=

 

 

 

 

 

 

 

 

 

dϕ =

(1 sin

 

ϕ )dϕ =

 

 

 

 

 

4

 

 

 

3

 

 

 

 

3

2

3

 

 

 

0

 

 

 

 

 

 

 

 

 

 

3 0

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

3

 

π

 

2

 

 

 

 

 

 

 

 

 

 

 

 

V =

 

 

a

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

3

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.4. Приложения двойных интегралов к задачам механики

Масса плоской пластинки переменной плотности.

Рассмотрим тонкую пластинку, расположенную на плоскости Оху и занимающую область D. Толщину этой пластинки считаем настолько малой, что изменением плотности по толщине её можно пренебречь.

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

Определённая таким образом поверхностная плотность будет зависеть только от положения данной точки, т.е. являться функцией её координат: = δ(x, y).

Рис.19

Если бы плотность была постоянной (δ = const), то масса всей пластинки равнялась бы Μ = δ S, где S – площадь пластинки. Найдём теперь массу неоднородной пластинки, считая, что её плотность является заданной функцией δ (x, y). Для этого разобьём область, занимаемую пластинкой, на частичные области σ1, σ2, …, σn с площадями σ1 σ2 , …, σn (рис. 19). Выбирая в каждой частичной области произвольную точку Pi (xi yi), будем считать, что плотность во всех точках частичной области постоянна и равна плотности δ(xi yi) в выбранной точке. Составим приближённое выражение для массы пластинки в виде интегральной суммы

n

M n = δ ( xi , yi ) σ i i =1

51

Для точного выражения массы следует найти предел суммы при условии n → ∞ и каждая

частичная область стягивается к точке. Тогда M = ∫∫ δ ( x , y ) d σ .

D

Статические моменты и центр тяжести пластинки.

Перейдём теперь к вычислению статических моментов рассматриваемой пластинки относительно осей координат. Для этого сосредоточим в точках Pi (xi, yi) массы соответствующих частичных областей и найдём статические моменты полученной системы материальных точек:

n

M x( n ) = yiδ ( xi , yi ) σ i , i =1

n

M y( n ) = xiδ ( xi , yi ) σ i . i =1

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

M x = ∫∫ yδ ( x, y)dδ , M у = ∫∫ xδ ( x, y )dδ

Находим координаты центра тяжести:

ξ

=

M y

=

∫∫ xδ ( x, y )dσ

 

M

 

 

∫∫ yδ ( x, y )d σ

D

η =

x

=

D

M

∫∫ δ ( x, y )dσ

M

 

∫∫ δ ( x, y )dσ

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

D

Если пластинка однородна, т.е. δ (x, у) = const, то формулы упрощаются :

ξ =

∫∫ xd σ

, η

=

∫∫ yd σ

,

D

D

 

 

 

 

 

S

 

 

S

 

где S – площадь пластинки.

Моменты инерции пластинки

Моментом инерции материальной точки Ρ с массой m относительно какой-либо оси называется произведение массы на квадрат расстояния точки Ρ от этой оси.

Метод составления выражений для моментов инерции пластинки относительно осей координат совершенно такой же, какой мы применяли для вычисления статических моментов. Приведём поэтому только окончательные результаты, считая, что δ(х, у) =1.

I x = ∫∫ y 2 dσ ,

I y = ∫∫ x 2 dσ (10)

D

D

Отметим еще, что интеграл ∫∫ xyd σ называется центробежным моментом инерции; он

D

обозначается Ιxy.

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

I0 = ∫∫(x2 + y 2 )dσ = I x + I y .

D

52