Министерство образования и науки Украины Харьковский национальный университет имени В. Н. Каразина Б. В. Кондратьев МЕТОДИЧЕСКИЕ УКАЗАНИЯ к решению задач по теории функций комплексного переменного для студентов радиофизического и физического факультетов Харьков – 2010 УДК 517.53/.54(075.8) ББК 22.161.5я73 М 54 Рецензенты: доктор физико-математических наук, профессор, заведующий кафедрой теоретической радиофизики радиофизического факультета Харьковского национального университета имени В. Н. Каразина Колчигин Н. Н.; доктор физико-математических наук, профессор, заведующий отделом радиофизической интроскопии Харьковского института радиофизики и электроники имени академика В. Я. Усикова Масалов С. А. Рекомендовано к печати Научно-методическим советом Харьковского национального университета имени В. Н. Каразина (протокол №6 от 19.06.09) М 54 Методические указания к решению задач по теории функций комплексного переменного для студентов радиофизического и физического факультетов / Сост. : Кондратьев Б. В. – Х. : ХНУ имени В. Н. Каразина, 2010. – 76 с. Методические указания имеют целью оказать помощь студентам радиофизического и физического факультетов при подготовке к практическим занятиям по теории функций комплексного переменного и ее приложениям (свойства аналитических функций, их отображения, интегралы и вычеты, операционные методы, плоское электростатическое поле, асимптотические методы и др.). Этот материал входит в первую часть курса "Методы математической физики" и излагается в четвертом семестре. В указаниях к каждой теме приводятся перечень основных вопросов, ссылки на литературу, решения типичных задач и задания для самостоятельной работы. УДК 517.53/.54(075.8) ББК 22.161.5я73 © Харьковский национальный университет имени В. Н. Каразина, 2010 © Кондратьев Б. В., сост., 2010 © Дончик И. Н., макет обложки, 2010 2 Оглавление ВСТУПЛЕНИЕ…………………………………………………………. 4 Тема 1. Комплексные числа и действия над ними. Функции комплексного переменного, их аналитичность…... 5 Тема 2. Конформные отображения, осуществляемые различными элементарными аналитическими функциям………………... 7 Тема 3. Интегрирование ФКП и степенные ряды……………………. 11 Тема 4. Ряды Лорана и особые точки ФКП…………………………... 14 Тема 5. Вычеты и их применения…………………………………….. 18 Тема 6. Операционное исчисление…………………………………… 25 Тема 7. Специальные функции и их свойства……………………….. 31 Тема 8. Теория плоского электростатического поля………………… 36 Тема 9. Основы асимптотических методов…………………………... 40 Задания для самостоятельного решения................................................ 57 СПИСОК ЛИТЕРАТУРЫ…………………………….......…………… 75 3 Вступление Теория функций комплексного переменного (ТФКП) как продолжает, так и расширяет идеи математического анализа функций действительного переменного. Обычные определения, известные из алгебры чисел и математического анализа функций действительного переменного, остаются почти без изменений, но их содержание меняется весьма существенным образом. Хорошо известно, что уже обычные простейшие операции над действительными числами могут вывести за пределы их области. И решения большинства алгебраических уравнений не могут быть выражены только обычными действительными числами. Поэтому приходится расширять область действительных чисел, а таким расширением этой области и является область комплексных чисел. Причем замечательным свойством комплексных чисел является тот факт, что основные математические операции над комплексными числами не выводят из области комплексных чисел. Однако, множество комплексных чисел не обладает свойством упорядоченности; не существует правил сравнения этих чисел понятий больше или меньше. Основное понятие комплексного анализа аналитическая функция. Это понятие позволяет доказать теоремы о существовании производных любого порядка от этих функций, о независимости интегралов от формы пути интегрирования. Позволяет сравнительно единообразно вычислять сложные интегралы с помощью вычетов и многое другое. Свойство аналитического продолжения различных математических соотношений позволяет единообразно перенести их в комплексную область. А построение графиков действительных функций заменяется конформными отображениями, осуществляемыми аналитическими функциями. На основе ТФКП построены все свойства специальных функций Эйлера, Бесселя, Лежандра и др. На той же основе построено операционное исчисление, с помощью которого можно упростить решение многих задач радиофизики и гидродинамики. А с помощью конформных отображений можно решать задачи электростатики и определять функцию Грина при решении задач математической физики. Часто возникает вопрос о приближенном вычислении сложных интегралов или решении дифференциальных уравнений при больших или малых значениях некоторого параметра. На такой вопрос могут дать ответ методы асимптотических оценок, хорошо разработанные в ТФКП. В предлагаемых методических указаниях приводятся образцы решения примеров по курсу ТФКП, решаемые со студентами на практических занятиях. 4 Тема 1 Комплексные числа и действия над ними. Функции комплексного переменного, их аналитичность Основные вопросы: Модуль и аргумент комплексного числа. Три формы комплексного числа. Условия дифференцируемости ФКП (условия Коши-Римана). Понятие аналитической ФКП в области. Построение аналитической функции по ее заданной части. Аналитичность и вычисление значений элементарных ФКП: степенной, показательной, логарифмической и др. Литература: 1, гл. 1, § 1-4, гл.3, § 1 (3); 2, п. 1-5, 8-10; 3, § 1, 4, 7. Задание: 4, № 1(4), 4(2,3), 23, 26, 28, 29, 68(4,6), 71, 74(1,5,6), 81(1,2,4,6), 112, 113, 137, 159, 160. Пример 1. Найти модуль и аргумент комплексного числа z  re i  1  i 3 . аргумент   Arg z  arg z  2 k , 2 где главное значение arg(1  i 3 ) =    arctg 3   . 3 Пример 2. Исследовать аналитичность функции w( z )  U  iV  z  Im z  xy  iy 2 . Так как мана Модуль r | z | 2 ; k Z , U V U  y,  2 y , x и y y x V  0 , то условия Коши- Риx U V V U   , выполняются только в точке z  0 , следоx y y x вательно, функция w( z ) нигде не аналитическая, но в нуле имеет производную w(0)  0 функция дифференцируемая. Пример 3. Построить аналитическую функцию w( z ) = U ( x, y )  iV ( x, y ) по заданной мнимой части: V ( x, y )  x 2  y 2  3 y . Найдём действительную часть 5 U ( x, y )  ( 0,0) ( x, y ) ( x, y ) V V dx  dy = y x  (0,0) (2 y  3)dx  2 xdy  3x  2 xy  C , тогда w( z )  U  iV  iz 2  3z  C . Пример 4. Найти все значения выражений: а) Вычислить Arg th(ln 2  i ) . Сначала найдём 4 th (ln 2  i )  4 Arg   1 ln 4   i e 2 1 ln 4   i e 2 1 1  1  4i 15  8i  , затем 1  4i 17 15  8i 8  arctg  2k , 17 15 k Z , где arctg    i i  Ln ( 2i )  exp i  ln 2  ( 4k  1)  = б) Вычислить (2i )  e 2      = cos(ln 2)  i sin(ln 2)  exp  ( 4k  1) при k  Z .  2  Пример 5. Найти все корни уравнений: iz  iz а) Уравнение cos z  2 преобразуем к виду e  e  4 , затем 8 1    0,155956 ...   . 15 6 e 2iz  4eiz  1  0 , откуда eiz  2  3 ; тогда z  Arc cos 2  i  Ln(2  3)  i  ln(2  3)  2 k при k  Z . Главное значение arccos 2   i  ln(2  3 ) при б) Уравнение k  0. 3e 2 z  2ie z  1  0 , 2ch z  sh z  i преобразуем к виду i i z откуда e  (1  2) , тогда z   Ln i   (4k  1) , 3 2 i z    Ln 3i   ln 3   (4k  1) при k  Z . 2 6 Тема 2 Конформные отображения, осуществляемые различными элементарными аналитическими функциями Основные вопросы: Геометрический смысл модуля и аргумента производной аналитической функции. Определение конформного отображения; необходимые и достаточные условия такого отображения. Принцип соответствия границ. Теорема Римана. Дробно-линейное отображение. Отображение на верхнюю полуплоскость луночки (двуугольника), плоскости с разрезом, полосы и полуполосы. Понятие римановой поверхности. Литература: 1, гл.1, §4(3), гл.3, §1(4), 2(3), гл.6, §1,2; 2, п.6, 24, 26, 31, 33; 3, §8, 34, 35. Задание: 4, № 43 (1,3,5), 117, 187, 210 (3), 213, 220(1,2), 229(1,3), 283(2,3), 287, 291-293, 306, 338(1,3,6), 350, 351, 356, 358, 359. Пример 6. Какая часть плоскости растягивается (сжимается) при отображении с помощью функции f ( z )  z  2 z ; как при этом ведёт себя бесконечно малая окрестность точки z0  1  i . Коэффициент растяжения равен K  f ( z )  2 z  1 , поэтому 1 2 2 2 внутренность круга z  1  ( x  1)  y  сжимается, а его внеш4 ность растягивается в K  f ( z 0 )  2 раза и поворачивается на угол    arg f ( z 0 )  arg(i )   . 2 2 Пример 7. Для функции f ( z )  z 2  1 при z  2 примем начальное значение Arg f ( 2)  0 . Считая, что Arg f ( z ) изменяется непрерывно, когда точка z делает полный оборот против часовой стрелки по окружности C  z  2 и возвращается в точку z  2 , найти приращение аргумента C Arg f ( z ). 1 i 1, 2 , тогда Arg f ( z )  (1   2 ) . При обОбозначим z  1  r1, 2  e 2 ходе по контуру C , внутри которого оказываются обе точки ветвления z  1 , каждый из углов 1 и  2 получит приращение 2 , поэтому  C Arg f ( z )  2 .   7 Пример 8. Круговую лунку (двуугольник) на рис.1 отобразить на верхнюю полуплоскость Im w  0 . На плоскости z будет z A  0 и zC  i , а координаты точек B и D найдём из системы уравнений 2  z  x2  y2  1   2  1  x 2  ( y  1) 2  1 z   , Рис. 1 1 1 3 , y  1 ; поэтому z B , D  ( 3  i ) . Дробно-линейная 2 2 2 z  zB 2z  3  i  выпрямляет дуги BCD и BAD и функция p ( z )  z  zD 2z  3  i отображает заданную лунку DABCD на сектор D1 BD2 в плоскости р откуда x   (рис.2), так как 4 pB  0 и pD   . 2 Поскольку i i 1 1 3 p A (0)   (1  i 3 )  e 3 и pC (i )   (1  i 3 )  e , то луч BD1 2 2 2 4 имеет направление arg p A   , а луч BD2 - направление arg pC   . 3 3 Величина угла, образованного границей лунки в точке B , не изме2   arg  arg  ABC p p . p ( z ) и равна няется при отображении A C 3 Рис. 2 Отображение q( p)  2  i pe 3 повернёт сектор D1 BD2 вокруг 2 его вершины B на угол   . На плоскости q луч BD2 совпадает 3 8 с действительной положительной полуосью; точки q A  2 i 3 e qC  1 . 3 плоскости w луч BD2 останется на месте, а луч BD1 совпадёт с действительной отрицательной полуосью. Таким образом, сектор D1 BD2 отобразится на верхнюю полуплоскость Im w  0 и wA,C  1 . Выраw через жая z, получим отображающую функцию 2 q  , то при отображении w  q 2 на arg  Так как arg qC  0 и A 3  2z  3  i   w( z )    2z  3  i    32 . Во всех плоскостях z , p, q , и w рассматри- ваемые области расположены слева от границы. Пример 9. Верхнюю полуплоскость с вырезанными полукругами (рис. 3), т.е. область Im z  0 и z  1  1 , отобразить на верхнюю полуплоскость Im w  0 . z Рис. 3 2 отобразит заданную область в полуполосу z П p  Re p  1, Im p  0 , т.е. границы области выпрямляются и  p B, F  1  i ; pC , E  1 ; p D  0 . Функция q ( p)    i  i ( p  1) Функция p ( z )    переводит полуполосу П P в полуполосу П q    0, 0    2 9 Рис. 4 (рис.4); при этом qB   2   i , qC   i , q D   2 i , qE  0 , qF   2 . Наконец, функция w  U  iV  ch q отобразит полуполосу П q на верхнюю полуплоскость Im w  0 ; причём wF , B  ch  2,50918 , 2  wE ,C  1 , wD  0 . Действительно, участок границы полуполосы EA1  0    ,   0 с помощью функции wEA1  ch   1 перейдёт на плоскости w в полупрямую EA1   U  1, V  0  ; участок CA2   0    ,     функция wCA2  ch  1 переведёт в полупрямую CA2   U  1, V  0; отрезок EC     0, 0      функция wEC  cos переведёт в отрезок EC   U  1, V  0 . Таким образом заданное отображение сделает функция  1 1 w( z )  ch i (  )   sin . z 2 z 10 Тема 3 Интегрирование ФКП и степенные ряды Основные вопросы: Методы вычисления контурных интегралов. Теорема и формула Коши. Степенные ряды, их радиус сходимости. Разложение аналитических функций в ряд Тейлора. Нули функций, их кратность. Свойство единственности аналитических функций. Литература: 1, гл.1, §5, 6, гл.2, §1-3, гл.3, §1(1,2), 2(1-3); 2, п.11-14, 17-20; 3, §2, 5, 9-12, 14. Задание: 4, №388 (1,2), 393 (2,3,5), 403-405, 412, 434, 470, 471, 506-509, 523 (1,2), 525.  Пример 10. Вычислить интеграл I   z L 1 3 dz в положительном наx  0  от той ветви правлении по левой полуокружности L  3  z  1, корня, для которой i  i . Для выделения заданной ветви, выберем нужное значение числа k из трёх возможных k  0,1 . Для этого в начальной точке контура интегрирования L запишем 3 i     i  exp   2 k   i  exp  i  , 3 2   2 1  1 3   k   2   2k   , значит k  1 .   или откуда 3 2 2 2 2  i Под знаком интеграла сделаем замену z ( )  e при  2   3 2 и с учётом требуемой ветви подынтегральной функции запишем 3 1  z 3 i  (  2 ) e 3 I   2  i (  2 ) e 3  i ei d  2 3 2 2i  i ie 3 e3 d 2    2 i ie 3 2 3 i  e3 2i  2 3 2  2 11 3 3  i i 3 3 i i 6 i 3  e  (e  e )   e  e 2 2 2 2    i i  e 6  e 6        i 3  3 e 6 2 5 i 3 e 6  3  6 2   3 3   3   cos  i sin   3  ( 3  i) . 2 6 6 4   cos  3e  dz 3 , где C самопереC z (1  z )   1 i  1   6 Пример 11. Вычислить интеграл I   секающийся замкнутый контур (”восьмёрка”), который обходит особую точку z0  0 в отрицательном направлении, а z1  1 в положительном. Деформируем контур C в ”гантель”, по перемычке которой интегралы взаимно уничтожаются и остаются интегралы по окружностям около особых точек C0   z    и C1   z  1    (здесь 0    1 ); поэтому dz  1  dz 1 I        . z  0 C1 z  ( z  1) 21 C0 1  z  Затем, воспользовавшись интегральной формулой Коши, получим 3  1 3 I  2 i   1 z     z 0 1 1    2!  z      4 i . z1   Пример 12. Функцию f ( z )  рестности точки z0  1 и найти область сходимости ряда. Воспользуемся формулой для суммы бесконечно убывающей геометрической прогрессии и получим  ( 1) k 2 1 f ( z)  1   1  2  k 1 ( z  1) k . 3 1  1 ( z  1) k 0 3 3 z разложить в ряд Тейлора в окz2 Важно учесть, что радиус сходимости ряда – это расстояние от точки z  2 раскладыразложения z0  1 до ближайшей к ней особой точки ~ 12 ваемой функции f ( z ) (в этой точке f (2)  f (2)   ); поэтому 2 R  z0  ~ z  3 , а область сходимости z  1  ( x  1) 2  y 2  9 . Пример 13. Разложить в ряд Тейлора в окрестности точки z0  i ту ветвь функции f ( z )  z 2 3 , для которой f (i )  i Представим функцию в виде 2 3  1 . 2 3 f ( z)  z 2 3  i  ( z  i )  i 2 3  1  i ( z  i )  3 . 2 Для выбора нужной ветви проще использовать первый множитель i 2 3  i  2 ki    e 2      2 3 e i  3 ( 4 k 1) при k  0,  1. Очевидно, нужная ветвь 2 i   1. Затем для функции получится при k  1 ; тогда i 3 k  1  e второго множителя можно принять главное значение и воспользоваться формулой разложения бинома   1  k   k!   (  1)(  2)...(  k  1) биномиальный коэффициент.    2i k  2 3 k   f ( z )   (  i )  ( z  i )   1  ( z  i)   Тогда получим  k  3 k 0    1 4i 7 2 3 ( z  i )  ( z  i )  ( z  i ) 4  ... 2 4 5 3 3 3 2   k (1  t )    k   t , где k 0     t 1 и Ряд сходится в круге z  i  x 2  ( y  1) 2  1 . Радиус сходимости ряда R  i  0  1 , где i центр круга разложения и 0 ближайшая особенность (точка ветвления третьего порядка). Пример 14. Определить порядки всех нулей функции f ( z )  ( z 2  4) z 5 . 13 Функция обращается в нуль при z1, 2  2i , но её производная в этих точках z1, 2 в 1 1 2 ( 3 z  20 )   0 , значит, точки 8 z6 z  2 i  2i являются простыми нулями. Для определения порядка нуля z3   z 1 t , бесконечности сделаем замену тогда f ( 2i )   f (1 t )   (t )  t 3 (4t 2  1) ; функция  (t ) имеет при t  0 нуль третьего порядка, значит, и функция f ( z ) имеет при третьего порядка. z3   так же нуль Пример 15. Проверка тождественности двух аналитических функций. Функции f1 ( z )  sin  плоскости переменного z , за исключением точки z  0 (там существенная особенность). Эти функции совпадают на последовательности z 2 и f 2 ( z )  sin  z аналитические во всей 1 (k  N , k  0) . Однако функции не тождественны, так как эта k последовательность имеет предельную точку z  0 , а функции f1 ( z ) и f 2 ( z ) в нуле не аналитические. zk  Тема 4 Ряды Лорана и особые точки ФКП Основные вопросы: Методы разложения функций в ряд Лорана. Разложение в окрестности бесконечности. Единственность разложения в ряд. Классификация особых точек аналитических функций; поведение функций в окрестности особых точек. Связь между нулями и полюсами. Литература: 1, гл.3, §2(5); гл.4, §1, 2; 2, п.21, 22, 25; 3, §17, 18. Задание: 4, №543, 545, 549, 551, 559, 561(1,2,4,8,11), 562(1-3,6), 566-571, 576-578, 580, 582-585, 601, 602. 1 Пример 16. Функцию f ( z )  при 0  a  b   разло( z  a)( z  b) жить в ряды Лорана в окрестности точек z1  a , z2  b , z3   и в кольце a  z  b ; определить области сходимости этих рядов. 14 Для разложения рациональных функций в ряды обычно бывает достаточно разложить их на простые дроби и воспользоваться формулой суммы бесконечно убывающей геометрической прогрессии; поэтому в окрестности точки z1  a получим   1 1 1  1  1  f ( z)      b  a  z  a b  a 1  z  a  b  a 2 ba  Этот ряд сходится при  za  ba .  k  1   k za  1 , т. е. внутри круга 0  z  a  b  a с ba выколотым центром; здесь R1  a  a  0 и R2  b  a  0 . Аналогично получается разложение в окрестности точки z 2  b k 1 f ( z)  b  a 2 z b   (1)  a  b  при 0  z  b  a  b .   k  1  k Разложение в окрестности бесконечности будет   1 1 1 1 1  1  bk  ak   f ( z)      ; b  a  z 1  b z 1  a  b  a k 1 z k 1  z z ряд сходится вне круга z  b , т. е. в области b  z   . Разложение внутри кольца 0  a  z  b   имеет вид:   1   z k 1  1 1 1 1  ak       k 1  f ( z)     k 1 ;  z a b  a  b 1 z 1   b  a k 0  b z   b z  здесь первая часть разложения сходится при a  1. z 15 z  1 , а вторая – при b По виду разложений заключаем, что в точках z1  a и z 2  b находятся простые полюсы функции f ( z ) , а в точке z3   - устранимая особенность. Вычеты функций в этих точках равны Re s f ( z )   Re s f ( z )  z a z b 1 , ba 3  1 z Re s f ( z )  0 . z  Пример 17. Функцию f ( z )  z  e разложить в ряды Лорана в окрестности всех её особых точек; определить области сходимости. При разложении в ряды Лорана элементарных аналитических функций обычно используют их известные разложения в ряды Тейлора. Предложенная функция имеет две особенности при z1  0 и z 2   ; её лорановское разложение будет: 1 1 1    z 3  1     ...   2 3 3! z  1! z 2! z  2 z z 1 1 1  z 3   ....    ... + 3! 1! 2! 4! z 5! z 2 Правильн. часть при z  0 Главн. часть при z  0 f ( z)  z e 3  1 z Главн. при z   Правильн. часть при z   Здесь разложения в окрестностях нуля и бесконечности по форме совпадают. Разложение сходится в области 0  z   . По виду разложения заключаем, что в точке z1  0 находится существенная особенность, а в точке z 2   полюс третьего порядка. Вычеты в этих точках равны: Re s f ( z )   Re s f ( z )   z 0 z  1 . 24 Пример 18. Найти особые точки функций и выяснить их характер. 1 3 имеет 3 . Так знаменатель z  z  z (1  z )(1  z ) zz простые нули при z1  0 , z 2  1 и z3  1 , то для функции f ( z ) это будут простые полюсы. Точка z0   является устранимой особенностью, так как f ()  0 . а) f ( z )  16 б) f ( z )  z e . Точка z0   будет существенно особой, так как f ()  0 и f ()  0 . в) f ( z )  z  cth z . Точка z0  0 устранимая особая, там f (0)  1. z k  ik  0 (k  Z , k  0) , в которых sh z k  0 и Точки 3 z ( sh z k )  ch z K  (1) k  0 будут простыми полюсами функции. В бесконечности z   находится неизолированная особенность – точка предельная для полюсов. г) 1 f ( z )  z 2  cos . Представим функцию в виде ряда Лорана z 1 1 1    ... , очевидно, она имеет при z1   полюс 2! 4! z 2 6! z 4 второго порядка и при z 2  0 существенную особенность. f ( z)  z 2  1 . Функция двузначна, она имеет точки ветвле1 z 1 ния второго порядка при z1  1 и z 2   . Кроме того, та ветвь функции, для которой  1  1 , имеет простой полюс в точке z3  2 ; д) f ( z )  действительно, числитель функции там постоянный, а её знаменатель (1  z  1) z  2  (1  1) 1 1  0 , но производная от знаменателя 1 1    0 . Для другой ветви, z 2 z 2 2 2 1 1 1 где 1  1 , точка z3  2 правильная, так как там существуют конечные (1  z  1)  1 2 z 1 и  значения f (2)  1 2 1 f (2)   . 8 Пример 19. Выяснить, допускают ли указанные функции f ( z ) (или их однозначные ветви) разложение в ряд в окрестности заданной точки z0 ; если разложение возможно, указать область сходимости ряда D. 1 при z0  0 . Точка нуль является изолированной z существенно особой точкой функции f ( z ) . Поэтому возможно лораа) f ( z )  z  sin новское разложение в кольцевой области D  0  z   .   17 б) f ( z )  ctg z при z0   . Точка бесконечность является неизолированной особой точкой, она предельная для полюсов zk   k (k  Z ) ; поэтому требуемое разложение невозможно. z0  1. Многозначная функция f ( z )  ln( z  1) при f k ( z )  Ln( z  1)  ln( z  1)  2 k i ( k  Z ) имеет при z0  1 точку в) ветвления логарифмического типа. Для выделения однозначной ветви f ( z )  ln( z  1) при k  0 нужно на плоскости z провести разрез из точки z0  1 в бесконечность. Поэтому в окрестности точки z0  1 нельзя найти неразрезанное кольцо вида D  0  z  1   , в котором было бы возможно разложение функции f ( z ) . 2 г) f ( z )  3 z ( z  1) при z1  0 , z 2  1 и z3   . Точки z1  0 и z 2  1 точки ветвления третьего порядка функции f ( z ) , поэтому в их окрестностях разложения в ряды невозможны. Бесконечность – простой  1 3 полюс функции f ( z )  z 1    1 , поэтому в окрестности точки  z z3   возможно лорановское разложение в области D  1  z  . 23 Тема 5 Вычеты и их применения Основные вопросы: Вычеты аналитических функций. Методы их вычисления в изолированных особых точках разного типа. Теорема Коши о вычетах. Приложение теории вычетов к вычислению определённых интегралов. Четыре основных леммы. Исследование нулей аналитических функций. Литература: 1, гл.5 §1,2; 2, 23, 24, 73-75; 3, §28-30. Задание: 4, № 621, 628, 643, 658, 659, 662, 673, 676, 678, 680, 683, 685, 692-694, 696-698, 702, 704, 714, 716, 733, 788, 793. z 9 dz 2 Пример 20. Вычислить контурный интеграл I r   5 Cr ( z  1) вдоль Cr   z  r  1. положительного направления обхода окружности 18    z k  5  1  exp i (2k  1)  при k  0, 1, 2, 3, 4 и устранимую  5  особенность в бесконечности (там f ()  0 ). Поэтому I r  0 при r  1 , а при r  1 получим: Подынтегральная функция f ( z ) на единичной окружности z  1 имеет пять полюсов второго порядка в точках I r  2 i   Re s f ( z )  2 i  Re s f ( z )  k 0 z  z k z  2 1  1   2 i  Re s  1  5    2 i . z  z  z      4 Пример 21. Вычислить интеграл от тригонометрической функции I ( )   tg (  i )d при действительном 0 3   0. Сделав замену z ( )  e турному интегралу I  i   3 0 2i ( 2 i ) и d  dz 2iz , перейдём к кон- e 2i ( i )  1 3 z 1 d   dz ;   e 2i ( i )  1 2 z ( z  1) C  2 1 здесь C  z  e   окружность, пробегаемая в положительном направлении. Подынтегральная функция f ( z )  z 1 имеет простые z ( z  1) полюсы в точках z1  0 и z 2  1 и устранимую особенность при z3   . Если число   0 , то внутрь контура C  z  e 2  1 попадает одна особая точка z1  0 ; если же   0 , то внутрь конту-    1 попадают две особенности z1  0 и z2  1 . ра C  z  e Поэтому при   0 получим:  2  I   3 2 C z 1 z 1 ...  3 i  Re s z  1  3 i   3 i ; z 0 z z  1 z 0 19 а при   0 будет z 1 1  1   3 z   3 i I     ...  3 i  Re s z  1  Re s . z 0 z  1 z  1  2 C z      Окончательно результат запишем в виде: I ( )  3 i  sign , I (0)  0 . Пример 22. Вычислить интеграл по полуоси I   a, b  0 . Выберем вспомогательный контурный sin ax  dx  x( x 2  b 2 ) 2 при 0 интеграл вида ~ I  e  dz 2 2 2 , где новая подынтегральная функция имеет простой C z( z  b ) полюс z1  0 , полюсы второго порядка z2,3  bi и существенно особую точку iaz z 4   . В верхней полуплоскости Im z  0 интегрирования вида C  ( , r )  C r  (r , )  C (рис. 5), где 0    b и b  r   . Тогда вспомогательный интеграл легко вычисляется по теореме о вычетах (все контуры интегрирования будем Рис. 5 обходить только против часовой стрелки в положительном направлении) выберем контур ~ I  2i  Re s z bi eiaz  2 2 2 z( z  b )  e iaz     2 i  2   z ( z  ib)    z bi i 2b 4 (ab  2) e  ab . ~ Затем распишем интеграл I по четырём участкам контура C ; это даст ~ r I    eiax dx eiaz dz eiax dx eiaz dz    2 2 2 2 2 2 2 2 2   z( z 2  b2 )2 . x ( x  b ) Cr z ( z  b )  r x ( x  b ) C  20 Сделав в третьем интеграле замену переменной x   x , после преобразований получим: sin ax  dx eia z dz eiaz dz  i  ab    2i    ( ab  2 ) e   . 2 2 2 2 2 2 2 2 2   2b 4  x( x  b ) Cr z ( z  b ) C z ( z  b ) r В этом равенстве перейдём к пределам   0 и r   , тогда sin ax  dx eiaz dz eiaz dz  lim   lim  2i     z ( z 2  b2 )2  2 2 2  0 r  z ( z 2  b 2 ) 2  x ( x b ) 0 Cr C   i (ab  2)e  ab . 4 2b  Здесь предел второго интеграла равен нулю по лемме Жордана, так как lim a0 и z  вычета равен: 1  0 , а предел третьего по лемме о половине z( z 2  b2 )2 e iaz eiaz dz i ( z 2  b2 )2 lim  i Re s      .  0  z ( z 2  b 2 ) 2 z 0 z 0 b4 C После простых преобразований получаем окончательно:  I sin ax  dx 1   1   ab   1  1  ab e  . 2 2 2 4 ( ) 2 x x b b     2  0 Из последнего выражения легко получить значения многих других нетривиальных интегралов, например: cos ax  dx   ab   2 Ia  ( ab  1 ) e , 2 2 4b 3 0 (x  b )   Ia a 0   dx   . 2 2 2 4b 3 0 (x  b ) Пример 23. Вычислить интеграл от алгебраической функции x p  dx I 2 при 0     и p  1 . x  2 x  cos   1 0  21 ~ I  Выберем вспомогательный интеграл вида z p  dz , 2 C z  2 z  cos   1 где новая подынтегральная функция в общем случае бесконечнозначная p (так как z  exp( p  Ln z )  exp( p  ln z  2 k pi ) при k  Z ) и имеет и z4   ); кроме того, каждая ветвь функции имеет простые полюсы i i ( 2  ) (здесь z1, 2  1 и 0     ), которые сопри z1  e и z 2  e логарифмические точки ветвления в нуле и бесконечности (при z3  0 ответствуют простым нулям знаменателя z 2  2 z  cos   1  ( z  z1 ) ( z  z2 ) . Выберем главную ветвь подынтегральной функции (когда k  0 и 0  arg z  2 ), для этого на плоскости комплексного переменного z проведём разрез по действительной положительной полуоси 0  Re z   и положим на верхнем берегу разреза arg z  0 . Выберем контур интегрирования C , проходящий по берегам этого разреза и дополненный двумя окружностями (рис.6); пусть C  Cr  ( r ,  ) arg z  2  C  ( , r ) arg z 0 , где 0    1  r   22 Рис. 6 Тогда вспомогательный интеграл I легко вычисляется по теореме о вычетах: ~  z p ( z  z2 ) z p (z  z 1 ) ~  I  2 i s  Res  Re  z  z1 zz2 z z  z z  1 2   z  z2  2 i  1 z1  z 2  2 i p p eip  eip ( 2  )  2 i  i  e  ei ( 2  ) sin p (   ) i p e . sin  ~ Распишем интеграл I по четырём участкам контура C и получим: r ~ 2p i I  (1  e )    x p  dx z p  dz  2  2  x  2 x cos  1 Cr z  2 z cos   1 z p  dz sin p (   ) i p   2 i   2 e .  z  2 z cos   1  sin C Здесь учтено, что на верхнем берегу разреза arg z  0 и z p  ( x  ei 0 ) p  x p  0 , arg z  2 а на нижнем и z p  ( x  e 2 i ) p  x p  e 2 p i (здесь x p  0 ). Затем перейдём к пределам   0 и r   , тогда (1  e 2 p i z p  dz z p  dz  lim   )  I  lim  2 r   z  2 z  cos   1  0  z 2  2 z  cos   1 Cr C  2 i  sin p (   )  p i e . sin  Здесь предел второго интеграла по лемме об оценке равен нулю, так как лемме об интегрируемой особенности тоже равен нулю, так как 23 z p 1 lim 2  lim z p 1  0 при p  1 ; предел третьего – по z  z  2 z cos   1 z  z p 1 lim 2  lim z p 1  0 при p  1. Следовательно, для z 0 z  2 z cos   1 z 0 параметра p можно указать интервал сходимости  1  p  1 . Наконец, после преобразований получим: x p  dx  sin p (   ) I 2   sin  p sin  0 x  2 x  cos   1 при 0     и p  1 . Из этого результата легко получить значения многих нетривиальных интегралов; например,  I p 0    sin  , I  0  , I    p и др.; после дифференцирования по параметру p поsin  p x p  ln x  dx  I p   2    x 2 x cos  1 0  sin  p        cos (  )  p    , sin   sin  p  sin  p   лучим: далее найдём I  p p 0  0 и т. п. 4 Пример 24. Сколько нулей функции f ( z )  z  5 z  1 находится в круге z  1 и в кольце 1  z  2 ? Положим f ( z )   ( z )   ( z )  z 4  5 z  1 . Если min  ( z )  max  ( z )  0 на окружности C   z  r, то по теореме Руше функции f ( z ) и  ( z ) имеют одинаковое число нулей 4 внутри круга z  r . Сначала выберем 1 ( z )  5 z и  1 ( z )  z  1 , тогда на окружности C1   z  1 4 получим min 1 ( z )  5 zC1 и max  1 ( z )  z 4  1 zC1 C1  ( z  1) C1  2 , поэтому внутри круга z  1 24 находится один корень функции f(z), так как там находится один корень функции 1 ( z ) . Затем выберем сти  2 ( z )  z 4 и  2 ( z )  5 z  1 , тогда на окружнополучим min  2 ( z )  z zC2 4 C2 C2   z  2  16 и max  2 ( z )  5 z  1 zC2 C2  (5 z  1) C2  11, поэтому внутри круга z  2 функция f ( z ) имеет четыре корня. Внутри кольца 1  z  2 функция f ( z ) имеет 4  1  3 корня; причём на окружностях z  1 и z  2 корней нет. Тема 6 Операционное исчисление Основные вопросы: Свойства функций оригиналов и изображений. Показатель роста. Функция Хевисайда. Интегралы Лапласа и Римана-Меллина. Теоремы подобия, смещения, дифференцирования и интегрирования, свёртки, разложения. Методы решения интегральнодифференциальных уравнений. Вычисление определённых интегралов. Расчёт радиоконтуров. Литература: 1, гл.8, §1, 2(1,3), 3(1); 2, п.79-81,83,84; 3, §47-49. Задание: 5, №3, 5(1,2), 7, 8(1-3), 13(1-2), 49(2,3,6), 50(1,2,4,10), 53(1), 96(2,4), 97(1,4,7), 100(1,3,4). Пример 25. Являются ли заданные функции оригиналами? e t f1 (t )   (t ) - да; t 1 et  (t ) - нет. f 2 (t )  t 1 Являются ли заданные функции изображениями? F2 ( p)  tg p p 2 - нет. at bt Пример 26. Найти изображение функции f (t )  (e  e ) t при a, b  0 и a  b . 25 F1 ( p)  th p p 2 - да; Функция f (t ) удовлетворяет всем трём условиям, накладываемым на функцию-оригинал; её показатель роста S 0  max(a, b) . Так как известно изображение более простой функции f 0 ( t )  e at  e bt  F 0 ( p )  по теореме об интегрировании  p b 1 f (t )  f 0 (t )  F ( p )   F0 ( q ) dq  ln . t pa p изображения. то 1 1  , p a p b изображения Полученная найдём функция F ( p ) удовлетворяет обоим условиям, накладываемым на функцииp2  a2 Пример 27. Найти оригинал для функции F ( p )  2 при ( p  a 2 )2 a0. Функция F ( p ) удовлетворяет обоим условиям для функцийизображений; для неё S 0  1 . Известен оригинал более простой функции F0 ( p)  p  f 0 (t )  cos at , 2 2 p a поэтому по теореме о дифференцировании изображения найдём  F ( p )   F0 ( p )  f (t )  t  f 0 (t )  t  cos at . Полученная функция f (t ) удовлетворяет всем трём условиям, накладываемым на оригиналы. Пример 28. Воспользовавшись формулой обращения (интеграл Римана-Меллина), найти оригинал f (t ) для функции-изображения F ( p)  Функция F ( p ) на плоскости p  S  i имеет две точки ветвления второго порядка при p  0 и p   ; её показатель роста S 0  0 , и можно выбрать путь интегрирования S1  1 . Для выделения главного значения функции F ( p ) проведём разрез    S  0 и положим на его верхнем берегу arg p   . 26 1 a e p p при a  0 Рис. 7 pt Рассмотрим вспомогательный интеграл I   F ( p )e dp по замкC 2 2 нутому контуру C  ( S1  i r  S1 , S1  i r  S1 )  2 2  Cr  ( r , ) arg p    C  ( , r ) arg p    Cr  при 0    r   (рис. 7) .По теореме о вычетах он равен нулю. На i верхнем берегу разреза p  Se ( S  0) , там F ( p )    i нижнем берегу p  Se ( S  0) и F ( p )   1 ia e S S ; на 1 ia e S S ; поэтому сумма S t  sin a S  интегралов по берегам разреза равна 2i   e r После перехода к пределам   0 и r   получим S1  i S1 i  dS . S  F ( p)  e pt dp  lim  0 r  Cr   F ( p)  e pt dp  lim  F ( p)  e pt dp   0 C  2i   e  st  sin a s  dS  0. S   Здесь пределы интегралов по контурам Cr   p  r , S  S1 ,   0   при r   обращаются в нуль по видоизменённой лемме Жордана, так 27   как показатель t  0 и F ()  0 . Предел интеграла по окружности при   0 после замены p   e легко вычисляется: pt C   p    i (     )  0 lim  F ( p )  e dp  i lim  exp(t  e i  a   e i 2 )  d  C   0  =  2 i . Интеграл по прямой Re p  S1  0 даёт искомый оригинал 2 i f (t ) , поэтому после преобразований получим ( S   ) 2 f (t )  1 S1  i 2 i S1 i  F ( p)  e dp  1  pt 2    e  t  0  2 sin a  d . Этот результат можно представить в виде (см. пример 33): 2   2  a  f (t )  1  Ф   e d .  , где Ф( )   2 t   0 Пример 29. Найти частное решение линейного обыкновенного дифференциального уравнения с постоянными коэффициентами и начальными условиями (задача Коши): y  4 y ( x)  2e 2 x ; y ( 0)  0 , y (0)  1 при 0  x   . Решение делится на три этапа: 1) По задаче для функции-оригинала y ( x) поставить задачу для функции-изображения Y ( p) : p 2  Y ( p )  p  y (0)  y (0)  4Y ( p )  2 или Y ( p )  ( p  4)  1  2 p2 2 . p2 p . ( p  2)( p 2  4) 2) Определить функцию-изображение: Y ( p)  28 Эта функция удовлетворяет обоим условиям для изображений; причём S0  3 . 3) По найденной функции-изображению Y ( p ) определяем функцию-оригинал y ( x) . Для этого методом неопределённых коэффициентов разложим изображение Y ( p ) на простые дроби p  1 1 2 ; Y ( p)    2  2  4 p2 p 4 p 4  откуда, зная операционные соответствия этих дробей, получаем 1 y ( x)  (e 2 x  sin 2 x  cos 2 x)  ( x) . 4 Этот результат можно также получить, воспользовавшись интегралом Римана-Меллина или теоремой об умножении изображений. Полученное решение удовлетворяет обоим условиям задачи y ( 0)  0 , y(0)  1. Пример 30. Операционные методы вычисления определённых интегралов. а) Формула обмена (равенство Парсеваля); если f (t )  F ( p)  и g (t )  G ( p ) , то  f ( x)  G ( x)dx   F ( x)  g ( x)dx . 0 0  Вычислим интеграл I    sin ax dx при a  0 ; x 0 Пусть f (t )  sin at  F ( p )   a p2  a2 и G ( p)   0 1  g (t )  1 , тогда p . a x  I 2 dx arctg 2 a 0 x a оригинал f (t )      2 б) Метод преобразования Лапласа под знаком интеграла. Примем за cos xt dx при a  0 ; тогда изображение 2 2 x a  0 29 F ( p )   f (t )e  pt  dt  0     0   1 p dx p 1    dx      2 2 2 2 2 2   2 2 2 2 p x x a p a 0 x a x p  1 x x 1 p   arctg arctg  2     2  p a p p a a    0 1 2a p  a ;  e a t .  возвращаясь к оригиналу, найдём F ( p )  f (t )   2a Пример 31. В электрической цепи последовательно соединены постоянные элементы: R активное сопротивление, L индуктивность, C ёмкость. Первоначально в цепи отсутствовали токи и заряды, но в момент времени t  0 включается ЭДС  (t )   0 (t ) , где  0  const ; требуется определить в цепи силу тока J (t ) при t  0 . С помощью правил Кирхгофа составим уравнение цепи: dJ 1 t R  J (t )  L    J ( )d   (t ) , dt C 0 J (0)  0 . Примем ток J (t ) за функцию-оригинал и перейдём к уравнению для изображений; если J (t )  I ( p) операторный ток и  (t )  E ( p) операторная ЭДС, получим уравнение для операторного тока:  0 1     ( 0 )  ( )  I ( p)  R pL L J E p   .   pC p   Выведем операторное сопротивление (импеданс) цепи Z ( p)  R  pL  2 R 1 L  0 - декремент  ( p   ) 2   2 , где   pC p 2L   1 R2  0 затухания;   CL 4 L2 частота свободных колебаний цепи. 30 0 L E ( p)  . Z ( p) ( p   ) 2   2 Функция I ( p ) обладает всеми свойствами изображения, причём S 0  0 . Теперь, воспользовавшись теоремой смещения, получим функJ (t ) t  0) цию-оригинал (силу тока при в виде Тогда операторный ток I ( p )  J (t )  0 L e  t sin  t   (t ) , т.е. в нашей цепи возбуждаются затухаю- щие колебания тока (и напряжения). Указание Для определения тока J (t ) в электрической цепи часто сразу записывают уравнение для операторного тока, где операторная ЭДС E ( P) находится по заданной  (t )  E ( p) , а импеданс всей цепи Z ( p) составляется из импедансов отдельных элементов цепи Z R  R , Z L  pL и Z C  1 pC по правилам Кирхгофа. Тема 7 Специальные функции и их свойства Основные вопросы: Гамма-функция Эйлера: определение, связь с факториалом, формулы смещения и дополнения, особые точки, значения при полуцелом аргументе, асимптотика. Бета-функция Эйлера: определение, связь с гамма-функцией. Интегралы вероятностей и Френеля; их свойства, представления в виде рядов. Функции Бесселя; представления в виде ряда и интеграла, уравнение, производящая функция, рекуррентные соотношения, функции Неймана и Ханкеля, нули, графики, функции с полуцелым индексом, асимптотики, интегралы от функций Бесселя. Функции Лежандра; формула Родрига, интегральные представления, уравнение, нули, графики, производящая функция, рекуррентные соотношения. Присоединенные функции Лежандра, их уравнение и свойства. Литература: 1, гл.8, §2(3,4); 2, п.16, 70(1,2,4), 80, 90, 94(1-3), 95 (5-7); 6, Дополнение II(1): §1(2,3), 3(1,2,), 5(2); Дополнение II(2):§1 (1-3,6,7), 2(1); 7, гл.13, §1, 2, 5-7, гл.16, §1, 3-7; 9, гл.4, §1-4. 31 Пример 32. Для  - и  - функций Эйлера ( )   e 0   x  1 x dx , B ( ,  )   x 1 (1  x)  1 dx  ( )  (  ) (   ) при Re  , Re   0 проверить следующие соотношения ( n, m  N ) :  ( n  1)  n! 0 1   (2n)  2 2 n 1  (n)   n   . 1  2  (2n)!!  B n  1,   . 2  (2n  1)!!  2  1   1 2 B(n, m)  (n  1)!(m  1)! . (n  m  2)!  2 B( ,  )  B (  , )  2   sin 0 2 1   cos   d     1 (1   )   d . 0  (1) n Re s ( )  .   n n! B( ,1   )   t  (1  t ) 1 dt  0   sin  при 0    1. 2 Указание. В интегралах сделать соответственно замены x  cos  ,   x (1  x) , t  x (1  x) . Пример 33. Проверить указанные ниже свойства интеграла вероят- ностей  ( z )  ций: 2  e 0 z  2 d при z C и связанных с ним функ- ( z )  ( z ) ; (0)  0 ;  ( )  1 . (1) k  z 2 k 1 2 z ( z )     ... при k k ( 2  1 )  !  k 0  2  z . 2        2   2  0 e  2 d  32 2    e  0  2 2   sin   d . Здесь разложить в ряд sin  , почленно проинтегрировать и воспользоваться значением  k  1 2 , тогда получится ряд для функции  ( z ) .   С функцией  ( z ) связаны интеграл ошибок Erf ( z )  и дополнительный интеграл ошибок  2 ( z ) Erfc ( z )   2 1  ( z )    e  d . 2  z С функцией  ( z ) связаны интегралы Френеля z z  i 1  2  2   C ( z )  i  S ( z )   cos z d  i  sin d .    2  2 2 2i  0 0  Эти функции нечётные и C ()  S ()  1; их можно разложить в ряды ( z  ) :  2 2  (1) k  z 4 k 1 2   z  ... . C z      k 0 (4k  1)  (2k )!     2 2  (1) k  Z 4 k 3 1 3 2  S z   z  ...     k 0 (4k  3)  (2k  1)! 3    Пример 34. Функция Бесселя J ( z ) и другие связанные с ними функции. Проверить, что функция J ( z )  является решением уравнения: k 0   (1)  z k k!(k   1)  2 2 k  при   R J  1   2 1 2   J     J ( z )  0 .  z z   Записать ряды для функций J n ( z ) при n  0, 1, 2 . Проверить, что J 0 (0)  1 ,   J 0 (0)  J n (0)  J n (0)  0 , J n ( z )  J n ( z )  (1) n  J n ( z ) и J n ( x)  1 при n  N . 33 Воспользовавшись рекуррентными формулами, доказать, что Проверить справедливость соотношений: 1    J 2 ( z)  J 0 ( z)  J 0 ( z)  2J 0 ( z)  J 0 ( z) . z H1 (1, 2 ) z  2 i   2 ( z)  e . z  2 J 1 ( z)   N 1 ( z)   cos z .   z 2 2 J 3 ( z)   N 3 ( z)  2  2  2 2  sin z    cos z  . z  z  3 Вычислить интегралы:  J 0 ( z )  zdz  z  J1 ( z ) ;  J 0 ( z )  z 2 dz  2 z 2  J 0 ( z )  ( z 3  4 z ) J1 ( z ) . 2 2   J n ( z )  zdz  z 2  J n 1 ( z )  J n ( z )  2nz  J n 1 ( z )  J n ( z )   2   ( z  n )  J n ( z)  z  J n ( z) . Доказать операционное соответствие: 1 2 2 2 2 2 2 f (t )  J 0 (at )  F ( p )  ( p  a ) , Ошибка! Объект не может быть создан из кодов полей редактирования.. 2 2  Воспользовавшись разложением в ряд для J (bx) , вычислить интеграл Ошибка! Объект не может быть создан из кодов полей редактирования.   ax  1  e  J (bx)  x dx  0 2 1  b    e 2a  2a     b2 4a . Проверить справедливость разложений cos( x  sin  )  J 0 ( x)  2   J 2 n ( x)  cos 2n , n 1 sin( x  sin  )  2   J 2 n 1 ( x)  sin(2n  1) . n 0  34 Для модифицированных функций I ( z )  e что I n ( z )   I n ( z ) при n  Z и I 1 ( z )  2  i  2  J (iz ) доказать, 2  sh z . z ции Pnk ( x) при  1  x  1 и n, k  Z 0 . Проверить, что полином Лежандра, записанный в виде формулы 1 dn 2 n Родрига Pn ( x)  n  n ( x  1) , является решением уравнения: 2 n! dx Пример 35. Полиномы Лежандра Pn ( x) и присоединённые функ-   (1  x 2 ) Pn  2 xPn  n(n  1)  Pn ( x)  0 . Воспользовавшись формулой Родрига, вычислить n  0, 1, 2, 3, 4. n Проверить, что Pn (1)  (1) , Pn ( x) при 1  Pn (1)  n (n  1) (1) n 1 , 2 P2 n 1 (0)  0 , Ошибка! Объект не может быть создан из кодов полей редактирования.,  P2 n (0)  0 , Pn ( x)  (1) n  Pn ( x) , Используя формулу Родрига, доказать, что Pn ( x)  1. (2n  1)   Pn ( x )dx  Pn1 ( x )  Pn1 ( x ) . Воспользовавшись производящей функцией F ( x,  )  (1  2 x   ) , 2  1 2 получить внутри круга   1 разложение:  1 2   (2n  1)  Pn ( x) n . 2 32 (1  2 x   ) n 0 Проверить, что (k ) присоединённая функция Лежандра Pnk ( x)  (1  k 2 2 x )  Pn ( x) является решением уравнения: 35 (1  x ) Pnk 2  k2     2 xPnk    n(n  1)  1  x 2    Pnk ( x)  0 .   n 2 2 x ) nk Показать, что Pnk ( x)  Pn ( x) при k  0 , Pnk ( x)  (1)  Pnk ( x) , , Pnk (1)  Pnn (1)  0 . Воспользовавшись формулой Родрига, вычислить Pnk ( x) при n, k  0, 1, 2, 3 и k  n . Pnk ( x)  0 при, Pnn ( x)  (2n  1)!!(1  Тема 8 Теория плоского электростатического поля Основные вопросы: Комплексный потенциал электростатического поля, силовая и эквипотенциальная функции. Определение по известному комплексному потенциалу напряжённости поля, работы и заряда, картины поля. Определение комплексного потенциала методом конформных отображений. Литература: 1, гл.7. §2(2); 2, п.46, 47 (3,4); 3, §39, 40; 6, Приложение к гл.4, п.5. Задание: 4, №1335, 1337, 1338, 1340, 1345, 1375. Пример 36. Исследовать структуру плоского электростатического поля по заданному комплексному потенциалу W ( z )  a  ln( z  z0 ) , где Сначала определим физический смысл размерной постоянной a . Пусть A работа, совершаемая пробным зарядом при движении по кон туру C , и Q поток вектора напряжённости E через контур C , создаваемый зарядами, расположенными внутри контура. Тогда a и z0 постоянные. A  iQ  i  C W ( z ) dz  ia   C dz  2 a , где z0  D и C  D z  z0 замкнутый контур, а D область внутри этого контура. В потенциальном поле всегда A C  0 . По теореме Гаусса Q  4 q , где q величина линейного заряда, создающего поток Q ; поэтому a  ia0 (Re a  0) и a0  2q . Таким образом, a  2iq и комплексный потенциал равен 36 W ( z )  2iq  ln( z  z0 ) . Исследуем картину поля в плоскости, перпендикулярной к наводящим линейным зарядам. Так как W ( z )  G ( x, y )  i ( x, y ) , где G ( x, y ) потенциальная функции поля, то, полагая силовая и  ( x, y ) получим y  y0 G ( x, y )  2q  2q  arctg x  x0 z  z0  rei , W  G  i  2q(  i ln r ) , и где  ( x, y )  2q ln r   q  ln ( x  x0 ) 2  ( y  y0 ) 2 . Поэтому уравнения силовых линий C1 будут G ( x, y )  C1 или y  y0  ( x  x0 )  tg лучи из точки 2q M ( x0 , y0 ) , а уравнения эквипотенциальных линий  ( x, y )  C 2 или  C2  ( x  x0 ) 2  ( y  y0 ) 2  exp  q   семейство окружностей с центром   линий взаимно ортогональны, так как ( grad G, grad )  0 . Полученная картина поля может создаваться только линейным зарядом q , расположенным в точке M ( x0 , y 0 ) ; заряд-изображение   q  находится в бесконечности. Комплексный потенциал поля W ( z ) должен иметь особенности в тех точках, где находятся заряды, создающие поле. Вектор напряжённости нашего поля равен в той же точке M ( x0 , y 0 ) . Семейства силовых и эквипотенциальных    E ( x, y )  E x ( x, y )  iE y ( x, y )   E ( x, y )  ei ( x , y )  i  W ( z )  2 q 2 q i  e при r  0 . r z Указания. Воспользовавшись принципом суперпозиции, комплексный потенциал поля, создаваемого несколькими зарядами ложенными в точках qk , в расповиде: W ( z )  2i   qk  ln( z  zk ) . Определение комплексного потенциала k zk , можно записать поля линейных зарядов над металлической поверхностью производится методом электростатических изображений. Напряжённость постоянного 37  магнитного поля H (наведенного постоянным линейным током I из   точки z0 ) равна H  i E , т.е. силовая и потенциальная функции меняются местами. Пример 37. Исследование поля диполя. Пусть линейный заряд  q расположен в начале координат z  0 , а заряд (  q ) в точке z    0 . Комплексный потенциал такой систе  мы зарядов равен W0 ( z )  2qi  ln 1   . Рассмотрим предельный  z случай для этой системы, когда   0 и q   так, что произведение 2q    p  const ; тогда получим комплексный потенциал точечного 1  p    lim 2q  ln1     i диполя W ( z )  i   с дипольным моментом 0 z z       p  0. W ( z )  G ( x, y )  i ( x, y ) , то уравнения силовых линий будут  py i G ( x, y )   p  Re  2  C1 или z x  y2 x 2   y  p 2C1    p 2C1  2 2 В поле диполя работа A C  0 и поток Q  4 (q  q)  0 . Так как семейство окружностей, которые в начале координат касаются оси абсцисс и имеют центры на оси ордиi нат. Уравнения эквипотенциальных линий  ( x, y )   p  Im  z   px  C2 или  x  p 2C2 2  y 2   p 2C2 2 это семейство ок2 2 x y ружностей, ортогональных предыдущему. При отображении функцией p W ( z )  i на плоскость W  G  i указанные семейства окружностей z переходят в семейства прямых G  C1 и   C2 ; т.е. образуют прямоугольную сетку.  Вектор напряжённости поля диполя равен: E ( x, y )  i W ( z )   p p 2i   e 2 2 z r при r  0 . 38 Аналогично исследуются поля мультиполей любого порядка. Лорановские члены разложения комплексного потенциала W ( z ) соответствуют мультипольному разложению системы зарядов, которые создают этот потенциал. Пример 38. Найти напряжённость электростатического поля над идеально проводящей незаряженной плоскостью с ортогональной к ней На рис. 8 изображено сечение рассматриваепластиной ширины a . мой системы плоскостью z , которая ортогональна к незаряженной плоскости AB1, 2 A и к пластине B1, 2 C . Рис. 8 отображает картину поля в плоскости z картину поля в плоскости  на  W , то искомая напряжённость равна Ez  EW  W ( z ) . Обычно выбирают такую новую плоскость, на которой силовые G ( x, y )  C1 и эквипотенциальные  ( x, y )  C2 линии образуют прямоугольную декартову сетку. Легко показать, что функция W ( z )  z  a отображает полуплоскость Im z  0 с разрезом 0  y  a на верхнюю полуплоскость 2 2  электростатического поля E z , наведенного некоторой системой  E проводников (зарядов), и известны W напряжённость поля в плоскости W  G  i и аналитическая функция W ( z ) , которая конформно Если в плоскости z  x  iy требуется определить напряжённость  Im w  0 (без разреза), где напряжённость поля известна EW  2 i  (поле над равномерно заряженной плоскостью); здесь   const поверхностная плотность зарядов. Очевидно, линии G  C1 и   C 2 39 образуют декартову сетку на плоскости W . Если в соотношении положить    0  const и исключить величину G , то уравнения эквипотенциальных линий в плоскости z W  G  i  z2  a2  a2    . На рис. 8 эти эквипотенполучаются в виде y ( x)   0 1  2 2   x  0  циали в плоскостях z и W изображены штриховыми линиями. 1 2   a2   Искомая напряжённость поля равна E z ( z )  2 i   1  z 2  .   Вычислим её значения в некоторых характерных  точках:  E z (ia )   напряжённость на конце острия, E z (0)  0 нет зарядов  внутри углов с вершинами B1 и B2 , E z ()  2i поле вдали от неод E B CB нородности 1 2 совпадает с полем W . 12 Тема 9 Основы асимптотических методов Основные вопросы: Асимптотические последовательности и ряды. Асимптотические методы вычисления определенных интегралов: Лемма Ватсона, методы Лапласа и стационарной фазы. Основы метода перевала. Получение ВКБЙ (WKBJ) – приближения при решении линейных дифференциальных уравнений второго порядка. Литература: 1, Прилож. 1; п. 1-3; 2, п. 76, 77; 3, § 41-46; 8, гл. 1, § 4. Задание: 4, № 776, 777, 779, 780, 782, 783, 786. Пример 39. Асимптотические оценки; асимптотические последовательности и ряды. а) Пусть lim z z  ( z)  A   для z  D и  ( z )  0 при z  z0   ; 0  ( z) тогда, если A  0 обозначают  ( z )  O  ( z )  при z  z0 , а если A  0 обозначают  ( z )  o ( z )  при z  z0 . Если A  1 , пишут  ( z ) ~  ( z ) при z  z0 (асимптотическое равенство). 40 Проверим оценки: 3 z 2  5 z 3  O( z 2 ) 5 z  sin 4 z  O( z 3 ) или  o( z 2 ) при z  0 ; (1  z )  1   z  o( z ) при   R и z  0 ; ln z ~ z  1 при z  1 ; 4 z 2  2 z 5  O( z 5 ) или  o( z ) при z   ; 6 или  o( z ) при z  0 ; e  z  o( z n ) при n  N и z   . б) Последовательность   k ( z ) называется асимптотической при z  z0 в области D , если  k 1 ( z )  o k ( z )  при z  z0 и k  N ; причем  k ( z )  0 при z  z0 . k Например,  k ( z )  ( z  z0 ) при z  z0 и  k ( z )  k при z z   образуют асимптотические последовательности. 1 в) Формальный ряд k 0  ck  k ( z )  называется асимптотическим раз ложением функции f ( z ) по асимптотической последовательности в области D и обозначается f ( z ) ~ n k 0  k ( z ) при z  z0 при  ckk ( z) z  z0 , если разность Rn ( z )  f ( z )   ckk ( z )  k 0 венное асимптотическое разложение функции f ( z ) при z  z0 , коэффициенты которого определяются по формулам  f ( z )  S n ( z )  o n ( z )  при z  z0 и n  N . Существует единст- ck  lim   f ( z )  Sk 1 ( z )   k ( z )  z  z0 Асимптотическое разложение при n   могут расходиться в каждой фиксированной точке z  z0  D , так как cn  0 при n   , но Sn ( z ) дает приближенное значение функции f ( z ) с точностью On1 ( Z )   on ( z )  при z  z0 и z  D . 41 при фиксированном n и z  z0 частичная сумма Часто члены асимптотического ряда функции f ( z ) ~ k 0  uk ( z )  при z  z0 сначала убывают до некоторого номера k  m( z0 ) , а потом начинают возрастать. Поэтому для вычисления значения функции следует брать частичную сумму Sm ( z0 ) и ошибка o m ( z0 )  будет при f ( z0 ) этом наименьшей. Значение номера k  m( z0 ) можно получить приближенно из соотношения u m 1 ( z0 ) u m ( z0 )  1 .  k ( z )   z  z0 k и т.п.) можно производить алгебраические действия, такие ряды можно почленно дифференцировать и интегрировать, если полученный ряд тоже оказывается асимптотическим. Легко проверить, что ряд Тейлора (или Лорана) функции f ( z ) в окрестности точки z  z0 является асимптотическим сходящимся рядом при z  z0 . Можно показать, что в асимптотическом разложении функции Над асимптотическими степенными рядами (когда f ( z )  e z в области по последовательности k ( z )  f ( z)  e z D   Re z  0  все коэффициенты ck  0 ; т.е. 0 0 ~ 0   ...  k ... это асимптотический нуль при z   . z z 1 zk при z Пример 40. Разложение в асимптотические ряды интегралов при большом значении параметра обычно производится методами интегрирования по частям или разложения в ряд подынтегральной функции, обобщением последнего метода является лемма Ватсона. а) Найдем асимптотическое разложение интеграла I ( )   e   x   dx при    . x Если этот интеграл проинтегрировать n  1 раз по частям, получим 42 I ( )  n 1   1!  2  2!  3  ...  (1) n  n!  n1   (1)  (n  1)!  e   x n dx k!    Rn ( )   n2 k 1 x k 0 (  )  Sn ( )  Rn ( ) . Найденное разложение интеграла I ( ) по асимптотической последовательности  k ( )  1 только в том случае, если выполняются достаточные условия оценки остаточного члена Rn ( )  I ( )  S n ( )  o n ( )  при    и n  N , т.е. для всех n должен существовать предел k при    будет асимптотическим рядом lim n  Rn ( )  0 . Сделаем оценку остаточного члена   Rn ( )  (n  1)! e     x dx dx (n  1)!  (n  1)! e   x n  2  n  2 , n2 x     1   2 n ( ) ( 1 )!    R n     O 2  при    и полученное n поэтому   разложение является асимптотическим рядом I ( )   e  x  dx k! ~  k 1 x k  0 (  ) при    . k 1 Так как отношение членов ряда ck 1 ck  k   , то полученный ряд расходится при всех значениях   0 . Чтобы обеспечить наибольшую точность при расчетах, ряд следует обрывать на члене с номером m  E ( ) , где E ( ) целая часть числа R.  1 x I (  )  f ( x )  x e dx  б) Лемма Ватсона: Если интеграл a     при  0  a  ,  и   0 0 абсолютно сходится при    , функ43 bx ция f ( x) при x  0 растет не быстрее экспоненты f ( x)  Me  xk (k ) ( M  0 и b  0 постоянные) и ее разложение f ( x)   f (0)  k! k 0 x  R   (   0 , R  a ), то сходится в некотором интервале  возможно разложение интеграла I ( ) в асимптотический ряд     k  k 1 (k )   I (  ) ~    f ( 0)     при    .    k 0 k!    1  в) Получим из леммы Ватсона при a   и тическое разложение для интеграла Лапласа     1 асимптопри    . I ( )   f (t )e 0   t dt ~ k 0  f ( k ) (0)  k 1  0  г) Найдем асимптотическое разложение интеграла I ( )   1    1 e   d  e    (1  t ) 12 t  1 2 e  t dt при    ; 1 2 здесь   1  t . Воспользуемся леммой Ватсона, считая  a   ,   ,  1 и f (t )  (1  t ) 12  1    1 k 0 k (2k  1)!! k 1 t при t  1 . (2k  2)!! После преобразований получим    3  (2k  1)!!2  I ( ) ~ e     k   2 k 1 (2k  2)!!(2 )   (3!!) 2 (5!!) 2    3 (1!!) 2     e   ...  2 3   2 4!!2 6!!(2 ) 8!!(2 )  здесь главный член разложения I ( ) ~ ченный ряд расходится при 44 3   e при    . Полу2  всех  0 так как ( 2k  3) 2 ck 1   при k   . Этот ряд следует обрывать на  ck ( 2 k  4)  2  члене с номером m  E ( ) . д) Найти асимптотическое разложение интеграла e t I ( )   2 dt при    . 0 t 1  Воспользуемся леммой Ватсона, считая a   ,     1 и f (t )  (1  t )   (t 2 ) K при t  1 . После преобразований получим 2 1 K 0  I ( ) ~ Этот k 0  (1)k  (2k )!  2 k 1  1   2!  3  4!  5  6!  7  ... при    .   0, так как ряд расходящийся при всех ck 1 (2k  1)(2k  2)    при k   . Поэтому ряд следует обck 2 1   m E  . рывать на члене с номером 2  Пример 41. Асимптотическое разложение функций-оригиналов f (t ) при t  0 . Часто оригиналы f (t ) получаются в виде неберущихся интегралов, для которых можно найти асимптотическое разложение по переменной t . Например, оригинал в примере № 28 может быть разложен в асимптотический ряд при t  0  0 . Рассмотрим соответствие F ( p)  1 a e p z 2 p  a  f (t )  1     2 t     при a  0 ; здесь  ( z )  2  0 Для функции  ( z ) найдем асимптотическое разложение при z   a   при t  0  0) . (z  2 t 45   e  d интеграл вероятностей (см. пример №33). Сначала запишем ( z)  1  Сделаем в 2    e  d  1  z  2 2   Erfc ( z ) . замену переменной интеграле Erfc ( z )   z 1  x  0 и получим новый интеграл  1  z 2   z 2 x2 Erfc ( z )  ze   e  (1  x) 2 dx , 2 0 1 к которому легко применить лемму Ватсона. Действительно, теперь   1 ,   2 , f ( x )  1   зований получим (2k  1)!!  ( x) k , x  1 ; после преобраk 1 ( 2k )!!  ( z) ~ 1  1 z  e z2  ( 2 k  1)!!    1   ( 2 z 2 ) k    k 1   1 1 z  e z2 при z   . Этот ряд расходится при всех z  0 , так как отношение его чле- ck 1 2k  1    при k   . При расчетах ряд следует обрынов 2 ck 2z вать на члене с номером m  E ( z ) . Главный член разложения для функции – оригинала f (t ) будет 2  a  2 t f (t )  1    e ~ 2 t  a   a2 4t при t  0  0 . Разложения оригинала f (t ) при t   непосредственно получаa  0. ется из тейлоровского разложения функции  ( z ) при z  2 t 46 Пример 42. Асимптотическое вычисление интегралов от функции с острым максимумом по методу Лапласа. а) Теорема: Имеет место асимптотическое представление интеграла I ( )   g ( x)  ef ( x ) dx  a b 2   1   g ( x0 )  ef ( x0 )  1  O   f ( x0 )     при x  x0  ( a   , b   ) , где 0    b  a и может быть разложена в ряд Тейлора в окрестности x  x0 ; функция g ( x ) непрерывна и изменяется медленнее, чем exp  f ( x )  . б) Для гамма–функции  ( z  1) при z  1 получим асимптотическую формулу Стирлинга при    (здесь    a  b   и Im f ( x )  0 ), если этот интеграл абсолютно сходится при   1 , функция f  x  имеет острый максимум ( z  1)   e  t dt  z 0  t z z 1   e z (ln x  x )  dx ; 0  Здесь t  zx , g ( x )  1 и f ( x )  ln x  x . Максимум функции f ( x ) достигается при x0  1 , при этом f (1)  f (1)  1 . Воспользовавшись формулой асимптотического представления интеграла, получим z ( z  1)  2 z    e z   1   1  O  при z  1 .  z   Формула Стирлинга дает результаты с хорошей точностью уже при z  4. 4  ( 4  1 )   ( 5 )  4 !  24  8      23,50617 ; здесь Например, e ошибка 2,05763%  2,06% . в) Получить асимптотическое представление интеграла 4 I ( )   x 0   x dx   e   x ln x dx при   1 ; 0 47  здесь g ( x)  1 и f ( x )   x  ln x . Максимум функции f ( x ) достигается 1 1 1 при x0  1 e , при этом f    и f    e . Асимптотическое предe e e ставление интеграла будет 2 e I ( )  e e    1   1  O  при   1 .     Пример 43. Асимптотическое вычисление интегралов от быстро осциллирующих функций по методу стационарной фазы (методу Стокса). а). Теорема. Имеет место асимптотическое представление интеграла I ( )   g ( x)  eif ( x ) dx  a b  2   1   g ( x0 )  ei f ( x0 )    1  O  ,  f ( x0 )       4  sign f ( x0 ) , при    (здесь    a  b   и Imf ( x )  0 ), если этот интеграл f ( x ) имеет экстремум при сходится при   1 . Функция x  x0  ( a   , b   ) , где 0    b  a , и может быть разложена в ряд Тейлора в окрестности x  x0 ; функция g ( x ) непрерывна и изменяется медленнее, чем expi  f ( x )  . Если пределы интегрирования a и b конечны и около них крайние полуволны функции Re g ( x )  e не погашаются, то иногда приходится учитывать значение интеграла I ( ) вблизи этих концов. б) Найти асимптотическое представление для функции Бесселя  1 J n ( z )    cos( z  sin t  nt )  dt  Re  ei ( zsin t  nt )  dt 2  0   i z f ( t )  1  при z  1 48 и nZ ; здесь g (t )  e  in t и f (t )  sin t . В интервале t   функция f (t ) имеет две стационарные точки t    i  2 ; при этом  n         g    e 2 , f       f     1 ,       . 4  2  2  2  2 Сумма вкладов обеих стационарных точек дает асимптотическое представление функции Бесселя J n ( z)  2     1   cos z  ( 2n  1)   O 3 2  при z  1 . 4 z   z    J n ( z )  ( 1) n J n ( z ) n n Далее, подставив в формулу для определения функции Неймана   Nn ( z)  при n  Z интегральные представления для функции Неймана в виде J  n (z) , получим асимптотику Nn ( z)  ем 2   1     sin  z  (2n  1)   O 3 2  при z  1 . z 4 z     Теперь легко находятся и асимптотики для функции Ханкеля, име i  z  ( 2 n 1)  2  1  (1, 2 )  H n ( z)  J n ( z)  i  Nn ( z)  e  4  O 3 2  z z  при z  1 . Полученные выше асимптотические выражения верны и при комплексных значениях аргумента z или индекса n , если z  1    z  n 2 0 .  в) Обобщением методов Лапласа и Стокса является метод перевала (седловой точки). Рассмотрим контурный интеграл I ( )   f ( z )  e L  ( z ) dz , где па- раметр   0 , а контур интегрирования L проводится в области анали49 тичности однозначных функций  ( z ) и f ( z ) при z  x  iy . Величина f ( z ) должна меняться медленнее экспоненты и не иметь особенностей на контуре L (иначе добавляются вычеты). Если при больших значениях   1 величина подынтегрального выражения быстро растет с увеличением функции u  Re  ( z )  0 и значение v  Im  ( z )  const , то это даст основной вклад в асимптотическое значение интеграла I ( ) . Так как  ( z )  u ( x, y )  i  v ( x, y ) , где функции u ( x, y ) и v ( x, y ) являются гармоническими и не имеют абсолютных экстремумов, то значение z  z0  x0  iy0 – корень уравнения  ( z0 )  0 является точкой минимакса (седловой точкой перевала на поверхности  ( z) ). Контур интегрирования L можно деформировать так, чтобы и далее вдоль направления наискорейшего спуска; он проходил именно через точку перевала z  z0 – частный экстремум функции u ( x, y ) т.е. направления v  Im  ( z )  const , так как (grad u , grad v )  0 . Вдоль такого контура окажется Re  ( z )   при z   . В окрестности седловой точки z  z0 сделаем разложение в ряд Тейлора L  ( z )   ( z0 )   ( z0 )  ( z  z0 ) 2  ... Здесь производные  ( z0 )  0 и  ( z0 )  0 , так как вдоль контура интегрирования L мнимая часть постоянна и мы движемся вниз от седловой точки. Сделаем замену переменной интегрирования 1 2  ( z )   ( z0 )   ( z0 )  ( z  z0 ) 2   t 2 2 , где z  z0  1 2  e i и   const . Тогда  ( z0 ) I ( )  f ( z 0 )  e где     e t 2 2  dz dt , dt dz dt 1  ei dz dt    ( z0 ) 50 . Теперь окончательный результат примет вид I ( )  f ( z 0 )  e   ( z0 )     ( z0 )  1  t 2 2 2 e dt  ei    2  f ( z0 )e ( z0 )  ei    ( z0 ) при   1 . Здесь интегрирование проводилось в бесконечных пределах, так как подынтегральная функция практически равна нулю при достаточно больших t и использован интеграл ошибок Гаусса. Угол  – аргумент контура, проходящего через седловую точку z  z0 . Это направление выбирается так, чтобы величина  была постоянной и значение Re  ( z ) – наибольшим. Важно отметить, что если экспоненциальный множитель подынтегральной функции записать в форме exp   ( z )  , то в методе наискорейшего спуска Лапласа путь интегрирования следует проводить вдоль кривой Im  ( z )  const , а в методе стационарной фазы Стокса – вдоль кривой Re  ( z )  const . Очевидно, эти направления взаимно ортогональны; причем контур в методе Лапласа проводят через перевал вниз в долины поверхности  ( z ) , а в методе Стокса – наоборот вверх вдоль вершин горной цепи этой поверхности. Пример 44. Асимптотический метод решения линейных дифференциальных уравнений второго порядка – метод ВКБЙ (Wentzel – Kramers – Brillouin –Jeffreys). Линейное дифференциальное уравнение второго порядка ~ y   p ( x )  ~ y   q( x)  ~ y ( x )  0 , где p ( x ) и q ( x) – заданные непрерывные функции, с помощью замены искомой функции  1  ~ y ( x)  y ( x)  exp   p ( x)dx  приводится к стандартному виду   2 y   I ( x )  y ( x )  0 , где инвариант уравнения I ( x)  q( x)  1 2 1 p ( x)  p ( x) 4 2 присущ только рассматриваемому уравнению и отличает его от всех других. 51 Пусть некоторое линейное дифференциальное уравнение удалось привести к виду y  2  U ( x )  y ( x )  0 С помощью замены искомой функции при    . 2 y ( x )  y ( x ) V ( x ) ( y   yV , y   yV  yV  ) получим нелинейное уравнение первого порядка типа Риккати V   V 2 ( x )  2  U ( x ) при    , решение которого будем искать методом разложения по малому параметру 1   1 в виде ряда Лорана 1 V ( x )  с1    с0  с1    с2  1  2  ... , где сn  сn ( x ) . Здесь разложение функции V ( x ) начинается с первой отрицательной степени параметра  , чтобы обе части уравнения Риккати были эквивалентны при    . Легко находятся значения  1 ( x )    с0  ( x )  с1  ( x)  V ( x)  с 2 1   ( x)   с2 1 2 2  ... , 1 1   V ( x )   c1    c0  c1  c2  2  ...       2 2  c    c0  c1 1 2 2 1 2  c2  1 2 1 4  ...  1  2c1c0    2c1c1  2c1c2  2 с0с1 1   2c1c3  1  2  ...    2 с0 с2 1  2  ...  2 с1с2  3  ... Приравняв в уравнении Риккати коэффициенты при одинаковых степенях  , получим 52 2 : с12  U ( x ),  1  2с1с0  0, 1 : с   с0 2  2с1с1  0, 0 : с0   2с1с2  2с0с1  0,  : с1   с12  2с0с2   2 : с2  2с1с3  0, 3 : 1 с1   U ( x )  1 2с1 с0   с   с0 с1   с0  2  2с 1   2с0с1  2с1 с2   с1   с12  2с0 с2 2с1 с3   с2   и т.д. В общем случае рекуррентная формула для коэффициентов сn ( x ) будет n      сm  сn  m  2с1 сn 1 ( x)    сn m 0   при n  0,1,2,... По этой рекуррентной формуле находятся все значения коэффицентов разложения, например: с1   U ( x ) , с0   U ( x ) 4U ( x ) , и т.д. с1   ( 4UU   5U 2 ) (32U 2 U ) В тех точках x  x0 , где U ( x0 )  0 , все выражения для коэффициентов теряют смысл и асимптотическое разложение невозможно. В общем случае нельзя сказать сходится ли ряд функции V ( x ) при n   ; можно только утверждать, что это формальное разложение будет асимптотическим при    . Вернемся к функции y ( x) ; из V  y  y получим y  ( x)  e  V ( x ) dx   1 1    exp    U dx  ln U   n   сn ( x)dx   4 n 1    U  1 4 e    U dx  1   exp  n   сn ( x ) dx  . n 1   53 (*) Наличие множителя x  x0 , где функция U ( x ) меняет знак, должен меняться и тип асимптотики – убывание, возрастание, колеблющееся поведение. При больших значениях  последний в формуле (*) экспоненциальный множитель близок к единице и не оказывает большого влияния на поведение функции y ( x ) . Постоянные множители везде опущены. Если получается функция U ( x,  ) , то решение тоже можно найти, использовав раз1 ложение U ( x,  )  U 0 ( x )  U 1 ( x )  ... Решения y ( x ) являются анали- exp     U ( x )  dx показывает, что в точках    тическими функциями в некоторой области, их можно неоднократно дифференцировать. Наличие в решении радикалов показывает, что нули и полюсы функции U ( x ) могут являться точками ветвления (точками поворота). Проведем более подробную оценку последнего экспоненциального множителя функции y ( x ) . Для функции V ( x ) имеем U  1 U 1   n сn ( x )    U   O  V ( x)   U  4U n 1  4U   при    . Последняя оценка равномерна по x , если сn ( x ) будут ограничен- ными функциями. Поэтому переходя к y  exp  U ( x ) dx , получим два линейно независимых решения, соответствующие двум разным знакам (+ и –) перед корнем U ; т.е. найдем   y ( x)  U   1 4  e  U dx 1  exp O    при    . 1 1 Экспонента с оценкой очевидно равна exp O   1  O  при        , что позволяет записать главный член асимптотики решения уравнения в окончательном виде y  ( x)  U  1 4 e     U dx   1   1  O        54 при    .  Теперь подробно выпишем асимптотику решения y ( x ) с точно1 стью до 2 ; имеем  V ( x)   U  U 1  ( 4UU   5U  2 ) 32U 2 U  O  4U   при    . После возвращения к функции y ( x ) получим  1 4 y ( x)  U  1 4    4UU   5U  2  exp       U  322U 2 U       1   1 dx O    2            U e   U dx  1 4UU   5U  2  1  dx  O 2   1    32U U      при    . Для детальной оценки порядков членов разложения удобно рассмотреть выражения с U ( x ) ; например, неравенство с1    с0 U  5 U   2 , т.е. достаточны оценки сведется к выражению 32 U U 4 U U  U U и U   . Аналогично находятся оценки для последующих U32 членов разложения. а) Применим полученные результаты к нахождению асимптотики решения уравнения Эйри d2y  t  y (t )  0 , dt 2 t 0; U (t )  t и t   . 23 С помощью замены t  x    0 получим стандартную форму уравнения d2y 2    x y ( x)  0 2 dx при x и   0 ,    . Запишем асимптотическое решение этого уравнения и вернемся к исходной переменной t , тогда 55 3   32    1 5 1  2   1  O y ( x)  4  exp   2 x  x     2 2   3 16   x          1  3t  4 e t 2 32  5   1 t  3   32    1   1  O  2    t   при t   . Здесь опущены произвольные постоянные множители. Если в уравнение Эйри принять t   , получим y (t )   e i  4 4 t e 2 3 2  5 3   i t  1 t  3  32    1    1  O 2   при t   .   t     Фазовый множитель e 4 здесь возникает за счет обхода точки ветвления t  0 ; при этом изменяется и тип решения с монотонного на колеблющийся. б) Легко проверить, что при    главный член асимптотического решения уравнения (типа уравнения Эрмита) i  y   ( a  2 x 2 )  y ( x )  0 , 2 a  const  0 , если выполняется условие ( x )  a , будет равен 1  2  x2   1   1  O    y (t )  e x      1 при    . в) Асимптотическое решение для уравнения вырожденной гипергеометрической функции xy   (b  x ) y   ay ( x )  0 , a и b  const  0 может быть представлено для одного из решений в виде y ( x )  x a b  e x при x   ; a y ( x)  x при x   . 56 Задания для самостоятельного решения Задания содержат основные типы примеров по теории аналитических функций (по теории функций комплексного переменного) и операционного исчисления. Образцы решения всех примеров приведены в первых шести разделах пособия. Каждый студент должен представить решение всех 16 примеров; номера пунктов в каждом примере соответствуют номеру, под которым фамилия студента стоит в списке группы. Пример 1. Вычислить все значения корня 1. z  1, n  3. 2. z  1, n  3. 3. z  i, n  3. 4. z  i, n  3. 5. z  1, n  4. 6. z  1, n  4. 7. z  i, n  4. 8. z  i, n  4. 9. z  i, n  2. 10. z  i, n  2. 11. z  1  i, n  2. 12. z  1  i, n  2. 13. z  1  i, n  2. 14. z  1  i, n  2. 15. z  1  i, n  3. z , если 16. z  1  i, n  3. n 17. z  1  i, n  3. 18. z  1  i, n  3. 19. z  1  i, n  4. 20. z  1  i, n  4. 21. z  1  i, n  4. 22. z  1  i, n  4. 23. z  1, n  5. 24. z  1, n  5. 25. z  i, n  5. 26. z  1  i, n  5. 27. z  1  i, n  5. 28. z  1  i, n  5. 29. z  1  i, n  5. 30. z  1, n  6. Пример 2. Выяснить геометрический смысл указанных соотношений. 1. 2. 3. Re(iz )  0. Im(iz )  1. Re( z  i )  1 / 2. 57 4. 5. 6. Im( z  i )  2. Re(2 z  3)  1. Im(iz  3i )  4. 7. 8. 9. Re(iz  3)  0. Im( z  2i )  2. 19. z  Im( z  i ). 20. arg( z  i )   / 3. 21. arg( z  i )   / 6. 22. arg( z  2i )  0. 23. arg( z  1)   . 24. arg( z  2i )   / 2. 25. arg( z  1)   / 2. 26. 0  arg( z  i )   / 2. 27. 0  iz  2  1 28. 1  z  3i  2. 29. 1  iz  2  2. 30. 1  i  z  z  i . Re(3  iz )  1. 10. Im(2  iz )  1. 11. z  3  iz . 12. z  i  z  1. 13. z  i  1  z . 14. z  i  z  i . 15. z  2i  iz . 16. z  Re z  1. 17. z  Im z  1. 18. z  Re( z  i ). Пример 3. Первоначальное значение Arg f ( z ) при z  a  0 принято равным нулю. Точка z делает один полный оборот против часовой стрелки по окружности z  z0  a и возвращается в точку z  a . Считая, что Arg f ( z ) изменяется непрерывно, найти значение Arg f ( z ) после такого оборота, если: 1. 2. 3. 4. 5. 6. 7. f ( z )  3 z  2 , a  3, z0  0. f ( z )  z  2 , a  1, f ( z )  3 z  2 , a  3, z 0  0. z0  2. f ( z )  4 z  2 , a  1, z0  2. f ( z )  z  2 , a  1, z0  2. f ( z )  5 z  2 , a  3, z0  2i. f ( z )  4 z  2 , a  3, z0  2i. 58 8. 9. f ( z )  3 z  2 , a  3, z0  2i. f ( z )  z 2  2 , a  3, z0  2i. 3 3 10. f ( z )  z 2  2 , a  3, z0  2i. z 0  0. z 0  0. 11. f ( z )  z 2  1, a  2, 12. f ( z )  z 2  2 z  3, a  2, 13. f ( z )  14. f ( z )  3 3 z 2  2 z  3, a  2, z0  0. z 2  2 z  3, a  2, z 0  0. 15. f ( z )  z 2  2 z  3, a  2, 16. f ( z )  17. f ( z )  18. f ( z )  3 19. f ( z )  3 z0  0. z 1 , a  2, z2 z0  0. z 1 , a  2, z0  0. z 1 z 1 , a  2, z 1 z 1 , a  2, z 1 z0  0. z0  0. 20. f ( z )  2 Ln z , a  2, z0  0. z0  0. 21. f ( z )  3Ln (2i / z ), a  2, 22. f ( z )  Ln z  Ln( z  1), a  2, 23. f ( z )  Ln z  Ln( z  1), a  2, 24. f ( z )  Ln z  3Ln( z  i ), a  2, 25. f ( z )  2 Ln (1  z ), a  2, 26. f ( z )  4 Ln (iz ), a  2, z 0  0. z 0  0. z0  0. z0  0. z 0  0. 59 27. f ( z )  3Ln 3i , a  2, z0  0. z 1 z 0  0. 1 i , a  3, z0  0. z z 0  0. 28. f ( z )  Ln(i z )  Ln ( z  i ), a  3, 29. f ( z )  3Ln 30. f ( z )  3Lnz  2 Ln ( z  1), a  3, Пример 4. Найти все корни уравнения. 1. sin z  2. 2. sin z  2. 3. sin z  2i. 4. sin z  2i. 5. cos z  i. 6. sin z  i. 7. sin z  1 / 2. 8. cos z  1 / 2. 9. sh z  i. 10. sh z  i / 2. 11. cos z  2. 12. cos z  2i. 13. sin z  i. 14. cos z  i. 15. ch z  1. 16. Ошибка! Объект не может быть создан из кодов полей редактирования. 17. ch z  1 / 2. 18. ch z  i. 19. tg z  2i. 20. ch z  i. 21. ctg z  2i. 22. th z  2i. 23. cth z  2i. 24. th z  2. 25. cth z  1 / 2. 26. sec z  3i. 27. sch z  3i. 28. cosec z  3i. 29. csch z  3i. 30. ch z  sh z  2. Пример 5. Найти все значения степени. 1. (i ) . 60 2i 2. i 2i . 3. (1)3i . 4. (i )3i . 5. (i ) 2i . 6. (1) 3i . 7. (1) 2i . 8. 1  2i 19. (1  i ) 2i . 20. (1  i ) 2i . 21. (1  i ) 2i . 1 i  22.   .  2 1  i  23.   .  2  1  i  24.   .  2   1 i  25.   .  2  26. i 2 i i i i . 9. (1  i ) 2i . 10. (1  i ) 2i . 11. (1  i ) 2i . 12. (1  i ) 2i . 13. 1 2. 14. 1  2 . 2 . 2 27. i  . 2 15. (1) 16. i . 17. 1 . 3i 3i . 28. (i ) 29. (i )  30. i i 2 . 2 . . 18. (1  i ) 2i . Пример 6. Определить, существует ли область D , в которой функция f(z) аналитическая (удовлетворяет условиям Коши-Римана); найти производную f ( z ) при z  D , если: 1. 2. 3. f ( z )  e z  2 z. f (z)  e z 4. f ( z )  sh z. f ( z )  cos z. f ( z )  ch z. f ( z )  sin 2 z.  2 z. 5. 6. 7. 61 f ( z )  sin z. 8. 9. f ( z )  sh 2 z. f ( z )  cos 2 z. 2 19. f ( z )  2 Ln z. 20. f ( z )  2 Ln z  z. 21. f ( z )  z 2  Re z. 22. f ( z )  z 2  Im z. 23. f ( z )  z  3z . 24. f ( z )  z  i z . 25. f ( z )  z / z. 26. f ( z )  z 2  2i z. 27. f ( z )  z 2  2i z. 28. f ( z )  z 2  2 z . 29. f ( z )  z 2  2 z. 30. f ( z )  z z . 10. f ( z )  ch z. 11. f ( z )  sin z . 12. f ( z )  sh z 2 . 13. f ( z )  cos z . 14. f ( z )  ch z 2 . 2 2 i 15. f ( z )  z  . z i 16. f ( z )  z  . z 2 i 17. f ( z )  z 3  . z 18. f ( z )  3i  1 . z2 Пример 7. Какая часть плоскости растягивается (сжимается) при конформном отображении, которое осуществляет функция f ( z ) . Найти также коэффициент растяжения к и угол поворота Ө в точке z0 , если: z 1 1. f ( z )  , z0  i. z i z 1 2. f ( z )  , z0  i. z 1 z i , 3. f ( z )  z  2i zi , 4. f ( z )  z  2i z0  1. z0  1. 5. f ( z )  6. f ( z )  z 1 , z 1 z0  0. z i , z0  0. z i z 0  1 . z0  i. 7. f ( z )  z 2  2i, 8. f ( z )  z 2  2i, 9. f ( z )  z 3  3i, z0  1. 62 10. f ( z )  z 3  3i, 11. f ( z )  2 / z , 12. f ( z )  2i / z , 13. f ( z )  ln z , 14. f ( z )  ln z , 15. f ( z )  e  z , 16. f ( z )  ei z , 17. f ( z )  e  2 z , 18. f ( z )  e 2 z , z0  i. 21. f ( z )  ln( z  1), 22. f ( z )  ln( z  i ), 23. f ( z )  2 z 3 , 24. f ( z )  3 z 3 , 25. f ( z )  2 z 3 , 26. f ( z )  3 z 3 , 27. f ( z )  ei z , 28. f ( z )  ei z , 29. f ( z )  e 3i z , 30. f ( z )  e 3i z , z0  i. z0  1. z0  1. z0  i. z0  2. z0  2i. z0  i. z0  i. z0  1  i. z0  1  i. z0  1  i. z0  1  i. z0  i. z0  2i. z0  2i. z0  2i. z0  i. z0  i. z0  i. 19. f ( z )  z 2  2 z , 20. f ( z )  z 2  2 z , z0  i. Пример 8. Конформно отобразить на верхнюю полуплоскость Im w  0 круговую лунку (двуугольник) вида: 1. z  1, 2. z  1, 3. z  1, 4. z  1, 5. z  1, 6. z  1, 7. z  1, 8. z  1, 9. z  1, z  1  1. z  1  1. z  i  1. z  i  1. z  1  1. z  1  1. z  i  1. z  i  1. z  1  1. 63 10. z  1, 11. z  1, 12. z  1, 13. z  1, 14. z  1, 15. z  1, 16. z  1, z  1  1. z  i  1. z  i  1. z  1  1. z  1  1. z  i  1. z  i  1. 17. z  1  1, Im z  0. 18. z  1  1, Im z  0. 19. z  i  1, Im z  0. 20. z  i  1, Im z  0. 21. z  1, Re z  0. 22. z  1, Re z  0. 23. z  1, Im z  0. 24. z  1, Im z  0. 25. z  1, Re z  0. 26. z  1, Im z  0. 27. z  2, 28. z  2, 29. z  2, 30. z  2, z  2  2. z  2  2. z  2  2. z  2  2. Пример 9. Конформно отобразить на верхнюю полуплоскость Im w  0 заданную односвязную область D c разрезом L, где: 1. D : Im z  0; L  0, 2i  2. D : Im z  0; L  0,3i . 3. D : Re z  0; L  2, . 4. D : Re z  0; L   ,2. 5. D  C ; L   1,1 . 6. D  C ; 7. D  C ; 16. 17. 18. 19. 20. 21. 22. D : z  1; L  1, 2. D : z  1; L   1,  . D : z  1; L    2,1 . D : z  1; L   i, 2i . D : z  1; L   1,  . D : z  1; L   i,i . D : z  1; Im z  0; L   0, i / 2 . D : z  1; Im z  0; L   i / 2 , i . D   0  x  1,    y   . D      x  , 0  y  1, . D    1  x  1,    y   . D   x  0, 0  y  1 . D   x  0,  1  y  0 . D   0  x  1, y  0 . L   i,i  . L   1  i, 1  i  . 8. D  C ; L   1  i,1  i . 9. D : Im z  0; L  z  e i , 0     / 2 . 10. 11. 12. 13.   23. 24. 25. 26. 27. 28. 29. 64 D : Im z  0; L  z  e i ,  / 2       D : z  1; L   0.5, 1. D : z  1; L   1,  0.5. D : z  1; L  0, 1. 14. 15. D : z  1; L   i / 2, i . D : z  1; L  0,  i  . 30. D : z  1  1, z  2  2, Imz>0. Пример 10. Используя интегральную формулу Коши и формулы для производных от нее, вычислить интегралы вида  f ( z )dz , если: z R ez 1. f ( z )  2 , R  3. z 4 2. f ( z )  3. f ( z )  11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. f ( z)  f ( z)  1 sin z , R  1. 2 z 1 sin z , R  1. z3 e z z2 1 z , R  2. e , R  2. z ( z  1) 1 f ( z )  sin 2 z , R  1. z f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  sin z , R  2. z ( z  1) sin z , R  2. z ( z  1) sin z , R  2. z 2 ( z  1) sin z , R  2. z ( z  i) sin z , R  2. z ( z  1) 2 sin z , R  2. z 2 ( z  i) sin z , R  2. 2 z ( z  i) ch z , R  2. z ( z  1) ez , R  2. 4. f ( z )  2 z ( z  1) ez , R  2. 5. f ( z )  3 z ( z  1) ez , R  2. 6. f ( z )  2 z ( z  1) ez , R  2. 7. f ( z )  z ( z  1)3 ez , R  2. 8. f ( z )  2 ( z  1) 2 ez , R  2. 9. f ( z )  2 2 ( z  1) 10. f ( z)  sin z , R  1. z 65 22. 23. 24. 25. 26. f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  ch z , R  2. z 2 ( z  i) ch z , R  2. 2 z ( z  i) ch z , R  2. 2 z sh z , R  2. z2 z , R  2. 2 2 ( z  1) 27. 28. 29. 30. f ( z)  f ( z)  f ( z)  f ( z)  z , R  2. ( z 2  1) 2 z , R  2. 2 z 1 z , R  2. z2 1 1 , R  2. 2 3 ( z  1) Пример 11. Определить нули и особые точки функции f ( z ) , если: 1. 2. 3. 4. 5. 6. 7. 8. 9. f ( z )  z  cos z . 10. f ( z )  11. f ( z )  12. f ( z )  13. f ( z )  1 f ( z )   cos z . z f ( z)  1 z2  cos z . z . 2 z 1 z 1 z . f ( z )  z  sin z . z . z 1 1 f ( z )   sin z . z f ( z)  f ( z)  sin z . z z sin z . 2 z z 1 . 2  sin z 1 14. f ( z )   tg 2 z. z 15. f ( z )  sin z . 16. f ( z )  f ( z)  z3 z . f ( z)  z . z 1 66 3 1 z e . 2 z 17. f ( z )  sin 18. f ( z )  z e 1 . z 1 1 z.  19. f ( z )  20. f ( z )  21. f ( z )  22. 1 2 z z e . 1 25. f ( z )  cosec z  . z 1 . 1  z 1 1 . 1 z  i   26. f ( z )   z  i   th z. 2    27. f ( z )   z    tg z. 2  28. f ( z )  z  cth z. 29. f ( z )  z  ctg z. 30. f ( z )  z  ctg 2 z. f ( z )  ctg 1 . z 23. f ( z )  e tg z . 24. f ( z )  z  cosec z . Пример 12. Функцию f ( z ) разложить в степенной ряд в окрестности точки z0 (или внутри кольца), определить область сходимости полученного рядя, если: 1 1 8. f ( z )  2 , z0  . , z0  . 1. f ( z )  z  1 2 ( z  1) 1 , z0  0. 9. f ( z )  2 1 2 , z  0 . 2. f ( z )  ( z  1) 0 ( z  1) 2 1 , z0  i. 10. f ( z )  2 1 2 3. f ( z )  , z 0  0. ( z  1) z ( z  1) 1 1 , z0  . 11. f ( z )  2 2 4. f ( z )  , z0  1. ( z  1 ) z ( z  1) 1 5. f ( z )  , z0  . z ( z  1) 6. f ( z )  12. f ( z )  z ( z  1) , 13. f ( z )  3 z 2 ( z  1) , 14. f ( z )  3 z ( z 2  1) , 15. f ( z )  ln z0  . z0  . z0  . 1 , z0  0. 2 z 1 1 , z0  i. 7. f ( z )  2 z 1 z 1 , z0  . z i 67 16. f ( z )  17. f ( z )  18. f ( z )  sin z , z0  0. z z sin z , z0  . z z cos z , z0  . z 25. f ( z )  4 z 3 ( z  1) , 26. z0  . f ( z)  0  z  2. z , ( z  1)( z  2) z 1 , z ( z  2) z2 , z ( z  1) 27. 28. f ( z)  f ( z)  0  z  2. 0  z  2. 19. f ( z )  3 ( z  1)( z  2)( z  3) , z0  . 20. f ( z )  z 4 e 1 / z , 21. f ( z )  z 3e 1 / z , 22. f ( z )  23. f ( z )  z 0  0. z0  . z 2  4z 29. f ( z )  сos , z 0  2. 2 ( z  2) z 2  4z 30. f ( z )  sin , z0  2. ( z  2) 2 1 , 0  z  1. z (1  z ) 1 , 0  z  1. 2 z (1  z ) z0  1. 24. f ( z )  z , z R Пример 13. Вычислить с помощью вычетов контурный интеграл  f ( z )dz , где: 1. f ( z )  z 3e 1 / z , R  1. 2. f ( z )  ze 1 / z , R  1. 3. f ( z )  cosec z , R   / 2. 4. f ( z )  sec z , R   . 5. f ( z )  tg z , R   . 6. f ( z )  ctg z , R  1. 7. f ( z )  cth z , R  1. 2 1 8. f ( z )  sin , R  1. z 1 9. f ( z )  z  cos , R  1. z 10. 11. 12. 68 1 f ( z )  sin 2 , R  1. z f ( z)  f ( z)  z z 1 2 , R  2. z , R  2. ( z  3)( z 5  1) 13. 14. 15. 16. 17. 18. 19. 20. 21. f ( z)  f ( z)  f ( z)  1 , R  1. z 6 ( z  2) 1 , R  1. 5 2 z ( z  2) z , R  2. 6 z 1 22. 23. 24. 25. 26. 27. 28. 29. 30. f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  f ( z)  1 , R  2. ( z 2  1) 2 1 , R  2. 2 2 ( z  1) sin z , R  4. 2 (z   ) cos z , R  3. ( z   / 2) 2 tg z , R  1. z2 th z , R  1. 2 z сtg z , R  1. z cth z , R  1. z 1  f ( z )  cos  cos , R  1. z   1 f ( z )  z  sin  sin , R  1. z  1  f ( z )  cos  z 3 cos , R  1. z  1  f ( z )  sin  z 3 sin , R  1. z  f ( z)  1 , R  2. (z  z3 ) 1 f ( z)  , R  2. 4 (z  z ) 1 f ( z )  z  tg , R  1. z Пример 14. Вычислить с помощью вычетов определенный инте- грал  f ( x)dx , где: a b 1. f ( x)  2. f ( x)  3. f ( x)  1 ; 2  сos x 1 ; 2  sin x 1 ; 2  сos x a  0, b  2 . a  0, b  2 . a  0, b  2 . 69 4. f ( x)  1 ; 2  sin x a  0, b  2 . a  0, b  2 . 5. f ( z )  e cos x  cos(2 x  sin x); 6. f ( x)  e cos x  cos(3 x  sin x); 7. f ( x)  tg( x  i / 3); 8. f ( x)  сtg( x  i / 3); 9. f ( x)  a  0, b  2 . a  0, b  5 . a  0, b  5 . 1 ; ( x 2  4) 2 a  0, b  . 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. x2 f ( x)  2 ; ( x  1) 2 f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  x cos x ; 2 x  2 x  10 x sin x ; x 2  4 x  20 sin x ; 2 x  2 x  10 x cos x ; x 2  4 x  20 cos x ; 2 x  5x  6 x sin x ; x2 1 cos x ; 2 x 1 sin x ; x sin 2 x ; x 3 a  0, b  .  a  b  .  a  b  .  a  b  .  a  b  .  a  b  . a  0, b  . a  0, b  . a  0, b  .  a  b  . 70 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  f ( x)  cos 2 x ; x 3  a  b  .  a  b  . a  0, b  . a  0, b  . a  0, b  . sin x ; ( x  1)( x 2  4) sin x ; 2 x ( x  4) cos x ; 2 2 ( x  1) 1  (cos 2 x  cos 4 x); x2 cos x ; 2 x 1 x sin x x2 1 ; a  0, b  . a  0, b  . 1 ; 2 x  ( x  1) 1 ; 2 x  ( x  1) ln x x2  1 ln 2 x x 1 2 a  0, b  . a  0, b  . f ( x)  3 f ( x)  f ( x)  ; a  0, b  . a  0, b  . ; Пример 15. Используя методы операционного исчисления, решить предложенные задачи Коши для обыкновенных дифференциальных уравнений (для систем таких уравнений), или вычислить определенные интегралы, или решить интегральные уравнения. Все параметры a,b>0 . 1. y  2 y  2 y ( x)  1; 2. y   6 y   9 y ( x)  2; y (0)  1, y(0)  2. y (0)  1, y(0)  3. 71 3. y  y  6 y ( x)  2; 4.  e ax sin bx dx. 0 y (0)  1, y(0)  0.  5. y  9 y ( x)  2  x; 6. y  4 y ( x)  4 x; 7. y  4 y ( x)  2 cos 2 x; 8.   sin ax 0 y (0)  0, y(0)  1. y (0)  1, y(0)  0. y (0)  0, y(0)  4. x dx. y (0)  0, y(0)  1. y (0)  y(0)  0. 9. y  y ( x)  x 3  6 x; 10. 11. 12. 13. 14. 15. 16. y  y ( x)  cos x  sin 2 x; y  4 y ( x)  4e 2 x ;  0 y (0)  1, y(0)  0. e ax sin x  cos bx dx. y (0)  y(0)  0. y  2 y  y ( x)  x  2; y  9 y ( x)  e x ;  0 y (0)  1, y(0)  0. y (0)  y(0)  0. y  4 y  4 y ( x)  2 sin x; e  ax sin bx dx. x y (0)  z (0)  0. y (0)  1, y(0)  0. 17.  dy  dx  3z  y;  dz   z  y  ex,  dx y  y ( x)  sin 2 x;  dy  dx  2 y  3 z;  dz   y  1,  dx 18. 19. y (0)  z (0)  0. 72 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30.  0  (cos ax  cos bx) dx . x y (0)  0, y(0)  0. y(0)  2. y(0)  2. y  y ( x)  2  e  x ; y  6 y  5 y ( x)  3 x; y  6 y  13 y ( x)  2 x;  x  sin x 0  y (0)  1,  y (0)  0,  x3 x 2 dx. dx. dx. 1  cos ax 0   sin ax  a  sin x 0  x2 x sh x   cos( x   )  y ( ) d . 0 sin x   sin( x   )  y ( ) d . 0 2 x y ( x)   e ( x  )  y ( ) d  cos 3 x. 0 x y ( x)   sin( x   )  y ( ) d  x 3 . 0 x Пример 16. Пользуясь теоремой Руше, найти количество лежащих внутри круга z  1 корней заданных уравнений: 1. 2. 3. 4. z  5 z  z  2  0. 2 z  z  3 z  z  8  0. z  2 z  z  8 z  2  0. z  7 z  1  0. 4 9 6 2 5 3 2 7 4 2 5. 6. 7. 8. 73 z 7  5 z 4  z 2  2  0. z 4  8 z  10  0. 2 z 5  z 3  3 z 2  z  8  0. z 6  6 z  1  0. 9. z 5  8 z  11  0. 20. z 5  6 z 2  z  1  0. 21. z 4  8 z 2  z  5  0. 22. z 4  5 z 2  3  0. 23. z 3  5 z 2  1  0. 24. z 3  6 z 2  z  1  0. 25. z 3  7 z 2  2 z  2  0. 26. z 5  6 z 2  12  0. 27. z 5  7 z 2  z  1  0. 28. z 5  8 z 2  z  4  0. 29. z 5  4 z  1  0. 30. z 5  5 z 2  1  0. 10. z 9  2 z 6  z 2  8 z  2  0. 11. z 4  7 z 2  z  1  0. 12. z 4  5 z 2  z  1  0. 13. z 8  4 z 2  1  0. 14. z 7  3 z 2  1  0. 15. z 7  5 z 2  1  0. 16. z 7  4 z 2  2  0. 17. z 5  5 z 2  z  1  0. 18. z 5  5 z 2  2 z  1  0. 19. z 5  4 z 2  1  0. 74 СПИСОК ЛИТЕРАТУРЫ 1. Свешников А.Г., Тихонов А.Н. Теория функций комплексной переменной. – М.: Наука, 1979. – 319 с. 2. Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. – М.: Наука, 1987. – 688 с. 3. Сидоров Ю.В., Федорюк М.В., Шабунин М.И. Лекции по теории функций комплексного переменного. – М.: Наука, 1982. – 477 с. 4. Волковыский Л.И., Лунц Г.Л., Араманович И.Г. Сборник задач по теории функций комплексного переменного. – М.: Наука, 1975. – 320 с. 5. Шелковников Ф.А., Такайшвили К.Г. Сборник упражнений по операционному исчислению. – М.: Высш. школа, 1978. – 184 с. 6. Тихонов А.Н., Самарский А.А. Уравнения математической физики. – М.: Наука, 1977. – 724 с. 7. Кошляков Н.С., Глинер Э.Б., Смирнов М.М. Уравнения в частных производных математической физики. – М.: Высш. школа, 1970. – 712 с. 8. Сборник задач по теории аналитических функций / Под ред. М.А. Евграфова. – М.: Наука, 1972. – 416 с. 9. Сборник задач и упражнений по специальным главам высшей математики / Под ред. Г.И. Кручковича. – М.: Высш. школа, 1970. – 512 с. 10. Лунц Г.Л., Эльсгольц Л.Э., Араманович И.Г. Функции комплексного переменного. Операционное исчисление. Устойчивость движения. – М.: Наука, 1968. – 416 с. 11. Шахно К.У. Элементы теории функций комплексной переменной и операционного исчисления. – Минск: Вышейш. школа, 1975. – 400 с. с илл. 12. Грищенко А.Е., Нагнибида Н.И., Настасиев П.П. Теория функций комплексного переменного (Решение задач). – К.: Вища шк., 1986. – 336 с. 75 Навчальне видання Кондратьєв Борис Володимирович МЕТОДИЧНІ ВКАЗІВКИ до розвʼязку задач з теорії функцій комплексного змінного для студентів радіофізичного та фізичного факультетів Коректор Л. Є. Стешенко Комп’ютерна верстка Н. Є. Пруднік Макет обкладинки І. М. Дончик Формат 60х84/16. Умов. друк. арк. 3,98. Наклад 100 прим. Зам. № 23/11 Видавець і виготовлювач Харківський національний університет імені В. Н. Каразіна, 61077, м. Харків, пл. Свободи, 4. Свідоцтво суб’єкта видавничої справи ДК № 3367 від 13.01.2009 Видавництво ХНУ імені В. Н. Каразіна Тел. 705-24-32 76