F =全(interpMatrix([1 1 1 1] 1 11日2))”;
P = F(1:6,一14);
Q = F (1:7, 1:16);
R = F (1:9, 1:20);
K = KronProd ({P, Q, R, 1},(1 2 3 4 4],[南、南、南3 3]);
重新= K *;