2.8. Аналитические функции
1. Рассмотрим последовательность
S0, S1, S2, … , Sm, ... (1)
тензоров одинаковой валентности. Как и при изучении «обычных» (числовых) последовательностей, каждой такой последовательности можно сопоставить ряд. Действительно, , пусть Т0 = S0, T1 = S1 — S0, T2 = S2 — S1 и т. д.; тогда (1) — последовательность «частичных сумм» (отрезков) ряда Т0 + Т1 + + Т2 +... . Обратно, каждому ряду T0 + T1 + T2 + … из данных тензоров одинаковой валентности можно сопоставить последовательность (1), полагая Sm = å{Тk: k = 0, 1, 2, …, m}. В результате вопрос о сумме ряда сводится к вопросу о пределе этой последовательности.
Пусть, например, (1) состоит из двухвалентных тензоров. Тогда с фиксированием базиса для E и, следовательно, для E Ä E имеем Sm = msijeiej. Естественно считать, что последовательность (1) сходится, когда сходится каждая из последовательностей 1sij, 2sij, 3sij, ... , образуемых компонентами тензоров (1). Тензор S = sijeiej с компонентами sij = lim {msi j : m à ¥} при этом называется пределом последовательности (1). С помощью известных предложений о пределах суммы и произведения на число сходящихся числовых последовательностей нетрудно показать, что будучи сходящейся в каком-либо базисе e1, e2, e3 Î E последовательность (1) сходится и в любом другом базисе, причем к одному и тому же предельному тензору S, который, таким образом, от выбора базиса не зависит. Но при таком «покомпонентном» определении предела последовательности (1) некоторые важные вопросы остаются в стороне. Более общий подход опирается на понятия нормы.
Напомним, что нормой (длиной) вектора х Î E называется число |х| = (х × х)1/2. Из этого определения и свойств метрической формы пространства E (которое по условию собственно евклидово) вытекает, что 1) |х| ³ 0, причем | х | = 0 <=> х = 0, 2) |cx| = |c||x|, 3) |x + y| £ |x| + |y| («неравенство треугольника»). Аналогично, нормой векторов пространства E Ä E— двухвалентных тензоров над E — называются значения отображения Т -> |Т|, Т Î E Ä E в R (поле скаляров), для которого 1) |Т|>0 и |T| = 0 <=> T = 0, 2) |cT| = |c||T|, 3) |Tu + Tv| £ |Tu| + |T| и, кроме того, 4) |Тu × Тv| £ |Тu| |Тv| при любых Т, Тu, Тv Î E Ä E (дополнительное по сравнению с имевшими место в предыдущем случае условие 4) обеспечивает непрерывность композиции Тu, Tv -> Tu × Tw в топологии, которую порождает норма на E Ä E.
В соответствии с изложенным в п. 2 2.6 для каждого Т Î E Ä E скалярный инвариант sp (T × T*) = tijtji — неотрицательное число, причем sp (T × T*) = 0 <-->Т = 0. Можно показать, что при |T| º (sp (T × T*))1/2 вместе с первым выполняются и все остальные из перечисленных выше условий 1) —4). Таким же образом обстоит дело и при отождествлении | Т | с | lм |, где lм — максимальное по модулю из собственных чисел тензора Т.
Последовательность (1) называется сходящейся, если существует такой тензор S Î E Ä E, что lim|Sm — S|=0 при m ® ¥. Если последовательность (1) сходится в смысле этого определения, то она сходится и покомпонентно, равно как и наоборот.
Выше мы считали, что (1) состоит из двухвалентных тензоров. Но с очевидными изменениями все это остается в силе и для последовательности из тензоров любой данной валентности.
2. Пусть теперь С0, С1, С2, ... —какие-либо тензоры валентности .соответственно 2, 4, 6 и т. д. С заданием такой бесконечной системы тензоров над E каждому Т Î E Ä E можно сопоставить ряд
С0 + С1:Т + С2:Т):Т + С3:Т):Т):Т + .. . (2)
(двоеточие здесь обозначает двухкратное „свертывание"; так
C1 : T = (cijkleiejekel) : (tpqepeq) = cijklpq(eiejekel × ep) × eq = cijkltlkeiej,
аналогично C2:T):T = cijkleiejekel и т. д.). По определению ряд (2) сходится, если сходится последовательность С0, С0 + С1 : Т, Со + С1 : T + C2 : Т): Т, ... , образуемая отрезками ряда. Пусть F(T)—сумма ряда, т. е. предел этой последовательности, и W Ì E Ä E — подмножество, образуемое всеми такими Т Î E Ä E, при которых этот предел существует. Иными словами, W — область сходимости ряда (2). Очевидно, что W всегда не пустое множество, ибо (2) заведомо сходится при Т = 0. Тем самым с заданием ряда (2) определяется тензорная функция — отображение по правилу
F: T -> F(T) = С0 + С1:Т + С2:Т):Т + С3:Т):Т):Т + .. ., T Î W, (3)
некоторого подмножества W Ì E Ä E в E Ä E. Свойства этой функции во многом аналогичны свойствам обычных (скалярных) аналитических функций. Факты же, в которых проявляется специфически тензорная природа аргумента и значений отображения (3), связаны с понятием группы симметрии.
Как уже упоминалось (п.1), норма на E Ä E допускает различные конкретизации. Но в любом случае значения нормы |Т| — некоторый скалярный инвариант тензоров T Î E Ä E. В соответствии с доказанным в п.1 2.7 это означает, что при А Î О(E) преобразование А2: Т -> А2(Т) = А* × Т × А всегда сохраняет норму: |А* × Т × А| = |Т| для каждого Т Î E Ä E и А Î О(E). Отсюда следует, что при любом Т Î W и любом А Î О(E) вместе с (2) сходится и ряд А* × С0 × А + A* × (C1 : Т) × А + ... . Сумма этого ряда есть образ A2(F(Т)) = A* × F(T) × A суммы F(T) исходного ряда при автоморфизме А2 пространства E Ä E, т. е.
A* × F(T) × A = A* × C0 × A + A* × (C1 : T) × A + . . . . (4)
Пусть G0, G1, G2, ... — группы симметрии тензоров соответственно С0, C1, С2, ... . В соответствии с определениями п. 2 2.7 A Î G0 <=> A2(C0) = A* × C0 × A = C0. Нетрудно видеть также, что A Î G1 <=> A* × (C1 : T) × A = C1 :(А* × Т × А) при любом T Î E Ä E, A Î G2 <=> A* × (C2 : T) : T) × A = C2 : (A* × T × A) : (A* × T × A) и т. д. Поэтому всегда, когда A Î G0 ⋂ G1 ⋂ G2 ⋂ . . . для каждого Т Î E, вместе с (2) и (4) сходится и ряд, получающийся заменой Т в (2) тензором А* × Т × А, причем F (А* × Т × А) = = A* × F(T) × A. С определенной оговоркой справедливо и обратное, а именно если для некоторого А Î О(E) из T Î W всегда вытекает А* × Т × А Î W и F(A* × T × A) = A* × F(T) × A, то необходимо A Î G0 ⋂ G1 ⋂ G2 ⋂ . . . = ∩{Gm: m = 0, 1, . . . , ¥}. Отсюда видно, во-первых, что группу симметрии функции (3) можно определить непосредственно по образцу определений п.3 2.7:
A Î GF, A Î O(E), F(A* × T × -A) = A* × F(T) × A, " T Î E Ä E. (5)
и, кроме того, что это равносильно определению, по которому
GF = G0 ⋂ G1 ⋂ G2 ⋂ . . . = ∩{Gm: m = 0, 1, . . . , ¥}.
3. Когда последовательность С0, С1 С2, ... состоит из произвольных данных тензоров (валентности 2, 4, 6 и т. д.), отображение, которое описанным образом определяется с помощью ряда (2), т. е. действует по правилу (3), — общего вида отображение из E Ä E в E Ä E, аналитическое в некоторой окрестности нулевого элемента этого пространства. Термин „аналитическое отображение" здесь вполне оправдывается, например, уже тем, что, как вытекает из (3), компоненты тензора F(T) Î E Ä E в любом базисе суть аналитические функции компонент тензора-аргумента в том же базисе для E Ä E. В соответствии с (6), с другой стороны, группой симметрии этой тензорной функции может служить любая из подгрупп полной ортогональной групп, представляющих группы симметрии тензоров четной валентности.
Однако обычно название „аналитическая" (тензорная функция), по крайней мере для функций интересующего нас сейчас вида (с аргументом и значениями — двухвалентными тензорами над E), употребляется в существенно более узком смысле. Именно, допустим, что (3)—изотропная функция, т. е. GF = O(E). Тогда и только тогда, как вытекает из (6), изотропными должны быть и все тензоры С0, С1, С2, ... —„коэффициенты" ряда (3). На основании доказанного в п. 4 2.7 это означает в свою очередь, что С0, C1, C2, ... в данном случае представляются линейными комбинациями изомеров тензоров соответственно I, II, III (= I Ä I Ä I) и т. д. В частности, для любого целого m ³ 0 возьмем в качестве Сm тот изомер тензора сm II. .. I (m – раз), для которого
Cm : T) : T) : . . . ) : T = cm T × T × . . . × T, T Î E Ä E,cm Î R (7)
(несложно проверить, что такой изомер существует при любом m ³ 0). С так выбранными „коэффициентами” Cm ряд в (3) записывается согласно (7) в виде
А(Е) = с0I + с1T + с2T × T + . . . , (8)
где с0, с1, с2, ... Î R.
4. Пусть z Î С — произвольное комплексное число. Как и в 2.5, ряду (8) естественно сопоставить скалярный степенной ряд
F(z) = c0 + c1z + c2z2 . . . . (8’)
Произвольному Т Î E Ä E при этом сопоставляется произвольное z Î С, а тензору I — число 1. Наоборот, каждому ряду (8') можно поставить в соответствие определенный степенной ряд с членами из E Ä E, а именно ряд (8), заменяя скалярные степени zm тензорами Т × Т × ... × Т. При том же „законе сопоставления" ряду
c0 + c1(z – b) + c2(z – b)2 + . . . (9)
соответствует ряд
c0I + c1(T – bI) + c2(T – bI) × (T – bI) + . . . . (9’)
Целесообразность такого сопоставления степенных рядов на комплексной плоскости и рядов с членами из E Ä E определяется в первую очередь тем обстоятельством, что хорошо изученные факты о сходимости ряда (9) позволяют сделать некоторые выводы и о сходимости соответствующего тензорного ряда. Основной из этих выводов состоит в следующем: ряд (9) сходится при | Т — bI | < r и расходится при | Т — bI | > r, где r — радиус круга сходимости ряда (9). Если воспользоваться одной из указанных в п. 1 конкретных норм в E Ä E, а именно |Т| = |lМ|, то это же самое можно сформулировать и так: ряд (9') сходится при каждом Т Î E Ä E, характеристические числа которого лежат внутри круга сходимости ряда (9), и расходится для таких Т Î E Ä E, хотя бы одно из характеристических чисел каждого из которых лежит вне этого круга. Для случая | Т — bI | = r ответ также зависит от свойств ряда (9); когда последний сходится (расходится) на границе круга \ z — b | £ r, соответственно при | Т — bI | = r сходится (расходится) и ряд (9').
Изучение тензорных степенных рядов в свою очередь связано с тем, что по аналогии с ситуацией в классическом анализе с помощью таких рядов определяется некоторый класс „элементарных" тензорных функций (отображений E Ä E в E Ä E).
Напомним, например, что по определению экспоненциальной функции z ® еz
ez = 1 + z + z2/2! + z3/3! + ... = l + å{zm/m!: m = 1, 2, . . . , ¥}. (10)
Этот ряд сходится при любом конечном значении z Î C, т. е. еz представляет собой „целую функцию". Целой функцией является и
az = 1 + (ln a) z + (ln a)2 z2 / 2! + (ln a)3z3 / 3! + . . . = 1 + å{(ln a)m zm / m! . (11)
Для функций ln z, (1 + z)a
ln z = z – 1 – (z – 1)2 / 2 + (z – 1)3 / 3 - . . . = å{(-1)m-1 (z – 1)m / m : m = 1, 2, . . . , ¥}, (12)
(1 + z)a = 1 + az + a (a - 1)z2 / 2! + a (a -1) (a - 2) z3 / 3! + . . . , (13)
причем, с некоторыми оговорками относительно значений z на границе круга сходимости, в первом случае |z — 1| < 1, во втором (при a > 0 и не целом) |zj < 1. В соответствии с этим „элементарные" тензорные функции еT, аT, ln Т, =(I + T)a определяются следующим образом:
eT = I + T + ½! T × T + 1/3! T × T × T + . . . , (10’)
aT = I + (ln a) T + (ln a)2 T × T / 2! + (ln a)3 T × T × T / 3! + . . . , (11’)
ln T = T – I – (T – I) × (T – I) / 2 + (T – I) × (T – I) × (T – I) / 3 - . . . , (12’)
(I + T)a = I + aT + a (a - 1) T × T / 2! + a (a - 1) (a - 2) T × T × T / 3! + . . . , (13’)
причем, в соответствии с изложенным выше критерием сходимости тензорных рядов, функции еT и аT определены своими рядами для любого Т Î E Ä E, функция In T — для | Т — I | < 1, а функция (I + T)a для |Т| < 1 (при не целом a > 0).
Существенно, что все эти функции — изотропные отображения E Ä E (или некоторого W Ì E Ä E) в E Ä E. Целесообразность такого определения основных тензорных функций связана с тем, что им при этом оказываются присущими многие (хотя и не все) свойства соответствующих функций числового -аргумента.
Так, для еT имеем
eT × е-T = I
и вообще для любых таких Т1, Т2 Î E Ä E, для которых T1 × T2 ¹ T2 × T1, будет
eT(1) × eT(2) = eT(1) + T(2)
Функции еT и ln T — взаимно обратные, т. е. ln eT = T, elnT=T и т. д. Функция еT, как будет видно далее (см. п. 6 ниже), позволяет связать ортогональные и кососимметрические тензоры.
5. Пусть Ps: T ® Ps(T) — „полиномиальное" изотропное отображение E Ä E в E Ä E, т. е.
Ps(T) = c0 I + c1 T + c2 T × T + . . . + cs T × T × . . . × T. (14)
С помощью формулы Гамильтона—Кэли (2.3) каждая степень Т × Т ×... × Т при целом m ³ 3 выражается через Т × Т, Т, I, и потому многочлен (14) при любом целом s сводится к многочлену, содержащему степени Т × Т ×... × Т не выше второй (и со скалярными коэффициентами, которые суть функции характеристических чисел или, что то же самое, — полной системы скалярных инвариантов тензора-аргумента Т). Но ряды (10')— (13') также содержат такие степени тензора-аргумента Т (или Т — bI), что дает основание думать, что и они в области сходимости сводятся к конечному тензорному полиному упомянутого типа. Однако, благодаря наличию в правых частях (10')—(13') как угодно больших степеней, проверить это с помощью непосредственно формулы Гамильтона—Кэли затруднительно. Более удобный путь — использование так называемой формулы Лагранжа—Сильвестра, которая одновременно доставляет и более общий, чем степенные ряды, алгоритм сопоставления функциям z ® f(z) тензорных функций.
Пусть z -> f{z) — произвольная комплекснозначная функция комплексного аргумента. Этой функции и каждому тензору Т из некоторого подмножества Wf Ì E Ä E сопоставим полином q(z) такой, что
q(lk) = f(lk), k = 1, 2, 3, (15)
(где по-прежнему l1, l2, lз — собственные числа тензора T Î Wf). Условимся, что при l1 ¹ l2 ¹ lз ¹l1 этот полином квадратичный: q(z) = bo + b1z + b2z2, для каждого Т, для которого два из чисел l1, l2, lз совпадают (но не равны третьему) будем брать q(z) = b0 + b1z и, наконец, при l1 = l2 = lз q(z) = b0. Иными словами, сопоставляя данной f(z) и каждому T Î Wf полином q(z) = b0 + b1 z + b2 z2, удовлетворяющий (15), для осесимметричного Т (с двумя совпадающими собственными числами) считаем b2 = 0, а для шарового Т — b2 = b1 = 0. При этом q(z) для каждого T Î Wf определяется из (15) единственным образом.
Пусть р(z) — минимальный многочлен тензора T Î Wf (2.5). Очевидно, что с приближением z к каждой из точек l1, l2, lз комплексной плоскости С числитель и знаменатель отношения
(f(z) – q(z) / p(z) (16)
стремятся к нулю. Будем считать, что при z -> lk для каждого k = 1, 2, 3 это отношение остается ограниченным (это условие, например, выполняется, когда f(z) в каждой из точек l1, l2, lз — регулярная функция). Тогда
f(z) = q(z) + p(z) w(z) (17)
где |w(z)| < ¥ при z = l1, l2, lз. Так как p(Т) = 0, формальная замена скалярного аргумента z в (17) тензорным Т приводит к определению тензорной функции f: Т à f(Т) Î E Ä E, соответствующей скалярной функции f: z à f(z), причем
f(T) = q(T) = b0 + b1T + b2T × T. (18)
Условие ограниченности отношения (16) при z ®lk фактически характеризует „область определения" Wf функции f(Т): для каждой данной f{z) и данного Т Î E Ä E (16) — вполне определенная функция от z и, перебирая различные Т Î E Ä E, можно определить подмножество Wf Ì E Ä E, на котором выполняется упомянутое условие ограниченности отношения (16).
Полином g(z), удовлетворяющий условиям (15), есть полином Лагранжа, который для рассмотренных выше случаев имеет вид:
a) q(z) = [(z - l2)(z - l3) / (l1 - l2)( l1 - l3) ] f(l1) + [(z - l3)(z - l1) / (l2 - l3)( l2 - l1) ] f(l2) +
+ [(z - l1)(z - l2) / (l3 - l1)( l3 - l2) ] f(l3) при l1 ¹ l2 ¹ lз;
б) q(z) = [(z - l2) / (l1 - l2)] f(l1) + [(z - l1) / (l2 - l1)] f(l2) при l2 = l3 ¹ l1;
в) q(z) = const = f(l1) (при l1=l2 = l3),
откуда и из определения (18)
f(T) = q(z) = [(T - l2I) × (T - l3I) / (l1 - l2)( l1 - l3) ] f(l1) + [(T - l3I) × (T - l1I) / (l2 - l3)( l2 - l1) ] f(l2) +
+ [(T - l1I) × (T - l2I) / (l3 - l1)( l3 - l2) ] f(l3) при l1 ¹ l2 ¹ lз; (19)
и соответственно для случаев б) и в).
Запишем (19) в виде
F(t) = - 1 / [(l1 - l2)( l2 - l3) ( l3 - l13) {[l2l3f(l1)(l2 - l3) + l3l1f(l2)(l3 - l1) + l1l2f(l3)(l1 - l2)] I –
[f(l1)( l22 + l32) + f(l2)( l32 + l12) f(l3)( l12 + l22)] T +
[f(l1)( l2 - l3) + f(l2)( l3 - l1) + f(l3)( l1 - l2)] T × T} (19’)
и обозначим a = l3 — l2, так что l3 = l2 + a. Внося это выражение в (19') и переходя к пределу при a -> 0, получим
f(T) = 1/(l1 - l2)2 {[l22f(l1) - l1l2(l1 - l2)f’(l2) – (2l1l2 - l12)f(l2)] I –
[2l2(f(l1) – f(l2) – (l12 - l22)f’(l2)] T –
[f’(l2)( l1 - l2) + f’(l2) – f(l1)] T × T}. (20)
Но при l2 = l3, как нетрудно проверить, справедливы выражения
l1 = 1/3 J1 ± 2 (1/9 J12 – 1/3 J2)1/2,
l2 = l3 = 1/3 J1 ± (1/9 J12 – 1/3 J2)1/2,
откуда
J1 = l1 + 2l2, (J12 – 3 J2)1/2 = ±(l1 - l2).
С учетом этого из (1) 2.5 и (20) получаем
f(T) = [(T - l2I)f(l1) – (T - l1I)f(l2)], (21)
т. е. путем указанного предельного перехода из (19) для случая l2 = l3 получается точно такое же выражение, как и из определения (18) при выборе соответствующего (случай б)) полинома q(z). Необходимо, однако, чтобы в точке z = l2 = l3 выполнялось условие \f{z)\ < ¥, иначе выражение (20) теряет смысл. Аналогичным образом обстоит дело и в случае l1 = l2 = l3; здесь (случай в)) f(Т) = f(l1)I, с другой стороны, это же выражение для f(Т) получается предельным переходом из (19) (или (19')) лишь если существуют и ограничены f’(l1) и f’’(l1). Это означает, что непрерывность тензорной функции f(Т) (эквивалентная, как легко видеть, непрерывности функции f(z)) недостаточна для непрерывности коэффициентов b0, b1, b2 в „трехчленном" представлении (18) функии f(Т). Точнее, коэффициенты b0, b1, b2 — непрерывные функции полной системы скалярных инвариантов тензора аргумента Т тогда и только тогда, когда тензорная функция f(Т) дважды дифференцируема (или, что то же самое, когда дважды дифференцируема соответствующая скалярная функция f(z)).
В заключение этого пункта подчеркнем еще раз, что изложенным выше способом (определение (18)) функциям f(z) ставятся в соответствие изотропные отображения подмножеств множества E Ä E в E Ä E.
6. Вернемся к изучению, экспоненциальной тензорной функции (п. 4): T -> eT, T Î E Ä E. Пусть S(E) Ì E Ä E - множество кососимметричных двухвалентных тензоров над E: T Î S(E) <=> T* = — T; O’(E), O’’(E) Ì O(E) — множества всех ортогональных преобразований пространства E соответственно первого и второго рода (п. 3, 2.6).
Покажем, что тензорная функция еT взаимно однозначно отображает множества S(E) и О'(E) (или S(E) и О"(E)) друг на друга.
Действительно, если Т Î S(E), то Т* = — Т и для тензора A = еT имеем А* = (еT)* = е-T = (еT)-1 = A-1, т. е. A* = А-1, откуда ((1) 2.6) А Î О(E). Но раз Т кососимметричен, то его собственные числа суть 0, ij, -ij, где j2 = J2 (п.2, 2.4). С учетом этого из формулы Лагранжа (19) получаем
A = eT = I + sin j / j T + (1 – cos j) / j2 T × T, (22)
причем собственные числа тензора А, согласно (15), суть 1, eij, e-ij и тем самым A = еT — ортогональное преобразование первого рода (п. 4, 2.6). Так как А × А = е2Т, из (22) получаем
A × A = I + (sin 2j / j) T + [(1 – cos 2j) / j2] T × T.
Исключая с помощью этого уравнения Т×Т из (22) и разрешая полученное уравнение относительно Т, приходим к выражению
T = -(j / 2sin j)[(1 + 2cos j) I - 2(1 + cos j) A + A × A] (23)
для функции, обратной к экспоненциальной (22), т. е. к выражению для логарифма Т = ln А (ортогонального преобразования первого рода: A Î О'(E); разумеется, (23) можно получить из формулы Лагранжа (19) точно так же, как и (22)). Осталось лишь проверить, что тензор Т, определенный равенством (23), кососимметричен, если А Î О'(E). Для этого воспользуемся диадным представлением (5) 2.6 тензоров A Î O'(E): А = аа + (bb + cc) cos j + (bc — cb) sin j, где, напомним, а — собственный вектор тензора А, соответствующий собственному значению +1, и векторы а, b, с образуют ортонормированный базис пространства E. Подставляя указанное диадное представление в (23), получаем, что T = j (cb — bc) и поэтому T* = j (bc — cb) = -T, т. е. Т Î S(E).
Таким образом, Т -> А = еT, Т Î S(E) есть отображение множества S(E) на О'(E), а А ->Т.= lnА, А Î О(E) —обратное этому отображение О(E) на S(E). Поскольку каждый тензор A Î O"(E) представим в виде А = — I ×А1 = — А1, A1 Î О'(E£) (п. 3, 2.6), очевидно, что Т -> -еT, Т Î S(E) —(взаимно однозначное) отображение множества S(E) на О"(E).