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（