N=32; M=N/2; temp=0.25; SM=zeros(M*M,N*N); for i=1:1:M for j=1:1:M SM((i-1)*M+j,(2*i-2)*N+2*j-1)=temp; SM((i-1)*M+j,(2*i-2)*N+2*j)=temp; SM((i-1)*M+j,(2*i-1)*N+2*j-1)=temp; SM((i-1)*M+j,(2*i-1)*N+2*j)=temp; end end save SampleMatrix32 SM N M;