for j=1:1:3
L=round(rand*3)+1;    
[A,G,B,H]=mmwgeneratechannelmatrix(L,36,0.1,32,64);
subplot(3,4,4*(j-1)+1)
imagesc(angle(A)*(180/pi))

subplot(3,4,4*(j-1)+2)
imagesc(abs(diag(G)))

subplot(3,4,4*(j-1)+3)
imagesc(angle(B')*(180/pi))

subplot(3,4,4*(j-1)+4)
imagesc(abs(H))
end
