Как найти тип особой точки

Изолированные особые точки функций и полюсы

Классификация особых точек

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

Исследование функции в особой точке z_0 определяется поведением ее в окрестности этой точки, т.е. исследованием limlimits_{zto z_0}f(z). Очевидно, имеют место три возможности:

а) limlimits_{zto z_0}f(z) не существует;
б) limlimits_{zto z_0}f(z) существует и равен конечному числу;
в) limlimits_{zto z_0}f(z) равен бесконечности.

Исследование пределов функции в комплексной области — задача более сложная, чем в действительной области, так как, согласно определению, переменная z стремится к z_0~(zto z_0) по любому направлению. Вычисление пределов в точках аналитичности не представляет интереса, так как в этих случаях limlimits_{zto z_0}f(z)=f(z_0).

Будем рассматривать limlimits_{zto z_0}f(z), где z_0 — особая точка.

Пример 4.1. Исследовать существование limlimits_{zto0}exp frac{1}{z},~ limlimits_{zto0} exp frac{1}{z^2} в случаях a) z=xinmathbb{R}; б) zinmathbb{C}.

Решение

a) В действительной области limlimits_{zto0}exp frac{1}{z} не существует, так как не равны односторонние пределы limlimits_{xto0+0}exp frac{1}{x}=infty,~ limlimits_{xto0-0}exp frac{1}{x}=0, но существует предел второй функции: limlimits_{xto0}exp frac{1}{x^2}=infty.

б) В комплексной области, очевидно, limlimits_{zto0}exp frac{1}{z} не существует, так как он не существует в частном случае z=x.

Но для второй функции полученного выше результата limlimits_{xto0}exp frac{1}{x^2}=infty не достаточно, так как рассмотрены только два направления на плоскости — по действительной положительной и действительной отрицательной полуосям.

Рассмотрим еще какое-нибудь направление, например по мнимой оси, т.е. z=iy,~ yto0colon

limlimits_{zto0}exp frac{1}{z^2}= limlimits_{zto0}exp frac{1}{(iy)^2}= limlimits_{zto0}exp frac{-1}{y^2}=0.

Сравнивая этот результат с полученным выше limlimits_{xto0}exp frac{1}{x^2}=infty, заключаем, что в комплексной области limlimits_{zto0}exp frac{1}{z^2} не существует.

Аналогично можно показать, что не существует limlimits_{zto0}exp frac{1}{z^4}, хотя limlimits_{zto0}exp frac{1}{z^4}=infty для случаев z=x и z=iy (по действительной и мнимой осям).

Эти простые примеры показывают, что исследование функции в особой точке с помощью limlimits_{zto z_0}f(z) может представлять большие сложности. Но, с другой стороны, в примере 3.36 при вычислении пределов функции в особых точках было использовано разложение функции в ряд.

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

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

В частности, конечная особая точка z_0inmathbb{C} является изолированной особой точкой функции f(z), если существует число r>0, такое, что в круге |z-z_0|<r эта точка- единственная особая точка f(z), а в проколотой окрестности, т.е. в 0<|z-z_0|<r функция f(z) аналитическая.

Бесконечно удаленная особая точка z_0=infty является изолированной особой точкой функции f(z), если существует число R>0, такое, что в области |z|>R эта точка — единственная особая точка f(z), а в кольце R<|z|<infty функция f(z) — аналитическая.

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

Пример 4.2. Исследовать поведение и вид ряда Лорана в окрестности особой точки z=0 функций:

а) frac{sin z}{z}; б) frac{sin z}{z^n},~n>1; в) exp frac{1}{z}.

Решение

Эти простые примеры показывают, что поведение функции в особой точке связано с видом главной части ряда Лорана: трем отмеченным выше случаям нахождения предела функции в точке z_0 соответствуют три различных случая вида главной части ряда Лорана в окрестности точки. В примере 4.2 исследовалась конечная особая точка. Такой же результат можно получить, рассматривая точку z=infty, например, для функций exp frac{1}{z},~ z^2+frac{1}{z} и e^z.


Типы особых точек функции

В зависимости от трех случаев поведения функции в особой точке (исследования limlimits_{zto z_0}f(z)) особые точки функций делят на три типа — производится классификация особых точек. В качестве определения типа особых точек можно выбрать либо поведение функции в особой точке, либо вид ряда Лорана. Выберем первый подход.

Изолированная особая точка z_0in overline{C} функции f(z) называется:

– устранимой особой точкой, если limlimits_{zto z_0}f(z) существует и конечен (4.1);
– полюсом, если limlimits_{zto z_0}f(z)=infty (4.2);
– существенно особой точкой, если limlimits_{zto z_0}f(z) не существует (4.3).

Замечание 4.1. Если в случае устранимой особой точки z_0 положить f(z_0)= limlimits_{zto z_0}f(z), то f(z) будет аналитической в O_{delta}(z_0) и точку z_0 можно считать правильной, т.е. не особой. В этом случае говорят, что в точке z_0 устранена особенность.

Пример 4.3. Определить тип особой точки z=0 для функций frac{sin z}{z};~~frac{sin z}{z^n},~n>1;~~exp frac{1}{z},~~ exp frac{1}{z^2}.

Решение

На основании результатов решения примеров 4.1, 4.2 заключаем, что z=0 является устранимой особой точкой функции frac{sin z}{z}; полюсом для frac{sin z}{z^n} при любом n>1; существенно особой точкой для функций exp frac{1}{z} и exp frac{1}{z^2}.

Пример 4.4. Определить тип особой точки z=infty для функций f_1(z)=z^nsin frac{1}{z} и f_2(z)=e^z.

Решение

Рассмотрим limlimits_{ztoinfty}f(z). Для удобства введем обозначение frac{1}{z}=xi. Для функции f_1(z) получим limlimits_{ztoinfty}f_1(z)= limlimits_{xito0} frac{sinxi}{xi^n}=infty~(n>1) (см. пример 4.2), поэтому z=infty является полюсом функции f_1(z)=z^nsin frac{1}{z}. Для функции f_2(z)=e^z точка z=infty является существенно особой, так как limlimits_{ztoinfty}e^z= limlimits_{xito0}exp frac{1}{xi} не существует (см. пример 4.1).

Пример 4.5. Найти все конечные особые точки функций: а) f_1=frac{sin z}{z^4+1} б) f_2(z)= frac{sin z}{sin z^{-1}} и определить их тип.

Решение

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

а) Так как числитель и знаменатель функции f_1(z) — функции аналитические, то ее особыми точками являются только нули знаменателя, т.е. корни уравнения z^4+1=0. Это четыре точки z_k=exp frac{(pi+2kpi)i}{4},~k=0,1,2,3, или в алгебраической форме: z_{1,2}= frac{pm1+i}{sqrt{2}},~ z_{3,4}= frac{pm1-i}{sqrt{2}}. Заметим, что точки расположены в вершинах квадрата, вписанного в окружность радиуса R=1 с центром в начале координат, и справедливы равенства z_2=iz,~ z_3=-z_1,~ z_4= overline{z}_1.

Очевидно, все точки изолированные и являются полюсами, так как limlimits_{zto z_k}f_1(z)=infty для любой точки z_k,~k=1,2,3,4.

б) Особыми точками функции f_2(z) являются нули знаменателя, т.е. точки, для которых frac{1}{z_k}=kpi или z_k=frac{1}{kpi},~ k=pm1,pm2,ldots, а также z=0 — особая точка знаменателя. Точки z_k являются полюсами, так как limlimits_{zto z_k}f_2(z)=infty. Точка z=0 — неизолированная особая точка функции, так как в любой ее окрестности |z|<r (r — любое число, r>0), кроме этой точки, расположено бесконечное множество особых точек вида z_k=frac{1}{kpi},~ k=pm1,pm2,ldots. Точку z=0 в таком случае называют предельной точкой полюсов z_k=frac{1}{kpi}, так как limlimits_{ktoinfty}z_k=0.


Теоремы Сохоцкого и Пикара

Для исследования поведения функции в существенно особой точке имеют место следующие две теоремы.

Теорема 4.1 (Сохоцкого). Если r_0 — существенно особая точка функции f(z), то для любого Ain overline{mathbb{C}} существует последовательность {z_n}, сходящаяся к точке z_0, такая, что limlimits_{ntoinfty}f(z_n)=A.

Теорема 4.2 (Пикара). В любой окрестности существенно особой точки функция f(z) принимает любое значение (причем бесконечное число раз) кроме, быть может, одного.

Пример 4.6. Исследовать поведение следующих функций в существенно особых точках, проиллюстрировать теоремы Сохоцкого и Пикара:

a) f_1(z)=exp frac{1}{z},~ f_2(z)= exp frac{1}{z^2},~ z_0=0; б) f_3(z)=e^z,~ z_0=infty.

Решение

В примерах 4.3 и 4.4 показано, что точки z_0=0 и z_0=infty являются существенно особыми точками соответствующих функций. Исследуем пределы функций.

а) Для иллюстрации теоремы Сохоцкого выбираем A=0 и A=infty. Используя результат примера 4.1, имеем limlimits_{zto0}exp frac{1}{z}=0, если z=x<0, и limlimits_{zto0}exp frac{1}{z}=infty, если z=x>0, то есть limlimits_{ntoinfty}f_1(z_n)=0для последовательности z_n=x_n, такой, что limlimits_{ntoinfty}x_n=0 и x_n<0, и limlimits_{ntoinfty}f_1(z_n)=0 для последовательности z_n=x_n, такой, что limlimits_{ntoinfty}x_n=0 и x_n>0.

Аналогично исследуем функцию f_2(z). Для числа A=0 выбираем z_n=iy_n, где limlimits_{ntoinfty}y_n=0 и тогда limlimits_{ntoinfty}f_2(z_n)=0, а для A=infty выбираем z_n= x_n, где limlimits_{ntoinfty}x_n=0 и тогда limlimits_{ntoinfty} f_2(z_n)=infty.

Справедливость теоремы Пикара для этих функций следует из рассмотрения уравнений exp frac{1}{z}=A,~ exp frac{1}{z^2}=A, которые, как известно, имеют бесконечное множество решений для любого Ain mathbb{C},~ Ane0.

Например, для функции f_1(z)= expfrac{1}{z} имеем expfrac{1}{z}=A. Отсюда получаем

frac{1}{z}= operatorname{Ln}A,quad frac{1}{z}= ln|A|+i(arg A+2kpi),quad k=0,pm1,pm2,ldots или z_k=frac{1}{ln|A|+i(arg A+2kpi)}.

В частности, функция exp frac{1}{z} в любой окрестности точки z_0=0 принимает значение A=1 бесконечное множество раз: в точках z_k=frac{-i}{2kpi},~ k=pm1,pm2,ldots (рис. 4.1).

Рис. 4.1.

б) Точка z=infty является существенно особой точкой функции e^z (пример 4.4). Обозначив z=frac{1}{xi}, можно повторить рассуждения предыдущего пункта для функции expfrac{1}{xi} и точки xi=0.


Ряд Лорана в окрестности особой точки

В предыдущем разделе на примере простых функций (см. пример 4.2) было высказано предположение, что вид ряда Лорана в окрестности особой точки зависит от типа особой точки и потому задача исследования функции в особой точке может быть сведена к исследованию соответствующего ряда Лорана . Подтверждением этого предположения в общем случае является доказательство соответствующих утверждений.

Утверждение 4.1

1. Для того чтобы особая точка функции была ее устранимой особой точкой, необходимо и достаточно, чтобы в разложении функции в ряд Лорана в окрестности этой точки отсутствовала главная часть. Это означает, что если z_0 — устранимая особая точка, то ряд Лорана функции f(z) имеет вид

f(z)= sum_{n=0}^{infty} c_n(z-z_0)^n,quad 0<|z-z_0|<r

(4.4)

для z_0 — конечной точки z_0inmathbb{C}, и (для z_0=infty)

f(z)= c_0+sum_{n=1}^{infty} frac{c_{-n}}{z^n},quad R<|z|<infty.

(4.5)

2. Для того чтобы особая точка функции была полюсом, необходимо и достаточно, чтобы главная часть ряда Лорана функции в окрестности этой точки содержала конечное число членов. Ряд Лорана функции f(z) в случае z_0 полюса имеет вид

f(z)=sum_{k=-n}^{infty} c_k(z-z_0)^k,quad 0<|z-z_0|<r,

(4.6)

если z_0inmathbb{C}, и (если z_0=infty)

f(z)=sum_{k=-infty}^{n}c_kz^k,quad R<|z|<infty,

(4.7)

3. Для того чтобы особая точка функции была ее существенно особой точкой, необходимо и достаточно, чтобы главная часть ряда Лорана функции в окрестности этой точки содержала бесконечное число членов. Ряд Лорана функции f(z) в случае z_0 — существенно особой точки имеет вид

f(z)= sum_{n=-infty}^{infty} c_n(z-z_0)^n,quad 0<|z-z_0|<r,

(4.8)

если z_0inmathbb{C}, и (если z_0=infty)

f(z)= sum_{n=-infty}^{infty}c_nz^n,quad R<|z|<infty.

(4.9)

Замечания 4.2

1. Номер старшего члена главной части ряда Лорана функции в ее разложении в окрестности полюса называется порядком полюса.

Так, точка z_0inmathbb{C} является полюсом порядка n (Pi(n)) функции f(z), если в разложении (4.6) c_{-n}ne0,~ c_k=0 при k<-n. Точка z_0=infty является полюсом порядка n (Pi(n)) функции f(z), если в разложении (4.7) c_nne0,~ c_k=0 при k>n.

2. Главная часть ряда Лорана в случае полюса порядка и записывается следующим образом:

а) в случае z_0inmathbb{C} в виде sum_{k=-n}^{-1} c_k(z-z_0)^n, или sum_{k=1}^{n} frac{c_{-k}}{(z-z_0)^k}, или, подробнее:

frac{c_{-n}}{(z-z_0)^n}+frac{c_{-n+1}}{(z-z_0)^{n-1}}+ ldots+ frac{c_{-1}}{z-z_0},quad c_{-n}ne0,;

(4.10)

б) в случае z_0=infty в виде sum_{k=1}^{n}c_nz^k, или sum_{k=-n}^{-1} frac{c_{-k}}{z^k} (см. (4.7)), или, подробнее:

c_ncdot z^n+ c_{n-1}cdot z^{n-1}+ldots+c_1cdot z,quad c_nne0.

(4.11)

3. Главная часть ряда Лорана в случае существенно особой точки записывается так:

а) в случае z_0inmathbb{C} в виде sum_{n=-infty}^{-1} c_n(z-z_0)^n, или sum_{n=1}^{infty} frac{c_{-n}}{(z-z_0)^n} (см.(4.8)), или, подробнее:

frac{c_{-1}}{z-z_0}+ frac{c_{-2}}{(z-z_0)^2}+ ldots+ frac{c_{-n}}{(z-z_0)^n}+ldots;

(4.12)

б) в случае z_0=infty в виде sum_{n=-infty}^{-1} frac{c_{-n}}{z^n} или sum_{n=1}^{infty}c_nz^n (см.(4.9)), или, подробнее:

c_1cdot z+c_2cdot z^2+ldots+c_ncdot z^n+ldots

(4.13)

Пример 4.7. Определить тип особых точек функций: а) f_1(z)=frac{z+2}{z^2-2z-3}; б) f_2(z)=frac{z+2}{(z+1)^2(z-3)}.

Решение

Особыми точками функций являются z_1=-1,~z_2=3,~z_3=infty. Чтобы определить тип особой точки, используем разложения функций в окрестности каждой точки, полученные в примерах 3.31 , 3.33 , 3.34.

a) f_1(z)= frac{-1}{4}cdot frac{1}{z+1}-sum_{n=0}^{infty} frac{5(z+1)^n}{4^{n+2}},~ 0<|z+1|<4. В главной части разложения — один член ряда: frac{-1}{4}cdot frac{1}{z+1}, здесь c_{-1}=frac{-1}{4}ne0, все c_n=0 для n<-1. Следовательно, в точке z=-1 — полюс первого порядка, т.е. простой полюс функции f_1(z).

Аналогично из разложения f_1(z)= frac{5}{4}cdot frac{1}{z-3}+ sum_{n=0}^{infty} frac{(-1)^{n+1}}{4^{n+2}}(z-3)^n,~ 0<|z-3|<4 получим такой же результат: точка z=3 — простой полюс функции f_1(z).

Разложение f_1(z)= sum_{n=1}^{infty} frac{(-1)^n+5cdot3^{n-1}}{4}cdot frac{1}{z^n},~ |z|>3. Функции в окрестности z=infty не содержит главной части — разложение имеет вид (4.5). Следовательно, точка z=infty — устранимая особая точка функции f_1(z).

б) Из разложения f_2(z)= frac{5}{16}cdot frac{1}{z-3}+ sum_{n=0}^{infty} frac{(-1)^{n+1}(n+6)}{4^{n+3}}(z-3)^n,~ 0<|z-3|<4 следует, что z=3 — простой полюс функции f_2(z).

Из разложения f_2(z)= frac{-5}{16}cdot frac{1}{z+1}+ frac{-1}{4}cdot frac{1}{(z+1)^2}+ sum_{n=0}^{infty} frac{-5(z+1)^n}{4^{n+3}},~ 0<|z+1|<4, где c_{-2}=frac{-1}{4}ne0 и все c_n=0 для n<-2, получаем, что z=-1 — полюс второго порядка функции f_2(z).

Разложение f_2(z) в окрестности z=infty и не содержит положительных степеней, в чем можно убедиться, проанализировав разложения элементарных дробей (см. пример 3.34). Поэтому z=infty — устранимая особая точка функции f_2(z).

Пример 4.8. Определить тип конечных особых точек для функций:

а) exp frac{1}{z-a},~ sin frac{1}{z-a},~ cos frac{1}{z-a}; б) frac{sin z}{z^n},~ frac{1-cos z}{z^3}.

Решение

а) Используем разложения функций по степеням (z-a)colon

begin{gathered}exp frac{1}{z-a}= 1+frac{1}{z-a}+ frac{1}{2!(z-a)^2}+ldots;qquad sinfrac{1}{z-a}= frac{1}{z-a}- frac{1}{3!(z-a)^3}+ldots;\ cosfrac{1}{z-a}= 1-frac{1}{2!(z-a)^2}+ldots end{gathered}

Убеждаемся, что для всех указанных функций точка z=a является существенно особой точкой, так как в разложениях главная часть содержит бесконечное число членов, т.е. имеется бесконечное число членов с отрицательными степенями (см. п.1 утверждения 4.1).

б) Запишем разложения функций по степеням zcolon

begin{aligned}frac{sin z}{z^n}&= frac{1}{z^n}! left(z-frac{z^3}{3!}+ldotsright)= frac{1}{z^{n-1}}-frac{1}{z^{n-3}cdot3!}+ldots;\ frac{1-cos z}{z^3}&= frac{1}{z^3}! left(1-left(1-frac{z^2}{2!}+frac{z^4}{4!}+ldotsright)right)= frac{1}{2!z}-frac{z}{4!}+ldots end{aligned}

Для первой функции при n=1 в разложении отсутствует главная часть — совокупность членов с отрицательными степенями. Следовательно, согласно п.1 утверждения 4.1, точка z=0 для frac{sin z}{z} является устранимой особой точкой.

При n>1 главная часть разложения содержит конечное число членов, поэтому точка z=0 для frac{sin z}{z^n} является полюсом (см. п.2 утверждения 4.1). Кроме того, так как при n>1 в разложении старшая отрицательная степень равна (n-1), то, согласно п. 1 замечаний 4.2, заключаем, что z=0 для frac{sin z}{z^n} при n>1 является полюсом порядка (n-1). Рассуждая аналогично, получаем, что z=0 является полюсом первого порядка — простым полюсом для функции frac{1-cos z}{z^3}.

Сравнивая разложения функций по степеням z в окрестности z_0=0 (формулы (4.4),(4.6),(4.8) при z_0=0) и z=infty (формулы (4.5), (4.7), (4.9)), можно сделать следующее заключение.

Утверждение 4.2

1. Чтобы z=infty была устранимой особой точкой функции f(z), необходимо и достаточно, чтобы точка xi=0 была устранимой (или не особой) для varphi(xi)= f!left(frac{1}{xi}right).

2. Чтобы z=infty была полюсом порядка n функции f(z), необходимо и достаточно, чтобы точка xi=0 была полюсом порядка n функции varphi(xi)= f!left(frac{1}{xi}right).

3. Чтобы z=infty была существенно особой точкой функции f(z), необходимо и достаточно, чтобы точка xi=0 была существенно особой точкой функции varphi(xi)= f!left(frac{1}{xi}right).

Замечание 4.3. Как и в случае конечной особой точки z_0, в которой функция не определена, но limlimits_{zto z_0}f(z)=0 (см. утверждение 3.5) , так и для z=0 в случае limlimits_{ztoinfty}f(z)=0, устранимую особую точку z=infty можно считать нулем функции f(z). Порядок нуля можно определить как порядок нуля функции varphi(xi) в точке xi=0.

Пример 4.9. Исследовать точку z=infty для функций: a) f(z)=frac{1}{z^2(z-2)}; б) f(z)=frac{z^2+1}{3z^2-2}; в) f(z)= zsinfrac{1}{z}.

Решение


Правила определения порядка полюса

Используя формулу (4.6) разложения функции в ряд в окрестности полюса, можно получить практически удобные правила определения порядка полюса, не требующие записи разложений в ряд в каждом конкретном случае.

Пусть z_0~(z_0in mathbb{C}) — полюс порядка n~(Pi(n)) функции f(z). Разложение (4.6), где главная часть имеет вид (4.10) , преобразуем следующим образом:

f(z)=frac{1}{(z-z_0)^n}bigl[c_{-n}+c_{-n+1}cdot (z-z_0)+ ldots+c_0(z-z_0)+ ldotsbigr] или f(z)=frac{varphi(z)}{(z-z_0)^n},

где varphi(z) — функция, аналитическая в точке z_0, как сумма степенного ряда, записанного в скобках, и varphi(z_0)=c_{-n}ne0.

Далее рассмотрим функцию F(z)=frac{1}{f(z)}, то есть F(z)=frac{(z-z_0)^n}{varphi(z)} или F(z)=(z-z_0)^n varphi_1(z), где varphi_1(z) — аналитическая в точке z_0 и varphi_1(z_0)ne0. Из этого, согласно утверждению 3.5, следует, что z_0 является нулем порядка n функции F(z). Можно доказать и обратное утверждение.

А именно, если функция представлена в виде f(z)=frac{varphi(z)}{(z-z_0)^n}, где varphi(z) — функция, аналитическая в точке z_0, и varphi(z_0)ne0, то z_0 — полюс порядка n функции f(z), а также, если z_0 — нуль порядка n функции F(z), то для функции f(z)=frac{1}{F(z)} эта точка является полюсом порядка n.

Кроме того, рассмотрим частное f(z)=frac{f_1(z)}{f_2(z)}, где точка z_0 является нулем порядка k для функции f_1(z) и нулем порядка m для функции f_2(z), то есть f(z)= frac{(z-z_0)^k varphi_1(z)}{(z-z_0)^m varphi_2(z)}. При m>k получаем f(z)= frac{varphi(z)}{(z-z_0)^{m-k}}, из чего, с учетом приведенных выше рассуждений, находим, что z_0 — полюс порядка (m-k). Заметим, что при m=k точка z_0 — устранимая особая точка; случай m<k рассмотрен ранее. Результаты приведенных рассуждений запишем в виде утверждения.

Утверждение 4.3

1. Для того чтобы точка z_0 была полюсом порядка n функции f(z), необходимо и достаточно, чтобы ее можно было записать в виде

f(z)= frac{varphi(z)}{(z-z_0)^n},quad varphi(z_0)ne0.

(4.14)

2. Для того чтобы точка z_0 была полюсом порядка n функции f(z), необходимо и достаточно, чтобы она была нулем порядка n функции frac{1}{f(z)}(связь нулей с полюсами).

3. Если точка z_0 является нулем порядка k функции f_1(z) и нулем порядка m функции f_2(z)~(m>k), то она — полюс порядка (m-k) для frac{f_1(z)}{f_2(z)}.

Пример 4.10. Определить порядок полюсов функций из примеров: а) 4.7 ; б) 4.8.

Решение


Замечания 4.4

1. Так как конечными особыми точками рациональной дроби frac{P_n(z)}{Q_m(z)} являются только нули знаменателя, то это либо полюсы, либо устранимые особые точки функции.

2. Такое же заключение можно сделать и для функции вида f(z)= frac{varphi(z)}{Q_m(z)}, где varphi(z) — аналитическая функция. При этом, используя определение устранимой особой точки (4.1) и правила определения порядка нуля и полюса (утверждения 3.5 и 4.3), можно сделать следующие выводы относительно особой точки z_0 — нуля порядка k~[0(k)] знаменателя:

а) z_0 — полюс порядка k функции f(z), если varphi(z_0)ne0;
б) z_0 — полюс порядка (k-n), если z_0 — нуль порядка n функции varphi(z) и k>n;
в) z_0 — устранимая особая точка функции f(z), если z_0 — нуль порядка k функции varphi(z);

г) z_0 — нуль порядка (n-k) функции f(z), если z_0 — нуль порядка n функции varphi(z) и n>k; при этом полагаем f(z_0)= limlimits_{zto z_0}f(z)=0.

▼ Примеры нахождения особых точек и определения их типа

Пример 4.11. Найти конечные особые точки следующих функций и определить их тип:

а) f(z)= frac{z^2-4z+3}{(z-1)^3(z+2)^2}; б) f(z)= frac{z+2i}{z^3+8i}.

Решение

Конечными особыми точками этих рациональных дробей являются нули знаменателя. Чтобы для каждой их этих точек определить, является ли она полюсом или устранимой особой точкой, нужно, согласно определению, найти предел функции в этой точке. В случае полюса, т.е. когда limlimits_{zto z_0}f(z)=infty, далее следует определить его порядок. Для этого используется утверждение 4.3.

Можно поступить иначе — согласно замечанию 4.4. Для этого нужно найти и нули числителя.

а) Особые точки функции z_1=1,~ z_2=-2. Для точки z_2=-2 можно применить формулу (4.14) и из f(z)= frac{varphi(z)}{(z+2)^2}, где varphi(z)= frac{z^2-4z+3}{(z-1)^3} и varphi(-2)ne0, получить, что эта точка — полюс второго порядка. Для точки z_1=1 формула (4.14) не применима, так как из varphi(z)= frac{z^2-4z+3}{(z+2)^2} имеем varphi(1)=0. Поступаем далее согласно замечанию 4.4. Раскладываем на множители числитель и записываем функцию

f(z)= frac{(z-1)(z-3)}{(z-1)^3(z+2)^2}= frac{z-3}{(z-1)^2(z+2)^2}= frac{varphi(z)}{(z-1)^2},quad varphi(1)ne0.

Получаем, что z=1 — полюс второго порядка для f(z).

б) Особые точки функции — корни уравнения z^3+8i=0, то есть z=sqrt[LARGE{3}]{-8i} или z_k= exp left[left(-frac{pi}{6}+ frac{2kpi}{3}right)!iright],~ k=0,1,2. Все эти точки: z_{1,3}=2! left(pmfrac{sqrt{3}}{2}-i frac{1}{2}right)!,~ z_2=2i — простые нули знаменателя, и так как числитель в этих точках не обращается в нуль, то они — простые полюсы функции f(z).

Пример 4.12. Найти конечные особые точки следующих функций и определить их тип:

а) f(z)= frac{z^6+z^4-z^2-1}{(z^2-4z+3)^2}; б) f(z)= frac{z-2i}{z^3+8i}.

Решение

Пример 4.13. Найти конечные особые точки функций и определить их тип:

а) f(z)=frac{cospi z}{(4z^2-1)(z^2+3)}; б) f(z)=frac{sinpi z}{z^4-1}.

Решение

Конечными особыми точками этих функций вида f(z)= frac{varphi(z)}{Q_4(z)}, где varphi(z) — аналитическая функция, являются только нули знаменателя.

а) Особые точки функции: z_1=frac{2}{2},~ z_2=-frac{1}{2},~ z_3=isqrt{3},~-i sqrt{3}. Точки z_3 и z_4 — простые полюсы, так как числитель в этих точках не обращается в нуль и функцию можно представить в виде f(z)= frac{varphi(z)}{z-z_0},~ varphi(z_0)ne0,~ z_0 — точка z_3 или z_4. В точках z=pmfrac{1}{2} числитель обращается в нуль. Очевидно, это простые нули числителя, и поэтому его можно записать в виде cospi z= (z-z_0)cdot varphi(z),~ varphi(z_0)ne0,~ z_0 — точка z_1 или z_2. Тогда для функции f(z) получаем

f(z)= frac{(z-z_0)cdot varphi(z)}{4 left(z-dfrac{1}{2}right)!! left(z+ dfrac{1}{2}right)! (z^2+3)},.

Так как limlimits_{zto z_k} f(z)neinfty для z_0=z_1=frac{1}{2} или z_0=z_2=-frac{1}{2}, то эти точки — устранимые особые точки функции f(z).

б) Особые точки функции: z_{1,2}=pm1,~ z_{3,4}=pm i. Точки z_3=i и z_4=-i — простые полюсы.

Для точек z_1=1 и z_2=-1 проводим рассуждения, как в предыдущем пункте, и находим, что они — устранимые особые точки f(z).

Пример 4.14. Определить тип особой точки z=0 для следующих функций: а) f(z)= frac{e^z-1-z}{e^{z^5}-1}; б) f(x)= frac{cos z^2-1}{z^3sin^2z^2}.

Решение. В точке z=0 и числитель, и знаменатель каждой из функций обращается в нуль. Определим порядок нуля в каждом случае и используем п.3 утверждения 4.4.

а) Из разложений по степеням z функций

begin{gathered} e^z-1-z&= left(1+z+frac{z^2}{2!}+ldotsright)-1-z= frac{z^2}{2!}+frac{z^3}{3!}+ldots\ e^{z^5}-1&= left(1+z^5+frac{z^{10}}{2!}+ldotsright)-1= z^5+frac{z^{10}}{2!}+ldots end{gathered}

находим, что z=0 — нуль второго порядка для числителя (0(2)) и нуль пятого порядка для знаменателя (0(5)). Следовательно, z=0 — полюс третьего порядка для функции f(z).

б) Используя правила определения порядка нуля, в частности, как и в предыдущем пункте, раскладывая функции в ряды по степеням z, находим, что z=0 является 0(2) для числителя и 0(7) для знаменателя. Следовательно, z=0 -полюс пятого порядка для f(z).

Пример 4.15. Найти конечные особые точки следующих функций и определить их тип:

а) f(z)= frac{e^{2+2i}}{z^4-16}; б) f(z)= frac{(z+pi) sin frac{pi z}{2}}{zsin^2z}.

Решение

Пример 4.16. Определить тип особой точки z=0 для следующих функций:

а) f(z)= frac{sin z-z}{cos z-1-dfrac{z^2}{2}}; б) f(z)= frac{sin z+z}{operatorname{ch}z-1-dfrac{z^2}{2}}; в) f(z)= frac{sin z-z}{cos z-1+dfrac{z^2}{2}}; г) f(z)= frac{sin z+z}{operatorname{ch}z-1+dfrac{z^2}{2}}.

Решение

Точка z=0 является нулем и знаменателя, и числителя для каждой из функций. Определим порядок нуля в каждом случае, используя правило определения порядка нуля (утверждение 3.5), в частности, раскладывая соответствующую функцию по степеням z.

а) Из разложений

begin{aligned}sin z-z&= left(z-frac{z^3}{3!}+ldotsright)-z=-frac{z^3}{3!}+ frac{z^5}{5!}+ldots\ cos z-1-frac{z^2}{2}= left(1-frac{z^2}{2!}+ frac{z^4}{4!}-ldotsright)-1-frac{z^2}{2}=-z^2+frac{z^4}{4!}-ldots end{aligned}

находим, что z=0 является 0(3) для числителя и 0(2) — для знаменателя, поэтому она — устранимая особая точка. Так как

limlimits_{zto0}f(z)= limlimits_{zto0} frac{z^3! left(-frac{1}{3}+ frac{z^2}{5!}+ ldotsright)}{z^2! left(-1+frac{z^2}{4!}-ldotsright)}=0,

то, полагая f(0)=0, можно считать, что z=0 — нуль для f(z), причем 0(1) (см. замечания 4.4).

б) Из разложений

sin z+z=2z-frac{z^3}{3!}+ldots и operatorname{ch}z-1-frac{z^2}{2}= left(1+frac{z^2}{2!}+ frac{z^4}{4!}+ldotsright)-1-frac{z^2}{2}= frac{z^4}{4!}+ frac{z^6}{6!}+ldots

находим, что z=0 является 0(1) для числителя и 0(4) — для знаменателя. Поэтому z=0 — полюс третьего порядка для f(z).

в) Как и в предыдущих пунктах, находим, что z=0 является 0(3) для числителя и 0(4) — для знаменателя. Поэтому z=0 — простой полюс для f(z).

г) Точка z=0 является простым нулем числителя, нулем второго порядка для знаменателя. Следовательно, это простой полюс для f(z).


Определение порядка полюса в бесконечно удаленной точке

Рассмотрим бесконечно удаленную точку. Тип особой точки можно определить, вычисляя limlimits_{ztoinfty}f(z) или раскладывая функцию в ряд Лорана (см. примеры 4.4, 4.7). Можно свести задачу к исследованию конечной точки xi=0,left(frac{1}{z}= xiright) (см. утверждение 4.2 и пример 4.9). В двух последних случаях определяется и порядок полюса.

Практически удобное правило определения порядка полюса z=infty можно получить, используя п. 2 утверждения 4.2 и правила определения порядка полюса в конечной точке (утверждение 4.3). Действительно, пусть z=infty-Pi(n) для функции f(z), тогда xi=0-Pi(n) для f! left(frac{1}{xi}right) и можно записать f! left(frac{1}{xi}right)= frac{F(xi)}{xi^n},~ F(xi)ne0 (см. (4.14)). Поэтому, обозначив f! left(frac{1}{z}right)= varphi(z), для f(z) получим

f(z)=z^ncdot varphi(z),qquad limlimits_{ztoinfty} varphi(z)ne0,qquad limlimits_{ztoinfty} varphi(z)neinfty.

(4.15)

Представление функции в виде (4.15) является необходимым и достаточным условием полюса порядка n функции f(z) в точке z=infty.

Замечание 4.5. Используя формулу (4.15), нетрудно убедиться, что если z=infty-Pi(m) для f_1(z) и Pi(k) для f_2(z), то z=infty — полюс порядка (m-k) для функции f(z)= frac{f_1(z)}{f_2(z)}.

Пример 4.17. Определить тип особой точки z=infty для функций: а) f(z)= z^2(z-2); б) f(z)=frac{z^5+z^2+1}{z^3-2z}.

Решение

Так как limlimits_{ztoinfty}f(z)=infty в обоих случаях, то z=infty для данных функций — полюс. Определим порядок полюса.

а) Точка z=infty является полюсом третьего порядка, в чем можно убедиться любым из следующих способов.

Первый способ. Разложение функции по степеням z имеет вид f(z)=z^3-2z,~c_3ne0, все c_n=0,~n>3, и по определению (см. формулы (4.7), (4.11)) заключаем, что z=infty-Pi(3).

Второй способ. Обозначим z=frac{1}{xi}, получим функцию f! left(frac{1}{xi}right)= frac{1-2xi}{xi^3}, для которой xi=0-Pi(3). Поэтому, согласно п. 2 утверждения 4.2, точка z=infty-Pi(3) для f(z).

Третий способ. Запишем функцию в виде f(z)=z^3! left(1-frac{2}{z}right) и, так как функция varphi(z)=left(1-frac{2}{z}right) — удовлетворяет условиям формулы (4.15), получим, что z=infty-Pi(3) для f(z).

б) Разложение функции в ряд по степеням z представляет некоторые трудности. Используем другие способы.

Первый способ. Обозначим z=frac{1}{xi}, получим f! left(frac{1}{xi}right)= frac{(xi^5+xi^3+1)cdotxi^3}{xi^5cdot(1-2xi^2)}, или f! left(frac{1}{xi}right)= frac{xi^5+xi^3+1}{xi^2cdot(1-2xi^2)}.

Поэтому xi=0 является Pi(2) для f!left(frac{1}{xi}right) и, следовательно, z=infty-Pi(2) для f(z).

Второй способ. Представим функцию в виде f(z)= frac{z^5cdot left(1+ dfrac{1}{z^3}+dfrac{1}{z^2}right)}{z^3cdot left(1-dfrac{2}{z^2}right)} или f(z)= z^2cdot varphi(z), где limlimits_{ztoinfty} varphi(z)=1, и, согласно формуле (4.15), z=infty-Pi(2) для f(z).

Третий способ. Используем замечание 4.5. Можно определить порядок полюса z=infty для дроби frac{f_1(z)}{f_2(z)}, зная соответствующие порядки полюсов числителя и знаменателя. Здесь, очевидно, z=infty-Pi(5) для числителя и Pi(3) — для знаменателя (см. формулы (4.7), (4.11)). Поэтому z=infty-Pi(2) для f(z).

Пример 4.18. Определить порядок полюса в точке z=infty для следующих функций: а) f(z)= frac{z^3}{z-1}; б) f(z)=z^2+1+frac{1}{z^2-1}.

Решение


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

Пусть z_0 — особая точка функций f_1(z) и f_2(z) и тип особой точки для каждой из функций известен. Требуется определить тип особой точки для функций f_1(z)pm f_2(z);~ f_1(z)cdot f_2(z);~ frac{f_1(z)}{f_2(z)}. Рассмотрим следующие случаи.

Первый случай. Пусть точка го является полюсом порядка m~(Pi(m)) для функции f_1(z) и полюсом порядка k~(Pi(k)) для функции f_2(z).

а) При исследовании суммы f(z)= f_1(z)+ f_2(z) воспользуемся формулой (4.14) (п.1 утверждения 4.3) и запишем слагаемые в виде

f_1(z)=frac{varphi_1(z)}{(z-z_0)^m},~~ f_2(z)=frac{varphi_2(z)}{(z-z_0)^k}, где varphi_1(z_0)ne0,~varphi_2(z_0)ne0.

При k=m для суммы f(z)= f_1(z)+ f_2(z) получаем f(z)= frac{varphi_1(z)+varphi_2(z)}{(z-z_0)^m} или f(z)=frac{varphi(z)}{(z-z_0)^m}, где varphi(z)= varphi_1(z)+ varphi_2(z). Если varphi(z_0)ne0, то z_0-Pi(m) для функции f(z). Однако для функций f_1(z),, f_2(z) может выполняться условие varphi_1(z_0)+ varphi_2(z_0)=0 и’ следовательно, varphi(z_0)=0. В этом случае формула (4.14) не применима и точка z_0 не будет полюсом порядка m для f(z). В соответствии с п.3 утверждения 4.3 порядок полюса будет меньше, чем m, и равен (m-p) в случае m>p, где p — порядок нуля функции varphi(z). Если p=m, то z_0 — устранимая особая точка для f(z).

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

б) Для исследования произведения f_1(z)cdot f_2(z) воспользуемся формулой связи нулей с полюсами (п.2 утверждения 4.3) и рассмотрим вспомогательные функции F_1(z)= frac{1}{f_1(z)},~ F_2(z)= frac{1}{f_2(z)}. Для первой из этих функций z_0-0(m), для второй соответственно z_0-0(k) . а поэтому для F(z)= F_1(z)cdot F_2(z) она будет 0(m+k). Согласно п.2 утверждения 4.3, z_0 является Pi(m+k) для f(z)= f_1(z)cdot f_2(z).

в) Аналогичные рассуждения для частного frac{f_1(z)}{f_2(z)} приводят к результату: при m>k точка z_0 является Pi(m-k) для frac{f_1(z)}{f_2(z)}.

Второй случай. Пусть точка z_0 является полюсом, устранимой особой точкой или не особой для f_1(z) и существенно особой для f_2(z). Так как limlimits_{zto z_0}f_2(z) не существует, то по свойству пределов он не существует для каждой из рассматриваемых комбинаций f_1(z)pm f_2(z);~ f_1(z)cdot f_2(z);~ frac{f_1(z)}{f_2(z)}. Следовательно, для каждой из них z_0 — существенно особая точка. Заметим, что для функции frac{1}{f_2(z)} эта точка является либо существенно особой точкой, либо не является изолированной особой точкой. Последнее проиллюстрировано в примере 4.5 для функции frac{1}{sin z^{-1}}.

Третий случай. Пусть z_0 — полюс порядка n для f_1(z) и устранимая особая точка для f_2(z). Разложения этих функций в ряд в окрестности z_0 имеют вид (4.6) и (4.4) соответственно.

а) При сложении рядов в общей области сходимости получится ряд, главную часть которого будет составлять главная часть ряда функции f_1(z). Следовательно, для f(z)= f_1(z)pm f_2(z) точка z_0 — полюс порядка n.

б) Аналогичные рассуждения приводят к заключению, что такой же результат получится и для f(z)= f_1(z)cdot f_2(z), если limlimits_{zto z_0}f_2(z)ne0.

Если limlimits_{zto z_0}f_2(z)=0 и z_0-0(p),~ p<n для функции f_2(z), то из равенства

f_1(z)cdot f_2(z)= frac{varphi_1(z)}{(z-z_0)^n}cdot varphi_2(z)cdot (z-z_0)^p= frac{varphi(z)}{(z-z_0)^{n-p}} заключаем, что z_0-Pi(n-p).

в) Для частного frac{f_1(z)}{f_2(z)} при условии limlimits_{zto z_0} f_2(z)ne0 из равенства f(z)= frac{f_1(z)}{f_2(z)}= frac{varphi_1(z)}{(z-z_0)^n f_2(z)}= frac{varphi(z)}{(z-z_0)^n} заключаем, что z_0-Pi(n) для f(z).

Если limlimits_{zto z_0}f_2(z)=0 и z_0-0(p) для f_2(z), то, используя условие кратного нуля, из равенства

f(z)= frac{f_1(z)}{f_2(z)}= frac{varphi_1(z)}{(z-z_0)^n varphi_2(z) (z-z_0)^p}= frac{varphi(z)}{(z-z_0)^{n+p}}

заключаем, что z_0 является Pi(n+p) для frac{f_1(z)}{f_2(z)}, где n — порядок полюса функции f_1(z),~p — порядок нуля функции f_2(z) в точке z_0.

Подводя итог, запишем следующее утверждение.
Утверждение 4.4

1. Пусть точка z_0 является Pi(m) для функции f_1(z) и Pi(k) для функции f_2(z). Тогда:

а) для f_1(z)pm f_2(z) она будет Pi(n),~ n leqslant max(m,k), а при n=0 — устранимой особой точкой;

б) для f_1(z)cdot f_2(z) она является Pi(n),~ n=m+k;

в) для frac{f_1(z)}{f_2(z)} она будет Pi(n),~ n=m-k.

2. Пусть z_0 — существенно особая точка для функции f_2(z) и устранимая особая точка или полюс для функции f_1(z). Тогда z_0 — существенно особая точка для f_1(z)pm f_2(z);~ f_1(z)cdot f_2(z);~ frac{f_1(z)}{f_2(z)}.

3. Пусть точка z_0 является Pi(n) для функции f_1(z) и устранимой особой точкой для функции f_2(z). Тогда:

а) для f_1(z)pm f_2(z) она будет Pi(n);

б) для f_1(z)cdot f_2(z) она является Pi(n), если limlimits_{zto z_0}f_2(z)ne0, и Pi(n-p), если limlimits_{zto z_0}f_2(z)=0 и p — порядок нуля f_2(z) в точке z_0;

в) для frac{f_1(z)}{f_2(z)} она будет Pi(n), если limlimits_{zto z_0}f_2(z)ne0, и Pi(n+p), если limlimits_{zto z_0}f_2(z)=0 и p — порядок нуля f_2(z) в точке z_0;

4. Если точка z_0-Pi(n) для varphi(z), то она существенно особая точка для сложной функции f(varphi(z)). В этом можно убедиться, рассматривая ряды для varphi(z) и f(varphi(z)) в окрестности z_0.

Пример 4.19. Определить тип особой точки z=0 для функции f(z), если f(z)= f_1(z)+ f_2(z), где f_1(z)= frac{1}{z^2}, а функция f_2(z) определяется следующим образом:

а) f_2(z)= frac{1}{z}-frac{2}{z^2}; б) f_2(z)= frac{1}{z}-frac{1}{z^2}; в) f_2(z)=1-frac{1}{z}.

Решение

Пример 4.20. Найти особые точки функции f(z)= frac{1}{e^z-1}-frac{1}{z}. Определить их тип.

Решение

Особыми точками функции являются особые точки первого слагаемого z=2kpi i,~ kinmathbb{Z}, особая точка второго слагаемого z=0 входит в это множество. Точки z_k=2kpi i,~kne0 являются простыми нулями знаменателя и поэтому простыми полюсами первой функции; для второго слагаемого эти точки не являются особыми. Поэтому точки z_k=2kpi i,~kne0 -простые полюсы f(z) (см. п. 3 “а” утверждения 4.4).

Точка z=0 — простой полюс и для первого, и для второго слагаемого. Для f(z) — это или простой полюс, или устранимая особая точка (см. п.1 “а” утверждения 4.4). Преобразуем разность в дробь: f(z)= frac{z-e^z+1}{(e^z-1)z}. Точка z=0 является нулем второго порядка и для числителя, и для знаменателя. Следовательно, это — устранимая особая точка, в чем можно убедиться, используя определение, т.е. находя limlimits_{zto0}f(z). Действительно,

limlimits_{zto0}f(z)= limlimits_{zto0} frac{-1-z-dfrac{z^2}{2!}+ldots+z+ 1}{left(-1-z-dfrac{z^2}{2!}+ldots-1 right)!z}= limlimits_{zto0} frac{-dfrac{z^2}{2!}-dfrac{z^3}{3!}-ldots}{z^2! left(1+dfrac{z}{2!}+ldotsright)}= limlimits_{zto0} frac{-z^2! left(dfrac{1}{2!}+ dfrac{z}{3!}+ldotsright)}{ z^2!left(1+dfrac{z}{2!}+ldotsright)}=-frac{1}{2},.

Точка z=infty для данной функции является неизолированной особой точкой, так как в любой ее окрестности |z|>R содержится бесконечное множество особых точек вида z_k=2kpi i. Эта точка- предельная точка полюсов. Заметим, что для знаменателя первого слагаемого функции она — существенно особая точка.

Пример 4.21. Найти особые точки следующих функций, определить их тип:

а) f(z)= frac{z-pi}{sin^2z}cos frac{1}{z-2i}+ frac{1}{z^6+1}; б) f(z)= frac{1}{z^2-1} sin frac{pi z}{2z+1}+ frac{1}{e^z+i}.

Решение

Обозначим f_1(z) — первое слагаемое, f_2(z) — второе слагаемое функции f(z), т.е. имеем f(z)=f_1(z)+f_2(z).

а) Для f_1(z) точка z=2i является существенно особой точкой, так как это существенно особая точка для cos frac{1}{z-2i} множителя этой функции. Поэтому она — существенно особая точка для f(z) (п. 2 утверждения 4.4).

Точки z_k=kpi,~kne1 — полюсы второго порядка функции f_1(z), так как ее можно записать в виде f_1(z)=frac{varphi(z)}{sin^2z}, где varphi(z_k)ne0, а для знаменателя эти точки — нули второго порядка . Так как для f_2(z) эти точки не особые, то z_k=kpi,~znepi — полюсы второго порядка для f(z) (п. 3 утверждения 4.4).

С помощью аналогичных рассуждений получаем, что z=pi — простой полюс для f(z).

Особыми точками f_2(z) являются корни уравнения z^6+1=0, то есть z_k=exp frac{(-pi+2kpi)i}{6},~ k=0,1,ldots,5. Все они — простые нули знаменателя- функции F(z), а потому — простые полюсы для f_2(z)= frac{1}{F(z)}. Так как эти точки не являются особыми для f_1(z), то для f(z) — это простые полюсы.

Точка z=infty — неизолированная особая точка f(z).

б) Точка z=frac{-1}{2} — полюс дроби frac{pi z}{2z+1} является существенно особой точкой для sinfrac{pi z}{2z+1} (п.4 утверждения 4.4), поэтому она — существенно особая точка для f_1(z) и, следовательно, для f(z).

Точка z=1 — простой полюс для f_1(z), так как можно записать f_2(z)= frac{varphi(z)}{z-1},~ varphi(z)ne0. Поскольку z=1 не является особой точкой для f_2(z), то она — простой полюс для f(z).

Точка z=-1 — устранимая особая точка для f_1(z), так как она — простой нуль и для числителя, и для знаменателя дроби frac{sindfrac{pi z}{2z+1}}{z^2-1}. Так как z=-1 не является особой точкой для f_2(z), то она — устранимая особая точка для f(z).

Особыми точками f_2(z) являются простые нули знаменателя — корни уравнения e^z+i=0, или e^z=-i, то есть z=operatorname{Ln}(-i). Все точки

z_k=ln|-i|+ibigl(arg(-i)+2kpibigr), или z_k=i! left(2kpi-frac{pi}{2}right)!,~ kin mathbb{Z}

являются простыми полюсами для f_2(z) и, следовательно, простыми полюсами для f(z).

Точка z=infty — неизолированная особая точка f(z).

Математический форум (помощь с решением задач, обсуждение вопросов по математике).

Кнопка "Поделиться"

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

Особые точки функций комплексного переменного

Опр.
Особой точкой функции

называется точка в которой

не определена или не дифференцируема.

Опр.
Особая точка называется изолированной,
если

такая ее окрестность, в которой нет
других особых точек.

Утв. Если

изолированная особая точка
,
то в окрестности
,

раскладывается в ряд Лорана.

Классификация
особых точек

Опр1.
Особая точка называется устранимой,
если в ряде Лорана в окрестности этой
точки отсутствует главная часть.

Опр2.
Изолированная особая точка называется
полюсом, если главная часть ряда Лорана
в окрестности этой точки имеет конечное
число членов:

Число N
называется кратностью (порядком полюса).

Утв.
Если

– полюс
,
то
.

Док-во:

Опр3.
Изолированная особая точка называется
существенно особой, если главная часть
разложения в ряд Лорана в окрестности
этой точки содержит бесконечное число
членов.

Лекция 8

Связь между нулем
и полюсом

Утв1.

имеет в точке

нуль порядка n

имеет в точке

полюс порядка n.

Док-во: {}

Утв2.

имеет существенно особую точку в точке

имеет в

неизолированную особую точку ИЛИ
существенно особую точку.

Пример.

;

;

;
;

Таким образом,
получаем не изолированную особую точку.

Утв3.
Если
,

,

,

то

имеет при:

1)mn
устранимую особую точку,

2)m>n
полюс порядка n-m.

Док-во: {для
2}

;

Теорема Сохоцкого.

Если
-существенно
особая точка функции
,
то
.

Док-во:

1)

а)

-сходится
при


сходится при

б) Предположим
противное:


ограничена в
окрестности точки
.

в)

при
(т.е.

ограничена в окрестности
).

г)В круге


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

д) Из б), в), г)
следует

ограничена на всей комплексной плоскости.


е)

ограничена на С,

аналитическая, по теореме Ляувилля
противоречие.

2)

а)
имеет
не изолированную особую точку.

б)
-изолированная
особая точка


имеет изолированную
особую точку в



имеет существенно особую точку

по Утв2

имеет существенно особую точку в

по
1)


Теорема доказана.

Особые точки в
бесконечности

Утв.
Если
-изолированная
особая точка
,
то

Док-во:

Пусть
.
Раскладываем

в окрестности нуля:

.

Вычеты

Опр.
-изолированная
особая точка.

называется вычетом, где

– коэффициент при -1 степени в разложении
ряда Лорана:

Основная теорема
о вычетах.

Если G
– односвязная область, Г – замкнутый
контур, Г ограничевает G,
G
содержит конечное число изолированных
особых точек

функции
,
то

.

Док-во:


Г

G

.

.

.

Окружит каждую
особую точку

окружностью

так, чтобы внутри

не было других особых точек, и чтобы

и

не пересекались(ij).

.

Вычисление
вычетов

1.

Утв.
Если

– устранимая особая точка
,
то
(Т.к.
главная часть ряда Лорана не содержит
ни одного члена
)

2.

а) Утв.
Если
-простой
полюс
(полюс
кратности 1), то
.

Док-во:

Пример.

,


имеет простой полюс.

б) Утв.
Если
,,,,
то
.

Док-во:

-полюс
I
порядка

3.

Утв.
Если
-полюс
порядка n
,
то
.

Док-во:

Переходим к

и делим на
:

Пример1.

;

Пример2.

Лекция 9

Опр.

– изолированная особая точка
,

,
где Г- замкнутый контур.

Утв.
Если
,
то
.

Док-во:

1)
.

2) С: {}

,

при

3)

4)

Теорема.
Если
-изолированная
особая точка, кроме

имеется конечное число особых точек,
то

Док-во:

Возьмем замкнутый
контур С, охватывающий все особые точки,
кроме

;

;

;

Логарифмический
вычет.

Опр.
Логарифмическим вычетом называется:

,
если С – замкнутый контур,

– аналитическая внутри С и на нем за
исключением конечного числа особых
точек, все особые точки лежат внутри С,
все особые точки – полюсы.

Утв1.
Если
,

нуль кратности

фунции
,
то
.

Док-во:

Для функции


полюс I
порядка.

.

Утв2.
Если
-полюс
кратности n
функции
,
то
.

Док-во:

.

Принцип аргумента.

Теорема.
Логарифмический вычет функции

относительно контура С равен приращению

аргумента

при обходе контура С, деленному на
,
равно разности между числом нулей М и
числом полюсов N
функции

в облости D,
ограниченной контуром С:

Док-во:

Z
W

z w

C

1)

2) Внутри С

будет иметь конечное число нулей, т.к.
она аналитическая в замкнутой области.
В силу Утв1 и Утв2 :

3)

Теорема Руше.

ЕСЛИ G
– односвязная область, С – замкнутый
контур, ограничивающий G,

и

аналитические в G
и на С,
на
С,
на
С,

сумма кратностей всех лежащих в G
нулей функции
,

сумма кратностей всех лежащих в G
нулей функции
+,
ТО
.

Док-во:

1)

2)

3)

w

Вектор из начала
координат в точку, при такой конфигурации
образа С, ни одного оборота не совершит.

.

4)

Пример.
Найти количество нулей, которые имеет
функция

в круге
.

,

при
:


имеет нуль кратности
5

w
имеет 5 нулей.

Утв.
Если
,
то

имеет n
корней.

Док-во:

С:

имеет
нуль кратности n,
т.о.

имеет n
нулей.

Теорема.
Если
,


аналитическая в G
,
то

– аналитическая.

Док-во:

1)

  1. аналогично
    доказываем

  2. Из пунктов 1) и 2)
    следует, что для F
    выполнены условия Коши-Римана,
    следовательно F
    аналитическая.

Лекция 10

Вычисление
несобственных интегралов с помощью
вычетов

Теорема.
Если


при x=z,
-изолированная
особая точка f(z),

имеет в

нуль не ниже II
порядка,

не имеет особых точек на действительной
оси,

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

Док-во:

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

Y

R

-R R x

.

Пример.
Найти интеграл:.

,
;

Операционное
исчисление

Опр.
Функция

называется оригиналом, если:

1)

определена при
,

и
являются
кусочно-непрерывными на любом конечном
интервале,

2) при

3).

Утв.
Если
-многочлен
степени n,

то
.

Док-во:

,
по
правилу Лопиталя

;.

Опр.

называется изображением, соответствующим
оригиналу f(t),
если F(p)
– интеграл Лапласа:

;

.

Теорема.
Если f(t)
оригинал, то

изображение
,

1)
сходится
в полуплоскости
,

2)
является
в полуплоскости

аналитической функцией от p.

Док-во:

1)

,
таким образом F(p)
сходится.

2) Аналитичность
следует из теоремы, доказанной в
предыдущей лекции.

След. Если F(p)
– изображение некоторого оригинала,
то

Зам. Если
,
то F(p)
сходится равномерно.

Свойства
преобразования Лапласа
:

  1. Линейность

  1. Однородность.

.

Док-во для 2:

Теорема о
дифференцировании оригинала.

Если f(t)
– оригинал,
-оригинал,
F(p)-изображение
f(t),
,

то
.

Док-во:

.

Следствие.
Если
-оригиналы,
то
.

Док-во:

далее по индукции.

Теорема о
дифференцировании изображения.

Если
,
то
.

Теорема об
интегрировании оригинала.

Если
,
то
.

Док-во:

1) Докажем, что
-оригинал.
а)
кусочная гладкость – по свойству
интеграла.

б)
,
t>0
–очевидно.

в)

2)
.

.

Лекция 11

Теорема об
интегрировании изображения

Если f(t)
– оригинал,

оригинал, то
.

Док-во:

.

,


.

Теорема о
запаздывание

Если
-оригинал,
,

то
.

Док-во:

.

Теорема смещения

Если
,
то
.

Таблица соответствий

1.
.

2.

3.

4.

5.

6.

.

7.

.

8.

9.

10.

11.

12.

13.

14.

15.

Опр.
Сверткой функций f
и g
называется

Утв.
Если
,
g(t)
– оригиналы, то f*g(t)
– оригинал.

Док-во:

Пункты 1) и 2) в
определении оригинала очевидно выполнены.
Докажем выполнение пункта 3).

,

,
где

Теорема о свертках.

Если f(t),
g(t)
– оригиналы,
,
,
то
.

Док-во:

.

Лемма Жордана

Лемма1.
Если f(z)
– аналитическая в верхней полуплоскости,
за исключением, быть может, конечного
числа точек,



полуокружность в верхней полуплоскости

.

Док-во:

;

.

Лекция12

Лемма2.
Если f(z)
– аналитическая в левой полуплоскости,
,

то
.

Док-во:

.

Лемма3.

Если f(z)
аналитическая,
,
то.


y

R


x

Док-во:

  1. Докажем, что
    .

.

2)Если

аналогично.

3)

по Лемме 2.

4) Из пунктов
1), 2), 3) следует
.

Лемма4.
Если f(z)
аналитическая
,
,


то

Докозательство
следует из Леммы3.

Теорема об
интеграле Фурье.

Если f(t)
кусочно непрерывна и кусочно дифференцируема
на R,
то
(сходится
абсолютно).

Теорема обращения
преобразования Лапласа.

Если f(t)
– оригинал,
,
то
.

Док-во:

;

;

;

Теорема разложения.

,
для

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

Док-во:

.

Пример.

;

;

.

Лекция13

Соседние файлы в папке Лекции и семинары

  • #
  • #

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

Определение 2. Точка называется изолированной особой точкой функции , если существует окрестность этой точки с исключенной точкой , в которой аналитична, кроме самой точки .

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

Определение 3. Точка называется устранимой особой точкой , если разложение ее в ряд Лорана в окрестности этой точки не содержит главной части.

Определение 4. Точка называется полюсом кратности N функции, если в разложении ее в ряд Лорана в окрестности точки главная часть содержит конечное число членов, причем младшим отличным от нуля коэффициентом является .

Определение 5. Точка называется существенно особой точкой функции , если главная часть ее разложения в ряд Лорана в окрестности этой точки содержит бесконечное число членов.

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

1) для того, чтобы точка была устранимой особой точкой функции , необходимо и достаточно, чтобы .

2) для того, чтобы точка была полюсом кратности N функции, необходимо и достаточно, чтобы , .

3) для того, чтобы точка была существенно особой точкой функции , необходимо и достаточно, чтобы .

Полезна следующая теорема. Для того, чтобы точка была полюсом порядка N функции, нужно, чтобы она была нулем N – го порядка функции (связь между нулями и полюсами).

Пример 1. Для функции особой точкой является . Имеем есть устранимая особая точка.

Пример 2. Для функции является особой точкой. Так как – это полюс. Так как для функции т. является нулем пятого порядка, то – полюс пятого порядка функции .

Пример 3. Для функции является особой точкой. Разложение в ряд Лорана: в главной части содержит бесконечное число членов: это существенно особая точка.

Пример 4. Найти все особые точки функции и определить их характер.

Решение. Особыми точками являются точка и точки, в которых знаменатель обращается в нуль. Имеем , откуда , причем эти точки являются нулями первого порядка. Следовательно, в точках , функция имеет простые полюса. Точка не является изолированной особой точкой, так как она является пределом полюсов: : это означает, что любая окрестность точки содержит бесконечное число особых точек .

Задачи для самостоятельного решения

У нижеследующих функций найти нули и определить их порядки:

132. . 133. . 134. . 135. . 136. .

137. .

Найти порядок нуля для следующих функций:

138. . 139. . 140. .

141. .

Определить характер особой точки для следующих функций:

142. . 143. . 144. .

Найти особые точки и определить их характер у следующих функций:

145. . 146. . 147. . 148. . 149. .

150. . 151. . 152. . 153. .

< Предыдущая   Следующая >

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

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

Особые точки векторных полей на плоскости[править | править код]

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

{displaystyle {dot {x}}=Ax},

где {displaystyle x=(x_{1},x_{2})} — точка на плоскости, A — матрица 2times 2. Очевидно, точка {displaystyle x=(0,0)} в случае невырожденной матрицы A является единственной особой точкой такого уравнения.

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

Тип собственных значений Собственные значения
на комплексной плоскости
Тип особой точки Тип фазовых траекторий Вид фазовых траекторий
Чисто мнимые Центр фазплоскость.svg Центр окружности, эллипсы Phase portrait center.svg
Комплексные с отрицательной действительной частью Устфокус фазплоскость.svg Устойчивый фокус Логарифмические спирали Phase Portrait Stable Focus.svg
Комплексные с положительной действительной частью Неустфокус фазплоскость.svg Неустойчивый фокус Логарифмические спирали Phase Portrait Unstable Focus.svg
Действительные отрицательные Устузел фазплоскость.svg Устойчивый узел параболы Phase Portrait Stable Node.svg
Действительные положительные Неустузел фазплоскость.svg Неустойчивый узел параболы Phase Portrait Unstable Node.svg
Действительные разных знаков Седло фазплоскость.svg Седло гиперболы Phase Portrait Sadle.svg

См. также[править | править код]

  • Теорема Пуанкаре о векторном поле

ЛЕКЦИЯ 4

Модели, описываемые системами двух автономных дифференциальных уравнений.

Фазовая плоскость. Фазовый портрет. Метод изоклин. Главные изоклины. Устойчивость стационарного состояния. Линейные системы. Типы особых точек: узел, седло, фокус, центр. Пример: химические реакции первого порядка.

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

(4.1)

P(x,y), Q(x,y) – непрерывные функции, определенные в некоторой области G евклидовой плоскости ( x,y ‑ декартовы координаты) и имеющие в этой области непрерывные производные порядка не ниже первого.

Область G может быть как неограниченной, так и ограниченной. Если переменные x, y имеют конкретный биологический смысл (концентрации веществ, численности видов) чаще всего область G представляет собой положительный квадрант правой полуплоскости:

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

Переменные x, y во времени изменяются в соответствии с системой уравнений (4.1), так что каждому состоянию системы соответствует пара значений переменных ( x, y) .

Изображающая точка на фазовой плоскости

Обратно, каждой паре переменных ( x, y) соответствует определенное состояние системы.

Рассмотрим плоскость с осями координат, на которых отложены значения переменных x,y. Каждая точка М этой плоскости соответствует определенному состоянию системы. Такая плоскость носит название фазовой плоскости и изображает совокупность всех состояний системы. Точка М(x,y) называется изображающей или представляющей точкой.

Пусть в начальный момент времени t=t0 координаты изображающей точки М0( x( t0) , y( t0)) . В каждый следующий момент времени t изображающая точка будет смещаться в соответствии с изменениями значений переменных x( t) , y( t) . Совокупность точек М( x( t) , y(t)) на фазовой плоскости, положение которых соответствует состояниям системы в процессе изменения во времени переменных x(t), y(t) согласно уравнениям (4.1), называется фазовой траекторией.

Совокупность фазовых траекторий при различных начальных значениях переменных дает легко обозримый “портрет” системы. Построение фазового портрета позволяет сделать выводы о характере изменений переменных x, y без знания аналитических решений исходной системы уравнений (4.1).

Для изображения фазового портрета необходимо построить векторное поле направлений траекторий системы в каждой точке фазовой плоскости. Задавая приращение D t>0, получим соответствующие приращения D x и D y из выражений:

Направление вектора dy/dx в точке ( x, y) зависит от знака функций P(x, y), Q(x, y) и может быть задано таблицей:

Задача построения векторного поля упрощается, если получить выражение для фазовых траекторий в аналитическом виде. Для этого разделим второе из уравнений системы (4.1) на первое:

. (4.2)

Решение этого уравнения y = y( x, c) , или в неявном виде F( x,y) =c, где с – постоянная интегрирования, дает семейство интегральных кривых уравнения (4.2) ‑ фазовых траекторий системы (4.1) на плоскости x, y.

Для построения фазового портрета пользуются методом изоклин – на фазовой плоскости наносят линии, которые пересекают интегральные кривые под одним определенным углом. Уравнение изоклин легко получить из (4.2). Положим

где А – определенная постоянная величина. Значение А представляет собой тангенс угла наклона касательной к фазовой траектории и может принимать значения от – ¥ до + ¥ . Подставляя вместо dy/dx в (4.2) величину А получим уравнение изоклин:

. (4.3)

Уравнение (4.3) определяет в каждой точке плоскости единственную касательную к соответствующей интегральной кривой за исключением точки, где P (x,y) = 0, Q ( x,y) = 0, в которой направление касательной становится неопределенным, так как при этом становится неопределенным значение производной:

.

Эта точка является точкой пересечения всех изоклин – особой точкой. В ней одновременно обращаются в нуль производные по времени переменных x и y.

Таким образом, в особой точке скорости изменения переменных равны нулю. Следовательно, особая точка дифференциальных уравнений фазовых траекторий (4.2) соответствует стационарному состоянию системы (4.1), а ее координаты – суть стационарные значения переменных x, y.

Особый интерес представляют главные изоклины:

dy/dx=0, P ( x,y) =0 – изоклина горизонтальных касательных и

dy/dx= ¥ , Q ( x,y) =0 – изоклина вертикальных касательных.

Построив главные изоклины и найдя точку их пересечения (x,y), координаты которой удовлетворяют условиям:

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

Система (4.1) обладает столькими стационарными состояниями, сколько точек пересечения главных изоклин имеется на фазовой плоскости.

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

Рис. 4.2. Пересечение главных изоклин на фазовой плоскости.

Таким образом, фазовые траектории системы – это проекции интегральных кривых в пространстве всех трех измерений x, y, t на плоскость x, y (рис.4.3).

Рис. 4.3. Траектории системы в пространстве ( x, y, t).

Если условия теоремы Коши выполнены, то через каждую точку пространства x, y, t проходит единственная интегральная кривая. То же справедливо, благодаря автономности, для фазовых траекторий: через каждую точку фазовой плоскости проходит единственная фазовая траектория.

Устойчивость стационарного состояния

Пусть система находится в состоянии равновесия.

Тогда изображающая точка находится в одной из особых точек системы, в которых по определению:

.

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

Состояние равновесия устойчиво, если для любой заданной области отклонений от состояния равновесия ( e ) можно указать область d ( e ) , окружающую состояние равновесия и обладающую тем свойством, что ни одна траектория, которая начинается внутри области d , никогда не достигнет границы e . (рис. 4.4)

Иллюстрация к определению устойчивости области e и d на плоскости ( x,y)

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

Рассмотрим систему двух линейных уравнений:

. (4.4)

Здесь a, b, c, d – константы, x, y ‑ декартовы координаты на фазовой плоскости.

Общее решение будем искать в виде:

. (4.5)

Подставим эти выражения в (4.4) и сократим на e l t :

(4.6)

Алгебраическая система уравнений (4.6) с неизвестными A, B имеет ненулевое решение лишь в том случае, если ее определитель, составленный из коэффициентов при неизвестных, равен нулю:

.

Раскрывая этот определитель, получим характеристическое уравнение системы:

. (4.7)

Решение этого уравнения дает значения показателя l 1,2 , при которых возможны ненулевые для A и B решения уравнения (4.6). Эти значения суть

. (4.8)

Если подкоренное выражение отрицательно, то l 1,2 комплексно сопряженные числа. Предположим, что оба корня уравнения (4.7) имеют отличные от нуля действительные части и что нет кратных корней. Тогда общее решение системы (4.4) можно представить в виде линейной комбинации экспонент с показателями l 1 , l 2 :

(4.9)

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

, (4.10)

допускающее более удобное представление на фазовой плоскости по сравнению с исходной системой (4.4). Введем новые координаты ξ , η по формулам:

(4.1)

Из курса линейной алгебры известно, что в случае неравенства нулю действительных частей l 1 , l 2 исходную систему (4.4) при помощи преобразований (4.11) всегда можно преобразовать к каноническому виду (4.10) и изучать ее поведение на фазовой плоскости ξ , η . Рассмотрим различные случаи, которые могут здесь представиться.

Корни λ 1 , λ 2 – действительны и одного знака

В этом случае коэффициенты преобразования действительны, мы переходим от действительной плоскости x,y к действительной плоскости ξ, η. Разделив второе из уравнений (4.10) на первое, получим :

. (4.12)

Интегрируя это уравнение, находим :

, где . (4.13)

Условимся понимать под λ 2 корень характеристического уравнения с большим модулем, что не нарушает общности нашего рассуждения. Тогда, поскольку в рассматриваемом случае корни λ 1 , λ 2 – действительны и одного знака, a >1 , и мы имеем дело с интегральными кривыми параболического типа.

Все интегральные кривые (кроме оси η, которой соответствует ) касаются в начале координат оси ξ, которая также является интегральной кривой уравнения (4.11). Начало координат является особой точкой.

Выясним теперь направление движений изображающей точки вдоль фазовых траекторий. Если λ 1 , λ 2 – отрицательны, то, как видно из уравнений (4.10), |ξ|, |η| убывают с течением времени. Изображающая точка приближается к началу координат, никогда, однако, не достигая его. В противном случае это противоречило бы теореме Коши, которая утверждает, что через каждую точку фазовой плоскости проходит лишь одна фазовая траектория.

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

Рис. 4.5. Особая точка типа узел на плоскости канонических координат ξ, η

Состояние равновесия типа узел при λ 1 , λ 2 0 устойчиво по Ляпунову, так как изображающая точка по всем интегральным кривым движется по направлению к началу координат. Это устойчивый узел. Если же λ 1 , λ 2 > 0, то |ξ|, |η| возрастают с течением времени и изображающая точка удаляется от начала координат. В этом случае особая точка – неустойчивый узел .

На фазовой плоскости x, y общий качественный характер поведения интегральных кривых сохранится, но касательные к интегральным кривым не будут совпадать с осями координат. Угол наклона этих касательных будет определяться соотношением коэффициентов α , β , γ , δ в уравнениях (4.11).

Корни λ 1 , λ 2 – действительны и разных знаков.

Преобразование от координат x,y к координатам ξ, η опять действительное. Уравнения для канонических переменных снова имеют вид (4.10), но теперь знаки λ 1 , λ 2 различны. Уравнение фазовых траекторий имеет вид :

где , (4.14)

Интегрируя (4.14), находим

(4.15)

Это уравнение определяет семейство кривых гиперболического типа, где обе оси координат – асимптоты (при a=1 мы имели бы семейство равнобочных гипербол) . Оси координат и в этом случае являются интегральными кривыми – это будут единственные интегральные кривые, проходящие через начало координат. Каждая из них состоит из трех фазовых траекторий : из двух движений к состоянию равновесия (или от состояния равновесия) и из состояния равновесия. Все остальные интегральные кривые – суть гиперболы, не проходящие через начало координат (рис. 4.6) Такая особая точка носит название «седло ». Линии уровня вблизи горной седловины ведут себя подобно фазовым траекториям в окрестности седла.

Рис. 4.6. Особая точка типа седло на плоскости канонических координат ξ , η

Рассмотрим характер движения изображающей точки по фазовым траекториям вблизи состояния равновесия. Пусть, например, λ 1 >0 , λ 2 . Тогда изображающая точка, помещенная на оси ξ, будет удаляться от начала координат, а помещенная на оси η – будет неограниченно приближаться к началу координат , не достигая его за конечное время . Где бы ни находилась изображающая точка в начальный момент (за исключением особой точки и точек на асимптоте η =0), она в конечном счете будет удаляться от состояния равновесия, даже если в начале она движется по одной из интегральных кривых по направлению к особой точке .

Очевидно, что особая точка типа седла всегда неустойчива . Только при специально выбранных начальных условиях на асимптоте η =0 система будет приближаться к состоянию равновесия. Однако это не противоречит утверждению о неустойчивости системы. Если считать , что все начальные состояния системы на фазовой плоскости равновероятны, то вероятность такого начального состояния, которое соответствует движению по направлению к особой точке, равна нулю. Поэтому всякое реальное движение будет удалять систему от состояния равновесия. Переходя обратно к координатам x,y, мы получим ту же качественную картину характера движения траекторий вокруг начала координат.

Пограничным между рассмотренными случаями узла и седла является случай, когда один из характеристических показателей, например λ 1 , обращается в нуль, что имеет место, когда определитель системы – выражение ad-bc=0 (см. формулу 4.8 ). В этом случае коэффициенты правых частей уравнений (4.4) пропорциональны друг другу :

и система имеет своими состояниями равновесия все точки прямой :

Остальные интегральные кривые представляют собой семейство параллельных прямых с угловым коэффициентом , по которым изображающие точки либо приближаются к состоянию равновесия, либо удаляются от него в зависимости от знака второго корня характеристического уравнения λ 2 = a+d. (Рис.4. 7 ) В этом случае координаты состояния равновесия зависят от начального значения переменных.

Рис. 4.7. Фазовый портрет системы, один из характеристических корней которой равен нулю, а второй отрицателен.

В этом случае при действительных x и y мы будем иметь комплексные сопряженные ξ , η ( 4.10) . Однако , вводя еще одно промежуточное преобразование, можно и в этом случае свести рассмотрение к действительному линейному однородному преобразованию. Положим :

(4.16)

где a,b, и u,v – действительные величины. Можно показать, что преобразование от x,y к u,v является при наших предположениях действительным, линейным, однородным с детерминантом, отличным от нуля. В силу уравнений (4.10, 4.16) имеем :

(4.17)

Разделив второе из уравнений на первое , получим :

которое легче интегрируется , если перейти к полярной системе координат ( r, φ ) . После подстановки получим , откуда :

. (4.18)

Таким образом, на фазовой плоскости u, v мы имеем дело с семейством логарифмических спиралей, каждая из которых имеет асимптотическую точку в начале координат. Особая точка, которая является асимптотической точкой всех интегральных кривых, имеющих вид спиралей , вложенных друг в друга, называется фокусом ( рис.4.8 ) .

Рис. 4.8. Фазовый портрет системы в окрестности особой точки типа фокус на плоскости координат u, v .

Рассмотрим характер движения изображающей точки по фазовым траекториям. Умножая первое из уравнений (4.17) на u , а второе на v и складывая , получаем :

где

Пусть a 1 0 ( a 1 = Re λ ) . Изображающая точка тогда непрерывно приближается к началу координат, не достигая его в конечное время. Это означает, что фазовые траектории представляют собой скручивающиеся спирали и соответствуют затухающим колебаниям переменных. Это – устойчивый фокус .

В случае устойчивого фокуса, как и в случае устойчивого узла, выполнено не только условие Ляпунова, но и более жесткое требование. Именно, при любых начальных отклонениях система по прошествии времени вернется как угодно близко к положению равновесия. Такая устойчивость, при которой начальные отклонения не только не нарастают, но затухают, стремясь к нулю, называют абсолютной устойчивостью .

Если в формуле (4.18) a1 >0 , то изображающая точка удаляется от начала координат, и мы имеем дело с неустойчивым фокусом . При переходе от плоскости u,v к фазовой плоскости x , y спирали также останутся спиралями, однако будут деформированы.

Рассмотрим теперь случай, когда a 1 =0 . Фазовыми траекториями на плоскости u, v будут окружности которым на плоскости x,y соответствуют эллипсы :

Таким образом, при a1 =0 через особую точку x= 0 , y=0 не проходит ни одна интегральная кривая. Такая изолированная особая точка, вблизи которой интегральные кривые представляют собой замкнутые кривые, в частности, эллипсы, вложенные друг в друга и охватывающие особую точку, называется центром.

Таким образом, возможны шесть типов состояния равновесия в зависимости от характера корней характеристического уравнения (4.7). Вид фазовых траекторий на плоскости x, y для этих шести случаев изображен на рис. 4.9.

Рис. 4.9. Типы фазовых портретов в окрестности стационарного состояния для системы линейных уравнений (4.4).

Пять типов состояния равновесия грубые, их характер не изменяется при достаточно малых изменениях правых частей уравнений (4.4). При этом малыми должны быть изменения не только правых частей, но и их производных первого порядка. Шестое состояние равновесия – центр – негрубое. При малых изменениях параметров правой части уравнений он переходит в устойчивый или неустойчивый фокус.

Бифуркационная диаграмма

. (4.11)

Тогда характеристическое уравнение запишется в виде:

. (4.12)

Рассмотрим плоскость с прямоугольными декартовыми координатами s , D и отметим на ней области, соответствующие тому или иному типу состояния равновесия, который определяется характером корней характеристического уравнения

. (4.13)

Условием устойчивости состояния равновесия будет наличие отрицательной действительной части у l 1 и l 2 . Необходимое и достаточное условие этого – выполнение неравенств s > 0, D > 0 . На диаграмме (4.15) этому условию соответствуют точки, расположенные в первой четверти плоскости параметров. Особая точка будет фокусом, если l 1 и l 2 комплексны. Этому условию соответствуют те точки плоскости, для которых , т.е. точки между двумя ветвями параболы s 2 = 4 D . Точки полуоси s = 0, D >0, соответствуют состояниям равновесия типа центр. Аналогично, l 1 и l 2 – действительны, но разных знаков, т.е. особая точка будет седлом, если D , и т.д. В итоге мы получим диаграмму разбиения плоскости параметров s , D , на области, соответствующие различным типам состояния равновесия.

Рис. 4.10. Бифуркационная диаграмма

для системы линейных уравнений 4.4

Если коэффициенты линейной системы a, b, c, d зависят от некоторого параметра, то при изменении этого параметра будут меняться и величины s , D . При переходе через границы характер фазового портрета качественно меняется. Поэтому такие границы называются бифуркационными – по разные стороны от границы система имеет два топологически различных фазовых портрета и, соответственно два разных типа поведения.

На диаграмме видно, как могут проходить такие изменения. Если исключить особые случаи – начало координат, – то легко видеть, что седло может переходить в узел, устойчивый или неустойчивый при пересечении оси ординат. Устойчивый узел может перейти либо в седло, либо в устойчивый фокус, и т.д. Отметим, что переходы устойчивый узел – устойчивый фокус и неустойчивый узел – неустойчивый фокус не являются бифуркационными, так как топология фазового пространства при этом не меняется. Более подробно мы поговорим о топологии фазового пространства и бифуркационных переходах в лекции 6.

При бифуркационных переходах меняется характер устойчивости особой точки. Например, устойчивый фокус через центр может переходить в неустойчивый фокус. Эта бифуркация называется бифуркацией Андронова-Хопфа по именам исследовавших ее ученых. При этой бифуркации в нелинейных системах происходит рождение предельного цикла, и система становится автоколебательной (см. лекцию 8).

Пример. Система линейных химических реакций

Вещество Х притекает извне с постоянной скоростью, превращается в вещество Y и со скоростью, пропорциональной концентрации вещества Y, выводится из сферы реакции. Все реакции имеют первый порядок, за исключением притока вещества извне, имеющего нулевой порядок. Схема реакций имеет вид:

(4.14)

и описывается системой уравнений:

(4.15)

Стационарные концентрации получим, приравняв правые части нулю:

. (4.16)

Рассмотрим фазовый портрет системы. Разделим второе уравнение системы (4.16) на первое. Получим:

. (4.17)

Уравнение (4.17) определяет поведение переменных на фазовой плоскости. Построим фазовый портрет этой системы. Сначала нарисуем главные изоклины на фазовой плоскости. Уравнение изоклины вертикальных касательных:

Уравнение изоклины горизонтальных касательных:

Особая точка (стационарное состояние) лежит на пересечении главных изоклин.

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

Если x=0, то .

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

Рассмотрим ось y=0 . В месте пересечения этой оси интегральными кривыми они описываются уравнением

.

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

при .

Затем при дальнейшем увеличении тангенс угла наклона уменьшается по абсолютной величине, оставаясь отрицательным и стремится к -1 при x ® ¥ . Зная направление касательных к интегральным кривым на главных изоклинах и на осях координат, легко построить всю картину фазовых траекторий.

Рис. 4.12. Фазовый портрет системы линейных химических реакций (4.15)

Особые решения дифференциальных уравнений

Решение дифференциального уравнения

называется особым , если в каждой его точке нарушается свойство единственности, т. е. если через каждую его точку кроме этого решения проходит и другое решение, имеющее в точке ту же касательную, что и решение , но не совпадающее с ним в сколь угодно малой окрестности . График особого решения будем называть особой интегральной кривой уравнения (1). Если функция и ее частные производные и непрерывны по всем аргументам , то любое особое решение уравнения (1) удовлетворяет также уравнению

Значит, чтобы отыскать особые решения (1), надо исключить из уравнений (1) и (2).

Полученное после исключения из (1) и (2) уравнение

Часто бывает так, что распадается на несколько ветвей . Тогда нужно установить, является ли каждая в отдельности ветвь решением уравнения (1), и если является, то будет ли оно особым решением, т.е. нарушается ли единственность в каждой его точке.

Пример 1. Найти особые решения дифференциального уравнения

а) Находим p-дискриминантную кривую. В данном случае и условие (2) принимает вид , отсюда . Подставляя это выражение для в уравнение (4), получаем

Кривая (5) есть p-дискриминантная кривая уравнения (4): она состоит из одной ветви — параболы.

б) Проверяем, является ли p-дискриминантная кривая решением заданного уравнения. Подставляя (5) и ее производную в (4), убеждаемся, что есть решение уравнения (4).

в) Проверяем, является ли решение (S) особым решением уравнения (4). Для этого найдем общее решение уравнения (4). Перепишем (4) в виде . Это уравнение Клеро. Его общее решение

Выпишем условие касания двух кривых и в точке с абсциссой :

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

Полагая , находим, что условия (7) принимают вид

Подставляя в первое из равенств (8), получаем или т.е. при первое равенство выполняется тождественно, так как есть абсцисса произвольной точки.

Итак, в каждой точке кривой (5) ее касается некоторая другая кривая семейства (6), а именно та, для которой . Значит, есть особое решение уравнения (4).

г) Геометрическое истолкование.
Общее решение уравнения (4) есть семейство прямых (6), а особое решение (5) является огибающей этого семейства прямых (рис. 19).

Огибающей семейства кривых

называется такая кривая, которая в каждой своей точке касается некоторой кривой семейства (9) и каждого отрезка которой касается бесконечное множество кривых из (9). Будем говорить, что кривые и касаются в точке , если они имеют в этой точке общую касательную.

Если (9) есть общий интеграл уравнения (1), то огибающая семейства кривых (9), если она существует, будет особой интегральной кривой этого уравнения. В самом деле, в точках огибающей значения совпадают со значениями для интегральной кривой, касающейся огибающей в точке , и, следовательно, в каждой точке огибающей значения удовлетворяют уравнению , т.е. огибающая является интегральной кривой.

Далее, в каждой точке огибающей нарушена единственность, так как через точки огибающей по одному направлению проходит, по крайней мере, две интегральные кривые: сама огибающая и касающаяся ее в рассматриваемой точке интегральная кривая семейства (9). Следовательно, огибающая является особой интегральной кривой.

Из курса математического анализа известно, что огибающая входит в состав C-дискриминантной кривой (коротко СДК), определяемой системой уравнений

Некоторая ветвь СДК заведомо будет огибающей, если на ней:

1) существуют ограниченные по модулю частные производные

где и — постоянные;

Замечание. Условия 1) и 2) лишь достаточны, а потому ветви СДК, на которых нарушено одно из этих условий, тоже могут быть огибающими.

Пример 2. Найти особые решения дифференциального уравнения

а) Находим C-дискриминантную кривую. Имеем , так что отсюда . Подставляя это значение в (14), получаем откуда

Это и есть C-дискриминантная кривая: она состоит из двух прямых и .

б) Непосредственной подстановкой убеждаемся, что каждая из ветвей СДК является решением уравнения (13).

в) Докажем, что каждое из решений (15) является особым решением уравнения (13). В самом деле, так как и , то на каждой ветви СДК имеем (предполагаем, что решение уравнения (13) рассматривается на отрезке

где — область допустимых значений .

Заметим, что на любой из ветвей СДК в области 0″ png;base64,iVBORw0KGgoAAAANSUhEUgAAADAAAAAQCAMAAABncAyDAAAAM1BMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADbQS4qAAAAEHRSTlMAMRDQiiHowAFBoWFRoLFx3eb7ogAAAMZJREFUKM+1UksSwyAIVUHAX+T+p602mTYkdqZd1AUL5fk+4NzfjiQvv/QXwkz++/6kyblOYfXmMd4vNxglaF//xu0KEeJZdVYXkDFUbhaSDCDqDtDhO3ASgOypGJbMyVh4A3A8bBpQq1URM1exAEcTUHaF4R5ZzFQXDE+FuDIfET4AiqZFe+PykiQHYIbb8rAgTsAM3lvTjvc5DCVeORANFjSxbhfOqn6ux5wPICRojOf2fJ81Uscj+bmEUc5q4jKCXucmPQAaYQaRCPmIUQAAAABJRU5ErkJggg==” />, так дх что выполняется одно из условий (12). Значит, условия (11) и (12) выполняются, а, следовательно, прямые (15) являются огибающими парабол (14).

Итак, установлено, что каждое из решений (15) есть особое решение.

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

Схема (16) означает, что уравнение p-дискриминантной кривой может распадаться на три уравнения:

1) — уравнение огибающей;

2) — уравнение геометрического места точек заострения (возврата);

3) — уравнение геометрического места точек прикосновения интегральных линий, причем множитель входит в в квадрате.

Схема (17) означает, что уравнение C-дискриминантной кривой может распадаться на три уравнения:

1) — уравнение огибающей;

2) — уравнение геометрического места узловых точек, причем множитель входит в в квадрате;

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

Не обязательно, чтобы для каждой задачи все составные части и фигурировали в соотношениях (16) и (17).

Из всех геометрических мест только огибающая есть особое решение дифференциального уравнения. Отыскание огибающей упрощается тем, что в схемы (16) и (17) она входит в первой степени.

В отношении других геометрических мест (точек заострения, узловых точек и точек прикосновения) требуется дополнительный анализ в каждом конкретном случае. То обстоятельство, что некоторый множитель входит в в квадрате (и совсем не входит в ) указывает на то, что здесь может быть геометрическое место точек прикосновения интегральных линий. Аналогично, если некоторый множитель входит в в квадрате (и совсем не входит в ), то здесь может быть геометрическое место узловых точек. Наконец, если множитель входит в в первой степени, а в — в третьей, то возможно наличие геометрического места точек заострения.

Пример 3. Найти особое решение дифференциального уравнения

Решение. Особое решение, если оно существует, определяется системой

где второе уравнение (19) получено из (18) дифференцированием его по . Исключив , получим p-дискриминантную кривую , которая распадается на две ветви

Подстановкой убеждаемся, что обе функции являются решениями уравнения (18).

Чтобы установить, являются ли решения (20) и (21) особыми или нет, найдем огибающую семейства

являющегося общим интегралом для (18).

Выпишем систему для определения C-дискриминантной кривой откуда, исключая , получаем , или и , что совпадает с (20) и (21). В силу того, что на линиях (20) и (21) условия (11) и (12) выполняются, заключаем, что линии и являются огибающими, а значит (20) и (21) есть особые решения заданного уравнения.

Интегральные кривые (22) суть параболы , а линии – огибающие этого семейства парабол (рис. 20).

Пример 4. Найти особые решения дифференциального уравнения

Решение. Дифференцируем (23) по

Исключая из (23) и (24), получим . Дискриминантная кривая есть ось ординат. Она не является интегральной кривой уравнения (23), но согласно схеме (16) может быть геометрическим местом точек прикосновения интегральных кривых.

Решениями уравнения (23) являются параболы и те гладкие кривые, которые можно составить из их частей (рис. 21).

Из чертежа видно, что прямая действительно есть геометрическое место точек прикосновения интегральных кривых уравнения (23).

Пример 5. Найти особые решения дифференциального уравнения

Решение. Найдем . Исключая из системы уравнений получаем

Преобразовав уравнение (25) к виду , находим его общий интеграл .

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

Итак, из (26) и (27) имеем

Множитель входит в p-дискриминант и в C-дискриминант в первой степени и дает огибающую, т. е. функция есть особое решение дифференциального уравнения (25). Непосредственной подстановкой убеждаемся, что действительно удовлетворяет уравнению.

Уравнение , входящее во второй степени в p-дискриминант и совсем не входящее в C-дискриминант, дает место точек прикосновения .

Наконец, уравнение , входящее в C-дискриминант во второй степени и совсем не входящее в p-дискриминант, дает место узловых точек (рис.22).

Пример 6. Найти особые решения дифференциального уравнения

а) Ищем p-дискриминантную кривую. Дифференцируя (28) по , получаем , откуда

Подставляя (29) в (28), найдем уравнение :

б) Ищем общий интеграл уравнения (28). Обозначив у’ через р, перепишем (28) в виде

Дифференцируя обе части (28) по и учитывая, что , будем иметь

Приравнивая нулю первый множитель , получаем (29), а соотношение дает

Исключая параметр из уравнений (31) и (32), найдем общее решение уравнения (28):

в) Находим C-дискриминантную кривую. Дифференцируя (33) по C, будем иметь

Подставляя (34) в (33), получаем уравнение .

Согласно символическим схемам (16) и (17) заключаем, что есть огибающая семейства полукубических парабол (33), а есть геометрическое место точек заострения (множитель входит в уравнение в кубе) (рис. 23). Подстановкой в уравнение (28) убеждаемся, что есть решение, а решением не является (при уравнение (28) не имеет смысла). Таким образом, решение есть особое (огибающая семейства интегральных линий).

[spoiler title=”источники:”]

http://mathhelpplanet.com/static.php?p=osobye-resheniya-differentsialnyh-uravnenii

[/spoiler]

Добавить комментарий