givens.m

使用QR分解方法计算矩阵特征值的matlab源码

function [c,s]=givens(A)
%Givens变换得到变换矩阵的参数
if A(2)==0
    c=1;
    s=0;
elseif abs(A(2))>=abs(A(1))
    t=A(1)/A(2);
    s=(1+t^2)^(-1/2);
    c=s*t;
else
    t=A(2)/A(1);
    c=(1+t^2)^(-1/2);
    s=c*t;
end
end