МАРТ_Часть1_Глава6
.pdfГлава 6. СПЕЦИАЛЬНЫЕ ФУНКЦИИ
6.1. Краткие сведения по теории функций комплексной перемен-
ной
Решение многих задач радиотехники, связанных, например, со спектральным анализом сигналов, изучением распространения электромагнитных волн в различных средах, исследованием пространственно-временной структуры электромагнитного поля в резонаторах, волноводах и излучающих системах, синтезом фильтров, этот список можно продолжить, приводит к использованию различных специальных функций. Иногда вместо термина специальные функции используют термин высшие трансцендентные функции, оставляя название элементарные функции за степенными, показательными, логарифмическими и тригонометрическими. Обычно в класс специальных включают гамма- и бета-функции, интегральную показательную функцию и родственные ей, интеграл ошибок и связанные с ним функции, дзета-функцию Римана, эллиптические интегралы, классические ортогональные полиномы и функции Бесселя. Для знакомства со специальными функциями нам потребуется освежить в памяти читателя основные понятия теории функций комплексного переменного. Более подробную информацию можно получить в многочисленных руководствах по данному разделу математики [15].
Говорят, что на множестве M точек плоскости z x jy задана функция w f z , если известно правило, по которому каждой точке z из M ставится в соответствие определённая точка или совокупность точек w u jv . В первом случае функция f (z) называется однозначной, а во втором – многозначной.
Обычно под множеством M понимают область D, для которой выполняются следующие условия:
–для каждой точки из D этому множеству принадлежит и достаточно малый круг с центром в этой точке (свойство открытости);
–любые две точки D можно соединить ломаной, все точки которой принадлежат D (свойство связанности).
Граничной точкой области называют такую точку, которая сама не принадлежит D, но в любом сколь угодно малом круге с центром в этой точке содержатся точки из D. Совокупность граничных точек называют грани-
96
цей области. Область D с присоединённой к ней границей называют замкнутой и обозначают D .
Как и для операторов и обычных функций множество M называется множеством, а в случае области – областью определения, а совокупность всех значений, которые f(z) принимает на M, – множеством (областью) значений.
Как видно из введённых обозначений z x jy и w u jv , задание функций комплексного переменного равносильно заданию функций двух действительных переменных x и y, а именно u u x, y и v v x, y . Например,
w u jv z 2 (x jy)2 x2 y 2 j2xy , |
|
|||
и по |
условию |
равенства |
комплексных чисел u x, y x2 y 2 , |
а |
v x, y 2xy . |
Для |
функции |
f z e z e x jy e x cos y j sin y |
и |
u(x, y) e x cos y , а v x, y e x sin y .
Важнейшими понятиями для функций комплексного переменного являются понятия предела, непрерывности и аналитической функции.
Функция f(z) имеет предел в точке z0 x0 jy0 , если она однозначна и определена в окрестности этой точки и существуют пределы
lim u(x, y) u0 и |
lim |
v(x, y) v0 , |
|
|
x x0 |
x x0 |
|
|
|
y y0 |
y y0 |
f z u0 jv0 |
|
|
что обозначают как |
lim |
w0 . |
||
|
z z0 |
|
|
|
Другое определение звучит так: условие |
lim f (z) w0 выполняется |
|||
|
|
|
|
z z0 |
тогда и только тогда, когда для 0 |
0, так что из неравенства |
|||
z z0 вытекает f (z) w0 .
Важно подчеркнуть, что если предел существует, то z может стремиться к
z0 по любому пути. Функция f (z) называется непрерывной в точке z0, если она определена в окрестности этой точки, включая саму эту точку z0, и
lim f (z) f z0 . Необходимым и достаточным условием непрерывности
z z0
f (z) в точке z0 является непрерывность функций u(x, y) и v(x, y) в точке (x0, y0). Функция называется непрерывной в области D, если она непрерывна в каждой точке этой области.
Для функции f (z) , непрерывной в замкнутой области D , справедливы следующие утверждения:
97
–f (z) ограничена в D , т. е. существует такая постоянная M, что для всех z из D : f (z) M ;
–f (z) достигает в D своего наибольшего и наименьшего значений,
т. е. z1, z2 |
|
, что для z |
|
|
|
|
f z1 |
|
|
|
|
f z |
|
|
|
|
|
f z2 |
|
|
|
|
f z |
|
; |
|||||||
D |
D |
|
|
|
|
и |
|
|
|
|
||||||||||||||||||||||
|
– f (z) равномерно непрерывна, |
т. |
е. |
для |
0 |
|
0 , зависящее |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
только от ε, так что для любой пары точек z1 |
и z2 из |
D из неравенства |
||||||||||||||||||||||||||||||
|
z1 z2 |
|
следует неравенство |
|
|
f z1 f z2 |
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
Условия дифференцируемости функции |
f z u x,y jv x,y дают так |
||||||||||||||||||||||||||||||
называемые условия Даламбера–Эйлера или Коши–Римана: функция f (z) , определённая в некоторой окрестности точки z и имеющая в этой окрестности дифференцируемые вещественную u(x, y) и мнимую v(x, y) части, является дифференцируемой, если u(x, y) и v(x, y) удовлетворяют условиям:
u x, y |
v x, y , |
u x, y |
v x, y . |
(6.1) |
x |
y |
y |
x |
|
Функция f (z) , дифференцируемая в каждой точке некоторой области, |
||||
называется аналитической в данной области. Например, функция |
f z e z |
|||
аналитична на всей плоскости комплексного переменного z. Её однозначность очевидна. Проверим выполнение условия (6.1):
e z e x jy e x cos y jsin y , т. е. u x, y ex cos y , v x, y ex sin y .
Тогда |
|
|
u x, y |
|
|
u x, y |
ex cos y ; |
ex sin y ; |
|
||
x |
|
|
y |
|
|
v x, y ex sin y ; |
v x, y ex cos y . |
(6.2) |
|||
x |
|
|
y |
|
|
Im z = y
z1
C3
C2
z0 C1
0 |
Re z = x |
Рис. 6.1
Из (6.2) видно, что условия Даламбера– Эйлера выполняются.
Интеграл от функции комплексного переменного задаётся как контурный интеграл вдоль некоторой ориентированной кривой C, на которой f (z) естественно полагается за-
данной: I f d , где буквой обозначена
C
комплексная переменная интегрирования.
98
Большой интерес представляют интегралы от регулярных или аналитических функций. Если функция f(z) аналитична в области, в которой лежит контур C, то интеграл f d зависит не от пути интегрирования, а лишь
C
от начальной z0 и конечной z1 точек (рис. 6.1). При этом естественно предполагается, что все рассматриваемые пути (С1, С2, С3, …) не выходят за пределы области D , в которой функция f (z) аналитична. Следствием этого факта является теорема Коши. Если f (z) регулярна в замкнутой области D с границей C, то интеграл по замкнутому контуру C равен нулю:
f ( ) d 0 .
C
Мощным аналитическим аппаратом является представление функций комплексного переменного с помощью степенных рядов. Справедлива следующая теорема.
Всякая функция f (z) , аналитическая в кольце K, r z a R , может быть представлена в этом кольце своим рядом Лорана:
|
f z |
|
|
|
z |
a n . |
|
|||
|
|
Cn |
(6.3) |
|||||||
|
|
|
|
|
|
n |
|
|
|
|
Коэффициенты Cn вычисляются по формуле |
|
|||||||||
Cn |
1 |
|
|
|
f d |
|
, |
n = 0, ±1, … , |
(6.4) |
|
2 j |
|
a n 1 |
||||||||
|
C |
ρ |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
где C – любой замкнутый контур, лежащий в кольце K и обходящий
точку
z = a один раз в положительном направлении (против часовой стрелки). Структура ряда Лорана зависит от поведения функции в окрестности точки z = a.
Если f (z) аналитична в точке z = a, т. е. аналитична в круге z a R ,
то ряд Лорана превращается в ряд Тейлора, для которого Cn= 0 |
при n = –1, |
|||||||||
–2, … . Если ряд (6.3) имеет вид |
|
|
|
|
|
|
|
|||
f z |
C m |
|
C m 1 |
... |
C 1 |
C |
0 |
C z a ... C |
n |
z a n ... |
|
|
|
||||||||
|
z a m |
|
z a m 1 |
|
z a |
|
1 |
|
||
|
|
|
|
|
|
|
|
|||
,
то говорят, что функция f (z) имеет в точке z a полюс порядка m.
99
Слагаемые |
C m |
... |
C 1 |
образуют главную часть ряда Лорана. |
|
z a m |
z a |
||||
|
|
|
Если главная часть ряда Лорана содержит бесконечное число членов, т. е. имеет вид
C
n n ,
n 1 z a
то точка z = a является существенно особой точкой функции f (z) . Среди коэффициентов ряда Лорана наибольший интерес представляет
C–1, называемый вычетом функции f (z) в точке z = a. Как видно из (6.4),
С 1 |
1 |
f d , |
|
|
|||
2 j |
|||
|
C |
||
|
|
т. е. интеграл по любому замкнутому контуру, охватывающему точку z = a и лежащему в кольце регулярности K равен 2 j C 1. Обобщением этого факта является основная теорема о вычетах.
Пусть C0 – простой замкнутый контур (это означает, что все точки контура проходятся один раз), внутри которого f (z) аналитична всюду за исключением конечного числа изолированных* особых точек а1, а2, а3, …, ап. Тогда
|
|
|
|
n |
|
|
z , ak , |
|
|
|
f d 2 j Re s f |
|
|
||||||
|
С0 |
|
k 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где Re s f z , ak – вычет функции f (z) |
в точке z = ak. |
|
|||||||
Рассмотрим вопрос об определении вычетов. Пусть |
f (z) |
имеет в точке |
|||||||
z = a полюс первого порядка, т. е. |
|
|
|
|
|
||||
f z |
C 1 |
C |
0 |
C z a ... C |
n |
z a n |
.... |
(6.5) |
|
|
|||||||||
|
z a |
|
1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
Умножая правую и левую части (6.5) на (z – a) и переходя к пределу
при z a , получим: Re s f z , a C 1 lim z a f z .
z a
Пусть f (z) имеет в точке z = a полюс кратности m:
* Напомним определение изолированной особой точки. Если f (z) регулярна при
0 z a R и не является регулярной в точке a, то точка a называется изолированной
особой точкой.
100
f z |
C m |
... |
C 1 |
C |
C z a .... |
|
|
||||
|
z a m |
|
z a |
0 |
1 |
|
|
|
|
Умножая правую и левую части на z a m , дифференцируя (m – 1) раз
по z и переходя затем к пределу при z 0, получим: |
|
|
|
|
||||||
|
|
1 |
|
|
m 1 |
|
m |
|||
Re s f z , a C 1 |
|
|
d |
|
f z z a |
|
||||
|
|
lim |
|
|
|
|
|
. |
||
m 1 ! |
|
|
|
m 1 |
|
|||||
|
|
z a |
|
d z |
|
|
||||
|
|
|
|
|
|
|
|
|
||
Вычисление большого числа несобственных интегралов базируется на лемме Жордана.
Если f z e j z F z , 0 и F z 0 при z по любому пути в верхней полуплоскости Im z 0 , и на действительной оси функция f (z) аналитична, а в верхней полуплоскости имеет конечное число изолированных особых точек а1, а2, а3, …, ап, то
|
n |
|
f x dx 2 j Re s f z , ak . |
|
k 1 |
Лемма Жордана может быть обобщена и на случай, когда функция
f (z) имеет на действительной оси конечное число полюсов первого порядка b1, b2, …, bm. В этом случае выполняется равенство
|
|
n |
|
|
|
1 |
m |
|
|
|
|
|
|
|
|
k |
|
|
k |
|
|
|
|
|
|
|
|
|||||
|
f x dx 2 j |
|
|
Re s f z , a |
|
|
|
|
Re s f z , b |
. |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
2 |
k 1 |
|
|
|
k 1 |
|
|
|
|
|||||
Рассмотрим применение леммы Жордана для вычисления интеграла
sin t dt , с которым мы столкнулись при изучении оператора Гильберта.
t
Яcно, что от ω зависит лишь знак интеграла, а не его величина, так как, умножив и разделив на ω и сделав замену переменной, мы получим интеграл
sin x dx . В то же время замена ω на –ω в силу нечётности синуса меняет
x
знак интеграла на противоположный. Применим лемму Жордана к функции e jz
z ( 1, f (z) 1 z , и условия леммы Жордана выполняются).
Интеграл
101
|
e |
jx |
|
|
|
cos x |
|
|
|
sin x |
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
1 |
m |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
k |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
dx |
|
|
|
|
|
|
dx j |
|
|
|
|
|
dx 2 j |
|
|
|
|
|
Re s f |
z , a |
|
|
|
|
|
|
Re s f z , b |
. |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
x |
|
|
x |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
k 1 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k 1 |
|
|
|
|
|
|
|
|||||||||||||||
Функция e jz z |
имеет единственную особую точку – полюс первого |
|
|
|||||||||||||||||||||||||||||||||||||||
порядка в точке z = 0 и вычет b |
lim |
z |
e jz |
|
1. Таким образом, |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
z 0 |
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
cos x |
dx j |
|
sin x |
dx j |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
и по условию равенства комплексных чисел |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
sin x |
|
|
|
|
|
|
|
sin t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
dx , а |
|
|
|
|
|
|
|
|
dt sign . |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
x |
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Завершим наш краткий экскурс в теорию функций комплексного пере- |
|
|||||||||||||||||||||||||||||||||||||||||
менного обсуждением понятия аналитического продолжения функции f (z) . |
|
|||||||||||||||||||||||||||||||||||||||||
Пусть однозначная функция f1(z) определена и аналитична всюду в об- |
|
|||||||||||||||||||||||||||||||||||||||||
ласти D1. Функция f2(z), определённая и аналитическая в области D2, являет- |
|
|||||||||||||||||||||||||||||||||||||||||
ся аналитическим продолжением функции f1(z), если существует открытая |
|
|||||||||||||||||||||||||||||||||||||||||
область D3, являющаяся пересечением областей D1 и D2, в которой функции |
|
|||||||||||||||||||||||||||||||||||||||||
f1(z) и f2(z) |
совпадают. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Аналитическое продолжение f2(z) определяется единственным образом |
|
|||||||||||||||||||||||||||||||||||||||||
по значениям f1(z) в области D3 |
и удовлетворяет каждому функциональному |
|
||||||||||||||||||||||||||||||||||||||||
уравнению, в частности каждому дифференциальному уравнению, которому |
|
|||||||||||||||||||||||||||||||||||||||||
удовлетворяет f1(z). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Стандартный метод построения аналитического продолжения базиру- |
|
|||||||||||||||||||||||||||||||||||||||||
ется на использовании ряда Тейлора, сходящегося в круге |
|
z a |
|
r . |
Для |
|
||||||||||||||||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||||||||||
каждой точки z = b внутри круга значения f(b), |
|
f (b), … известны и можно |
|
|||||||||||||||||||||||||||||||||||||||
записать разложения f(z) в ряд Тейлора в окрестности точки z = b. Этот но- |
|
|||||||||||||||||||||||||||||||||||||||||
вый степенной ряд сходится внутри круга |
|
|
|
z b |
|
r , часть которого может |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
лежать вне первого круга. Таким образом, мы получаем аналитическое про- |
|
|||||||||||||||||||||||||||||||||||||||||
должение f(z) в часть второго круга, лежащую вне первого. Этот процесс мо- |
|
|||||||||||||||||||||||||||||||||||||||||
жет быть продолжен. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
102
6.2. Гамма-функция
Изучение специальных функций во всех руководствах начинается с гамма-функции. Это связано с тем, что знакомство с её свойствами создаёт необходимую базу для изучения других специальных функций.
Рассмотрим хорошо знакомый из курса математического анализа инте-
грал I t ne t dt , где n – целое число. Интегрируя по частям, легко показать,
0
что I = n!, где n! 1 2 3 ... n и 0! 1. Если теперь рассматривать интеграл
|
t t z 1dt , где |
|
z e |
z x j y – комплексное число, причём x > 0, то при |
|
0 |
|
|
произвольных z он не берётся (не выражается через элементарные функции). Этот интеграл, как функция параметра z, называется гамма-функцией
Эйлера. При Re z 0 z является аналитической функцией и может быть аналитически продолжена на всю комплексную плоскость. z аналитична на всей комплексной плоскости за исключением точек 0, –1, –2, …, –n, в которых она имеет полюса первого порядка.
Интегрируя z по частям, получим основное функциональное соотношение z z 1 z 1 . Если z = n, то последовательное применение этого соотношения даст нам связь между гамма-функцией и факториаломn n 1 !. Таким образом, гамма-функция обобщает понятие факториала для произвольных значений аргумента.
Важную роль при различных преобразованиях и вычислениях выражений, содержащих гамма-функции, играют следующие функциональные соот-
ношения: |
|
|
|
|
|
|
|
|
|
|
|
|
22z 1 |
|
1 |
|
|
|
|
|
|
|
|
|
|||||
z 1 z |
|
, |
z z |
|
|
2z . |
|||
sin z |
|
||||||||
|
|
|
|
2 |
|
|
|
|
|
Часто бывает необходимо знать значение гамма-функции при полуцелом аргументе z n 1 2 . Пользуясь первым из записанных соотношений,
при z 1 2 получим: 2 1
2 или 1 2 
, а на основе функционального соотношения n 1 2 n 1 2 n 1 2 получим окончательно:
103
n 1 2 2n 1 !! |
|
|
2n 1 !! 1 3 5 ... |
, где |
|||
2n |
|
|
|
Рассмотрим асимптотическое представление гамма-функции, считая z = x вещественной переменной, причем x >> 1. При решении этой задачи будем пользоваться методом Лапласа, являющимся частным случаем метода перевала. Рас-
|
|
p x dx , где p – |
смотрим интеграл I p |
f x e |
|
|
|
|
2n 1 .
(x)
x exp(–p1 (x))
exp(–p2 (x))
фиксированное положительное число, называемое |
|
p2 > p1 |
в дальнейшем большим параметром; f (x) и φ(x) – |
|
x |
непрерывные функции. Предположим вначале, |
x0 |
|
б |
|
|
что функция φ(x) имеет один минимум в точке x0 |
Рис. 6.2 |
|
(рис. 6.2). |
|
|
При больших значениях p существенным является поведение функции φ(x) в малой окрестности точки минимума x0. Разлагая функцию φ(x) в ряд Тейлора в окрестности точки x = x0 и удерживая первые три члена, получим:
x x |
x0 x x 2 . |
|
||||||
|
|
|
0 |
|
2! |
0 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
0 , а |
|
Заметим, что в точке x0 функция φ(x) имеет минимум, поэтому x0 |
||||||||
x0 0. Считая далее, что функция |
f (x) |
мало меняется в области, в кото- |
||||||
рой отлична от нуля функция e p x |
(чем больше p, тем уже эта область и |
|||||||
справедливее это предположение), и полагая |
f x f x0 , получим: |
|
||||||
|
|
|
|
|
|
|||
I p f x e p x dx f x0 e |
p x0 e p x0 x x0 2 dx . |
|
||||||
|
|
|
|
|
|
|||
После замены переменной x x0 |
t и вычисления интеграла получим |
|||||||
окончательно: |
|
|
|
|
|
|||
|
|
|
|
|
|
|
||
I p |
|
2 |
|
f x |
e p x0 . |
(6.6) |
||
p x0 |
||||||||
|
|
0 |
|
|
||||
|
|
|
|
|
||||
Полученная формула без труда обобщается и на случай, когда функция |
||||||||
φ(x) имеет несколько минимумов в области интегрирования. В этом случае
|
|
|
|
n |
f xk |
|
|
|
|
|
|
|
|
|
|
I p |
2 |
|
|
|
|
p x |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
e |
|
k |
|
, |
(6.7) |
||
|
|
|
|
|
|
|
|
||||||||
p |
|
x |
k |
|
|
|
|
||||||||
|
|
|
|
k 1 |
|
|
|
|
|
|
|
|
|
||
104
где х1, х2, …, хп – точки минимума функции φ(x). Вывод формул (6.6), (6.7) базировался на предположении, что функция f (x) не меняется в малой окрестности точки x0. Более точный результат можно получить, если функ-
цию f (x) разложить в ряд Тейлора в окрестности точки x0, т. е.
|
|
I p |
|
k 0
|
|
|
|
|
|
|
|
x0 |
x x0 |
|
|
|
k x |
|
|
|
|
|
p x0 |
|
|
|
|
f |
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
|
||||||
x x |
|
k e |
|
|
|
|
|||||
|
0 |
|
|
0 |
|
|
|
|
dx . |
||
|
k! |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Меняя местами операции суммирования и интегрирования и выполняя
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
замену переменной |
|
|
p x0 |
|
x x |
|
t , получим: |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
2 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
f k x |
|
e p x0 2 k 1 2 |
|
t |
2 |
|
|
|
||||||||||||
I p |
|
0 |
|
|
|
|
|
|
|
|
t k e |
dt |
, |
(6.8) |
||||||||||||
|
|
|
k 1 |
2 |
|
|||||||||||||||||||||
|
|
|
k 0 |
|
|
|
k! p x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
или с учётом того, что интеграл в (6.8) при нечётных k равен нулю, |
|
|||||||||||||||||||||||||
|
|
|
f |
2k x |
e p x0 2 2k 1 2 |
|
|
|
|
|
2 |
|
|
|||||||||||||
I p |
|
|
0 |
|
|
|
|
|
|
|
2 t 2k e t |
|
dt . |
(6.9) |
||||||||||||
|
|
|
|
|
2k 1 2 |
|
|
|||||||||||||||||||
|
k 0 |
|
2k ! p x0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Интеграл, входящий в (6.9), равен |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2k e |
t 2 dt |
|
2k |
1 !! |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
t |
|
, |
|
|
|
|
6.10) |
||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
2k |
|
|
|
|
||||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
в чем можно убедиться, дифференцируя k раз по левую и правую части ра-
|
|
|
|
|
|
|
|
|
|
|
|
|
|
венства e t 2 dt |
1 |
|
|
|
и полагая затем = 1. |
|
2 |
|
|
|
|||
0 |
|
|
|
|
||
|
|
|
|
|
|
Подставляя (6.10) в (6.9) и учитывая, что 2k! 2k !! 2k 1 !!, где (2k)!! – произведение последовательных чётных чисел от 2 до 2k включительно, получим окончательно:
|
|
|
|
|
|
|
|
|
|
|
2k x0 |
|
|
|
I p |
2 |
|
|
|
p x |
|
|
f |
|
|
||||
|
|
|
e |
|
0 |
|
|
|
|
|
. |
(6.11) |
||
p x0 |
|
|
|
2k !! p x0 k |
||||||||||
|
|
|
|
|
|
k 0 |
|
|
||||||
Естественно, что формула (6.11) может быть обобщена и на случай, когда функция φ(x) имеет несколько минимумов. При интегрировании по конечному или полубесконечному интервалу и использовании полученных формул возникают дополнительные ошибки, связанные с тем, что максимум функции
e p x может располагаться достаточно близко к краям интервала интегриро-
105
