Метод аналитического решения (в радикалах) алгебраического уравнения n-ой степени с возвратом к корням исходного уравнения. Собственные значения для нахождения функций от матриц. Устойчивость решений линейных дифференциальных и разностных уравнений.
Лейбница Э.В. фон Чирнгауз (1651 - 1708) опубликовал в журнале "Acta Eruditorum" метод преобразования алгебраического уравнения в уравнение той же степени с меньшим числом членов. Пусть q = (-Cn) ** (1/n), тогда уравнение относительно переменной yi легко может быть решено yi = q* (COS (2* (i - 1) *PI/ n) j*SIN (2* (i - 1) *PI/ n), где q - арифметический корень n - ой степени из числа (-Cn), i - порядковый номер корня уравнения, i = 1, n; Исаак Ньютон (1643 - 1727) после безуспешных попыток точно решить уравнение пятой степени разработал приближенный метод численного определения действительного корня алгебраического уравнения произвольной степени, получивший его имя и используемый до сих пор (так называемый метод касательных Ньютона). Эйлер установил, что уравнения второй, третьей, четвертой степеней сводятся к уравнениям первой, второй и третьей степеней, которые он назвал "разрешающими уравнениями", резольвентами. Уравнение (2) может быть получено умножением исходного уравнения (1) на уравнение для корней, взятых с обратным знаком.Метод основан на последовательном получении общего алгебраического уравнения относительно квадратов независимой переменной и его Решении с последующим возвратом к корням исходного уравнения. Для решения уравнений разработанным Методом не требуется знания специальных разделов Высшей Алгебры: теорий групп Абеля, Галуа, Ли и пр. и специальной математической терминологии: полей, колец, идеалов, изоморфизмов и т.д., нужно лишь умение решать квадратные уравнения и извлекать корни n - ой степени из комплексного числа.
Введение
Проблема решения в радикалах алгебраического уравнения произвольной степени, так называемого Векового уравнения, интересовала математиков всех времен и народов. Удача Тартальи и Феррари в решении уравнений третьей и четвертой степеней внесла надежду на успехи в этом направлении и далее. Однако Решения долгое время найти не удавалось / 1/. Могу с уверенностью сказать, что все Великие математики, в течение последних пятисот лет, занимались решением уравнений высших степеней. Уравнение пятой степени решали Ньютон, Лейбниц, Лагранж, Эйлер, Гаусс, Тэйлор, Абель, Галуа, Пуанкаре, Клейн, Гильберт и многие другие (Список можно было бы еще долго продолжать). В справочниках по высшей Математике сказано, что НЕ СУЩЕСТВУЕТ решения в радикалах алгебраических уравнений выше четвертой степени / 2/. Казалось бы, не существует и решать не надо! Однако в Технике очень важно выбирать параметры Систем в соответствие с принципами Оптимальности, чтобы Объекты, описываемые системами дифференциальных или разностных уравнений, удовлетворяли заданному Критерию качества (например, минимуму потребляемой Энергии или максимальному быстродействию).
Для пояснения дальнейших рассуждений введем систему условных обозначений.
УСЛОВНЫЕ ОБОЗНАЧЕНИЯ: * - знак умножения, ** - знак возведения в степень, ABS (x) - абсолютная величина комплексной переменной x, Re x, Im x - действительная и мнимая величины комплексной переменной x соответственно, Mod x, Fi x - модуль и угол комплексной переменной x соответственно, SIN (x), COS (x) - тригонометрические функции sinx и cosx, ARCTAN (Im x, Re x) - обратная тригонометрическая функция arctg ( (Im x) / (Re x)).
SQRT (x) - операция извлечения квадратного корня из действительного числа x.
PI = 3.141592653589793 - число ?.
В 1683 году друг Г.В. Лейбница Э.В. фон Чирнгауз (1651 - 1708) опубликовал в журнале "Acta Eruditorum" метод преобразования алгебраического уравнения в уравнение той же степени с меньшим числом членов.
Чирнгауз из уравнения
(x**n) A1* (x** (n - 1)) A2* (x** (n - 2)) … An = 0, и уравнения с неопределенными коэффициентами y = B1* (x** (n - 2)) B2* (x** (n - 3)) … Bn-1, исключал x. Он полагал, что в полученном уравнении
(y**n) C1* (y** (n - 1)) C2* (y** (n - 2)) … Cn = 0, можно будет подобрать коэффициенты Bi, от которых зависят Ci, так, что все коэффициенты Ci, кроме одного, обратятся в нуль. Тогда последнее уравнение примет вид
( y**n) Cn = 0, и исходное уравнение относительно переменной x будет разрешимо в радикалах.
Отметим, что в общем случае коэффициент Cn может быть комплексной величиной, для которой, в соответствие с теорией функций комплексного переменного, существуют понятия модуля и угла вектора на комплексной плоскости. Для упрощения рассуждений будем полагать коэффициент Cn действительной величиной ( (-Cn) > 0)
Пусть q = (-Cn) ** (1/n), тогда уравнение относительно переменной yi легко может быть решено yi = q* (COS (2* (i - 1) *PI/ n) j*SIN (2* (i - 1) *PI/ n), где q - арифметический корень n - ой степени из числа (-Cn), i - порядковый номер корня уравнения, i = 1, n;
j - квадратный корень из ( - 1), мнимая величина.
Выражение COS (2* (i - 1) *PI/ n) j*SIN (2* (2* (i - 1) *PI/ n) задает корни уравнения
Последнее представляет собой выражение для суммы n членов геометрической прогрессии с основанием x.
Чирнгаузу удалось решить таким образом уравнение при n = 3, но в общем случае прием к цели не приводил. Лейбниц, которому Чирнгауз сообщил письмом в 1677 году идею метода, заметил, что ничего не получается даже для уравнения пятой степени.
Исаак Ньютон (1643 - 1727) после безуспешных попыток точно решить уравнение пятой степени разработал приближенный метод численного определения действительного корня алгебраического уравнения произвольной степени, получивший его имя и используемый до сих пор (так называемый метод касательных Ньютона). Суть метода заключается в следующем: Предположим, что действительный корень заданного алгебраического уравнения y1 находится в интервале (a, b).
Вычисляют значение алгебраической функции F (a) или F (b), ( F (a) = (a**n) A1* (a** (n - 1)) A2* (a** (n - 2)) … An ), записывают уравнение касательной в этой точке и определяют точку пересечения касательной с осью абсцисс, которой присваивают новое значение a или b.
Процесс вычислений выполняют до тех пор, пока не будет достигнута требуемая степень точности вычислений EPS (y1 = a или y1 = b в зависимости от того с какой стороны (слева или справа) решено приблизиться к корню y1).
Метод всегда сходится, но НИЧЕГО не говорит об оптимальных значениях коэффициентов уравнения, которые непосредственно связаны с параметрами Систем.
Следующий этап развития теории решения уравнений связан с творчеством Леонарда Эйлера (1707 - 1783), который, как и все предшественники, считал возможным решение уравнений любой степени.
Эйлер установил, что уравнения второй, третьей, четвертой степеней сводятся к уравнениям первой, второй и третьей степеней, которые он назвал "разрешающими уравнениями", резольвентами.
Для приведенного уравнения четвертой степени (x**4) B2* (x**2) B3*x B4 = 0, он рекомендовал подстановку x = (A** (1/ 4)) (B** (1/ 4)) (C** (1/ 4)).
Тем самым он открыл ДРУГОЙ способ решения уравнения четвертой степени, отличный от решения Феррари.
Эйлер полагал, что приведенное уравнение n-ой степени
(x**n) B2* (x** (n - 2)) B3* (x** (n - 3)) … Bn = 0, может быть решено с помощью подстановки x = (A** (1/ n)) (B** (1/ n)) … (G** (1/ n)), где число слагаемых равно (n - 1). Им использовались и другие подстановки. Однако уравнение выше четвертой степени Эйлеру решить не удалось.
При доказательстве невозможности решения уравнения пятой степени Н.Х. Абель (1802 - 1829) опирался на предложенную Эйлером подстановку x = w A* ( (v** (1/ 5)) B* ( (v** (2/ 5)) C* ( (v** (3/ 5)) D* ( (v** (4/ 5)), применив опыт великого Математика в своей работе.
Феликсом Клейном (1849 - 1925) написана монография / 3/, в которой наиболее полно показана сложность нахождения точного решения уравнения пятой степени. Книга содержит 336 страниц текста, а решения - нет! Оговорюсь сразу, что я вовсе не собираюсь принижать вклад Великих математиков в Науку, напротив, преклоняюсь перед их Волей и Настойчивостью при решении столь сложной Задачи. Они, как все лучшие представители Человечества, опережали свое Время. При отсутствии средств вычислительной техники все попытки были обречены: не было не только персональных компьютеров, но даже простых калькуляторов. Точность вычислений на логарифмической линейке для этой цели оставляла желать лучшего.
Мне удалось решить алгебраическое уравнение n - ой степени в радикалах, но Решение это - приближенное и требует вычислений с высокой степенью точности. За все надо платить, бесплатно НИЧЕГО не дается! Для определения корней уравнения не требуется знание интервала, где алгебраическая функция меняет свой знак (интервала нахождения действительного корня), что отличает разработанный Метод решения от численных методов расчета. Для определения корней уравнения не требуется знания теорий групп Абеля, Галуа, Ли и пр. и применения специальной математической терминологии: колец, полей, идеалов, изоморфизмов и т.д. Для решения алгебраического уравнения n - ой степени нужно только умение решать квадратные уравнения и извлекать корни из комплексного числа. Корни могут быть определены с любой степенью точности, если мощность персонального компьютера позволяет избежать влияния погрешностей округления на вычисления.
Отметим также, что с Решением Векового уравнения решаются Проблемы собственных значений при вычислении Функций от Матриц и Устойчивости решений линейных дифференциальных и разностных уравнений, описывающих движения сложных технических Объектов с постоянной и переменной структурой (например, вентильных преобразователей). В любом учебнике по Теории Автоматического Управления / 4/ можно прочитать: Решение линейного дифференциального уравнения устойчиво, если все корни характеристического уравнения лежат в левой полуплоскости комплексной плоскости корней. Решение разностного уравнения устойчиво, если все корни характеристического уравнения находятся внутри круга единичного радиуса на комплексной плоскости с центром в начале координат.
Оптимальное управление Системами требует отдельного рассмотрения. Скажу лишь, что Оптимальные параметры Систем могут быть достигнуты на Границе устойчивости.
Ниже приводятся СУТЬ метода Решения алгебраических уравнений и конкретные Примеры определения корней уравнений с третьей по восьмую степень включительно, доказывающие ПРАВИЛЬНОСТЬ полученных результатов и уже изложенные автором в других работах / 5, 6/.
ПОСТАНОВКА ЗАДАЧИ
Общий вид алгебраического уравнения n - ой степени
(x**n) A1* (x** (n-1)) A2* (x** (n-2)) … A (n-1) *x An = 0, (1) где n - порядок алгебраического уравнения, ___
Ai - коэффициенты уравнения, любые действительные числа, i = 1,n.
Случай комплексных коэффициентов уравнения в данной работе не рассматривается.
Поскольку Вычисления на персональном Компьютере обладают конечной точностью, целесообразно уравнение (1) нормировать по старшему коэффициенту An, чтобы не происходило переполнения разрядной сетки. Нормирующий коэффициент RCN = (ABS (An)) ** (1/n). Если n - нечетная величина, знак абсолютной величины обычно опускают. Вычисления на персональном Компьютере всегда ведутся с определенной степенью точности EPS, которая задает Критерий окончания Счета.
Критерий окончания Счета: Если алгебраическая функция, заданная уравнением (1), при вычисленном значении корня xi меньше величины ABS (EPS*An), то вычисления названного корня прекращают. Далее понижают порядок исходного уравнения до величины (n - 1), если корень xi - действительный, или до величины (n - 2), если xi принадлежит паре комплексно - сопряженных корней. Вся процедура повторяется сначала для полученного уравнения более низкого порядка до тех пор, пока не будут найдены все корни исходного уравнения (1). Если возможности Компьютера не достаточны, следует понизить степень точности EPS (в ущерб точности вычисления корней) или приобрести более мощную Персоналку. (Персоналка - персональная вычислительная Машина для каждого Пользователя)
Очевидно, что чем мощнее Компьютер, тем больше возможностей для решения уравнений более высоких Степеней n.
ЛОГИКА РАССУЖДЕНИЙ.
В общем случае, корни алгебраического уравнения отличаются друг от друга по величине. Следовательно, ВСЕГДА можно выделить в Решении наибольший по модулю (доминирующий) и наименьший корни. (Уместно оговориться сразу, что наименьший по модулю корень будет доминирующим в уравнении, обратном данному).
Попробуем последовательно возводить корни в квадрат и сравнивать их по величине между собой. После нескольких таких операций легко убедиться, что все корни уравнения для квадратов относительно переменной xc = (x** (2**J)) - ничтожно малы, кроме доминирующего корня xc1.
ВСЕ коэффициенты уравнения, кроме первых двух, будут стремиться к нулю и, следовательно, ими можно пренебречь. Тогда корень xc1 может быть найден из квадратного уравнения, а корень исходного алгебраического уравнения определится выражением x1 = (xc1** (1/ (2**J))).
Зачастую, при обеспечении заданной степени точности EPS, раньше вычисляется доминирующий корень обратного уравнения, поэтому РЕКОМЕНДУЕТСЯ определять доминирующие корни как прямого, так и обратного, уравнений.
При этом удается минимизировать затраты машинного времени и, следовательно, добиться максимальной скорости вычислений.
Уравнение (1) является частным случаем другого алгебраического уравнения n - ой степени для переменной xc = (x** (2**J)), где J - шаг преобразования, J = 1,m, m и n - любые натуральные числа.
(хс**n) B1* (xc** (n-1)) B2* (xc** (n-2)) … B (n-1) *xc Bn = 0, (2) где
Уравнение (2) может быть получено умножением исходного уравнения (1) на уравнение для корней, взятых с обратным знаком. Например, для случая n = 3 это выглядит следующим образом: ( (x**3) A1* (x**2) A2*x A3) * ( (x**3) - A1* (x**2) A2*x - A3) = 0.
Тогда относительно переменной xc = (x**2) получают уравнение (2) при J = 1
Пусть L = (2**J) - величина степени корня xc1 на J -ом шаге преобразования, xc1 = (x1**L).
Как уже отмечалось выше, на определенном шаге преобразований J все коэффициенты уравнения (2), кроме первых двух B1 и B2, становятся пренебрежительно малы и их можно отбросить. Тогда корень xc1 может быть найден из квадратного уравнения, получаемого путем отбрасывания ничтожно малых старших коэффициентов. (Не следует забывать, что исходное уравнение (1) уже нормировано по старшему коэффициенту An).
Если алгебраическая Функция при вычисленном значении корня x1 F (x1) не удовлетворяет Критерию окончания Счета, переходят к следующему шагу преобразования (J присваивают значение J 1) до тех пор, пока не будет достигнута требуемая точность вычислений EPS.
Уместно отметить, что величина xc1 может быть как действительной, так и комплексной величиной. При вычислении корня x1 следует подвергать Проверке ВСЕ КОРНИ степени L из переменной xc1: Если xc1 - комплексная величина (общий случай), тогда PI = 3.141592653589793, I2 = 1, L
Теорема: Для любого алгебраического уравнения при заданной степени точности EPS всегда существует такая величина J, при которой корень квадратного уравнения (3) совпадает с одним из корней исходного уравнения (1).
При выборе формулы расчета следует помнить, что Если I1 = 1 или I1 = 2, то вычисление xc1 осуществляется по формуле (3) для прямого уравнения (2).
Если I1 = 3 или I1 = 4, то вычисление xc1 происходит по формуле (3) для уравнения, обратного уравнению (2).
Теорема может быть доказана с помощью Метода Математической Индукции.
В заключение отметим, что в работе / 5/ коэффициенты квадратного уравнения (3) определены несколько иначе, однако корни исходного алгебраического уравнения (1) вычисляются с той же степенью точности EPS. Ввиду того, что коэффициенты Аі алгебраического уравнения (1) являются независимыми переменными, но возможны и ЧАСТНЫЕ СЛУЧАИ, указать величину J заранее не представляется возможным. Программы, используемые для проверочных расчетов, составлены автором на алгоритмическом языке FORTRAN - 90 и доказали свою высокую Эффективность.
Проверка всегда позволяет избежать Ошибок.
ПРОВЕРКА.
Дано алгебраическое уравнение третьей степени
(x**3) - 11* (x**2) - 10*x 200 = 0.
Решение: Степень точности EPS = 0.00001.
Нормирующий коэффициент для исходного уравнения RC3 = 5,8480.
Нормирующий коэффициент для исходного уравнения RC8 = 8,7261.
Коэффициент выбора формулы расчета I1 = 1.
I2 = 16.
Порядковый номер преобразования J = 5.
Корни x7, x8 - комплексно-сопряженные
Re x7 = - 15,000; Im x7 = 4,0002;
Re x8 = - 15,000; Im x8 = - 4,0002;
Корни x5, x6 - комплексно-сопряженные
Re x5 = - 2,0026; Im x5 = 2,9975;
Re x6 = - 2,0026; Im x6 = - 2,9975;
Корни x3, x4 - комплексно-сопряженные
Re x3 = - 0,9999; Im x3 = 12,000;
Re x4 = - 0,9999; Im x4 = - 12,000;
Корни x1, x2 - комплексно-сопряженные
Re x1 = - 6,9976; Im x1 = 4,9993;
Re x2 = - 6,9976; Im x2 = - 4,9993.
Вывод
Предложен Метод приближенного решения алгебраического уравнения n-ой степени в радикалах, характеризующийся простотой и доступностью для практического применения.
Метод основан на последовательном получении общего алгебраического уравнения относительно квадратов независимой переменной и его Решении с последующим возвратом к корням исходного уравнения.
Для решения уравнений разработанным Методом не требуется знания специальных разделов Высшей Алгебры: теорий групп Абеля, Галуа, Ли и пр. и специальной математической терминологии: полей, колец, идеалов, изоморфизмов и т.д., нужно лишь умение решать квадратные уравнения и извлекать корни n - ой степени из комплексного числа.
Разработанный Метод решения может быть использован при проведении оптимизационных расчетов и определении Оптимальных параметров сложных технических Систем, часть которых может быть достигнута на Границе устойчивости.
На конкретных примерах доказана ПРАВИЛЬНОСТЬ разработанного Метода и приведены Примеры решения алгебраических уравнений с третьей по восьмую степень включительно.
Решение может быть проверено Студентами, обладающими математическими знаниями в объеме институтского курса и имеющими навыки программирования на языках высокого уровня.
Список литературы
1. В.А. Никифоровский. В мире уравнений - Москва, Издательство "Наука", (Серия "История науки и техники") АКАДЕМИЯ НАУК СССР, 1987. - 176 с.
2. И.Н. Бронштейн, К.А. Семендяев. Справочник по математике для инженеров и учащихся ВТУЗОВ. - Москва, "Наука", Главная редакция физико-математической литературы, 1980. - 976 с., ил.
3. Ф. Клейн. Лекции об икосаэдре и решении уравнений пятой степени: Пер. с нем. / Под ред. А.Н. Тюрина. - Москва, "Наука", Главная редакция физико-математической литературы, 1989. - 336 с.
4. Справочник по теории автоматического управления / Под ред. А.А. Красовского. - Москва, "Наука", Главная редакция физико-математической литературы, 1987. - 712 с.
5. В.А. Будников. Классическая Алгебра. - Новосибирск, Типография ООО "ЮГУС - ПРИНТ", 2008. - 16 с.
6. В.А. Будников. Метод решения алгебраических уравнений. Решение Векового уравнения. - СТАТЬИ ДЕПОНИРОВАНЫ в "СИБКОПИРАЙТ", № 2480 от 02.09.08., Новосибирск, 2008. - 21 с.
Вы можете ЗАГРУЗИТЬ и ПОВЫСИТЬ уникальность своей работы
