A=[0 1;-2 -3];
B=[1;2];
P=[1 1;1 -1];
Ahat= P*A*inv(P)
Bhat=P*B
disp(Ahat,"A^=")
disp(Bhat,"B^=")