通过matlab实现正交表
论文：正交表的构造方法及Matlab 实现 河南师范大学数学与信息科学学院( 453007) 庞善起 鹿姗姗
ps:因为我是照着上述论文进行修改的，所以就当做非原创吧，尊重作者 ，侵权联系必删
%%生成L4(2 3)的正交表为例
----------------------------------------------------------------------------------------------------------------------
clc
clear all
t=2;u=2 %Lt^u(t)
p=t;
V=[]
K=[]
L=[]
D=[]
L1=[]
L2=[]
G=transpose(0:(t-1))
i=1
while(i<=u)
R=kron(kron(ones(t^(i-1),1),G),ones(t^(u-i),1))
V=[V,R]
i=i+1
end
V
U=rref(V')
[m,n]=find(U==1)
[s,t]=find(U((m+1):u,n)<=1);
B=U(:,unique(n(t)))'
L=V*B'
L=mod(L,p)
L=L+ones(size(L,1),size(L,2))
----------------------------------------------------------------------------------------------------------------------
运行结果如下：（有时间再补齐资料）