clear all; close all; clc;
v=zeros(1,8);
v(2:2:end)=1;
A=[1:8;1:2:15;2:2:16;v;2*v;2+2*v];

Q=orth(A);
A_hat=(A'*Q)';
