用Matlab算也可以，写一个函数就搞定了，M文件内容如下：

function [X,V]=fun(B,L,Q)
   P=inv(Q);
   NBB=Rot90(B)*P*B;
   W=Rot90(B)*P*L;
   X=inv(NBB)*W
   V=B*X-L
end

赋好B,L,Q的值后调用fun(B,L,Q)就行了