可以使用MATLAB自帶的函數 rowexch 生成任意因素和水平的正交表,參考代碼如下:
clear,clc %% 生成正交實驗矩陣 nfactors = 5; nruns = 32; [dRE,X] = rowexch(nfactors,nruns,'interaction','tries',10); dRE(dRE==-1)=0; sortrows(dRE)
程序默認生成的是由-1和1組成的正交表,這里將其轉換為0和1,其他復雜需求可查閱參考資料。
參考資料:
https://ww2.mathworks.cn/help/stats/rowexch.html?searchHighlight=rowexch&s_tid=srchtitle
https://blog.csdn.net/stereohomology/article/details/21938943#opennewwindow