matlab实现基于Household变换的QR分解程序

function [Q_ymt R_ymt]=QR_Household（A）
rownum_A = size（A1）;
colnum_A = size（A2）;
rank = min（rownum_Acolnum_A）;
t=eye（rank）;  %在循环中求Q时的中间变量矩阵t
for i=1:rank-1
    a=zeros（rank+1-i1）;
    for k=1:rank+1-i
        a（k1）=A（k+i-1i）;
    end
    delta = norm（a）*sign（a（11））;
    I=eye（