Предложен интегро-алгоритмический метод аппроксимации и итерационной коррекции для вычисления с высокой точностью логарифмов матрицы. Метод основан на применении линейных многошаговых формул численного интегрирования разностного типа, а также разностно-дифференциальных формул Обрешкова с учетом высших производных. В нем обобщены известные алгоритмы логарифмирования на основе формул Паде. В результате итерационных коррекций повышается их порядок и точность. Предложенные соотношения и программные решения позволяют определять необходимые параметры для организации процессов вычисления логарифмов матрицы с произвольно заданной высокой точностью.
Запропоновано інтегро-алгоритмічний метод апроксимації та ітераційної корекції для обчислення з високою точністю логарифмів матриці. Метод базований на застосуванні лінійних багатокрокових формул чисельного інтегрування різницевого типу, а також різницево-диференціальних формул Обрешкова з урахуванням вищих похідних. В ньому узагальнено відомі алгоритми логарифмування на основі формул Паде. В результаті використання ітераційної корекції підвищується їх порядок і точність. Запропоновані співвідношення і програмні рішення дозволяють визначати необхідні параметри для організації процесів обчислення логарифмів матриці з довільно заданою високою точністю.
The integro-algorithmic method of approximation and iteration correction for high-accuracy computation ofmatrix logarithms is proposed. Themethod is based on the use of linearmultistep formulas of numerical integration of the difference type as well as the Obreshkov difference-differential formulas with allowance for higher derivatives. Due to iterations in this case there is not a necessity to choose a high-fidelity primary approximating formula and a basic criterion is a receipt of high-rate of convergence. A method also summarizes the well-known algorithms of taking the logarithm based on the Pade formulas and increases their order and accuracy due to the additional use of iterative correction. The proposed relations and program solutions permit determining necessary parameters for organizing the processes of matrix logarithms computation with arbitrary preset high accuracy.
