方法1
de2bi(0:3^8-1,[],3)'+1
方法2
[x1,x2,x3,x4,x5,x6,x7,x8]=ndgrid(1:3);[x1(:),x2(:),x3(:),x4(:),x5(:),x6(:),x7(:),x8(:)]'
其中第二句也可以用如下语句代替
reshape(cat(9,x1,x2,x3,x4,x5,x6,x7,x8),[],8)'