function [ A ] = generate_matrix_new(y1, y2) A = []; l = length(y1); temp1 = y1; temp2 = y2; for i = 1:l t1 = circshift(temp1, i-1); t2 = circshift(temp2, i-1); if i - 1 > 0 t1(1:i - 1,1) = 0; end if i - 1 > 0 t2(1:i - 1,1) = 0; end A = [A,t1,t2]; end end