矩阵的扩展

rushuidetian

如何将3*3的矩阵扩展成9*9的？
例如：[a1  b1  c1
           a2  b2  c2
           a3  b3  c3]扩展成如下矩阵：

[a1 b1 c1 0   0   0   0   0   0
0   0   0   a1 b1 c1 0   0   0
0   0   0   0   0   0   a1 b1 c1
a2  b2 c2 0   0   0   0   0   0
0    0   0   a2 b2 c2 0   0   0
0    0   0   0   0   0   a2 b2 c2
a3  b3 c3 0   0   0   0   0   0
0    0   0   a3 b3 c3 0   0   0
0    0   0   0   0   0   a3 b3 c3]

最好有命令，因为所处理的问题维数较大！

rushuidetian

怎么没人指导呢？是不是没有这种现成的扩展命令？

wwbeyondww

想不到特别好的方法，以下code可以参考一下:
x=[1 2 3;4 5 6;7 8 9];
a=[1 0 0 ; 0 1 0; 0 0 1];
X=[kron(x,a(:,1)) kron(x,a(:,2)) kron(x,a(:,3))];

rushuidetian

就是这个思想，我自己再想想！谢谢帮忙！:@)

rushuidetian

如何将3*3的矩阵扩展成9*9的？
例如：[a1  b1  c1
             a2  b2  c2
             a3  b3  c3]扩展成如下矩阵：

[a1  b1  c1  b1  0    0       c1    0   0
a2  b2  c2  0    b1  0       0      c1  0
a3  b3  c3  0    0    b1     0      0   c1
a2  0     0   a1  b1  c1     c2    0   0
0    a2   0   a2  b2  c2      0     c2   0
0    0     a2  a3  b3  c3     0     0     c2
a3  0     0    b3  0     0      a1   b1  c1     
0    a3   0      0   b3  0      a2   b2  c2      
0    0     a3    0   0    b3    a3   b3  c3]

最好有命令，因为所处理的问题维数较大

wwbeyondww

clear;
a=[1 2 3; 4 5 6; 7 8 9];
b=[1 0 0; 0 1 0; 0 0 1];
c=tril(a,-1)+triu(a,1);
X=kron(b,a)+kron(c,b);

强烈要求你论文致谢我，:lol

rushuidetian

呵呵，谢谢啦！希望我这个算法能早日全部实现！