user3553714
user3553714

Reputation: 3

How do I reduce a symmetric matrix to a tridiagonal matrix with the same eigenvalues?

I have a symmetric matrix found from the following commands:

a=rand(n);
A=triu(a)+triu(a,1)'

I was wondering how I reduce this symmetric matrix into a tridiagonal matrix?

Any help is appreciated. Thank you very much.

Upvotes: 0

Views: 1405

Answers (1)

riklund
riklund

Reputation: 1071

What you want is to find the upper Hessenberg form of your matrix. For a symmetric matrix, this is tridiagonal.

Use the command hess to do this:

B=hess(A); 

Upvotes: 1

Related Questions