《矩阵论引论》田振际——P30-1.4

clear
clc
%《矩阵论引论》田振际 王永铎...
%P30 第4题
format rat
syms a
a1=[1 0 0 0]';
a2=[0 1 0 0]';
a3=[0 0 1 0]';
a4=[0 0 0 1]';
disp('解：基I:')
I=[a1 a2 a3 a4]
b1=[2 1 -1 1]';
b2=[0 3 1 0]';
b3=[5 3 2 1]';
b4=[6 6 1 3]';
disp('基II:')
II=[b1 b2 b3 b4]
disp('（1）由基（I)到基（II)的过渡矩阵：')
%设II=I*P
P=inv(I)*II
disp('（2）向量在基（II）下的坐标：')
A=eye(4);
B=inv(P)*A
disp('（3）对两组基由相同坐标的非零向量：')
%I*X=II*X  (II-I)X=0
Z=a*null([II-I],'r')             %求齐次线性方程组通解
disp('其中a≠0')