function?[N_QN_RT]?=?LLL（QRM）
%UNtitleD5?Summary?of?this?function?goes?here
%???Detailed?explanation?goes?here
%?H=QRHT=N_Q*N_RT=I
N_Q=Q;
N_R=R;
T=eye（M）;%or?T=eye（size（Q））
k=2;
c=3/4;??%1/4while?k????for?p=k-1:-1:1
????????temp=N_R（pk）/N_R（pp）;
????????temp=round（temp）;
????????%if?temp>=2
????????%????temp=2;
????????%elseif?temp<=-2
????????%????temp=-2;
????????%else?temp=temp;
????????%end
????????if?temp~=0
????????????N_R（1:pk）=N_R（1:pk）-temp*N_R（1:pp）;
????????????T（:k）=T（:k）-temp*T（:p）;
????????