a=[1 0;0 2];b=[1;3];syms x1 x2;y=a*[x1;x2]+b;s=solve(y(1),y(2),'x1','x2');s.x1 %w2值s.x2 %w2值
操作如下：a=[1 0;0 2];b=[-1; -3];x=inv(a)*b
