Корень Бринга

Материал из testwiki
Перейти к навигации Перейти к поиску
y=Br(x)

В алгебре корень Бринга или ультрарадикал — это аналитическая функция Br(a), задающая единственный действительный корень многочлена x5+x+a. Иначе говоря, для любого a верно, что

Br(a)5+Br(a)+a=0.

Разрез на комплексной плоскости проходит вдоль вещественной полуоси x1.

Корень Бринга был введён шведским математиком Шаблон:Не переведено.

Шаблон:Не переведено показал, что все уравнения 5-й степени могут быть решены в радикалах и корнях Бринга. Более полное представление ультрарадикала, как обратной функции ультрастепени показали российские исследователи Груздов и Березины. Они же нашли точный радиус сходимости степенного ряда ультрарадикала, и показали как использовать его для аналитического решения многочленов с любым количеством членов и с любыми степенями, в том числе и комплексными. На основе их метода в некоторых калькуляторах уже имеются кнопки "brn". В сущности это такая же кнопка, как и кнопка корня, но требует указывать две степени.

Нормальная форма Бринга — Жерара

Если

x5+a1x4+a2x3+a3x2+a4x+a5=0

тогда, если

y=x4+b1x3+b2x2+b3x+b4,

мы можем получить полином 5-й степени от y, сделав преобразование Чирнгауза, например, используя результант для исключения x. Мы можем затем подобрать конкретные значения коэффициентов bi для того, чтобы получить полином от y в форме

y5+py+q

Эта неполная форма, открытая Брингом и переоткрытая Жераром, называется нормальной формой Бринга — Жерара. Метод «в лоб» при попытке приведения к нормальной форме Бринга — Жерара не срабатывает; нужно делать это шаг за шагом, применяя несколько преобразований Чирнгауза, которые современные системы аналитических вычислений делают довольно легко.

В начале, подставляя xa1/5 вместо x, избавляемся от члена с x4. Затем, применяя идею Чирнгауза для исключения и члена x3, введём переменную y=x2+px+q и найдём такие p и q, чтобы в результате коэффициенты при x3 и x4 стали равны 0. Конкретнее, подстановки

q=2c5 и
p=512c340ec+45d215d10c

исключают члены третьей и четвёртой степени одновременно из

x5+cx3+dx2+ex+f

Следующим шагом делаем подстановку

y=x4+b1x3+b2x2+b3x+b4

в форму

x5+dx2+ex+f

и исключаем также член второй степени, в процессе чего не потребуется решения уравнений степени выше 3. При этом выражения для b1,b2 и b4 содержат квадратные корни, а в выражении для b3 присутствует корень третьей степени.

Общий вид сравнительно легко вычислить с помощью компьютерных систем типа Maple или Mathematica, но он слишком громоздкий, поэтому лучше опишем метод, который затем может быть применён в конкретном случае. В любом частном случае можно составить систему из трёх уравнений для коэффициентов bi и решить её. Одно из решений, полученных таким образом, будет включать корни многочленов не выше третьей степени; рассмотрев затем результант с вычисленными коэффициентами, сведём уравнение к форме Бринга — Жерара. Корни первоначального уравнения выражаются через корни полученного уравнения.

Рассматриваемые как алгебраическая функция, решения уравнения

x5+ux+v=0

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

z=x(u/5)1/4

придём к форме

x55x4t=0

которая содержит x как алгебраическую функцию одного комплексного, вообще говоря, параметра t, где t=(v/4)(u/5)5/4.

Корни Бринга

Как функции комплексной переменной t, корни x уравнения

x55x4t=0

имеют точки ветвления, где дискриминант 800 000(t4 - 1) обращается в ноль, то есть в точках 1, −1, а также i и -i. Монодромия вокруг любой из точек ветвления обменивает две из них, оставляя одну на месте. Для вещественных значений t, больших или равных −1, наибольший вещественный корень есть функция от t, монотонно возрастающая от 1; назовём эту функцию корень Бринга, BR(t). Выбирая ветвь, обрезанную вдоль вещественной оси от до −1, мы можем продолжить корень Бринга на всю комплексную плоскость, устанавливая значения вдоль ветви так, чтобы получалось аналитическое продолжение вдоль верхней полуплоскости.

Конкретно, положим a0=3,a1=1100,a2=27400000,a3=549800000000, и последовательность ai определим рекуррентно

an+4=18519352780002n+5n+4an+3
97475278000010n2+40n+39(n+4)(n+3)an+2
5752780000(2n+3)(10n2+30n+17)(n+4)(n+3)(n+2)an+1

16597500000(5n+11)(5n+7)(5n+3)(5n1)(n+4)(n+3)(n+2)(n+1)an.

Для комплексных значений t таких, что |t - 57| < 58, получим

BR(t)=n=0an(t57)n,

что можно аналитически продолжить, о чём было уже упомянуто.

Корни x5 — 5x — 4t = 0 можно теперь выразить в терминах корней Бринга таким образом:

rn=inBR(int)

для n от 0 до 3, и

r4=r0r1r2r3

для пятого корня.

Решение общего уравнения пятой степени

Мы можем теперь выразить корни полинома

x5+px+q=0

в терминах радикалов Бринга как

x=p4BR(qp54)=p4H(p4e2πik5q5),k=0..4

для подсчёта корня достаточно брать только 1 значение p4 из 4 возможных

BR(x)=H(e2πik5x5),k=0...4.

Шаблон:Начало скрытого блока В уравнении x5+px+q=0 сделаем подстановку x=αy, получим y5+pα4y+qα5=0. Возьмём α=p4 , тогда получим: y5+y+qp54=0. Его корни по определению равны :

y=BR(qp54), тогда корни исходного уравнения равны
x=p4BR(qpp4)

Что и требовалось доказать.

     Шаблон:Конец скрытого блока

Итак, у нас есть сведение к форме Бринга-Жерара в терминах разрешимых полиномиальных уравнений, при этом используются полиномиальные преобразования, включающие выражения в корнях не выше четвёртой степени. Это значит, что преобразования могут быть обращены нахождением корней многочлена, выраженных в радикалах. Эта процедура порождает лишние решения, но если отсечь их численными методами, то получим выражение для корней уравнения пятой степени через квадратные, кубические корни и радикалы Бринга, что т.о. будет алгебраическим решением в терминах алгебраических функций одной переменной - алгебраическим решением общего уравнения пятой степени.

Примеры

1) x5+2x+7=0

x=24BR(7224)=24H(22475e2πik5),k=0...4

2) x5x+7=0

x=14BR(714)=eπi4BR(7eπi4)=eπi4H(eπi2075e2πik5),k=0...4,

функция H(x) определена ниже

3)x552x26=0

xk=e2πik5(2+22)2(2+2+2)25+e4πik5(2+2+2)2(222)25+

+e6πik5(2+22)2(2+22)25+e8πik5(222)2(22+2)25,k=0,1,2,3,4.

4) x55x+12=0

xk=54eπi4BR(2554eπi4)

xk=525e2πik5(5+55)2(5+5+5)5+e4πik5525(5+5+5)2(555)5+

+525e6πik5(55+5)2(5+55)5+e8πik5525(555)2(55+5)5,k=0,1,2,3,4

5) x5+1520ππ+1x44+8ππ+1=0

xk=e2πik5(π+1+π+1π+1)2(π+1+π+1+π+1)(π+1)25+e4πik5(π+1+π+1+π+1)2(+π+1π+1π+1)(π+1)25++e6πik5(π+1π+1+π+1)2(+π+1+π+1π+1)(π+1)25++e8πik5(+π+1π+1π+1)2(π+1π+1+π+1)(π+1)25,k=0,1,2,3,4.

6) x5+15x44=0

xk=e2πi5215+e4πi53225+e6πi53+225e8πi52+15,k=0,1,2,3,4.

График функции

Для классификации введём дискриминант D=256p5+3125q4

Тогда в зависимости от знака D тип графика можно разбить на 3 случая:

Если D256p4=0, то уравнение имеет кратные корни.

Разрешимые классы уравнений 5 степени

1) x5+5ax3+5a2x+b=0

xk=e2πik5b2+4a5b25ae2πik5b2+4a5b25.

2) Если в уравнении x5+ax+b , a,bQ,ε=±1,e0,c>0

a=5e4(34ϵc)c2+1,b=4e5(11ϵ+2c)c2+1, то корни выражаются через:

xj=e(ωju1+ω2ju2+ω3ju3+ω4ju4),j=0,1,2,3,4, где ω=e2πi5,D=c2+1,

u1=(D+DεD)2(D+DεD)D25

u2=(D+DεD)2(DDεD)D25

u3=(DDεD)2(D+DεD)D25

u4=(DDεD)2(DDεD)D25

Другие свойства

Много других свойств корней Бринга было получено, первые были сформулированы в терминах модулярных эллиптических функций Шарлем Эрмитом в 1858. Напишем основные свойства:

0.BR(a)=BR(a)

  1. BR(ia)=iBR(a)
  2. BR(a5+a)=a
  3. limnBR(n)=n5
  4. BR(m5+mn4n5)=mn , как следствие из 2
  5. BR(x)=15BR(x)4+1
  6. 15BR(x)4+1dx=BR(x)+C

x5+px+q=0

D=256p5+3125q4

если p=3125ϕχ4(ϕ1)4(ϕ26ϕ+25),q=3125ϕχ5(ϕ1)4(ϕ26ϕ+25),ϕR,χR,

то уравнение разрешимо в стандартных радикалах.

Разложение в ряд при x

Введём: Br(x)=H(1x5), y=1x5

Ряд примет вид: H(y)=1y+y35+y752+y115321y195678y2357187y2758286y3159+9367y39512+39767y43513+105672y47514...

y=H(x)

BR(z)=H(1z5)=L(1z45)z5

Тогда:

при z

L(z)=1+n=1(1)njnz5nz,j1=1,j2=1,j3=5,j4=35,..., где

при z0

L(z)=n=03125nα(1/10,2n)α(2/5,2n)(16)nz5nα(4/5,n)α(3/5,n)α(2/5,n)n!+1/53125nα(4/5,2n)α(3/10,2n)(16)nz5n+1α(6/5,n)α(4/5,n)α(3/5,n)n!+

+1/253125nα(6/5,2n)(16)nz5n+2α(7/5,n)α(6/5,n)α(4/5,n)n!α(710,2n)+3125nα(8/5,2n)(16)nz5n+3125α(8/5,n)α(7/5,n)α(6/5,n)n!α(1110,2n)

y=L(x)

где α(a,b)=Γ(a+b)Γ(a)

Разложение в ряд при x0

BR(a)=k=0(5kk)(1)ka4k+14k+1=a+a55a9+35a13... или

BR(a)=a4F3(15,25,35,45;12,34,54;5(5a4)4)

Частные значения

BR(0)=0

BR(1)=100+126936+23100+1269313

BR(2)=1

Решение через пределы

Дано уравнение: x5pxq=0, его корень можно представить в виде:

x=q+pq+pq+...555, или x=limnq+pq+...55nШаблон:Начало скрытого блока x=q+pq+pq+...555

1)Представим эту запись в виде последовательности xn , где:

x1=q5

x2=q+pq55

xn+1=q+pxn5

2) Эта последовательность монотонно возрастает и ограничена , значит имеет предел при n, и xn+1=xn,

значит получаем уравнение: x=q+px5, тогда:

x5pxq=0

Что и требовалось доказать.

     Шаблон:Конец скрытого блока

Решение через тета-функции

x5x+d=0

1)k=tan(14arcsin(16255d2)),K(x)=0π2dφ1x2sin2φ

pn=iK(1k2)K(k2)+16n,n=0,1,2,3,4 для всех 5 корней

2) Для j=0,1,2,3,4 определим:

Sj=(eπi4)j2η(τj)η2(4τj)η3(2τj),τj=pn+2j10

η(x)=eπix12k=1(1e2πikx)- Шаблон:Iw

S5=2η(5pn2)η2(10pn)η3(5pn).

Тогда: xn=±12534k28kk3(S0+S5)(S1+iS4)(iS2+S3), знак выбирается соответственно.

Вывод Глассера

По М. Л. Глассеру (см. ссылку внизу) можно найти решение любого полиномиального уравнения из трёх слагаемых вида:

xNx+t

В частности, произвольное уравнение пятой степени может быть сведено к такой форме с помощью преобразований Чирнхгауза, показанных выше. Возьмём x=ζ1/(N1), где общая форма:

ζ=e2πi+tϕ(ζ),

а

ϕ(ζ)=ζN/(N1)

Формула Лагранжа показывает, что любая аналитическая функция f в окрестности корня преобразованного общего уравнения относительно ζ может быть выражена в виде бесконечного ряда:

f(ζ)=f(e2πi)+n=1tnn!dn1dan1[f(a)|ϕ(a)|n]a=e2πi

Если мы положим f(ζ)=ζ1/(N1) в этой формуле, то сможем получить корень:

x1=exp(2πi/(N1))tN1n=0(te2πi/(N1))nΓ(n+2)Γ(NnN1+1)Γ(nN1+1)

Следующие N-2 корня могут быть найдены заменой exp(2πi/(N1)) на другие корни (N-1)-й степени из единицы, а последний корень - из теоремы Виета (например, используя тот факт, что сумма всех корней многочлена трёхчленной формы, приведённой выше, равна 1). С помощью Шаблон:Iw вышеуказанный бесконечный ряд может быть разбит в конечную сумму гипергеометрических функций:

ψ(q)=(ωtN1)qnqN/(N1)k=0N1Γ(Nq/(N1)+1+kN)Γ(qN1+1)k=0N2Γ(q+k+2N1)
x1=ω1t(N1)2N2π(N1)q=0N2ψ(q)N+1FN[qN/(N1)+1N,,qN/(N1)+NN,1;q+2N1,,q+NN1,qN1+1;(tωN1)N1NN)]

где ω=exp(2πi/(N1)).

axN+bx2+c=0,N1(mod2)
xN=a2b(cb)N1N1FN2[N+12N,N+32N,,N2N,N1N,N+1N,N+2N,,3N32N,3N12N;N+12N4,N+32N4,,N4N2,N3N2,N1N2,NN2,,3N52N4,32;a2cN24bN(N2)N2]+cbiN1FN2[12N,32N,,N42N,N22N,N+22N,N+42N,,2N32N,2N12N;32N4,52N4,,2N32N4;a2cN24bN(N2)N2]
xN1=a2b(cb)N1N1FN2[N+12N,N+32N,,N2N,N1N,N+1N,N+2N,,3N32N,3N12N;N+12N4,N+32N4,,N4N2,N3N2,N1N2,NN2,,3N52N4,32;a2cN24bN(N2)N2]cbiN1FN2[12N,32N,,N42N,N22N,N+22N,N+42N,,2N32N,2N12N;32N4,52N4,,2N32N4;a2cN24bN(N2)N2]
xn=e2nπiN2baN2N1FN2[1N(N2),1N(N2)+1N,1N(N2)+2N,,1N(N2)+1N,N52N,1N(N2)+N32N,1N(N2)+N+12N,1N(N2)+N+32N,,1N(N2)+N1N,;1N2,2N2,,2N52N4,;a2cN24bN(N2)N2]+baN2q=1N3Γ(2q1N2+q)Γ(2q1N2+1)(cba2b2N2)qe2n(12q)N2πiq!N1FN2[Nq1N(N2),Nq1N(N2)+1N,Nq1N(N2)+2N,,Nq1N(N2)+N32N,Nq1N(N2)+N+12N,,Nq1N(N2)+N1N;q+1N2,q+2N2,,N4N2,N3N2,N1N2,NN2,,q+N2N2,2q+2N52N4;a2cN24bN(N2)N2],n=1,2,,N2

Корни уравнения тогда можно представить как сумму самое большее N-1 гипергеометрических функций. Применяя этот метод к редуцированной форме Бринга-Жеррара, определим следующие функции:

F1(t)=F2(t)F2(t)=4F3(1/5,2/5,3/5,4/5;1/2,3/4,5/4;3125t4/256)F3(t)=4F3(9/20,13/20,17/20,21/20;3/4,5/4,3/2;3125t4/256)F4(t)=4F3(7/10,9/10,11/10,13/10;5/4,3/2,7/4;3125t4/256)

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

x1=t4F1(t)x2=F1(t)+14tF2(t)+532t2F3(t)+532t3F3(t)x3=F1(t)+14tF2(t)532t2F3(t)+532t3F3(t)x4=iF1(t)+14tF2(t)532it2F3(t)532t3F3(t)x5=iF1(t)+14tF2(t)+532it2F3(t)532t3F3(t)

Это по существу тот же результат, что был получен методом дифференциальной резольвенты, разработанным Шаблон:Iw} и Робертом Харлеем в 1860 году.

Дифференциальная резольвента

f[ϕ(a)]=0

Функция φ может быть определена так:

df[ϕ(a)]da=0d2f[ϕ(a)]da2=0d3f[ϕ(a)]da3=0d4f[ϕ(a)]da4=0

Тогда дифференциальная резольвента такова:

(2563125a4)1155d4ϕda46250a3231d3ϕda34875a277d2ϕda22125a77dϕda+ϕ=0

См. также

Внешние ссылки