givens.m

来自「使用QR分解方法计算矩阵特征值的matlab源码」· M 代码 · 共 15 行

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