Halley's method for the matrix sector function

Cetin Kaya Koc, Bertan Bakkaloglu

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

The matrix n-sector function is a generalization of the matrix sign function; it can be used to determine the number of eigenvalues of a matrix in a specific sector of the complex plane and to extract the eigenpairs belonging to this sector without explicitly computing the eigenvalues. It is known that Newton's method, which can be used for computing the matrix sign function, is not globally convergent for the matrix sector function. The only existing algorithm for computing the matrix sector function is based on the continued fraction expansion approximation to the principal nth root of an arbitrary complex matrix. In this paper, we introduce a new algorithm based on Halley's generalized iteration formula for solving nonlinear equations. It is shown that the iteration has good error propagation properties and high accuracy. Finally, we give two application examples and summarize the results of our numerical experiments comparing Newton's, the continued fraction, and Halley's method.

Original languageEnglish (US)
Pages (from-to)944-949
Number of pages6
JournalIEEE Transactions on Automatic Control
Volume40
Issue number5
DOIs
StatePublished - May 1995
Externally publishedYes

Fingerprint

Newton-Raphson method
Nonlinear equations
Experiments

ASJC Scopus subject areas

  • Control and Systems Engineering
  • Electrical and Electronic Engineering

Cite this

Halley's method for the matrix sector function. / Kaya Koc, Cetin; Bakkaloglu, Bertan.

In: IEEE Transactions on Automatic Control, Vol. 40, No. 5, 05.1995, p. 944-949.

Research output: Contribution to journalArticle

@article{852b3003524840fa9951fa080e36fd48,
title = "Halley's method for the matrix sector function",
abstract = "The matrix n-sector function is a generalization of the matrix sign function; it can be used to determine the number of eigenvalues of a matrix in a specific sector of the complex plane and to extract the eigenpairs belonging to this sector without explicitly computing the eigenvalues. It is known that Newton's method, which can be used for computing the matrix sign function, is not globally convergent for the matrix sector function. The only existing algorithm for computing the matrix sector function is based on the continued fraction expansion approximation to the principal nth root of an arbitrary complex matrix. In this paper, we introduce a new algorithm based on Halley's generalized iteration formula for solving nonlinear equations. It is shown that the iteration has good error propagation properties and high accuracy. Finally, we give two application examples and summarize the results of our numerical experiments comparing Newton's, the continued fraction, and Halley's method.",
author = "{Kaya Koc}, Cetin and Bertan Bakkaloglu",
year = "1995",
month = "5",
doi = "10.1109/9.384237",
language = "English (US)",
volume = "40",
pages = "944--949",
journal = "IEEE Transactions on Automatic Control",
issn = "0018-9286",
publisher = "Institute of Electrical and Electronics Engineers Inc.",
number = "5",

}

TY - JOUR

T1 - Halley's method for the matrix sector function

AU - Kaya Koc, Cetin

AU - Bakkaloglu, Bertan

PY - 1995/5

Y1 - 1995/5

N2 - The matrix n-sector function is a generalization of the matrix sign function; it can be used to determine the number of eigenvalues of a matrix in a specific sector of the complex plane and to extract the eigenpairs belonging to this sector without explicitly computing the eigenvalues. It is known that Newton's method, which can be used for computing the matrix sign function, is not globally convergent for the matrix sector function. The only existing algorithm for computing the matrix sector function is based on the continued fraction expansion approximation to the principal nth root of an arbitrary complex matrix. In this paper, we introduce a new algorithm based on Halley's generalized iteration formula for solving nonlinear equations. It is shown that the iteration has good error propagation properties and high accuracy. Finally, we give two application examples and summarize the results of our numerical experiments comparing Newton's, the continued fraction, and Halley's method.

AB - The matrix n-sector function is a generalization of the matrix sign function; it can be used to determine the number of eigenvalues of a matrix in a specific sector of the complex plane and to extract the eigenpairs belonging to this sector without explicitly computing the eigenvalues. It is known that Newton's method, which can be used for computing the matrix sign function, is not globally convergent for the matrix sector function. The only existing algorithm for computing the matrix sector function is based on the continued fraction expansion approximation to the principal nth root of an arbitrary complex matrix. In this paper, we introduce a new algorithm based on Halley's generalized iteration formula for solving nonlinear equations. It is shown that the iteration has good error propagation properties and high accuracy. Finally, we give two application examples and summarize the results of our numerical experiments comparing Newton's, the continued fraction, and Halley's method.

UR - http://www.scopus.com/inward/record.url?scp=0029304771&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0029304771&partnerID=8YFLogxK

U2 - 10.1109/9.384237

DO - 10.1109/9.384237

M3 - Article

VL - 40

SP - 944

EP - 949

JO - IEEE Transactions on Automatic Control

JF - IEEE Transactions on Automatic Control

SN - 0018-9286

IS - 5

ER -