如何求

a(1,1)=[ 1 0 0    
         0 0 0 
         0 0 0] 
a(2,1)=[ 0 0 0
         1 0 0
         0 0 0]
a(3,1)=[ 0 0 0
         0 0 0
         1 0 0]

依次类 a(2,1),a(3,2),a(3,3)......

b=eye(5)
G 为3*5的矢量
如何 求

K(1,1)=G'*kron(a(1,1),b)*G

不能用循环,因为矩阵维数很大,循环太多不行!

谢谢