通过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))

 

----------------------------------------------------------------------------------------------------------------------

运行结果如下：（有时间再补齐资料）