[bohr1982]

这么长时间一直没有回复真是对不起！我琢磨出处理方法,和大家分享一下！
题设：
a = [a1 b1 a4 b4]'
b = [a2 b2 a3 b3]'
c = [b1 b2 b3 b4]'
d = [a1 a2 a3 a4]'
a = M1*b
求： 一个M2使c = M2*d

解：
[b1 0 0 b4]' = [0 1 0 0;0 0 0 0;0 0 0 0;0 0 0 1]*a
_________ = [0 1 0 0;0 0 0 0;0 0 0 0;0 0 0 1]*M1*b
[0 b2 b3 0]' = [0 0 0 0;0 1 0 0;0 0 0 1;0 0 0 0]*b
将上面两式相加，得
[b1 b2 b3 b4] = {[0 1 0 0;0 0 0 0;0 0 0 0;0 0 0 1]*M1+[0 0 0 0;0 1 0 0;0 0 0 1;0 0 0 0]}*b
同样的处理方法得到
[a1 a2 a3 a4] = {[1 0 0 0;0 0 0 0;0 0 0 0;0 0 1 0]*M1+[0 0 0 0;0 1 0 0;0 0 0 1;0 0 0 0]}*b
将上两式想除得到
M2 = {[0 1 0 0;0 0 0 0;0 0 0 0;0 0 0 1]*M1+[0 0 0 0;0 1 0 0;0 0 0 1;0 0 0 0]}/{[1 0 0 0;0 0 0 0;0 0 0 0;0 0 1 0]*M1+[0 0 0 0;0 1 0 0;0 0 0 1;0 0 0 0]}