How to find Eigen Values without using eig() function ???

Hello

I want to calculate Eigen Values without using eig() function. I am actually generating a HDL code of this code, but when i use eig() function, I get an error of not supported function.

Please help me with this

[John D'Errico]

You seriously DON'T want to write eig yourself. Figure out why you can't use eig, rather than writing it yourself.

[Mischa Kim]

John, just out of curiosity, what exactly is wrong with Saad trying to write his/her own algorithm if he/she wishes or decided to take this path?

[John D'Errico]

If Saad lacks understanding about the basic eigenvalue algorithms that he needs to ask questions like this, then wanting to write eig from scratch is silly. There is a huge difference between professionally written code and whatever mess Saad could cobble together, so that his NEXT question would probably be, "Why does my eigenvalue code fail in this crazy way?" And that question will re-appear every time Saad would trip over some new (to him) problematic case.

Sorry but resolve the problem, don't replace high quality code with amateur code. If you do re-write eig, expect amateur results.

Answer 1

[Walter Roberson]

Saad, to generate HDL you would have to convert your algorithm to use the Fixed Point Toolbox. Have you tried that? HDL Coder would refuse to generate for floating point.

http://www.mathworks.com/help/releases/R2012b/hdlcoder/gs/example-generating-hdl-code-from-matlab-code.html

Answer 2

[KAJING TAYENG]

function [eigVal,eigVec]=spec_calculation(A)

s = size(A);

if s(1)~=s(2)

error('Error: Input must be square.')

end

I = eye(length(A));

%

syms x

eq1 = det(A-I*x) == 0 ;

eigVal = double(solve(eq1,x));

%

eigVec = zeros(s);

for i = 1:length(A)

syms y

eq2 = (A-eigVal(i)*I)*y == 0;

eigVec(:,i) = double(solve(eq2,y));

end

end