Криволинейный интеграл

Вычисление криволинейного интеграла второго рода. Примеры.

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

Пусть кривая LL задана параметрическими уравнениями L:[x=x(t),y=y(t),z=z(t);t_0⩽t⩽t_kL:\left[{ \begin{array} { l } x=x(t), \\ y=y(t), \\ z=z(t); \\ \end{array} }\right.\quad t\_0 \leqslant t\leqslant t\_k, где x(t),;y(t),;z(t)x(t),;y(t),;z(t)- непрерывно дифференцируемые функции, и пусть точкам A_i(x_i,y_i,z_i)A\_i (x\_i ,y\_i ,z\_i ), которые задают разбиение кривой, соответствуют значения параметра t_it\_i, т.е. x_i=x(t_i),;y_i=y(t_i),;z_i=z(t_i)x\_i =x(t\_i ),;y\_i =y(t\_i ),;z\_i =z(t\_i ).

Тогда по теореме Лагранжа существуют такие точки tˉ_i;(i=1,2,…,n)\bar { t } \_i ;(i=1,2,\ldots ,n), что Δx_i=x_i−x_i−1=x′(tˉ_i)(t_i−t_i−1)\Delta x\_i =x\_i -x\_ { i-1 } = { x } '(\bar { t } \_i )(t\_i -t\_ { i-1 } ).

Выберем точки M_iM\_i, получающиеся при этих значениях параметра: M_i(x_i,y_i,z_i)=M_i(x(tˉ_i),y(tˉ_i),z(tˉ_i))M\_i (x\_i ,y\_i ,z\_i )=M\_i (x(\bar { t } \_i ),y(\bar { t } \_i ),z(\bar { t } \_i )).

Тогда интегральная сумма для криволинейного интеграла ∑_i=1nP(x_i,y_i,z_i)⋅Δx_i\sum\limits\_ { i=1 } ^n { P(x\_i ,y\_i ,z\_i )\cdot \Delta x\_i } будет равна интегральной сумме для определенного интеграла ∑_i=1nP(x(tˉ_i),y(tˉ_i),z(tˉ_i))⋅x′(tˉ_i)⋅Δt_i\sum\limits\_ { i=1 } ^n { P(x(\bar { t } \_i ),y(\bar { t } \_i ),z(\bar { t } \_i ))\cdot { x } '(\bar { t } \_i )\cdot \Delta t\_i }.

Так как Δl_i→0⇔Δt_i→0\Delta l\_i \to 0\Leftrightarrow \Delta t\_i \to 0, то, переходя к пределу при max⁡_i=1,2,…,nΔl_i→0,; или, ; что ; тоже ; самое ,;max⁡_i=1,2,…,nΔt_i→0\mathop { \max } \limits\_ { i=1,2,\ldots ,n } \Delta l\_i \to 0,;\text{ или, } ;\text{ что } ;\text{ тоже } ;\text{ самое } ,;\mathop { \max } \limits\_ { i=1,2,\ldots ,n } \Delta t\_i \to 0 в равенстве ∑_i=1nP(x_i,y_i,z_i)⋅Δx_i=∑_i=1nP(x(tˉ_i),y(tˉ_i),z(tˉ_i))⋅x′(tˉ_i)⋅Δt_i\sum\limits\_ { i=1 } ^n { P(x\_i ,y\_i ,z\_i )\cdot \Delta x\_i } =\sum\limits\_ { i=1 } ^n { P(x(\bar { t } \_i ),y(\bar { t } \_i ),z(\bar { t } \_i ))\cdot { x } '(\bar { t } \_i )\cdot \Delta t\_i }, получим

∫_LP(x,y,z)⋅dx=∫_t_0t_kP(x(t),y(t),z(t))⋅x′(t)⋅dt.\int\limits\_L { P(x,y,z)\cdot dx } =\int\limits\_ { t\_0 } ^ { t\_k } { P(x(t),y(t),z(t))\cdot { x } '(t)\cdot dt } .

Аналогично доказываются формулы для интегралов по другим координатам. Окончательно ∫_LP(x,y,z)⋅dx+Q(x,y,z)⋅dy+R(x,y,z)⋅dz==∫_t_0t_k[P(x(t),y(t),z(t))⋅x′(t)+Q(x(t),y(t),z(t))⋅y′(t)+R(x(t),y(t),z(t))⋅z′(t)]dt.\begin{array} { l } \int\limits\_L { P(x,y,z)\cdot dx+Q(x,y,z)\cdot dy+R(x,y,z)\cdot dz } = \\ =\int\limits\_ { t\_0 } ^ { t\_k } { \left[ { P(x(t),y(t),z(t))\cdot { x } '(t)+Q(x(t),y(t),z(t))\cdot { y } '(t)+R(x(t),y(t),z(t))\cdot { z } '(t) }\right]dt } . \\ \end{array}

Таким образом, вычисление криволинейного интеграла второго рода ни чем не отличается от вычисления интеграла первого и сводится к интегрированию по параметру. Направление интегрирования определяется условиями задачи.

Пример 1

Найти ∫_AB∪ydx−xdy+zdz\int\limits\_ { \mathop { AB } \limits^\cup } { ydx-xdy+zdz }, где AB∪\mathop { AB } \limits^\cup - виток винтовой линии  x = a ⋅cos⁡ t , y = a ⋅sin⁡ t , z = a ⋅ t ,0⩽ t ⩽\mathbf { \textit { x } } =\mathbf { \textit { a } } \cdot \cos \mathbf { \textit { t } } , \mathbf { \textit { y } } =\mathbf { \textit { a } } \cdot \sin \mathbf { \textit { t } } , \mathbf { \textit { z } } =\mathbf { \textit { a } } \cdot \mathbf { \textit { t } } , 0 \leqslant\mathbf { \textit { t } } \leqslant 2π\pi.

Решение:

∫_AB∪ydx−xdy+zdz=∫_02π(asin⁡t⋅(−asin⁡t)−acos⁡t⋅acos⁡t+at⋅a)dt=a2∫_02π(t−1)dt=2a2(π2−π)\int\limits\_ { \mathop { AB } \limits^\cup } { ydx-xdy+zdz } =\int\limits\_0^ { 2\pi } { (a\sin t\cdot (-a\sin t)-a\cos t\cdot a\cos t+at\cdot a)dt } =a^2\int\limits\_0^ { 2\pi } { (t-1)dt } =2a^2(\pi ^2-\pi )

Пусть плоская кривая задана в декартовой системе координат уравнением  y = y ( x )\mathbf { \textit { y } } =\mathbf { \textit { y } } (\mathbf { \textit { x } } ),  A ( a \mathbf { \textit { A } } (\mathbf { \textit { a } }, y ( a ))\mathbf { \textit { y } } (\mathbf { \textit { a } } )),  B ( b \mathbf { \textit { B } } (\mathbf { \textit { b } }, y ( b ))\mathbf { \textit { y } } (\mathbf { \textit { b } } )).

Тогда ∫_AB∪Pdx+Qdy=∫_ab(P(x,y(x))+Q(x,y(x))y′(x))dx.\int\limits\_ { \mathop { AB } \limits^\cup } { Pdx+Qdy } =\int\limits\_a^b { (P(x,y(x))+Q(x,y(x)) { y } '(x))dx } . > Пример 2

Найти ∫_AB∪(x+y)dx+y2dy\int\limits\_ { \mathop { AB } \limits^\cup } { (x+y)dx+y^2dy } по кривой [y=0,;0⩽x⩽1y=ln⁡x,;1⩽x⩽e\left[{ \begin{array} { l } y=0,;0\leqslant x\leqslant 1 \\ y=\ln x,;1\leqslant x\leqslant e \\ \end{array} }\right..

Решение:

Пользуясь свойством аддитивности, разобьем интеграл на сумму двух:

∫_01(x+o)dx+∫_1e(x+ln⁡x+ln⁡2x⋅1x)dx=x22∣_01+(x22+xln⁡x−x+ln⁡3x3)∣_1e=e22+43\int\limits\_0^1 { (x+o)dx+\int\limits\_1^e { (x+\ln x+ } } \ln ^2x\cdot \frac { 1 } { x } )dx=\left. { \frac { x^2 } { 2 } }\right|\_0^1 + \left. { (\frac { x^2 } { 2 } +x\ln x-x+ \frac { \ln ^3x } { 3 } ) }\right|\_1^e =\frac { e^2 } { 2 } +\frac { 4 } { 3 }

Пример 3

Найти ∮_cx(1−y)dx+xdy\oint\limits\_c { x(1-y)dx+xdy }, где  C \mathbf { \textit { C } } - окружность, проходимая в отрицательном направлении { по часовой стрелке } .

Решение:

Параметризуем окружность  x =2cos⁡ t , y =2sin⁡ t ,0⩽ t ⩽2π\mathbf { \textit { x } } =2\cos \mathbf { \textit { t } } , \mathbf { \textit { y } } =2\sin \mathbf { \textit { t } } , 0 \leqslant \mathbf { \textit { t } } \leqslant 2\pi. Движению по часовой стрелке соответствует уменьшение параметра  t \mathbf { \textit { t } }, поэтому интегрируем от 2π2\pi до 0:

∮_Сx(1−y)dx+xdy=∫_2π0(2cos⁡t(1−2sin⁡t)(−2sin⁡t)+2cos⁡t⋅2cos⁡t)dt=∫_2π0(−4cos⁡tsin⁡t++8cos⁡t⋅sin⁡2t+4cos⁡2t)dt=(−2sin⁡2t+83sin⁡3t+2t+sin⁡2t)∣_2π0=−4π.\begin{array} { l } \oint\limits\_С { x(1-y)dx+xdy=\int\limits\_ { 2\pi } ^0 { (2\cos t(1-2\sin t)(-2\sin t)+2\cos t\cdot 2\cos t)dt=\int\limits\_ { 2\pi } ^0 { (-4\cos t\sin t } } } + \\ +8\cos t\cdot \sin ^2t+4\cos ^2t)dt=\left. { (-2\sin ^2t+\frac { 8 } { 3 } \sin ^3t+2t+\sin 2t) }\right|\_ { 2\pi } ^0 =-4\pi . \\ \end{array}

Пример 4

Вычислить ∫_L2xydx+(x2+y2)dy\int\limits\_L { 2xydx+(x^2+y^2)dy } по каждому из путей, соединяющих точки  O (0,0)\mathbf { \textit { O } } (0,0) и  A (2,4)\mathbf { \textit { A } } (2,4), и изображённых на рисунке:

  1. Ломаная  OBA \mathbf { \textit { OBA } }, состоящая из прямолинейных отрезков;

  2. Ломаная  OCA \mathbf { \textit { OCA } };

  3. Прямолинейный отрезок  OA \mathbf { \textit { OA } };

  4. Параболы  ODA \mathbf { \textit { ODA } } и  OEA \mathbf { \textit { OEA } }, симметричные относительно координатных осей.

Решение:

1). ∫_OBA2xydx+(x2+y2)dy=∫_OB2xydx+(x2+y2)dy++∫_BA2xydx+(x2+y2)dy\int\limits\_ { OBA } { 2xydx+(x^2+y^2)dy } =\int\limits\_ { OB } { 2xydx+(x^2+y^2)dy } + \quad +\int\limits\_ { BA } { 2xydx+(x^2+y^2)dy } { по свойству аддитивности } . На  ОВ \mathbf { \textit { ОВ } } в качестве параметра естественно выбрать переменную  х \mathbf { \textit { х } }, при этом y=0,;dy=0y=0,;dy=0, поэтому ∫_OB2xydx+(x2+y2)dy=0\int\limits\_ { OB } { 2xydx+(x^2+y^2)dy } =0.

На  ВА x=2,;dx=0,;0⩽y⩽4\mathbf { \textit { ВА } } \quad x=2,;dx=0,;0\leqslant y\leqslant 4, поэтому ∫_BA2xydx+(x2+y2)dy=∫_04(22+y2)dy=(4y+y33)∣_04=3713\int\limits\_ { BA } { 2xydx+(x^2+y^2)dy } =\int\limits\_0^4 { \left( { 2^2+y^2 }\right)dy } =\left. { \left( { 4y+\frac { y^3 } { 3 } }\right) }\right|\_0^4 =37\frac { 1 } { 3 }.

Окончательно, I_1=∫_OBA2xydx+(x2+y2)dy=3713I\_1 =\int\limits\_ { OBA } { 2xydx+(x^2+y^2)dy } =37\frac { 1 } { 3 }.

2). ∫_OCA2xydx+(x2+y2)dy=∫_OC2xydx+(x2+y2)dy+∫_CA2xydx+(x2+y2)dy\int\limits\_ { OCA } { 2xydx+(x^2+y^2)dy } =\int\limits\_ { OC } { 2xydx+(x^2+y^2)dy } +\int\limits\_ { CA } { 2xydx+(x^2+y^2)dy }.

∫_OC2xydx+(x2+y2)dy=∣x=0,;dx=0,;0⩽y⩽4∣=∫_04y2dy=643=2113.∫_CA2xydx+(x2+y2)dy=∣y=4,;dy=0,;0⩽x⩽2∣=∫_028xdx=16.\int\limits\_ { OC } { 2xydx+(x^2+y^2)dy } =\left| { x=0,;dx=0,;0\leqslant y\leqslant 4 }\right|=\int\limits\_0^4 { y^2dy } =\frac { 64 } { 3 } =21\frac { 1 } { 3 } . \\ \int\limits\_ { CA } { 2xydx+(x^2+y^2)dy } =\left| { y=4,;dy=0,;0\leqslant x\leqslant 2 }\right|=\int\limits\_0^2 { 8xdx } =16. Окончательно, I_2=∫_OСA2xydx+(x2+y2)dy=3713I\_2 =\int\limits\_ { OСA } { 2xydx+(x^2+y^2)dy } =37\frac { 1 } { 3 }.

3). На прямолинейном отрезке  OA y=2x\mathbf { \textit { OA } } \quad y=2x, поэтому ∫_OA2xydx+(x2+y2)dy=∫_02(2x⋅2x⋅1+(x2+4x2)⋅2)dx=∫_0214x2dx=1123=3713\int\limits\_ { OA } { 2xydx+(x^2+y^2)dy } =\int\limits\_0^2 { \left( { 2x\cdot 2x\cdot 1+(x^2+4x^2)\cdot 2 }\right)dx } =\int\limits\_0^2 { 14x^2dx } =\frac { 112 } { 3 } =37\frac { 1 } { 3 }.

4). Уравнение параболы  OEA \mathbf { \textit { OEA } } имеет вид x=ky2x=ky^2, значение коэффициента  k \mathbf { \textit { k } } найдём, подставляя в это уравнение координаты точки  А \mathbf { \textit { А } }:

2=k42→k=1/8→x=y2/8; x′_y=y/42=k4^2\rightarrow k=1/8\rightarrow x=y^2/8;~x'\_y=y/4, поэтому

∫_OEA2xydx+(x2+y2)dy=∫_04[(2⋅y28⋅y)⋅y4+(y464+y2)]dy=∫_04(5y464+y2)dy=3713\int\limits\_ { OEA } { 2xydx+(x^2+y^2)dy } =\int\limits\_0^4 { \left[ { \left( { 2\cdot \frac { y^2 } { 8 } \cdot y }\right)\cdot \frac { y } { 4 } +\left( { \frac { y^4 } { 64 } +y^2 }\right) }\right]dy } =\int\limits\_0^4 { \left( { \frac { 5y^4 } { 64 } +y^2 }\right)dy } =37\frac { 1 } { 3 }.

5). Совершенно также убеждаемся, что интеграл по параболе  ODA \mathbf { \textit { ODA } } имеет то же значение.

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

Убедимся в том, что это не так, на примере интеграла ∫_L2xydx−(x2+y2)dy\int\limits\_L { 2xydx-(x^2+y^2)dy }: I_1=∫_OBA2xydx−(x2+y2)dy=−∫_04(22+y2)dy=−3713I\_1 =\int\limits\_ { OBA } { 2xydx-(x^2+y^2)dy } =-\int\limits\_0^4 { (2^2+y^2)dy } =-37\frac { 1 } { 3 }; I_2=∫_OCA2xydx+(x2−y2)dy=∫_04y2dy−∫_028xdx=643−16=513≠I_1I\_2 =\int\limits\_ { OCA } { 2xydx+(x^2-y^2)dy } =\int\limits\_0^4 { y^2dy } -\int\limits\_0^2 { 8xdx } =\frac { 64 } { 3 } -16=5\frac { 1 } { 3 } \ne I\_1.

Следующие разделы будут посвящены ответу на поставленный вопрос: при каких условиях криволинейный интеграл второго рода ∫_СPdx+Qdy\int\limits\_С { Pdx+Qdy } не зависит от формы пути, соединяющего начальную и конечную точки, и определяется только этими точками? В трёхмерном случае этот вопрос будет изучаться в теории поля.

Другие статьи из категории

Похожие статьи

Все статьи →

Криволинейный интеграл

Свойства криволинейного интеграла второго рода

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

Криволинейный интеграл

Формула Грина

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

Криволинейный интеграл

Условия независимости криволинейного интеграла от пути интегрирования

Независимость криволинейного интеграла от пути интегрирования является ключевым свойством потенциальных полей и имеет важное значение в…

Криволинейный интеграл

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

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

Криволинейный интеграл

Выражение площади плоской области через криволинейный интеграл

Криволинейный интеграл предоставляет мощный инструмент для вычисления площади плоской области без необходимости непосредственного…

Нужен расчёт?

Пришлите исходные данные — уточним состав работы и ответим по срокам.

Связаться