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
3 个评论
John D'Errico
2014-4-18
You seriously DON'T want to write eig yourself. Figure out why you can't use eig, rather than writing it yourself.
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
0 个投票
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.
KAJING TAYENG
2020-11-4
0 个投票
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
类别
在 帮助中心 和 File Exchange 中查找有关 Linear Algebra 的更多信息
2020-11-4
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
