reduce_hess_tril.m

MATLAB科学计算与工程分析源代码源程序4

function [L,P]=reduce_hess_tril(H)
% 此函数用来将下Hessenberg矩阵化为下三角矩阵L
% 输出参数P为变换矩阵，即:H=L*P
% 编写日期：2007-5-13

[m,n]=size(H);
if m~=n
    error('输入的矩阵不是方阵！');
else
    P=eye(n);
    for i=1:n-1
        x=H(i,i:i+1);
        [G,y]=planerot(x');
        H(i,i:i+1)=y';
        H(i+1:n,i:i+1)=H(i+1:n,i:i+1)*G';
        P(:,i:i+1)=P(:,i:i+1)*G';
    end
    L=H;
end