How to find Eigen Values without using eig() function ???
18 次查看(过去 30 天)
显示 更早的评论
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
3 个评论
Mischa Kim
2014-4-18
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
2014-4-29
编辑:John D'Errico
2014-4-29
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.
回答(2 个)
Walter Roberson
2014-4-29
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.
0 个评论
KAJING TAYENG
2020-11-4
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
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!